• 대한조선학회논문집
• /
• 제39권1호
• /
• pp.8-15
• /
• 2002

B-spline을 이용하여 물체의 형상과 포텐셜을 표현함으로써, 저차경계요소법의 단점들을 극복하고 수치계산의 정도를 높이기 위한 고차 패널법을 개발하였다. 물체표면과 자유표면에 법선 다이폴과 쏘스를 분포시켰으며, 자유표면 및 방사조건을 만족시키기 위해 상류차분식을 사용하는 대신 B-spline의 기저함수를 미분하여 선형화된 자유표면 경계조건에 직접 적용하였다. 이 방법을 적용함으로써 Dawson 방법에서 문제가 되었던 수치감쇠 문제를 피할 수 있었다. 수치계산 프로그램을 검증하기 위해 2차원 원주주위의 유동계산과 날개면 주위의 유동해석을 수행하였으며, B-spline 기저 고차패널법에 의한 수치계산 결과가 저차패널법에 비해 빠른 수렴성과 정확성을 보였고 계산에 필요한 패널 수가 현저히 줄어드는 대단히 만족스러운 결과를 얻었다.

• 수산해양기술연구
• /
• 제17권2호
• /
• pp.109-113
• /
• 1981

Most existing theories on ship waves and wave resistance are based on the perturbation of the flow field by a small pararr.eter which specifies the slenderness of the ship hull. Since however, ship hulls in practice are neither so slender nor thin enough to secure the validity of the linearized theory, the agreen:ent between the theoretical prediction and the experimental result is not generally satisfactory. The author pointed out that the contribution by the non-linear term in the free surface condition can be represented by sorr.e source distribution over the still water plane. This paper leads to a forrr.ula for the wave resistance of not slender ships at low Froude nurr.bers. and deals with the asynptotic expression. As a nurr.erical example, the wave resistance of Wigley model is calculated, and the result is compared with experimental values. It is concluded that the wave resistance coefficient varies in the rate of Fn6 at low speed limit in general. A comparison with the result derived from the linearized free surface condition shows that the non-linearity of the free surface is irr portant at low speed.

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

• Ocean Systems Engineering
• /
• 제6권3호
• /
• pp.245-264
• /
• 2016

The iterative boundary element method (IBEM) developed originally before for cavitating two-dimensional (2-D) and three-dimensional (3-D) hydrofoils moving under free surface is modified and applied to the case of 2-D (two-dimensional) airfoils and 3-D (three-dimensional) wings over water. The calculation of the steady-state flow characteristics of an inviscid, incompressible fluid past 2-D airfoils and 3-D wings above free water surface is of practical importance for air-assisted marine vehicles such as some racing boats including catamarans with hydrofoils and WIG (Wing-In-Ground) effect crafts. In the present paper, the effects of free surface both on 2-D airfoils and 3-D wings moving steadily over free water surface are investigated in detail. The iterative numerical method (IBEM) based on the Green's theorem allows separating the airfoil or wing problems and the free surface problem. Both the 2-D airfoil surface (or 3-D wing surface) and the free surface are modeled with constant strength dipole and constant strength source panels. While the kinematic boundary condition is applied on the airfoil surface or on the wing surface, the linearized kinematic-dynamic combined condition is applied on the free surface. The source strengths on the free surface are expressed in terms of perturbation potential by applying the linearized free surface conditions. No radiation condition is enforced for downstream boundary in 2-D airfoil and 3-D wing cases and transverse boundaries in only 3-D wing case. The method is first applied to 2-D NACA0004 airfoil with angle of attack of four degrees to validate the method. The effects of height of 2-D airfoil from free surface and Froude number on lift and drag coefficients are investigated. The method is also applied to NACA0015 airfoil for another validation with experiments in case of ground effect. The lift coefficient with different clearance values are compared with those of experiments. The numerical method is then applied to NACA0012 airfoil with the angle of attack of five degrees and the effects of Froude number and clearance on the lift and drag coefficients are discussed. The method is lastly applied to a rectangular 3-D wing and the effects of Froude number on wing performance have been investigated. The numerical results for wing moving under free surface have also been compared with those of the same wing moving above free surface. It has been found that the free surface can affect the wing performance significantly.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.381-384
• /
• 2002

