• 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 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.31 no.2
• /
• pp.45-56
• /
• 1994

A computationally efficient numerical method based on potential flow is developed for time-domain simulation of the nonlinear free-surface flows around a 2-dimensional hydrofoil. This numerical method, namely, spectral/boundary-element method, is a mixed one of the high-order spectral method and the boundary-element method in time-domain. The high-order spectral method is used to calculate the nonlinear evolution of free-surface, and the boundary-element method is used to calculate the effects of the hydrofoil and the shed vortex. As application examples, nonlinear free-surface flows around a 2-dimensional hydrofoil which starts from the rest and translates near the free-surface with or without harmonic oscillations are calculated. Nonlinear/unsteady results of free-surface waves and hydrodynamic farces are shown and discussed. Particularly, the results of steady-state which are obtained as a special case of the present unsteady solution are compared with others' calculated and experimental results, and good agreements are observed.

• 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 Korean Society of Fisheries and Ocean Technology
• /
• v.39 no.1
• /
• pp.27-32
• /
• 2003

Under the assumption of potential flow, free-surface flow around a hydrofoil is calculated by the high-order spectra1!boundary-integral method, This method is one of the most efficient numerical methods by which the nonlinear interactions between hydrofoil and free-surface can be simulated in time-domain. In this method. the wave potential which represents the nonlinear evolution of free-surface is solved by the high-order spectral method and the body potential which provides the effects of hydrofoil and shed vortex is solved by the boundary-integral method. The calculated free-surface profiles which are generated by a uniformly translating hydrofoil are compared with other experimental results. And they show relatively good agreements each other. As another example, free-surface flow generated by a heaving and translating hydrofoil is calculated and discussed.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.27-37
• /
• 1991

In this paper the nonlinear free surface flows for an axisymmetric submerged body oscillating beneath the free surface are solved and the forces acting on the body are calculated. A boundary integral method is applied to solve the axisymmetric boundary value problem and the Runge-Kutta 4-th order method is used for the time stepping of the free surface location. The nonlinear forces acting on the axisymmetric body are computed and compared with published results.

• 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 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.21 no.6
• /
• pp.53-58
• /
• 2007

A particle method, recognized as one of gridless methods, has been developed to investigate non-linear free-surface motions interacting with structures. This method is more feasible and effective than conventional grid methods for solving flow fieldswith complicated boundary shapes. The method consists of particle interaction models representing pressure gradient, diffusion, incompressibility, and the free-surface boundary conditions without grids. In the present study, broken dam problems with various viscosity values are simulated to validate the developed method.

• 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.