• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.24 no.4
• /
• pp.277-286
• /
• 2012

In the numerical simulation of wave fields using a multi-phase flow model that considers simultaneous flows of materials with different states such as gas, liquid and solid, there is need of an accurate representation of the interface separating the fluids. We adopted an algebraic interface capturing method called tangent of hyperbola for interface-capturing(THINC) method for the capture of the free-surface in computations of multi-phase flow simulations instead of geometrical-type methods such a volume of fluid(VOF) method. The THINC method uses a hyperbolic tangent functions to represent the surface, and compute the numerical flux for the fluid fraction functions. One of the remarkable advantages of THINC method is its easy applicability to incorporate various numerical codes based on Navier-Stokes solver because it does not require the extra geometric reconstruction needed in most of VOF-type methods. Several tests were carried out in order to investigate the advection of interfaces and to verify the applicability of the THINC method to wave fields based on the one-field model for immiscible two-phase flows (TWOPM). The numerical results revealed that the THINC method is able to track the interface between air and water separating the fluids although its algorithm is fairly simple.

• Journal of Ocean Engineering and Technology
• /
• v.27 no.2
• /
• pp.39-46
• /
• 2013

In this study, the impact loads on tank walls by sloshing phenomena and on a tall structure in a three-dimensional rectangular tank were predicted using multiphase flow simulations. The solver was based on the CIP/CCUP (Constraint interpolation CIP/CIP combined unified procedure) method, and the THINC-WLIC (Tangent hyperbola for interface capturing-weighted line interface calculation) scheme was used to capture the air-water interface. For the convection terms of the Navier-Stokes equations, the USCIP (Unsplit semi-lagrangian CIP) method was adopted. The results of simulations were compared with those of experiments. Overall, the comparisons were reasonably good.

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

In the present paper, a numerical method based on the constraint interpolation profile (CIP) method is applied for simulating two-dimensional violent sloshing problems. The free surface boundary value problem is considered as a multiphase problem which includes water and air. A stationary Cartesian grid system is adopted, and an interface capturing method is used to trace the shape of free surface profile. The CIP combined unified procedure (CCUP) scheme is applied for flow solver, and the tangent of hyperbola for interface capturing (THINC) scheme is used for interface capturing. Numerical simulations have been carried out for partially-filled 2D tanks under forced sway and roll motions at various filling depths and frequencies. The computational results are compared with experiments and/or the other numerical results to validate the present numerical method.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.1-10
• /
• 2011

A new computational program, which is based on the CIP/CCUP(Constraint Interpolation Profile/CIP Combined Unified Procedure) method, has been developed to numerically analyse sloshing phenomena dealt as multiphase-flow problems. For the convection terms of Navier-Stokes equations, the RCIP(Rational function CIP) method was adopted and the THINC-WLIC(Tangent of Hyperbola for Interface Capturing-Weighted Line Interface Calculation) method was used to capture the air/water interface. To validate the present numerical method, two-dimensional dam-breaking and sloshing problems in a rectangular tank were solved by the developed method in a stationary Cartesian grid system. In the case of sloshing problems, simulations by using a improved MPS(Moving Particle Simulation) method, which is named as PNU-MPS(Pusan National University-MPS), were also carried out. The computational results are compared with those of experiments and most of the comparisons are reasonably good.

• International Journal of Ocean System Engineering
• /
• v.1 no.3
• /
• pp.136-143
• /
• 2011

In this paper, three-dimensional free-surface flows are simulated by using two different numerical methods, the constrained interpolation profile (CIP)-based and finite volume (FV)-based methods. In the CIP-based method, the governing equations are solved on stationary staggered Cartesian grids by a finite difference method, and an immersed boundary technique is applied to deal with wave-body interactions. In the FV-based method, the governing equations are solved by applying collocated finite volume discretization, and body-fitted meshes are used. A free-surface boundary is considered as the interface of the multi-phase flow with air and water, and a volumeof-fluid (VOF) approach is applied to trace the free surface. Among many variations of the VOF-type method, the tangent of hyperbola for interface capturing (THINC) and the compressive interface capturing scheme for arbitrary meshes (CICSAM) techniques are used in the CIP-based method and FV-based method, respectively. Numerical simulations have been carried out for dam-breaking and wave-body interaction problems. The computational results of the two methods are compared with experimental data and their differences are observed.

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

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.6
• /
• pp.461-468
• /
• 2012

In this study, a Cartesian-grid method based on finite volume approach is applied to simulate the ship motions in large amplitude waves. Fractional step method is applied for pressure-velocity coupling and TVD limiter is used to interpolate the cell face value for the discretization of convective term. Water, air, and solid phases are identified by using the concept of volume-fraction function for each phase. In order to capture the interface between air and water, the tangent of hyperbola for interface capturing (THINC) scheme is used with weighed line interface calculation (WLIC) method which considers multidimensional information. The volume fraction of solid body embedded in the Cartesian grid system is calculated using a level-set based algorithm, and the body boundary condition is imposed by a volume weighted formula. Numerical simulations for the two-dimensional barge type model and Wigley hull in linear waves have been carried out to validate the newly developed code. To demonstrate the applicability for highly nonlinear wave-body interactions such as green water on the deck, numerical analysis on the large-amplitude motion of S175 containership is conducted and all computational results are compared with experimental data.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.2
• /
• pp.79-87
• /
• 2013

In this paper, an Euler equation solver based on a Cartesian-grid method and non-uniform staggered grid system is applied to predict the ship motion response and added resistance in waves. Water, air, and solid domains are identified by a volume-fraction function for each phase and in each cell. For capturing the interface between air and water, the tangent of hyperbola for interface capturing (THINC) scheme is used with a weighed line interface calculation (WLIC) method. The volume fraction of solid body embedded in a Cartesian-grid system is calculated by a level-set based algorithm, and the body boundary condition is imposed by volume weighted formula. Added resistance is calculated by direct pressure integration on the ship surface. Numerical simulations for a Wigley III hull and an S175 containership in regular waves have been carried out to validate the newly developed code, and the ship motion responses and added resistances are compared with experimental data. For S175 containership, grid convergence test has been conducted to investigate the sensitivity of grid spacing on the motion responses and added resistances.

• Journal of Ocean Engineering and Technology
• /
• v.29 no.2
• /
• pp.154-162
• /
• 2015

In this study, experimental and numerical methods were applied to observe sloshing impact phenomena. A two-dimensional rectangular tank filled with water and air was considered with a specific excitation condition that induced a hydrodynamic impact without an air pocket at the top corner of the tank. High-speed cameras and a pressure measurement system were synchronized, and a particle image velocimetry (PIV) technique was applied to measure the velocity field and corresponding pressure. The experimental condition was implemented in a numerical computation to solve incompressible two-phase flows using a Cartesian-grid method. The discretized solution was obtained using the finite difference and constraint-interpolation-profile (CIP) methods, which adopt a fractional step scheme for coupling the pressure and velocity. The tangent of the hyperbola for interface capturing (THINC) scheme was used with the weighed line interface calculation (WLIC) method to capture the interface between the air and water. The calculated impact pressures and velocity fields were compared with experimental data, and the relationship between the local velocity and pressure was investigated based on the computational results.