The free surface flow problem has been one of the most interesting and challenging topic in the area of the naval ship hydrodynamics and ocean engineering field. The problem has been treated mainly in the scope of the potential theory and its governing equation is well known Laplace equation. But in general, the exact solution to the problem is very difficult to obtain because of the nonlinearlity of the free surface boundary condition. Thus the linearized free surface problem has been treated often in the past. But as the computational power increases, there is a growing trend to solve the fully nonlinear free surface problem numerically. In the present study, a time-dependent finite element method is developed to solve the problem. The initial-boundary problem is formulated and replaced by an equivalent variational formulation. Specifically, the computations are made for a highly nonlinear flow phenomena behind a transom stern ship and a vertical strut piercing the free surface.

• International Journal of Ocean System Engineering
• /
• 제1권1호
• /
• pp.28-31
• /
• 2011

An attempt was made to compute the free surface deformation due to the impact of a water droplet. The Cauchy Poisson, i.e. the initial value problem, was solved with the kinematic and dynamic free surface boundary conditions linearized. The zero order Hankel transformation and Laplace transform were applied to the related equations. The initial condition for the free surface profile was derived from a captured video image. The effect of the surface tension was not significant with the water mass used in this investigation. The computed and observed free surface deformations were compared.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.156-161
• /
• 1998

This describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows of fluid behaviour with free-surface. The elliptic differential equations governing the flows have been linearized by means of finite-difference approximations, and the resulting equations have been solved via a fully-implicit iterative method. The free-surface is defined by the motion of a set of marker particles and interface behaviour was investigated by way of a 'Lagrangian' technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions or experimental results from the literature. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• International Journal of Ocean System Engineering
• /
• 제2권1호
• /
• pp.32-36
• /
• 2012

This paper presents a comparison of potential and viscous computational codes for the water entry problem. A po-tential code was developed which adopted the boundary element method to solve the problem. A nonlinear free surface boundary condition was integrated to find new locations of free surface. The dynamic boundary condition was simplified by taking constant potential values for every time steps. The simplified dynamic boundary condition was applied in the new position of the free surface not at the mean level, which is the usual practice for linearized theory. The commercial code FLUENT was used to solve the water entry problem from the viscosity point of view. The movement of the air-liquid interface is traced by distribution of the volume fraction of water in a computational cell. The pressure coefficients were compared with each other, while experimental results published by other researchers were also examined. The characteristics of each method were discussed to clarify merits and limitations when they were applied to the water entry problems.

• 대한조선학회지
• /
• 제27권4호
• /
• pp.32-42
• /
• 1990

선체표면 및 자유표면에 Rankine Source를 분포하는 방법에 의하여 선체 주위의 유동의 수치계산을 수행하였다. 선체표면 및 자유표면은 사각형 Panel들로 표시되며 자유표면 조건은 이중모형 흐름에 의해 선형화 되어 C.W. Dawson의 유한차분법에 따라 교란없는 자유표면에 적용되었다. Wigley 선형 및 Series 60, $C_B=0.6$ 선형에 대한 Fixed Condition에서의 조파저항, 선측파고, 압력분포 및 Trim & Sinkage 등을 계산하였으며 계산된 결과는 국내외 수조에서의 계측치와 비교하였다. 또한, 선체표면과 자유표면의 Panel 분할조건 및 자유표면의 설정영역의 변화에 따른 계산치의 영향도 아울러 조사하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.191-192
• /
• 2003

A numerical method for a wave diffraction problem in three-dimensional channels is developed. The physical models are various shapes of channel connected to the open sea. When a ship or an offshore structure is moored in various configurations of channel connected to an open sea, the prediction of the hydrodynamic force exerting on the moored ship could be important for the prediction of its motion. It is assumed that the fluid is inviscid and incompressible and its motion is irrotational. From the continuity equation, the Laplace equation can be obtained as the governing equation. The surface tension at free surface is neglected, and wave amplitude is assumed to be small compared to the wave length. Then the free surface condition can be linearized. The numerical method used here is the localized finite element method based on a variational formulation