• Journal of computational fluids engineering
• /
• v.13 no.1
• /
• pp.28-34
• /
• 2008

A particle method recognized as one of the gridless methods has been developed to investigate the nonlinear free-surface motions interacting to the structures. The method is more feasible and effective than convectional grid methods for solving the non-linear free-surface motion with complicated boundary shapes. The right-handed side of the governing equations for incompressible fluid, which includes gradient, viscous and external force terms, can be replaced by the particle interaction models. In the present study, the developed method is applied to the dam-broken problem on dried- and wet-floor and its adequacy will be discussed by the comparison with the experimental results.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.258-263
• /
• 2007

A particle method recognized as one of the gridless methods has been developed to investigate the nonlinear free-surface motions interacting to the structures. The method is more feasible and effective than convectional grid methods for solving the non-linear free-surface motion with complicated boundary shapes. The right-handed side of the governing equations for incompressible fluid, which includes gradient, viscous and external force terms, can be replaced by the particle interaction models. In the present study, the developed method is applied to the dam-broken problem on dried- and wet-floor and its adequacy will be discussed by the comparison with the experimental results.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.1
• /
• pp.11-19
• /
• 2010

A code is developed to simulate incompressible free surface flows using the Roe's flux-difference splitting scheme. An interface of two fluids is considered as a moving contact discontinuity. The continuities of pressure and normal velocity across the interface are enforced by the conservation law in the integral sense. The fluxes are computed using the Roe's flux-difference splitting scheme for two incompressible fluids. The interface can be identified based on the computed density distribution. However, no additional treatment is required along the interface during the whole computations. Complicated time evolution of the interface including topological change can be captured without any difficulties. The developed code is applied to simulate the Rayleigh-Taylor instability of two incompressible fluids in the density ratio of 7.2:1 and the broken dam problem of water-air. The present results are compared with other available results and good agreements are achieved for the both cases.

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

• Journal of Ocean Engineering and Technology
• /
• v.21 no.6
• /
• pp.53-58
• /
• 2007

A particle method, recognized as one of gridless methods, has been developed to investigate non-linear free-surface motions interacting with structures. This method is more feasible and effective than conventional grid methods for solving flow fieldswith complicated boundary shapes. The method consists of particle interaction models representing pressure gradient, diffusion, incompressibility, and the free-surface boundary conditions without grids. In the present study, broken dam problems with various viscosity values are simulated to validate the developed method.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.6
• /
• pp.388-396
• /
• 2017

A Lagrangian approach based computational fluid dynamics (CFD) was used to simulate large and/or sharp deformations and fragmentations of interfaces, including free surfaces, through tracing each particle with physical quantities. According to the concept of the particle-based CFD method, it is possible to apply it to both fluid particles and solid particles such as sand, gravel, and rock. However, the presence of more than two different phases in the same domain can make it complicated to calculate the interaction between different phases. In order to solve multiphase problems, particle interaction models for multiphase problems, including surface tension, buoyancy-correction, and interface boundary condition models, were newly adopted into the moving particle semi-implicit (MPS) method. The newly developed MPS method was used to simulate a typical validation problem involving dam breaking. Because the soil and other particles, excluding the water, may have different viscosities, various viscosity coefficients were applied in the simulations for validation. The newly developed and validated MPS method was used to simulate the mobile beds induced by broken dam flows. The effects of the viscosity on soil particles were also investigated.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.165-169
• /
• 2009

The effect of the Dirichlet boundary condition for the redistance equation of level set method on the solutionof sloshing problem is investigated by adopting four Dirichlet boundary conditions. For the solution of the incompressible Navier-Stokes equations, P1P1 four-step fractional finite element method is employed and a least-square finite element method is used for the solutions of the two hyperbolic type equations of level set method; advection and redistance equation. ALE (Arbitrary Lagrangian Eulerian) method is used to deal with a moving computational domain. It has been shown that the free surface motion in a sloshing tank is strongly dependent on the type of the Dirichlet boundary condition and the results of broken dam and sloshing problems using various Dirichlet boundary conditions are discussed and compared with the existing experimental results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.9
• /
• pp.1230-1243
• /
• 1997

A numerical technique for simulating incompressible viscous flow with free surface is presented. The flow field is obtained by penalty finite element formulation. In this work, a modified volume of fluid (VOF) method which is compatible with 4-node element is proposed to track the moving free surface. This scheme can be applied to irregular mesh system, and can be easily extended to three dimensional geometries. Numerical analyses were done for two benchmark examples, namely the broken dam problem and the solitary wave propagation problem. The numerical results were in close agreement with the existing data. Illustrative examples were studied to show the effectiveness of the proposed numerical scheme.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.12
• /
• pp.1555-1569
• /
• 2000

Numerical simulation of fluid flow with moving free surface has been carried out. For the free surface flow, a VOF(Volume of Fluid)-based algorithm utilizing a fixed grid system has been investigated. In order to reduce numerical smearing at the free surface represented on a fixed grid system, a new free surface tracking algorithm based on the donor-acceptor scheme has been presented. Novel features of the proposed algorithm are characterized as two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed algorithm can be easily implemented in any irregular non-uniform grid systems that are usual in finite element method (FEM). Moreover, the proposed algorithm can be extended and applied to the 3-D free surface flow problem without additional efforts. For computation of unsteady incompressible flow, a finite element approximation based on the explicit fractional step method has been adopted. In addition, the SUPG(streamline upwind/Petrov-Galerkin) method has been implemented to deal with convection dominated flows. Combination of the proposed free surface tracking scheme and explicit fractional step formulation resulted in an efficient solution algorithm. Validity of the present solution algorithm was demonstrated from its application to the broken dam and the solitary wave propagation problems.