• Journal of the Korean GEO-environmental Society
• /
• v.4 no.1
• /
• pp.5-10
• /
• 2003

The solids and the fluids in porous media have a relative velocity to each other. Due to physically and chemically different material properties and their relative velocity, the behavior of saturated porous media is extremely complicated. Thus, in order to describe and clarify the deformation behavior of saturated porous media, constitutive models for deformation of porous media coupling several effects need to be developed in frame of Arbitrary Lagrangian Eulerian(ALE) description. The aim of ALE formulations is to maximize the advantages of Lagrangian and Eulerian elements, and to minimize the disadvantages. Therefore, this method is appropriate for the analysis of porous media which are considered for the behavior of the solids and the fluids. For this reason, mass balance equations for saturated porous media are derived here in ALE description frames. ALE formulations of mass conservation for the solid phase and the fluid phase are expressed. Then, linear momentum balance equation for porous media as multiphase media is expressed.

• Journal of Korea Water Resources Association
• /
• v.32 no.2
• /
• pp.99-110
• /
• 1999

The method of moments, a Lagrangian scheme, considers the zeroth, first, and second moments of the grid cell spatial distributions of the concentration and then advects the concentration by maintaining conservation of the moments. The reasonable inital description of the first and second moments as well as the mean concentration, the zeroth moments, in grid element is important in the method of moments. In this study, the description methods of each initial moment are reviewed, and the method of moments is extended to overcome the restrictions of Courant number. Its performance is compared with those of available Eulerian and Lagrangian schemes. As the results, the method is successfully extended to overcome the stability restriction and is an accurate scheme for the advection simulation of concentration distribution, especially of which the gradient is steep. In addition, the method is very promising scheme in terms of computational efficiency when the mixing is confined in a relatively small region to the entire domain in two-dimensional problem.

• Journal of the Korea Computer Graphics Society
• /
• v.25 no.5
• /
• pp.21-30
• /
• 2019

We propose a physical simulation method based on the shallow water equation(SWE) to represent water surface effectively. In this paper, the water which can be represented has a much larger width compared to the depth does not have a large vertical direction flow. In order to calculate the water flow efficiently, we start with the shallow water equation as the governing equation, which is a simplified version of the Navier-Stokes equation. In order to numerically calculate the advection term of the SWE, we introduce a new conservtive USCIP(CUSCIP) method which improves the Constrained Interpolation Profile (CIP) method to preserve the physical quantity while increasing the numerical accuracy. The proposed method is based on Kim et. al.'s Unsplit Semi-lagrangian CIP[9], and calculates advection term with additional constraints on term that consider integral values. The experimental results show that the CUSCIP method is robust to the loss of physical quantity due to numerical dissipation, which improves wave detail and persistence.

• Water for future
• /
• v.28 no.5
• /
• pp.151-162
• /
• 1995

A two-dimensional stream tube dispersion model is developed to simulate accurately transverse dispersion processes of pollutants in natural channels. Two distinct features of the stream tube dispersion model derived herein are that it employs the transverse cumulative discharge as an independent variable replacing the transverse distance and that it is developed in a natural coordinate system which follows the general direction of the channel flow. In the model studied, Eulerian-Lagrangian method is used to solve the stream tube dispersion equation. The stream tube dispersion equation is decoupled into two components by the operator-splitting approach; one is governing advection and the other is governing dispersion. The advection equation has been solved using the method of characteristics and the results are interpolated onto Eulerian grid on which the dispersion equation is solved by centered difference method. In solving the advection equation, cubic spline interpolating polynomials is used. In the present study, the results of the application of this model to a natural channel are compared with a steady-state flow measurements. Simulation results are in good accordance with measured data.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.12
• /
• pp.2306-2314
• /
• 1992

Using the Lagrangian equation, nonlinear vibration analysis of laminated hybrid composite plates is carried out. The effects of stacking sequences, aspect ratios, number of modes, number of layers and various elastic properties on nonlinear vibration are investigated. The presence of bending-extension coupling in antisymmetric plates yields a second power term in addition to a cubic nonlinear term in governing differential equation of motion. In the other symmetric case, this second term vanishes. The fundamental frequency of analytic results are compared with that of ABAQUS FEM analysis. For nonlinear vibration of antisymmetric unimaterial plate, the result of reference is presented for comparison with this result.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.186-191
• /
• 2004

