• 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.30 no.3
• /
• pp.29-40
• /
• 1993

The fundamental characteristics of nonlinear free-surface waves generated by a shallowly submerged 3-dimensional hydrofoil are investigated. The fluid is assumed inviscid, incompressible and its motion irrotational. The surface tension on the free-surface is neglected. The hydrofoil is represented by a horseshoe vortex system whose shape is assumed fixed. Also the strengths of vortices are assumed given. The exact problem for the wave potential due to the horseshoe vortex system is formulated by the variational principle based on the classical Hamilton's principle. The localized finite element method is used in the numerical computations. In order to increase the numerical efficiency, an intermediate nonlinear-to-linear transition buffer subdomain for a smooth matching is introduced between the fully nonlinear computation subdomain and the truncated linear infinite subdomain. Also used is the modal analysis to reduce the computation tome drastically. The effect of inflow velocity, submergence depth of the hydrofoil and the shape of circulation distribution on the wave profiles are thoroughly examined. Especially it was possible to investigate the nonlinear influence of the free vortex on the free vortex. The nonlinear free-surface effect on the induced forces on the hydrofoil is also investigated.

• 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 the Korean Society for Marine Environment & Energy
• /
• v.7 no.4
• /
• pp.216-223
• /
• 2004

A nonlinear sloshing problem is numerically simulated. During excessive sloshing the sloshinginduced impact load can cause a critical damage on the tank structure. A three-dimensional free-surface flow in a tank is formulated in the scope of potential flow theory. The exact nonlinear free-surface condition is satisfied numerically. A finite-element method based on Hamiltons principle is employed as a numerical scheme. The problem is treated as an initial-value problem. The computations are made through an iterative method at each time step. The hydrodynamic loading on the pillar in the tank is computed.

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

Two-dimensional nonlinear free-surface wave problems are analyzed with consideration of viscosity. Navier-Stokes equation and continuity equation are solved by the application of Finite Analytic Method, and MAC scheme is used far the treatment of free surface. Surface tension effect is also considered and laminar flow is assumed. The free-surface waves in shallow water, the flows around a vortex-pair with free surface and the wave ahead of a rectangular body are simulated to test the present numerical scheme. In the shallow water problem, viscous effect due to the friction on the bottom is observed. In the second problem, the approach of a vortex-pair to the free surface is simulated to examine the interaction of vortex-pair with the free surface. In the third problem, the wave ahead of a semi-infinite floating body is simulated.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.1
• /
• pp.67-81
• /
• 2000

An accurate and efficient numerical method for two-dimensional nonlinear radiation problem has been developed. The wave motion due to a moving body is described by the assumption of ideal fluid flow, and the governing Laplace equation can be effectively solved by the higher-order boundary element method with the help of the GMRES (Generalized Minimal RESidual) algorithm. The intersection or corner problem is resolved by utilizing the so-called discontinuous elements. The implicit trapezoidal rule is used in updating solutions at new time steps by considering stability and accuracy. Traveling waves caused by the oscillating body are absorbed downstream by the damping zone technique. It is demonstrated that the present method for time marching and radiation condition works efficiently for nonlinear radiation problem. To avoid the numerical instability enhanced by the local gathering of grid points, the regriding technique is employed so that all the grids on the free surface may be distributed with an equal distance. This makes it possible to reduce time interval and improve numerical stability. Special attention is paid to the local flow around the body during time integration. The nonlinear radiation force is calculated by the "acceleration potential technique". Present results show good agreement with other numerical computations and experiments.

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

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

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.3
• /
• pp.80-89
• /
• 1992

This paper deals with the open boundary problems, and two numerical schemes are used for the implementation of open boundary condition. One is to add the artificial damping term to dynamic free-surface boundary condition. Determination of suitable damping coefficient and the damping cone is the most important in this scheme. The other scheme is a modified Orlanski's method. This will be useful for the problems with unidirectional waves. A few typical free-surface wave problems are modeled for the numerical test. Method of solution is fundamental source-distribution method and the fully nonlinear boundary conditions are applied. The computed results are compared with those of others for the proof of practicality of these schemes.

• 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.41 no.4
• /
• pp.38-44
• /
• 2004

The motion response analysis of a submerged elliptic cylinder in waves is presented and the elliptic cylinder is a simplification of the section of submarine in this paper. The method is based on boundary integral method and two-dimensional 3 degree motions are calculated in regular harmonic waves. The fully nonlinear free surface boundary condition is assumed in an numerical domain and this solution is matched along an assumed boundary as a linear solution composed of transient Green function, The large amplitude motions of an elliptic cylinder are directly simulated and effects of wave frequency, wave amplitude and the distance from buoyancy center to gravity center are discussed.