• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.334-336
• /
• 2008

A novel adaptive mesh refinement (AMR) strategy based on the Moment-of-Fluid (MOF) method for volume-tracking dynamic interface computation is presented. The Moment-of-Fluid method is a new interface reconstruction and volume advection method using volume fraction as well as material centroid. The mesh refinement is performed based on the error indicator, the deviation of the actual centroid obtained by interface reconstruction from the reference centroid given by moment advection process. Using the AMR-MOF method, the accuracy of volume-tracking computation with evolving interfaces is improved significantly compared to other published results.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.334-336
• /
• 2008

A novel adaptive mesh refinement (AMR) strategy based on the Moment-of-Fluid (MOF) method for volume-tracking dynamic interface computation is presented. The Moment-of-Fluid method is a new interface reconstruction and volume advection method using volume fraction as well as material centroid. The mesh refinement is performed based on the error indicator, the deviation of the actual centroid obtained by interface reconstruction from the reference centroid given by moment advection process. Using the AMR-MOF method, the accuracy of volume-tracking computation with evolving interfaces is improved significantly compared to other published results.

• Journal of computational fluids engineering
• /
• v.14 no.1
• /
• pp.18-23
• /
• 2009

A novel adaptive mesh refinement(AMR) strategy based on the Moment-of-Fluid(MOF) method for volume-tracking dynamic interface computation is presented. The Moment-of-Fluid method is a new interface reconstruction and volume advection method using volume fraction as well as material centroid. The adaptive mesh refinement is performed based on the error indicator, the deviation of the actual centroid obtained by interface reconstruction from the reference centroids given by moment advection process. Using the AMR-MOF method, the accuracy of volume-tracking computation with evolving interfaces is improved significantly compared to other published results.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.3
• /
• pp.279-290
• /
• 2010

In this paper, a VOF method called RHRIC (refined high resolution intertace capturing) is introduced for solving the motion of the free surface and applied to the simulation of the advection of rigid interiaces of different shapes and a 20 dam-break problem, which are typical benchmark test cases. The numerical results for the interface advection cases are compared to the analytic solutions, while the available experimental data and other numerical results of various free surface methods for the dam-break problem are provided for the validation of the proposed VOF method. The same simulations were also carried out using the original HRIC scheme and a modified HRIC scheme called MHRIC for comparison. Although the RHRIC uses a simple order scheme, a basis of the original HRIC scheme, lower than the third-order ULTIMATE-QUICKEST used by the MHRIC, it provides an improved accuracy over the two previous HRIC methods.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.109-114
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.675-682
• /
• 2010

The flood impact pressure acting on a vertical wall resulting from a dam-breaking problem is simulated using a navier-Stokes(N-S) solver. The N-S solver uses Eulerian Finite Volume Method(FVM) along with Volume Of Fluid(VOF) method for 2-D incompressible free surface flows. A Split Lagrangian Advection(SLA) scheme for VOF method is implemented in this paper. The SLA scheme is developed based on an algorithm of Piecewise Linear Interface Calculation(PLIC). The coupling between the continuity and momentum equations is affected by using a well-known Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm. Several two-dimensional numerical simulations of the dam-breaking problem are presented to validate the accuracy and demonstrate the capability of the present algorithm. The significance of the time step and grid resolution are also discussed. The computational results are compared with experimental data and with computations by other numerical methods. The results showed a favorable agreement of water impact pressure as well as the global fluid motion.

• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.25-31
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP method. The result using improved CIP method shows the better result of the experiments, comparison with result of original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.568-571
• /
• 2008

The computation of moving interface by the level set method typically requires reinitializations of level set function. An inaccurate estimation of level set function ${\phi}$ results in incorrect free-surface capturing and thus errors such as mass gain/loss. Therefore, accurate and robust reinitialization process is essential to the free-surface flows. In the present paper, we pursue further development of the reinitialization process, which evaluates directly level set function ${\phi}$ using a normal vector in the interface without solving the re-distancing equation of hyperbolic type. The Taylor-Galerkin approximation and P1P1splitting FEM are adopted to discretize advection equation of the level set function and the Navier-Stokes equation, respectively. Advection equation of free surface and re-initialization process are validated with benchmark problems, i.e., a broken dam flow and time-reversed single vortex flow. The simulation results are in good agreement with the existing results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.10
• /
• pp.754-760
• /
• 2008

Computation of moving interface by the level set method typically requires the reinitialization of level set function. An inaccurate estimation of level set function $\phi$ results in incorrect free-surface capturing and thus errors such as mass gain/loss. Therefore, an accurate and robust reinitialization process is essential to the simulation of free-surface flows. In the present paper, we pursue further development of the reinitialization process, which evaluates level set function directly using a normal vector on the interface without solving there-distancing equation of hyperbolic type. The Taylor-Galerkin approximation and P1P1 splitting/SUPG (Streamline Upwind Petrov-Galerkin) FEM are adopted to discretize advection equation of the level set function and the incompressible Navier-Stokes equation, respectively. Advection equation and re-initialization process of free surface capturing are validated with benchmark problems, i.e., a broken dam flow and timereversed single vortex flow. The simulation results are in good agreement with the existing results.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2008.05a
• /
• pp.235-237
• /
• 2008

A unified hydrocode, ExLO, in which Largrangian, ALE and Eulerian solvers are incorporated into a single framework, has recently been developed in Korea. It is based on the three dimensional explicit finite element method and written in C++. ExLO is mainly designed for the calculation of structural responses to highly transient loading conditions, such as high-speed impacts, high-speed machining, high speed forming and explosions. In this paper the numerical schemes are described. Some improvements of the material interface and advection scheme are included. Details and issues of the momentum advection scheme are provided. In this paper the modeling capability of ExLO has been described for two extreme loading events; high-speed impacts and explosions. Numerical predictions are in good agreement with the existing experimental data. Specific applications of the code are discussed in a separate paper in this journal. Eventually ExLO will be providing an optimum simulation environment to engineering problems including the fluid-structure interaction problems, since it allows regions of a problem to be modeled with Lagrangian, ALE or Eulerian schemes in a single framework.