• Journal of computational fluids engineering
• /
• v.17 no.2
• /
• pp.35-41
• /
• 2012

This paper presents incompressible Navier-Stokes solution algorithm for 2D Free-surface flow problems on the Cartesian mesh, which was implemented to run on Graphics Processing Units(GPU). The INS solver utilizes the variable arrangement on the Cartesian mesh, Finite Volume discretization along Constrained Interpolation Profile-Conservative Semi-Lagrangian(CIP-CSL). Solution procedure of incompressible Navier-Stokes equations for free-surface flow takes considerable amount of computation time and memory space even in modern multi-core computing architecture based on Central Processing Units(CPUs). By the recent development of computer architecture technology, Graphics Processing Unit(GPU)'s scientific computing performance outperforms that of CPU's. This paper focus on the utilization of GPU's high performance computing capability, and presents an efficient solution algorithm for free surface flow simulation. The performance of the GPU implementations with double precision accuracy is compared to that of the CPU code using an representative free-surface flow problem, namely. dam-break problem.

• Journal of computational fluids engineering
• /
• v.19 no.4
• /
• pp.20-28
• /
• 2014

An efficient solution algorithm for simulating free surface problem is presented. Navier-Stokes equations for variable density incompressible flow are employed as the governing equation on Cartesian meshes. In order to describe the free surface motion efficiently, VOF(Volume Of Fluid) method utilizing THINC(Tangent of Hyperbola for Interface Capturing) scheme is employed. The most time-consuming part of the current free surface flow simulations is the solution step of the linear system, derived by the pressure Poisson equation. To solve a pressure Poisson equation efficiently, the PCG(Preconditioned Conjugate Gradient) method is utilized. This study showed that the proper application of the preconditioner is the key for the efficient solution of the free surface flow when its pressure Poisson equation is solved by the CG method. To demonstrate the efficiency of the current approach, we compared the convergence histories of different algorithms for solving the pressure Poisson equation.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.3
• /
• pp.933-944
• /
• 1996

In this paper, two new techniques, the pattern filling and the refined flow field regeneration, based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible viscous flow with free surfaces. The gorerning equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The penalty and Newton-Raphson methods are used effectively for finite element formulation. The flow front surface and the volume inflow rate are calculated using the pattern filling technique to select an adequate pattern among five filling patterns at each quadrilateral control volume. By the refined flow field regeneration technique, the new flow field which renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. Using the new thchniques to be developed, the dam-breaking problem has been analyzed to predict flow phenomenon of fluid and the predicted front positions versus time have been compared with the reported experimental result.

• Journal of Korea Water Resources Association
• /
• v.44 no.6
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• 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.

• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.228-230
• /
• 2009

The moment-of-fluid (MOF) method is a new volume-tracking method that accurately treats evolving material interfaces. Based on the moment data (volume and centroid) for each material, the material interfaces are reconstructed with second-order spatial accuracy in a strictly conservative manner. The MOF method is coupled with a stabilized finite element incompressible Navier-Stokes solver for two fluids, namely water and air. The effectiveness of the MOF method is demonstrated with a free-surface dam-break problem.

• Journal of Ocean Engineering and Technology
• /
• v.14 no.3
• /
• pp.35-40
• /
• 2000

SMAC method is, one of the numerical simulation techniques, modified from the original MAC for the time-dependent variation of fluid flows. The Navier-Stokes equation for incompressible time-dependent viscous flow is applied and, also marker particles which move with the fluid are used. Two-dimensional numerical computations of fluid flow are carried out in a virtual water tank. This paper has shown very well the movements of marker particles using SMAC method.

• Ocean Systems Engineering
• /
• v.1 no.3
• /
• pp.207-222
• /
• 2011

A stabilized incompressible smoothed particles hydrodynamics (ISPH) method is utilized to simulate free falling rigid body into water domain. Both of rigid body and fluid domain are modeled by SPH formulation. The proposed source term in the pressure Poisson equation contains two terms; divergence of velocity and density invariance. The density invariance term is multiplied by a relaxed parameter for stabilization. In addition, large eddy simulation with Smagorinsky model has been introduced to include the eddy viscosity effect. The improved method is applied to simulate both of free falling vessels with different materials and water entry-exit of horizontal circular cylinder. The applicability and efficiency of improved method is tested by the comparisons with reference experimental results.