• Journal of Soil and Groundwater Environment
• /
• v.11 no.1
• /
• pp.37-44
• /
• 2006

This study summarizes advantages and disadvantages of numerical methods and compares ELLAM and LEZOOMPC to develop an efficient numerical modeling technique on contaminant transport. Eulerian-Lagrangian method and Eulerian method are commonly used numerical techniques. However Eulerian-Lagrangian method does not conserve mass globally and fails to treat boundary in a straightforward manner. Also, Eulerian method has restrictions on the size of Courant number and mesh Peclet number because of time truncation error. ELLAM (Eulerian Lagrangian Localized Adjoint Method) which has been popularly used for past 10 years in numerical modeling, is known for overcoming these numerical problems of Eulerian-Lagrangian method and Eulerian method. However, this study investigates advantages and disadvantages of ELLAM and suggests a change for the better. To figure out the disadvantages of ELLAM, the results of ELLAM, LEZOOMPC (Lagrangian-Eulerian ZOOMing Peak and valley Capturing), and visual MODFLOW are compared for four examples having different mesh Peclet numbers. The result of ELLAM generates numerical oscillation at infinite of mesh Peclet number, but that of LEZOOMPC yields accurate simulations. The simulation results suggest that the numerical error of ELLAM could be alleviated by adopting some schemes in LEZOOMPC. In other words, the numerical model which combines ELLAM with backward particle tracking, forward particle tracking, adaptively local zooming, and peak/valley capturing of LEZOOMPC can be developed for not only overcoming the numerical error of ELLAM, but also keeping the numerical advantage of ELLAM.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.381-388
• /
• 2006

A porous medium is composed of solids, fluids, and gas which have different physical and chemical properties. In addition, these constituents have a relative velocity between each other. So far, in order to analyze porous media using finite element method, Lagrangian or Eulerian method has been used. However, the numerical analyses for porous media have a defect that the methods do not describe the movements of constituents. In this paper, numerical analysis for unsaturated porous media was performed in frame of ALE method which has advantages of Lagrangian and Eulerian. Namely, the Lagrangian description was used in solid phase, and the Eulerian description was used in fluid or gas phase in a porous medium Then the relationship between each other was controlled by the convective term in ALE method. Finally, the numerical results of ALE were compared with tile results of Lagrangian analysis.

• Proceedings of the Korea Water Resources Association Conference
• /
• 1993.07a
• /
• pp.151-157
• /
• 1993

• Journal of Environmental Science International
• /
• v.11 no.9
• /
• pp.907-914
• /
• 2002

Various Eulerian-Lagrangian models for the one-dimensional longitudinal dispersion equation in nonuniform flow were studied comparatively. In the models studied, the transport equation was decoupled into two component parts by the operator-splitting approach; one part is governing advection and the other is governing dispersion. The advection equation has been solved by using the method of characteristics following fluid particles along the characteristic line and the results were interpolated onto an Eulerian grid on which the dispersion equation was solved by Crank-Nicholson type finite difference method. In the solution of the advection equation, Lagrange fifth, cubic spline, Hermite third and fifth interpolating polynomials were tested by numerical experiment and theoretical error analysis. Among these, Hermite interpolating polynomials are generally superior to Lagrange and cubic spline interpolating polynomials in reducing both dissipation and dispersion errors.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.6 no.1
• /
• pp.1-15
• /
• 2002

Advection-dominated transport problems possess difficulties in the design of numerical methods for solving them. Because of the hyperbolic nature of advective transport, many characteristic numerical methods have been developed such as the classical characteristic method, the Eulerian-Lagrangian method, the transport diffusion method, the modified method of characteristics, the operator splitting method, the Eulerian-Lagrangian localized adjoint method, the characteristic mixed method, and the Eulerian-Lagrangian mixed discontinuous method. In this paper relationships among these characteristic methods are examined. In particular, we show that these sometimes diverse methods can be given a unified formulation. This paper focuses on characteristic finite element methods. Similar examination can be presented for characteristic finite difference methods.

• Nuclear Engineering and Technology
• /
• v.55 no.12
• /
• pp.4395-4407
• /
• 2023

Bubble columns are widely encountered in several industries, especially in the field of nuclear safety. The Eulerian-Eulerian and the Eulerian-Lagrangian methods are commonly used to investigate bubble columns. Eulerian approaches require additional tasks such as strict volume conservation at the interface and a predefined well-structured grid. In contrast, the Lagrangian approach can be easily implemented. Hence, we introduce a fully Lagrangian approach for the simulation of bubble columns using the discrete bubble model (DBM) and moving particle semi-implicit (MPS) methods. Additionally, we propose a rigorous method to estimate the volume fraction accurately, and verified it through experimental data and analytical results. The MPS method was compared with the experimental data of Dambreak. The DBM was verified by analyzing the terminal velocity of a single bubble for each bubble size. It agreed with the analytical results for each of the four drag correlations. Additionally, the improved method for calculating the volume fraction showed agreement with the Ergun equation for the pressure drop in a packed bed. The implemented MPS-DBM was used to simulate the bubble column, and the results were compared with the experimental results. We demonstrated that the MPS-DBM was in quantitative agreement with the experimental data.

• Journal of computational fluids engineering
• /
• v.6 no.3
• /
• pp.19-26
• /
• 2001

An Eulerian-Lagrangian method, so called immersed boundary method, is used for analysing viscous flow around arbitrary bodies, where governing equations are discretized on a regular grid by using a finite volume method. To improve the accuracy of flow near body boundaries, a second-order accurate interpolation scheme is used and a level-set based grid deformation method is presented to construct the adaptive grids around body boundaries. The present scheme is used to simulate steady flow around a semicircular cylinder mounted on the bottom of flow domain and calculated results are validated by results of a body fitted grid method. Finally, present method is applied to a complex flow around multi body and the usefulness is checked by investigating calculated results.

• Journal of the Korean Geotechnical Society
• /
• v.17 no.6
• /
• pp.85-98
• /
• 2001

Cone penetration was analyzed by arbitrary Lagangian-Eulerian(ALE) method. In order to simulate full penetration, steady state analyses were performed using ABAQUS/Explicit, which models upward flow of soil layers. In the analysis of homogeneous layer it was found that the paths and the strain of soil particles were consistent with the result of the strain path method and that the ultimate resistance were reasonably evaluated. The cone penetration through different soil layers was also analyzed and that showed the transfer of cone resistance. The steady state ALE analysis could perform full penetration through the layered soils.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.43-48
• /
• 2003

When the Eulerian-Lagrangian method is used to analyze the particle laden two-phase flow, a large number of particles should be used to obtain statistically meaningful solutions. Then it takes too much time to track the particles and to average the particle properties in the numerical analysis of two-phase flow. The purpose of this paper is to reduce the computation time by means of a set of particle gird separate to the flow grid. Particle motion equation here is the simplified B-B-O equation, which is integrated to get the particle trajectories. Particle turbulent dispersion, wall collision, and wall roughness effects are considered but the two-way coupling effects between gas and particles are neglected. Particle laden 2-D channel flow is solved and it is shown that the computational efficiency is indeed improved by using the current method

• Transactions of Materials Processing
• /
• v.4 no.1
• /
• pp.69-78
• /
• 1995

The arbitrary Lagrangian-Eulerian(ALE) finite element analysis is applied to the axisymmetric hot extrusion through continuous dies. In order to simulate hot forming problems, an ALE scheme for temperature analysis is proposed. The computed results are compared with experimental results as with those by pure Lagrangian method. In the present study mesh control is accomplished by the use of isoparametric mapping of quadrilaterals.