The reactive flowfield of the transverse injecting combustor has been studied using Euler-Lagrange method in order to develop an efficient solution procedure for the understanding of liquid spray combustion in the transverse injecting combustor which has been widely used in ramjets and turbojet afterburners. The unsteady two-dimensional gas-phase equations have been represented in Eulerian coordinates and the liquid-phase equations have been formulated in Lagrangian coordinates. The gas-phase equations based on the conservation of mass, momentum, and energy have been supplemented by combustion. The vaporization model takes into account the transient effects associated with the droplet heating and the liquid-phase internal circulation. The droplet trajectories have been determined by the integration of the Lagrangian equation in the flow field obtained from the separate calculation without considering the iterative effect between liquid and gas phases. The reported droplet trajectories had been found to deviate from the initial conical path toward the flow direction in the very end of its lifetime when the droplet size had become small due to evaporation. The integration scheme has been based on the TEACH algorithm for gas-phase equation, the second order Runge-Kutta method for liquid-phase equations and the linear interpolation between the two coordinate systems. The calculation results has shown that the characteristics of the droplet penetration and recirculation have been strongly influenced by the interaction between gas and liquid phases in such a way that most of the vaporization process has been confined to the wake region of the injector, thereby improving the flame stabilization properties of the flowfield.

• Journal of Mechanical Science and Technology
• /
• v.17 no.3
• /
• pp.440-448
• /
• 2003

In the present paper, numerical simulations of buoyant plume dispersion in a neutral and stable atmospheric boundary layer have been carride out. A Lagrangian Stochastic Model (LSM) with a Non-Linear Eddy Viscosity Model (NLEVM) for turbulence is used to generate a Reynolds stress field as an input condition of dispersion simulation. A modified plume-rise equation is included in dispersion simulation in order to consider momentum effect in an initial stage of plume rise resulting in an improved prediction by comparing with the experimental data. The LSM is validated by comparing with the prediction of an Eulerian Dispersion Model (EDM) and by the measured results of vertical profiles of mean concentration in the downstream of an elevated source in an atmospheric boundary layer. The LSM predicts accurate results especially in the vicinity of the source where the EDM underestimates the peak concentration by 40% due to inherent limitations of gradient diffusion theory. As a verification study, the LSM simulation of buoyant plume dispersions under a neutral and stable atmospheric condition is compared with a wind-tunnel experiment, which shows good qualitative agreements.

• Structural Engineering and Mechanics
• /
• v.47 no.5
• /
• pp.599-622
• /
• 2013

Seismic response of two dimensional liquid tanks is numerically simulated using fully nonlinear velocity potential theory. Galerkin-weighted-residual based finite element method is used for solving the governing Laplace equation with fully nonlinear free surface boundary conditions and also for velocity recovery. Based on mixed Eulerian-Lagrangian (MEL) method, fourth order explicit Runge-Kutta scheme is used for time integration of free surface boundary conditions. A cubic-spline fitted regridding technique is used at every time step to eliminate possible numerical instabilities on account of Lagrangian node induced mesh distortion. An artificial surface damping term is used which mimics the viscosity induced damping and brings in numerical stability. Four earthquake motions have been suitably selected to study the effect of frequency content on the dynamic response of tank-liquid system. The nonlinear seismic response vis-a-vis linear response of rectangular liquid tank has been studied. The impulsive and convective components of hydrodynamic forces, e.g., base shear, overturning base moment and pressure distribution on tank-wall are quantified. It is observed that the convective response of tank-liquid system is very much sensitive to the frequency content of the ground motion. Such sensitivity is more pronounced in shallow tanks.

• Journal of Korean Society for Atmospheric Environment
• /
• v.15 no.6
• /
• pp.757-765
• /
• 1999

Lagrangian particle dispersion model(LPDM) is an effective tool to calculate the dispersion from a point source since it dose not induce numerical diffusion errors in solving the pollutant dispersion equation. Fictitious particles are released to the atmosphere from the emission source and they are then transported by the mean velocity and diffused by the turbulent eddy motion in the LPDM. The concentration distribution from the dispersed particles in the calculation domain are finally estimated by applying a particle count method or a Gaussian kernel method. The two methods for calculating concentration profiles were compared each other and tested against the analytic solution and the tracer experiment to find the strength and weakness of each method and to choose computationally time saving method for the LPDM. The calculated concentrations from the particle count method was heavily dependent on the number of the particles released at the emission source. It requires lots fo particle emission to reach the converged concentration field. And resulting concentrations were also dependent on the size of numerical grid. The concentration field by the Gaussian kernel method, however, converged with a low particle emission rate at the source and was in good agreement with the analytic solution and the tracer experiment. The results showed that Gaussian kernel method was more effective method to calculate the concentrations in the LPDM.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.1
• /
• pp.83-106
• /
• 2016

This paper reviews and compares three different methods for modeling incompressible and immiscible ternary fluid flows: the immersed boundary, level set, and phase-field methods. The immersed boundary method represents the moving interface by tracking the Lagrangian particles. In the level set method, an interface is defined implicitly by using the signed distance function, and its evolution is governed by a transport equation. In the phase-field method, the advective Cahn-Hilliard equation is used as the evolution equation, and its order parameter also implicitly defines an interface. Each method has its merits and demerits. We perform the several simulations under different conditions to examine the merits and demerits of each method. Based on the results, we determine the most suitable method depending on the specific modeling needs of different situations.