• Journal of the Institute of Electronics Engineers of Korea TC
• /
• v.39 no.4
• /
• pp.201-208
• /
• 2002

In this paper, we present a new formulation using time-domain electric field integral equation (TD-EFIE) to obtain transient scattering response from arbitrarily shaped three-dimensional conducting bodies. The time derivative of the magnetic vector potential is approximated with a central finite difference and the scalar potential is time averaged by dividing it into two terms. This approach with an implicit method using central difference results in accurate and more stable transient scattering responses from conducting objects. Detailed mathematical steps are included and several numerical results are presented and compared with the inverse discrete Fourier transform (IDFT) of the frequency-domain solution.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.19-26
• /
• 1999

This paper describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows with free-surface. The Navier-Stokes equations governing the flows have been discretized by means of finite-difference approximations, and the resulting equations have been solved via the SIMPLE-C algorithm. The free-surface is defined by the motion of a set of marker particles and the interface behaviour was investigated by means of a "Lagrangian" technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.1
• /
• pp.51-60
• /
• 1994

A one-dimensional numerical model has been developed through a Lagragian formulation, which is applicable to advection and diffusion of dissolved matters in storages. The study was conducted to the pollutants released into a canal network in Burnt Store Isles, Punta Gorda, Florida. USA. The hydrodynamic model was developed by using an implicit finite difference scheme. In the computational domain, the network system consists of prismatic channels and storages. The finger canals and small tributaries ale treated as storages in the simplified flow network The numerical results show relatively good agreement with field experiments.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.6
• /
• pp.1163-1170
• /
• 1992

A computer program was developed to analyze the two-dimensional unsteady incompressible viscous flow behind a rectangular bluff body between two parallel plates. The Peaceman-Rachford alternating direction implicit numerical method and Wachspress parameter were adopted to solve the governing equations in vorticity-transport and stream function formulation. The steady state flow and the vortex flow behind a rectangular bluff body in a chemical were investigated for Reynolds numbers of 200 and 500. The vortex shedding was generated by a physical pertubation numerically imposed at the center of the flow field for a short time. It was observed that the perturbed flow became periodic after a transient period.

• The Transactions of the Korean Institute of Electrical Engineers C
• /
• v.51 no.11
• /
• pp.551-558
• /
• 2002

A time-domain combined field integral equation (CFIE) is presented to obtain the transient scattering response from arbitrarily shaped three-dimensional conducting bodies. This formulation is based on a linear combination of the time-domain electric field integral equation (EFIE) with the magnetic field integral equation (MFIE). The time derivative of the magnetic vector potential in EFIE is approximated using a central finite difference approximation and the scalar potential is averaged over time. The time-domain CFIE approach produces results that are accurate and stable when solving for transient scattering responses from conducting objects. The incident spectrum of the field may contain frequency components, which correspond to the internal resonance of the structure. For the numerical solution, we consider both the explicit and implicit scheme and use two different kinds of Gaussian pulses, which may contain frequencies corresponding to the internal resonance. Numerical results for the EFIE, MFIE, and CFIE are presented and compared with those obtained from the inverse discrete Fourier transform (IDFT) of the frequency-domain CFIE solution.

• Structural Engineering and Mechanics
• /
• v.58 no.6
• /
• pp.1021-1044
• /
• 2016

Different versions of a damage index (DI) along with a formulation to find the number of cycles at failure due to fatigue, applicable to reinforced concrete (RC) structures are presented. These are based on an energetic analysis method and applicable to both global and local levels. The required data can be found either from the numerical simulation of structures or from the experimental tests. A computer program has been developed to simulate numerically the nonlinear behavior of RC columns under cyclic loading. The proposed DI gives a regular distribution of structural damages up to failure and is validated by the results of the tests carried out on RC columns subjected to cyclic loading. In general, the local and global damage indices give approximately similar results, while each of them has its own advantages. The advantage of the implicit version of DI is that, it allows the comparison of the results with those of the monotonic loading case, while the explicit version makes it possible to estimate the number of loading cycles at failure due to fatigue, and the advantage of the simplified version is that; the monotonic loading data is not needed for the cyclic loading case.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.7 s.100
• /
• pp.836-842
• /
• 2005

Nonlinear dynamic characteristics of active piezolaminated plates are investigated with respect to the thermopiezoelastic behaviors. For largely deformed structures with small strain, the incremental total Lagrangian formulation is presented based on the virtual work principles. A multi-field layer-wise finite shell element is proposed for assuring high accuracy and non-linearity of displacement, electric and thermal fields. For dynamic consideration of thermopiezoelastic snap-through phenomena, the implicit Newmark's scheme with the Newton-Raphson iteration is implemented for the transient response of various piezolaminated models with symmetric or eccentric active layers. The bifurcate thermal buckling of symmetric structural models is first investigated and the characteristics of piezoelectric active responses are studied for finding snap-through piezoelectric potentials and the load-path tracking map. The thermoelastic stable and unstable postbuckling, thermopiezoelastic snap-through phenomena with several attractors are proved using the nonlinear time responses for various initial conditions and damping loss factors. Present results show that thermopiezoelastic snap-through phenomena can result in the difficulty of buckling and postbuckling control of intelligent structures.

• Journal of Mechanical Science and Technology
• /
• v.19 no.spc1
• /
• pp.388-394
• /
• 2005

Recently, the authors have developed a method for real-time dynamics of multibody systems, which combines a semi-recursive formulation to derive the equations of motion in dependent relative coordinates, along with an augmented Lagrangian technique to impose the loop closure conditions. The following numerical integration procedures, which can be grouped into the so-called structural integrators, were tested : trapezoidal rule, Newmark dissipative schemes, HHT rule, and the Generalized-${\alpha}$ family. It was shown that, for large multi body systems, Newmark dissipative was the best election since, provided that the adequate parameters were chosen, excellent behavior was achieved in terms of efficiency and robustness with acceptable levels of accuracy. In the present paper, the performance of the described method in combination with another group of integrators, the Implicit Runge-Kutta family (IRK), is analyzed. The purpose is to clarify which kind of IRK algorithms can be more suitable for real-time applications, and to see whether they can be competitive with the already tested structural family of integrators. The final objective of the work is to provide some practical criteria for those interested in achieving real-time performance for large and complex multibody systems.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.7 no.5
• /
• pp.53-58
• /
• 1998

This paper concentrates on the development of a computational design program to determine nozzle size in water jet, combing the numerical optimization technique with the flow analysis code. To achieve the above objective, a two-dimensional model was developed for investigating the fluid flow in water jet and calculating the velocity and pressure distributions. The mathematical formulation as a standard ${k}-\varepsilon$ model was solved employing a general thermo fluid-mechanics computer program, PHOENICS code, which is based on the Semi-Implicit Method Pressure Linked Equations(SIMPLE) algorithm. The developed code was applied to water jet design to determine the nozzle size, and investigated the effect of the change of nozzle location. Calculated results showed that the flow pattern is not changed as the change of nozzle location.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.5
• /
• pp.787-797
• /
• 1986

An implicit approach is employed to obtain a general boundary integral formulation of axisymmetric elastic problems in terms of a pair of singular integral equations. The corresponding kernel functions from the solutions of Navier's equation are derived by applying a three dimensional integral and a direct axisymmetrical approach. A numerical discretization schem including the evaluation of Cauchy principal values of the singular integral is described. Finally the typical axisymmetric elastic models are analyzed, i.e. the hollow sphere, the constant thickness and the V-notched round bar.