• 한국해안·해양공학회논문집
• /
• 제24권4호
• /
• pp.277-286
• /
• 2012

기체와 액체의 유동 및 고체의 변형을 동시에 고려할 수 있는 혼상류모델을 이용하여 파동장을 해석하는 경우, 서로 다른 비혼합의 유체 경계면의 시간변형을 고정도로 추적하는 것이 대단히 중요하다. 본 연구에서는 경계면의 추적에 있어서 VOF(Volume of Fluid)법으로 대표되는 경계면 형상의 재구축이 필요한 Geometrical-type의 경계추적법의 대신에 Algebraic-type의 경계추적법인 THINC(Tangent of Hyperbola for INterface-Capturing)법을 적용하였다. THINC법은 경계면에 대한 형상의 구축이 필요하지 않으므로 VOF법에 비해 비교적 간단한 알고리즘을 가지며, 기존의 Navier-Stokes solver에로 적용성이 용이한 장점을 갖는다. 본 연구에서는 THINC법의 기본적인 이류특성을 고찰하고, 혼상류수치모델인 TWOPM(one-field Model for immiscible TWO-Phase flows)과 결합한 수치모델을 파동장에 적용하여 비혼합 혼상류에서 경계면의 추적능을 검토하였다. 그 결과, 혼상류의 경계면 추적에 있어서 상대적으로 간단한 알고리즘의 THINC법이 기존의 VOF법과 유사한 정도를 갖는 해석법이라는 것을 확인할 수 있었고, 따라서 향후 기포의 연행을 동반하는 쇄파 및 쇄파력의 해석 등에 그의 적용성이 기대된다.

• 한국해양공학회지
• /
• 제27권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.

• 대한조선학회논문집
• /
• 제47권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.

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

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

• 대한조선학회논문집
• /
• 제49권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.

• 대한조선학회논문집
• /
• 제50권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.

• 한국해양공학회지
• /
• 제29권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.