• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.19-26
• /
• 1999

This paper describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows with free-surface. The Navier-Stokes equations governing the flows have been discretized by means of finite-difference approximations, and the resulting equations have been solved via the SIMPLE-C algorithm. The free-surface is defined by the motion of a set of marker particles and the interface behaviour was investigated by means of a "Lagrangian" technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.6 s.144
• /
• pp.583-592
• /
• 2005

The turbulent free surface flow around a self-propelled KRISO 138K LNG Carrier is numerically simulated using the finite volume based multi-block RANS code, WAVIS developed at HRISO. The realizable k-$\varepsilon$ turbulence model with a wail function is employed for the turbulence closure. The free surface is captured with the Level-Set method and body forces are used to model the effects of a propeller without resolving the detail blade flow. In order to obtain an accurate free surface solution and stable convergence, the computations are executed with a proper fine grid refinement around the free surface and with an adoption of implicit discretization scheme for the Level-Set formulation. The computed velocity vectors at the several stations and wave patterns show a good agreement with the experimental results measured at the KRISO towing tank.

• International Journal of Ocean System Engineering
• /
• v.2 no.1
• /
• pp.32-36
• /
• 2012

This paper presents a comparison of potential and viscous computational codes for the water entry problem. A po-tential code was developed which adopted the boundary element method to solve the problem. A nonlinear free surface boundary condition was integrated to find new locations of free surface. The dynamic boundary condition was simplified by taking constant potential values for every time steps. The simplified dynamic boundary condition was applied in the new position of the free surface not at the mean level, which is the usual practice for linearized theory. The commercial code FLUENT was used to solve the water entry problem from the viscosity point of view. The movement of the air-liquid interface is traced by distribution of the volume fraction of water in a computational cell. The pressure coefficients were compared with each other, while experimental results published by other researchers were also examined. The characteristics of each method were discussed to clarify merits and limitations when they were applied to the water entry problems.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.13 no.1
• /
• pp.260-277
• /
• 2021

The appropriate radiation conditions of ship motion problem with advancing speed in frequency domain are investigated from a theoretical and practical point of view. From extensive numerical experiments that have been conducted for evaluation of the relevant radiation conditions, a hybrid radiation technique is proposed in which the Sommerfeld radiation condition and the free surface damping are mixed. Based on the comparison with the results of the translating and pulsating Green function method, the optimal damping factor of the hybrid radiation technique is selected, and the observed limitations of the proposed hybrid radiation technique are discussed, along with its accuracy obtained from the numerical solutions. Comparative studies of the forward-speed seakeeping prediction methods available confirm that the results of applying the hybrid radiation technique are relatively similar to those obtained from the translating and pulsating Green function method. This confirmation is made in comparisons with the results of solely applying either the free surface damping, or the Sommerfeld radiation condition. By applying the proposed hybrid radiation technique, the wave patterns, hydrodynamic coefficients, and motion responses of the Wigley III hull are finally calculated, and compared with those of model tests. It is found that, in comparison with the model test results, the three-dimensional Rankine source method adopting the proposed hybrid radiation technique is more robust in terms of accuracy and numerical stability, as well as in obtaining the forward speed seakeeping solution.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.4
• /
• pp.440-451
• /
• 1985

This study is concerned with the analysis of cutting force system acting on ball-nose end mill in three-dimensional surface machining process. Conical tipped circular cutting edge element model and free surface machining process types are proposed to apply oblique cutting theory, and then derived equations are used for numerical approach of cutting force curves by matrix method. This approach has a good agreement with experimental results both in magnitude and shape within the range of 15 percent, which was conformed on 6061-T6 aluminum workpiece having twofold curvatured surface. From the cutting load variation to edge location, it is confirmed that circular cutting edge shapes has a better cutting ability than that of straight and both have a singularity near a tool point. It is also verified that what kind of machining condition is recommendable for three-dimensional machining process in connection with deflection of the cutter to workpiece and tool point wearing or system stability.

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

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.54-57
• /
• 2011

In this paper we present an algorithm about how to simulate two dimensional dam breaking with lattice Boltzmann method (LBM). LBM considers a typical volume element of fluid to be composed of a collection of particles that represented by a particle velocity distribution function for each fluid component at each grid point. We use the modified Lattice Boltzmann Method for incompressible fluid. This paper will represent detailed information on single phase flow which considers only the water instead of both air and water. Interface treatment and conservation of mass are the most important things in simulating free surface where the Interface is treated by mass exchange with the water region. We consider the surface tension on the interface and also bounce back boundary condition for the treatment of solid obstacles. We will compare the results of the simulation with some methods and experimental results.

• Journal of Ocean Engineering and Technology
• /
• v.4 no.1
• /
• pp.62-70
• /
• 1990

Cauchy의 적분공식을 복소속도(complex velocity)에 적용하여 포텐시얼 유동을 해석하는 복소경계요소법이 개발되었다. 이 결과로 얻어지는 적분방정식은 경계면에서의 접선속도(tangential velocity)와 법선속도(normal velocity)의 함수로 주어진다. 자유수면에서의 접선속도의 시간변화(evolution of tangential velocity)를 수식화하기 위하여 새로운 비선형 동역학적 자유수면경계조건(nonlinear dynamic free surface boundary condition)을 유도하였다. 복소포텐시얼 대신 복소속도를 이용하는 이 방법은 유장내의 특이점(field singularity)을 용이하게 고려할 수 있으며, 수치미분없이 직접 경계면에서의 유속을 해로서 구하게 된다. 그러나 자유수면이 존재하는 문제의 경우에는, 자유수면에서의 동역학적 경계조건을 만족 시키기 위한 계산과정에 접선 벡타의 변화량을 추정하는 것이 포함되게 되어, 계산과정이 다소 복잡하게 된다.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.90-100
• /
• 2006

A digital wave tank (DWT) simulation technique has been developed by authors to investigate the interactions of fully nonlinear waves with 3D marine structures. A finite-difference/volume method and a modified marker-and-cell (MAC) algorithm have been used, which are based on the Navier-Stokes (NS) and continuity equations. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique or the Level-Set (LS) technique developed for one or two fluid layers. In this paper, some applications for various engineering problems with free-surface are introduced and discussed. It includes numerical simulation of marine environments by simulation equipments, fully nonlinear wave motions around offshore structures, nonlinear ship waves, ship motions in waves and marine flow simulation with free-surface. From the presented simulations, it seems that the developed DWT simulation technique can handle various engineering problems with free-surface and reliably predict hydrodynamic features due to the fully-nonlinear wave motions interacting with such marine structures.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.1
• /
• pp.15-23
• /
• 1998

Characteristics of turbulence beneath the free surface around a blunt bow are numerically-investigated. Three dimensional wavier-Stokes and continuity equations are solved for the simulations.. The Large Eddy Simulation(LES) with the external disturbance is performed to simulate the turbulent free surface flow called sub-breaking wave. The result shows that the free surface fluctuates beyond a certain critical condition and the characteristics of the fluctuation are similar to the turbulent boundary layer flow around a solid body.