• Journal of Ocean Engineering and Technology
• /
• v.12 no.1
• /
• pp.85-98
• /
• 1998

본 논문의 수치해법은 경계치문제를 풀기 위하여 코시이론(Cauchy's theorem)을 사용하였다. 경계치문제는 완전한 물체표면조건과 자유표면조건을 만족시키는 초기치문제로 귀결된다. 현 수치해법에서 무한영역은 수치계산 영역인 비선형 영역과 선형 자유표면조건을 만족하는 선형영역으로 나누어진다. 선형영역의 해는 과도 그린(Green)함수를 사용하여 정합조건을 부과함으로써, 수치계산은 비선형 영역에서만 수행된다. 본 논문에서 저자는 수치계산 영역에서 코시이론을 사용하여 적분방정식을 도출하였고, 무한영역의 해는 정합면에서 과도 그린함수를 사용하여 표현하였다. 본 수치계산에서 자유표면에 요소 재분배법을 적용함으로써 쇄파현상에 대해서도 안정적인 수치해석을 할 수 있었다. 본 논문에서 개발된 수치방법을 적용한 문제는 다음과 같다. 첫째는 자유표면에서 실린더가 강제동요하는 경우에 자유표면형상과 힘을 계산하여 이전의 실험치 및 계산치와 비교하였다. 두번째로는 실린더가 자유수면하에서 일정한 속도로 항주하는 경우에는 조파저항과 양력을 계산하여 고차 스펙트럴법과 비교하였다.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.1
• /
• pp.41-49
• /
• 1997

Computational results from Navier-Stokes equations are presented for the Stokes-wave/flat-plate boundary-layer and wake for small wave steepness(Ak=0.01), including exact and approximate treatments of the viscous free-surface boundary conditions. The macro-scale flow indicate that the variations of the external-flow pressure gradients cause acceleration or deceleration of the streamwise velocity component and alternating direction of the cross flow. Remarkably, the wake displays a greater response, i.e., a bias with regard to favorable as compared to adverse pressure gradients. The micro-scale flow indicates that the free-surface boundary conditions have a profound influence over the boundary layer and near/intermediate wake. Order-of-magnitude estimates are conformed to the computational results. And appreciable errors are introduced through approximations to the free-surface boundary conditions.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.1
• /
• pp.1-8
• /
• 2001

The computations of the turbulent flow around the ship models with the free-surface effects were carried out. Incompressible Reynolds-Averaged Navier-Stokes equations were solved by using an explicit finite-difference method with the nonstaggered grid system. The method employed second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration. For the turbulence closure, a modified Baldwin-Lomax model was exploited. The location of the free surface was determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and a free-surface conforming grid was generated at each time step so that one of the grid boundary surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition was applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method developed in the present study, the computations were carried out for beth Wigley and Series 60 $C_B=0.6$ ship model and the computational results showed good agreements with the experimental data.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.54-65
• /
• 1993

In computing the wave resistance numerically, satisfying the boundary condition(BC) on the body surface is not so difficult, and then what form of the BC on the free surface(FS) be used is a crucial question. To shed some light on this, we examine the various BC's on the FS, namely, the Poisson's[1], the Ogilvie's[2] and the Dawson's[3] BC, using the same panel method for submerged bodies in two-dimension. We also show the performance of the Poisson's BC for a submerged sphere and the Wigley hull. It seems that we are still in need of a theory which gives a BC on the FS more accurate than those tested, and more practically applicable than the exact nonlinear BC.

• Journal of Korean Port Research
• /
• v.12 no.1
• /
• pp.87-93
• /
• 1998

For the calculation of the free-surface elevation, a new finite difference scheme is studied where the third derivative term for the wave elevation is artificially added in the Eulerian expression of the free-surface boundary condition. The paper presents a comparative analysis with simulations performed by the classical MAC method. More schematic computations are carried out by changing the submergence-depth and angle-of-attack. The present study shows that this new method is very efficient for the simulation of free-surface elevation around the trailing edge.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.112-122
• /
• 1993

A potential-based panel method is presented for the analysis of a partially or supercavitating two-dimensional hydrofoil at a finite submergence beneath a free surface, treating without approximation the effects of the finite Froude number and the hydrostatic pressure. Free surface sources and normal dipoles are distributed on the foil and cavity surfaces, their strength being determined by satisfying the kinematic and dynamic boundary conditions on the foil-cavity boundary. The cavity surface is determined iteratively as a part of the solution. Numerical results show that the wave profile is altered significantly due to the presence of the cavity. The buoyancy effect due to the hydrostatic pressure, which has usually been neglected in most of the cavitating flow analysis, is found playing an important role, especially for the supercavitating hydrofoil; the gravity field increases the cavity size in shallow submergence, but decreases it when deeply submerged, while the lift reduces at all submergence depth.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.1
• /
• pp.8-15
• /
• 2002

A two-dimensional higher order panel method using B-splines has been developed to overcome the disadvantages of the low order panel method and to obtain more accurate solution. The sources and the normal dipoles are distributed on both the body and the free surface. Instead of applying the upwind finite difference schemes to satisfy the linearized free surface and the radiation condition, the derivatives of the basis functions of the B-splines are directly applied to the linearized free surface condition. Numerical damping in the Dawson's method are avoided in the Present computations. In order to validate the present method, numerical computations are carried out for a submerged cylinder and a two-dimensional hydrofoil steadily moving beneath a free surface. The numerical results show that fast convergence and better accuracies have been achieved by the present method.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.3
• /
• pp.56-72
• /
• 1990

The hydrodynamic aspects of bow flare impact are discussed and the assumptions implicit in typical flare impact calculations are examined. Calculations based upon a pressure release free surface boundary condition are presented. While this simplified boundary condition eliminates the possibility of important factors such as splash-up and above-surface wetting, it has the significant advantage of being computationally simple. Both three-dimensional and two-dimensional results are compared with experiments. Errors resulting from the approximate free surface condition are examined. Also included in the two-dimensional and three-dimensional analysis are longitudinal distributions of the vertical bending moment and vertical shear force.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.1
• /
• pp.33-42
• /
• 2003

This paper deals with the FEM analysis of nonlinear sloshing of incompressible, invicid and irrotational flow in two dimensional rectangular tank. We use laplace equation based on potential theory as governing equation. For large amplitude sloshing motion, kinematic and dynamic free surface conditions derived from Bernoulli equation are applied. This problem is solved by FEM using 9-node elements. For the time integration and accurate velocity calculation, we introduce predictor-corrector time marching scheme and least square method. Also, numerical stability in tracking of free surface is obtained by direct calculation of free surface location to time variation. Numerical results of sloshing induced by harmonic excitations, while comparing with those of linear theory and references, prove the accuracy and stability. After verification of our program, we analyze sloshing response characteristics to the fluid height and the excitation amplitude.

• Journal of Ocean Engineering and Technology
• /
• v.12 no.1
• /
• pp.99-106
• /
• 1998

3차원 파수조에서 완전 비선형파를 시뮬레이션하기 위하여 우선 랜킨 소스를 기저로한 적분방정식을 고차경계요소법을 이용하여 이산화하였다. 그리고 방사경계조건은 파흡수 비치와 포텐셜 스트레칭 기법을 이용하여 모델링하였으며, 비선형 자유표면과 경계조건식은 고차 예측 및 보정 기법을 이용하여 시간 적분하였다. 파흡수 비치는 파의 진행방향 특성에 따라 수조내에 다양하게 배치할 수 있으며 비칭서 흡수가 덜된 파는 수조의 길이 방향 끝단에서 포텐셜 스트레칭 기법에 의하여 반사없이 진행하도록 하였다. 수치실험 결과 일-에너지 보존법칙과 모멘텀-임펄스 보존 법칙이 만족됨으로써 본 수치기법의 효용성이 검증되었다.