• International Journal of Naval Architecture and Ocean Engineering
• /
• 제2권3호
• /
• pp.127-131
• /
• 2010

A two-layer fluid with free surface is simulated in the time domain by a two-dimensional potential-based Numerical Wave Tank (NWT). The developed NWT is based on the boundary element method and a leap-frog time integration scheme. A whole domain scheme including interaction terms between two layers is applied to solve the boundary integral equation. The time histories of surface elevations on both fluid layers in the respective wave modes are verified with analytic results. The amplitude ratios of upper to lower elevation for various density ratios and water depths are also compared.

• 한국해양공학회지
• /
• 제20권3호
• /
• pp.1-6
• /
• 2006

Modulation theory describes propagation of surface waves with deep wave number and frequency modulation. Locally spectrally narrow wave packet can have accumulated large scale frequency shift of carrier wave during propagation. Some important nonlinear modulation effects, such as negative frequencies, phase kinks, crest pairing, etc., often observed experimentally at long fetch propagation of finite amplitude surface wave trains, are reproduced by the proposed theory. The presented model permits also to analyze the appropriately short surface wave packets and modulation periods. Solutions show the wave phase kinks to arise on areas' of relatively small free surface displacement in complete accordance with the experiments.

• 대한조선학회논문집
• /
• 제51권4호
• /
• pp.265-273
• /
• 2014

Wave-body interaction is simulated using a developed code based on the flux-difference splitting scheme for immiscible and incompressible fluids and the hybrid Cartesian/immersed boundary method. A free surface is captured as a moving contact discontinuity within a fluid domain and an approximated Riemann solver is used to estimate the inviscid flux across the discontinuity. Immersed boundary nodes are identified inside an instantaneous fluid domain near a moving body, then dependent variables are reconstructed at those immersed boundary nodes based on interpolation along local normal lines to the boundary. Free surface flows around an oscillating cylinder are simulated and the computed wave elevations are compared with other reported results. The generation of a solitary wave by a moving wave-maker is simulated and the time histories of wave elevations at two different points are compared with other results. The developed code is applied to simulate body motion of an elastically mounted circular cylinder as a solitary wave passes the body. The force acting on an elastically mounted cylinder is compared with the force acting on a fixed cylinder. Grid independency of the computed body motion is established based on a comparison of results using three different-size grids.

• Ocean Systems Engineering
• /
• 제2권2호
• /
• pp.147-158
• /
• 2012

A time-domain simulation of a land-based Oscillating Water Column (OWC) with various irregular waves as a form of PM spectrum is performed by using a two-dimensional fully nonlinear numerical wave tank (NWT) based on the potential theory, mixed Eulerian-Lagrangian (MEL) approach, and boundary element method. The nonlinear free-surface condition inside the OWC chamber was specially devised to describe both the pneumatic effect of the time-varying pressure and the viscous energy loss due to water column motions. The quadratic models for pneumatic pressure and viscous loss are applied to the air and free surface inside the chamber, and their numerical results are compared with those with equivalent linear ones. Various wave spectra are applied to the OWC system to predict the efficiency of wave-energy take-off for various wave conditions. The cases of regular and irregular waves are also compared.

• 대한조선학회논문집
• /
• 제46권3호
• /
• pp.219-231
• /
• 2009

Numerical simulations have been carried out for a 138K LNG Carrier (KLNG) model ship with free surface, using WAVIS 2.0 and Fluent 6.3.26 with various $y^+$ values and different grid densities. Level-set method for free surface capturing was adopted in WAVIS, while VOF has been used in Fluent. The calculated results were compared with the experiment data. Resistance coefficient, wave pattern, wave profile along the hull surface, axial velocity contours and transverse vectors had been analyzed. When the first $y^+$ value was fixed at 60, the simulation results from both WAVIS and Fluent were improved as the number of grids increased. The convergence time of WAVIS was much shorter than that of Fluent. Furthermore, WAVIS predicted the velocity field and the wave profile along the hull surface better than Fluent. However, Fluent gave better wave patterns.

• 한국해양공학회지
• /
• 제24권1호
• /
• pp.1-9
• /
• 2010

