• Selected Papers of The Society of Naval Architects of Korea
• /
• 제2권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
• /
• 제78권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.

• 대한토목학회논문집
• /
• 제28권6B호
• /
• pp.691-696
• /
• 2008

본 연구에서는 자유수면 흐름에 적용할 수 있는 연직방향에 대해 좌표변환된 3차원 동수압 모형을 제시하였다. 제시한 모형은 자유수면과 동수압의 해석을 위하여, 2중 예측-수정(double predictor-corrector)방법을 적용하였다. 본 연구에서는 정확한 동역학적 경계조건(자유수면에서의 압력은 0인 조건)을 적용하는 방법을 검토하였고, 이 경계조건은 기존에 개발된 모형에 미소한 수정을 통하여 적용 가능함을 보여주었다. 본 연구에서 제시한 모형과 기존 모형의 계산결과를 비교하였을 때 동역학적 경계조건의 정확한 적용이 매우 중요함을 알 수 있다.

• Ocean Systems Engineering
• /
• 제7권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.

• 한국방재학회 논문집
• /
• 제10권1호
• /
• pp.103-109
• /
• 2010

본 연구에서는 자유수면 흐름에 적용할 수 있는 연직방향에 대해 좌표변환된 3차원 동수압 모형을 제시하였다. 제시한 모형은 자유수면과 동수압의 해석을 위하여, 2중 예측-수정(double predictor-corrector)방법을 적용하였다. 본 연구에서는 정확한 동역학적 경계조건(자유수면에서의 압력은 0인 조건)을 적용하는 방법을 검토하였고, 제시한 모형을 이용한 수치모의 결과를 해석해와 비교하여 본 연구에서 제시한 모형의 우수성을 검증하였다.

• 한국해양공학회지
• /
• 제26권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.

• 대한조선학회지
• /
• 제18권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.

• 대한조선학회지
• /
• 제19권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
• /
• 제2권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년도 춘계 학술대회논문집
• /
• 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.