• 한국지하수토양환경학회지:지하수토양환경
• /
• 제11권1호
• /
• pp.37-44
• /
• 2006

본 연구는 오염물 거동에 대한 수치해석을 위해 보편적으로 사용되고 있는 수치 방법들의 장단점을 총괄적으로 나타내고, 효율적인 수치모델링 기법 개발을 위해 ELLAM과 LEZOOMPC를 비교분석하였다. 지하수 분야에서 가장 많이 사용되는 수치 방법은 Eulerian-Lagrangian 방식과 Eulerian 방식인데, Eulerian-Lagrangian 방식은 수치영역 내에서 일반적으로 질량을 보존하지 못하고, 경계조건을 체계적으로 처리하지 못하는 한계를 갖고 있다. 반면에 Eulerian 빙식은 시간 및 공간 절삭 오차로 인해서 시간 간격 및 격자 크기를 극히 줄여야 하는 제약을 갖고 있다. 최근 10 년간 지하수 분야에서 크게 대두되고 있는 수치기법인 ELLAM(Eulerian Lagrangian Localized Adjoint Method)은 Eulerian-Lagrangian 방식과 Eulerian 방식에서 나타나는 수치 제약점이나 한계점을 동시에 해결하는 수치기법으로 알려져 왔다. 그러나 본 연구에서는 ELLAM의 장단점을 파악하고 보완점을 제안한다. ELLAM의 단점을 파악하기 위해, mesh Peclet number가 다른 예제들을 설정하고, 그 예제들에 대한 ELLAM, LEZOOMPC(Lagrangian-Eulerian ZOOMing Peak and valley Capturing)와 visual MODFLOW의 수치결과들을 해석해와 비교하였다. Mesh Peclet number가 무한대일 때 ELLAM의 수치결과는 수치진동으로 인해 해석해와 일치하지 않았으나, LEZOOMPC의 수치 결과는 해석해와 일치했다. 위의 결과는 ELLAM의 수치오차가 LEZOOMPC의 특성을 이용하여 개선 및 보완될 수 있는 가능성을 시사해 준다. 따라서 ELLAM에 LEZOOMPC의 후향 입지추적, 전향 입지추적, 선택적 국부 격자 세립화 과정과 최고/최저 농도점 이동 추적 과정을 결합하면 ELLAM의 수치적 장점을 유지하면서 mesh Peclet number에 제약을 받지 않는 효율적인 수치모델링 기법을 개발할 수 있을 것으로 판단된다.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2006년도 정기 학술대회 논문집
• /
• 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.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 1993년도 수공학연구발표회논문집
• /
• pp.151-157
• /
• 1993

• 한국환경과학회지
• /
• 제11권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
• /
• 제6권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
• /
• 제55권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.

• 한국전산유체공학회지
• /
• 제6권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.

• 한국지반공학회논문집
• /
• 제17권6호
• /
• pp.85-98
• /
• 2001

Arbitrary Lagrangian-Eulerian(ALE) 기법에 의거하여 콘 관입문제를 해석하였다. 완전한 관입을 모의하기 위하여 ABAQUS/Explicit을 이용하여 지반의 상향 유동을 모델링하는 정상상태해석(steady state analysis)을 수행하였다. 단일 지층의 해석에서는 흙 입자의 유동 경로와 변형률이 strain path method와 일관된 결과를 나타내고 있음을 확인하였고 극한저항을 합리적으로 계산할 수 있었다. 상이한 지층에 콘을 관입하는 경우에 대해서도 콘 저항이 전이하는 경향을 해석할 수 있었다. 따라서 ALE 해법의 정상상태해석으로 층상 지층에 대한 완전한 관입을 해석할 수 있었다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 추계 학술대회논문집
• /
• 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

• 소성∙가공
• /
• 제4권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.