• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.4
• /
• pp.813-824
• /
• 1988

Use of a plasma arc as the source of energy for penetration welding of thin plates gives rise to a cylindrical hole surrounded by the molten metal. Material moves from the front to the rear of the hole by flowing around the hole as the workpiece is translated relatively to the arc. Based on the finite difference method, three different computer models have been proposed for the steady state, two dimensional heat and mass flow during the plasma arc welding. In the formulation energy equation was derived by the energy blance method through the cell control volume, and all the governing equations derived for the fixed coordinates was translated for the moving coordinate system. The driving force for fluid flow being considered was only electromagnetic force. The calculated and measured molten poon and HAZ width were compared and better agreement was obtained for the models considering the keyhole effect.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.29-37
• /
• 1997

An implicit viscous turbulent flow solver is developed for two-dimensional geometries on unstructured triangular meshes. The flux terms are discretized based on a cell-centered finite-volume formulation with the Roe's flux-difference splitting. The solution is advanced in time using an implicit backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with the Gauss-Seidel relaxation scheme. The effect of turbulence effects is approximated with a standard $k-{\varepsilon}$ two-equation model which is solved separately from the mean flow equations using the same backward-Euler time integration scheme. The triangular meshes are generated using an advancing-front/layer technique. Validations are made for flows over the NACA0012 airfoil and the Douglas 3-element airfoil. Good agreements are obtained between the numerical results and the experiment.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.4
• /
• pp.493-499
• /
• 2007

The Taylor expansion expresses a differentiable function and its coefficients provide good approximations for the given function and its derivatives. In this study, m-th order Taylor Polynomial is constructed and the coefficients are computed by the Moving Least Squares method. The coefficients are applied to the governing partial differential equation for solid problems including crack problems. The discrete system of difference equations are set up based on the concept of point collocation. The developed method effectively overcomes the shortcomings of the finite difference method which is dependent of the grid structure and has no approximation function, and the Galerkin-based meshfree method which involves time-consuming integration of weak form and differentiation of the shape function and cumbersome treatment of essential boundary.

• 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.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.19 no.4
• /
• pp.427-435
• /
• 2008

In this paper, a stable marching-on in time(MOT) method with 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 derivatives in the EFIE and MFIE are approximated using a central finite difference scheme and other terms are averaged over time. This time domain CFIE approach produces results that are accurate and stable when solving for transient scattering responses from conducting objects. Numerical results with the proposed MOT scheme are presented and compared with those obtained from the conventional method and the inverse discrete Fourier transform(IDFT) of the frequency domain CFIE solution.

• Economic and Environmental Geology
• /
• v.27 no.3
• /
• pp.281-287
• /
• 1994

The present paper deals with numerical modeling of Love channel waves in transversely isotropic elastic medium. First, an explicit finite-difference scheme of second order approximation is formulated with the wave equation of SH particle displacement in transversely isotropic medium. Since it is a heterogeneous formulation, it should enable efficient modeling of complex model structures without additional treatment of the internal boundary matching. With a model of isotropic coal seam embedded in high velocity host rock, seismograms are synthesized and tutn out to be essentially identical with published ones of Korn and $St{\ddot{o}}ckl$. Next, anisotropic coal seams are investigated. It is found that the horizontal velocity of the seam appears to play a major role of determining the group velocity of Love channel waves. The group velocity increases with the increase of the horizontal velocity or vice versa. However, further study will be needed to exploit fully Love channel waves for the determination of lithology, stratification, fracture in sedimentary rocks, for instance, for hydrocarbon exploration and development.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.1
• /
• pp.141-152
• /
• 1993

A generalized method is presented for shape design sensitivity analysis of axisymmetric thermal conducting solids. The shape sensitivity formula of a general performance functional arising in shape optimal design problem is derived using the material derivative concept and the adjoint variable method. The method for deriving the formula is based on standard axisymmetric boundary integral equation formulation. It is then applied to obtain the sensitivity formulas for temperature and heat flux constraints imposed over a small segment of the boundary. To show the accuracy of the sensitivity analysis, numerical implementations are done for three examples. Sensitivities calculated by the presented method are compared with analytic sensitivities for two examples with analytic solutions, and compared with sensitivies by finite difference for a cooling fin example.

• Proceedings of the Korean Institute of Information and Commucation Sciences Conference
• /
• v.9 no.1
• /
• pp.363-366
• /
• 2005

A micro flow sensor on silicon substrate allows the fabrication of small components where many different functions can be integrated so that the functionality of the sensors can be increased. Further more, due to the small size of the elements the sensors can be quite fast. A thermal mass flow sensor measures the asymmetry of temperature profile around the heater which is modulated by the fluid flow. In normal, a mass flow sensor is composed of a central heater and a pair of temperature sensing elements around the heater. A new 2-D wide range micro flow sensor structure with three pairs of temperature sensors and a central heater was proposed and numerically simulated by the Finite difference formulation to confirm the feasibility of the flow sensor structure.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.267-277
• /
• 1993

A finite element method(FEM) is presented and applied to the two-dimensional bottom turbulent boundary layer. The time-dependent incompressible motion of a viscous fluid is formulated by using the well-known Navier-Stokes equations and vorticity equation in terms of the velocity and pressure fields. The general numerical formulation is based on Galerkin method and solved by introducing the mixing length theory of Prandtl for eddy kinematic viscosity of a turbulent flow field. Numerical computations of the transport of sediment on an arbitrary sea-bed due to wave motion in the turbulent boundary layer are carried out. The results obtained by the FEM made clear the difference in characteristic features between the boundary layer due to oscillatory flow and the boundary layer due to wave motion.

• Ocean Systems Engineering
• /
• v.10 no.4
• /
• pp.399-414
• /
• 2020

A floating bridge is an innovative solution for deep-water and long-distance crossing. This paper presents a curved floating bridge's dynamic behaviors under the wind, wave, and current loads. Since the present curved bridge need not have mooring lines, its deep-water application can be more straightforward than conventional straight floating bridges with mooring lines. We solve the coupled interaction among the bridge girders, pontoons, and columns in the time-domain and to consider various load combinations to evaluate each force's contribution to overall dynamic responses. Discrete pontoons are uniformly spaced, and the pontoon's hydrodynamic coefficients and excitation forces are computed in the frequency domain by using the potential-theory-based 3D diffraction/radiation program. In the successive time-domain simulation, the Cummins equation is used for solving the pontoon's dynamics, and the bridge girders and columns are modeled by the beam theory and finite element formulation. Then, all the components are fully coupled to solve the fully-coupled equation of motion. Subsequently, the wet natural frequencies for various bending modes are identified. Then, the time histories and spectra of the girder's dynamic responses are presented and systematically analyzed. The second-order difference-frequency wave force and slowly-varying wind force may significantly affect the girder's lateral responses through resonance if the bridge's lateral bending stiffness is not sufficient. On the other hand, the first-order wave-frequency forces play a crucial role in the vertical responses.