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

• 대한조선학회논문집
• /
• 제42권6호
• /
• 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
• /
• 제2권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
• /
• 제13권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.

• 대한기계학회논문집
• /
• 제9권4호
• /
• pp.440-451
• /
• 1985

본 논문에서는 곡면절삭력계의 해석을 통하여 절삭성과 내마멸성이 우수한 공 구기하형상과 곡면절삭과정에서 절삭조건에 따라 변화하는 가공정밀도 및 안정성을 규 명하고자 한다.

• 한국수자원학회논문집
• /
• 제44권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년 춘계학술대회논문집
• /
• 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.

• 한국해양공학회지
• /
• 제4권1호
• /
• pp.62-70
• /
• 1990

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

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

• 대한조선학회논문집
• /
• 제35권1호
• /
• pp.15-23
• /
• 1998

뭉툭한 스트럿에 의하여 생성되는 자유 표면의 난류 현상에 대한 기본 특성을 수치해석적 방법으로 연구하였다. 지배 방정식으로는 3차원 Navier-Stokes 방정식과 연속 방정식을 사용하였으며, 이들 지배 방정식은 유한 차분법으로 이산화 하였다. 물체앞 자유 표면에서의 난류 유동을 모사하기 위하여 임의의 작은 외부 교란을 도입하여 Large Eddy Simulation을 수행하였다. 물체앞 자유 표면은 어떠한 속도 이상에서 격렬하게 진동하는데 그 기본 특성은 물체 주위의 난류 흐름과 유사함을 수치 계산으로 보였다.