• Selected Papers of The Society of Naval Architects of Korea
• /
• v.2 no.1
• /
• pp.63-78
• /
• 1994

This paper deals with the treatment of the open boundary in two-dimensional free-surface wave problems. Two numerical schemes are investigated for the implementation of the open boundary condition. One is to add the artificial damping term to the dynamic free-surface boundary condition, in which the determination of suitable damping coefficient and the damping zone is the most important. The other is a modified Orlanski's method, which is known to be very useful for the uni-directional waves. Using these two schemes, numerical tests have been conducted for a few typical free-surface wave problems. To obtain the numerical solution of the free-surface boundary value problem, the fundamental source-distribution method is used and the fully nonlinear free-surface boundary conditions are applied. The computed results are presented in comparison with those of others for the proof of practicality of these two schemes.

• Structural Engineering and Mechanics
• /
• v.78 no.4
• /
• pp.497-518
• /
• 2021

In this study, a new method for treating the wall boundary in smoothed particle hydrodynamics (SPH) is proposed to simulate free surface flows effectively. Unlike conventional methods of wall boundary treatment through boundary particles, in the proposed method, the wall boundary condition is directly imposed by adding boundary truncation terms to the mass and momentum conservation equations. Thus, boundary particles are not used in boundary modeling. Doing so, the wall boundary condition is accurately imposed, boundary modeling is simplified, and computation is made efficient without losing stability in SPH. Performance of the proposed method is demonstrated through several numerical examples: dam break, dam break with a wedge, sloshing, inclined bed, cross-lever rotation, pulsating tank and sloshing with a flexible baffle. These results are compared with available experimental results, analytical solutions, and results obtained using the boundary particle method.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6B
• /
• pp.691-696
• /
• 2008

In this paper, a three-dimensional non-hydrostatic pressure model for free surface flows using a normalized vertical coordinate system is presented. To strongly couple the free surface and non-hydrostatic pressure in the momentum equations, a double predictor-corrector method is employed. This research is especially focused on implementing the dynamic boundary condition (a zero pressure condition) at the free surface. This boundary condition can be specified accurately with a small modification to existing models. Numerical results with and without this modification clearly show that a precise implementation of the dynamic boundary condition is paramountly important.

• Ocean Systems Engineering
• /
• v.7 no.4
• /
• pp.371-398
• /
• 2017

Calm water wave resistance plays a very important role in ship hull design. Numerical methods are meaningful for this reason. In this study, two prevailing methods, the Neumann-Kelvin and the Rankine source method, were implemented and compared. The Neumann-Kelvin method assumes linearized free surface boundary condition and only needs to mesh the hull surface. The Rankine source method considers nonlinear free surface boundary condition and meshes both the ship hull surface and free surface. Both methods were implemented and the wave resistance of a Wigley III and three Series 60(Cb=0.6, 0.7, 0.8) hulls were analyzed. The results were compared with experimental results and the merits of both numerical techniques were quantified. Based on the results, it is concluded that the Rankine source method is more accurate in the calculation of the wave-making resistance. Using the Neumann-Kelvin method, it is found to be easier to model the hull and can be used for slender ships to solve problems like wave current coupling calculation.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.10 no.1
• /
• pp.103-109
• /
• 2010

In this study, a three-dimensional non-hydrostatic pressure model based on a normalized vertical coordinate system for free surface flows is presented. To strongly couple the free surface and non-hydrostatic pressure with the momentum equations, a double predictor-corrector method is employed. The study is especially focused on implementing the dynamic boundary condition (a zero pressure condition) at the free surface with ignoring of the atmospheric pressure. It is shown that the boundary condition can be specified easily with a slight modification to existing models.

• Journal of Ocean Engineering and Technology
• /
• v.26 no.5
• /
• pp.55-62
• /
• 2012

In this paper, the rigid lid boundary condition is applied to simulate the influence of floating structures such as ships or pontoons, and the pressure term in both the momentum equations and continuity equation are modified. The pressure of a floating structure under the free surface is dependent on the draft of the structure, generally called a ship. If the free surface is covered by a floating structure, the free surface cannot move freely. The water level should be fixed, using a rigid lid boundary condition. This boundary condition is implemented by reducing the storage area of the grid cell with a factor between zero and one. The numerical model developed by Hong (2009) is verified through a comparison with experimental results, and the influence of the reduction factor is investigated using the verified numerical model.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.18 no.4
• /
• pp.1-11
• /
• 1981

A flow problem around a ship hull with the nonlinear free surface boundary condition has been considered within the potential flow assumption. The Green's function based on the hull boundary condition is constructed numerically and used to satisfy the free surface boundary condition. This singularity to be distributed ideally on the undulating free surface is put actually, for practical reasons, on the flat water surface with the assumption of linear variation of velocities between the two positions. The surfaces of singularity distribution are approximated by Hess and Smith type quadrilaterals. The radiation condition is only crudely satisfied and this produced one of the major difficulties arising in the present way of attacking the problem.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.19 no.2
• /
• pp.13-18
• /
• 1982

The method of Rankine source distribution is emerging as a powerful yet simple alternative for the solution of complicated free surface problems. But it has been uncertain whether the radiation condition could be satisfied exactly by distributing the simple sources on the free surface only. In this paper, it is proved rigorously that the Rankine sources, whose intensities are varying sinusoidally along the axis satisfying the free surface boundary condition, generate the radiation waves both in the infinite and finite-depth flows. A formula is derived to give the upper and lower bounds of the errors in the induced velocity computation that will be introduced by truncating the extent of source distribution on the free surface. Since the truncation is inevitable in the numerical analysis, this formula may be used as a criterion to limit the position of the field points, where velocity computation is made, away from the truncation boundary. A typical analysis shows that the maximum error will be 3.4 percent of the exact induced velocity when the field point is on the free surface two wave lengths away from the truncation boundary.

• Journal of Advanced Research in Ocean Engineering
• /
• v.2 no.1
• /
• pp.1-11
• /
• 2016

The objective of this paper is to present an analytical solution in deep water waves and verify the validity of the theory (Shin, 2015). Hence this is a follow-up to Shin (2015). Instead of a variational approach, another approach was considered for a more accurate assessment in this study. The products of two coefficients were not neglected in this study. The two wave profiles from the KFSBC and DFSBC were evaluated at N discrete points on the free-surface, and the combination coefficients were determined for when the two curves pass the discrete points. Thus, the solution satisfies the differential equation (DE), bottom boundary condition (BBC), and the kinematic free surface boundary condition (KFSBC) exactly. The error in the dynamic free surface boundary condition (DFSBC) is less than 0.003%. The wave theory was simplified based on the assumption tanh $D{\approx}1$ in this paper. Unlike the perturbation method, the results are possible for steep waves and can be calculated without iteration. The result is very simple compared to the 5th Stokes' theory. Stokes' breaking-wave criterion has been checked in this study.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.145-155
• /
• 1998

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.