• Journal of the Korean Institute of Navigation
• /
• v.17 no.3
• /
• pp.85-92
• /
• 1993

In connection to the design of high speed vessels, the numerical simulation is carried out to make clear the property of flows and breaking phenomena around the catamaran. It is because the bradking phenome-non is closely related to the free-surface turbulent flow. The free-surface wave and transverse velocity vectors are calculated around the twin and demi hull of the catamaran. Computed results are applied to detect the appearance of sub-breaking waves around the hull. The critical condition for their appearance is studied at two Froude numbers of 0.45 and 0.95. The nu-merical analysis shows that the breaking is more serious near the twin hull rather the demi hull. To simu-late the flows, the Navier-Stokes solver is invoked with a free-surface. The computation is made only in half a domain because it is symmetric in the shape.

• Journal of Mechanical Science and Technology
• /
• v.14 no.7
• /
• pp.789-795
• /
• 2000

Free-surface flows with an arbitrary deformation, induced by a submerged hydrofoil, are simulated numerically, considering two-fluid flows of both water and air. The computation is performed by a finite volume method using unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell wise local mesh refinement. The integration in space is of second order, based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels. The linear equations are solved by conjugate gradient type solvers, and the non-linearity of equations is accounted for through Picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations, the continuity equation, the conservation equation of one species, and the equations for two turbulence quantities. Finally, a comparison is quantitatively made at the same speed between the computation and experiment in which the grid sensitivity is numerically checked.

• 한국방재학회:학술대회논문집
• /
• 2007.02a
• /
• pp.181-185
• /
• 2007

A free-surface correction(FSC) method is presented to solve the 3-D shallow water equations. Using the mode splitting process, FSC method can simulate shallow water flows under the hydrostatic assumption. For the hydrostatic pressure calculation, the momentum equations are firstly discretized using a semi-implicit scheme over the vertical direction leading to the tri-diagonal matrix systems. A semi-implicit scheme has been adopted to reduce the numerical instability caused by relatively small vertical length scale compare to horizontal one. and, as the free surface correction step the final horizontal velocity fields are corrected after the final surface elevations are obtained. Finally, the vertical final velocity fields can be calculated from the continuity equation. The numerical model is applied to the calculation of the simulation of flow fields in a rectangular open channel with the tidal influence. The comparisons with the analytical solutions show overall good agreements between the numerical results and analytical solutions.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.1
• /
• pp.104-114
• /
• 1993

In this paper, effects of free surface on viscous flow is investigated. Continuity equation coupled with Navier-Stokes equations are solved numerically by using an artificial compressibility method[1, 2]. The body-fitted generalized curvilinear coordinate system is employed to deal with arbitrary body shape. The IAF scheme with finite difference method is used to solve the equations, and a diagonal algorithm is applied to time-varying Jacobian matrices for the computational economics. Free surface shape is obtained by applying zero pressure condition to still water surface at each time step. A numerical test is made for larminar flow around a circular cylinder vertically piercing the free surface. Computed flow patterns are largely affected by the existance of free surface in low Reynolds number flows treated in this paper. Free surface causes viscous pressure drag to vary much in depth direction in accordance with the variations of flow pattern.

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

A volume capturing method using unstructured grid system for numerical analysis of multiphase flows is introduced in the present paper. This method uses an interface capturing method (CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The novelty of CICSAM lies in the adaptive combination of high resolution discretization scheme which ensures the preservation of the sharpness and shape of the interface while retaining boundedness of the field, and no explicit interface reconstruction which is perceived to be difficult to implement on unstructured grid system. Several typical test cases for multiphase flows are presented, which are simulated by an in-house solution code(PowerCFD). This code employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with CICSAM. It is found that the present method simulates efficiently and accurately complex free surface flows such as multiphase flows.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.8
• /
• pp.1142-1149
• /
• 2003

The present study describes a numerical analysis for simulation of the sloshing of flows with free-surface which contained in a rectangular tank moving in harmonic or pitching motion. The VOF function, representing the volume fraction of a cell occupied by the fluid, is calculated for each cells, which gives the location of the free-surface filling any some fraction of cells with fluid. The time-dependent changes of free-surface height are used for visualization subject to several conditions such as fluid height, horizontal acceleration, sinusoidal motion, and viscosity. The free-surface heights were used for comparing wall-force, which is caused by sloshing of flows. Damping effects by baffles were extensively investigated for various conditions in terms of baffle shape and position.

• Journal of computational fluids engineering
• /
• v.20 no.1
• /
• pp.57-64
• /
• 2015

This paper describes a volume of fluid(VOF) method, mRHRIC for the simulation of free surface flows around ship hulls and provides its validation against benchmark test cases. The VOF method is developed on the basis of RHRIC method developed by Park et al. that uses high resolution differencing schemes to algebraically preserve both the sharpness of interface and the boundedness of volume fraction. A finite volume method is used to solve the governing equations, while the realizable ${\kappa}-{\varepsilon}$ model is used for turbulence closure. The present numerical results of the resistance performance tests for DTMB5415 and KCS hull forms show a good agreement with available experimental data and those of other free surface methods.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.27-37
• /
• 1991

In this paper the nonlinear free surface flows for an axisymmetric submerged body oscillating beneath the free surface are solved and the forces acting on the body are calculated. A boundary integral method is applied to solve the axisymmetric boundary value problem and the Runge-Kutta 4-th order method is used for the time stepping of the free surface location. The nonlinear forces acting on the axisymmetric body are computed and compared with published results.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2558-2561
• /
• 2008

In the present study, a parallel Taylor-Galerkin/level set based two-phase flow code was developed using finite element discretization and domain decomposition method based on MPI (Message Passing Interface). The proposed method can be utilized for the analysis of a large scale free surface problem in a complex geometry due to the feature of FEM and domain decomposition method. Four-step fractional step method was used for the solution of the incompressible Navier-Stokes equations and Taylor-Galerkin method was adopted for the discretization of hyperbolic type redistancing and advection equations. A Parallel ILU(0) type preconditioner was chosen to accelerate the convergence of a conjugate gradient type iterative solvers. From the present parallel numerical experiments, it has been shown that the proposed method is applicable to the simulation of large scale free surface flows.

• Journal of Ship and Ocean Technology
• /
• v.8 no.3
• /
• pp.8-17
• /
• 2004

A numerical method is developed for computing the free surface flows around a transom stern of a ship at a high Froude number. At high speed, the flow may be detached from the flat transom stern. In the limit of the high Froude number, the problem becomes a planning problem. In the present study, we make the finite-element computations for a transom stern flows around a wedge-shaped floating ship. The numerical method is based on the Hamilton's principle. The problem is formulated as an initial value problem with nonlinear free surface conditions. In the numerical procedures, the domain was discretized into a set of finite elements and the numerical quadrature was used for the functional equation. The time integrations of the nonlinear free surface condition are made iteratively at each time step. A set of large algebraic equations is solved by GMRES(Generalized Minimal RESidual, Saad and Schultz 1986) method which is proven very efficient. The computed results are compared with previous numerical results obtained by others.