A series of experiments employing particle image velocimetry (PIV) technique was conducted to produce benchmark wave kinematics data for regular waves having four different wave slopes in 2-D wave tank. Water velocities and accelerations near the free surface of regular waves were computed from image pair obtained by PIV systems. With the measured wave velocity field, the wave accelerations were computed using a centered finite difference scheme. Both local and convective components of the total accelerations are obtained from experimental data. With increasing the wave slope, the horizontal velocity and the vertical accelerations near the wave crest obtained by PIV technique became larger than theoretical results, which are well-known phenomena of the wave nonlinearity. It is noted that the relative magnitude of convective acceleration to the local acceleration became larger with increasing wave slope.

• 한국전산유체공학회지
• /
• 제2권2호
• /
• pp.73-79
• /
• 1997

VOF 방법에 의한 자유표면 유동계산의 정확성을 개선하기 위해 이차정도 모델을 개발하였다. 개발된 이차원 모델의 정확성을 비교하기 위하여 여러 가지 크기의 원형 및 Solitary wave형상의 자유표면 유동을 통하여 기존에 개발된 두 가지의 일차정도 모델과 비교하였다. 비교결과 반경이 큰 원과 같이 곡률이 작은 형상의 경우에는 일차정도 모델도 비교적 정확한 결과를 보여주고 있으나 작은 반경의 원형이나 Solitary wave와 같이 곡률이 큰 형상의 경우 일차정도 모델은 많은 오차를 보여주는 반면에 이차정도 모델은 어느 경우에나 매우 정확한 결과를 보여준다.

• 대한조선학회논문집
• /
• 제47권3호
• /
• pp.291-305
• /
• 2010

A three-dimensional time-domain calculation method is of crucial importance in prediction of the motions and wave loads of a ship advancing in a severe irregular sea. The exact solution of the free surface wave-ship interaction problem is very complicated because of the essentially nonlinear boundary conditions. In this paper, an approximate body nonlinear approach based on the three-dimensional time-domain forward-speed free-surface Green function has been presented. The Froude-Krylov force and the hydrostatic restoring force are calculated over the instantaneous wetted surface of the ship while the forces due to the radiation and scattering potentials over the mean wetted surface. The time-domain radiation and scattering potentials have been obtained from a time invariant kernel of integral equations for the potentials which are discretized according to the second-order boundary element method (Hong and Hong 2008). The diffraction impulse-response functions of the Wigley seakeeping model advancing in transient head waves at various Froude numbers have been presented. A simulation of coupled heave-pitch motion of a long rectangular barge advancing in regular head waves of large amplitude has been carried out. Comparisons between the linear and the approximate body nonlinear numerical results of motions and wave loads of the barge at a nonzero Froude number have been made.

• 대한조선학회논문집
• /
• 제43권6호
• /
• pp.638-646
• /
• 2006

In this paper a new hull-panel generation algorithm named 'Panel Cutting Method' was developed to solve the flow phenomena around a ship advancing on the free surface with a constant speed. In this algorithm the non-linearity of the free surface boundary conditions was taken into account using the iterative method and the raised panel was used at each iteration step. Numerical calculations were performed to investigate the validity of the developed algorithm using the series $60(C_B=0.60)$ hull The wave resistance coefficients, the wave patterns and the wave heights were compared between the computed and the experimental results at Fn=0.25 and 0.316. The comparison showed good agreement between computation and experiment.

• Ocean Systems Engineering
• /
• 제4권2호
• /
• pp.83-97
• /
• 2014

Wave propagation in a three-dimensional (3D) fully nonlinear numerical wave tank (NWT) is studied based on velocity potential theory. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing algorithm, B-spline, is applied to eliminate the possible saw-tooth instabilities. The artificial wave speed employed in MTF (multi-transmitting formula) approach is investigated for fully nonlinear wave problem. The numerical results from incorporating the damping zone (DZ), MTF and MTF coupled DZ (MTF+DZ) methods as radiation condition are compared with analytical solution. An effective MTF+DZ method is finally adopted to simulate the 3D linear wave, second-order wave and irregular wave propagation. It is shown that the MTF+DZ method can be used for simulating fully nonlinear wave propagation very efficiently.