• Journal of applied mathematics & informatics
• /
• v.26 no.1_2
• /
• pp.45-60
• /
• 2008

Here we consider the evaluation of the the dynamic component of the second order force due to wave diffraction by a circular cylinder analytically and numerically. The cylinder is fixed, vertical, surface piercing in water of finite uniform depth. The formulation of the wave-structure interaction is based on the assumption of a homogeneous, ideal, incompressible, and inviscid fluid. The nonlinearity in the wave-structure interaction problem arises from the free surface boundary conditions, namely, dynamic and kinematic free surface boundary conditions. We expand the velocity potential and free surface elevation functions in terms of a small parameter and then consider the second order diffraction problem. After deriving the pressure using Bernoulli's equation, we obtain the analytical expression for the dynamic component of the second order force on the cylinder by integrating the pressure over the wetted surface. The computation of the dynamic force component requires only the first order velocity potential. Numerical results for the dynamic force component are presented.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6B
• /
• pp.691-696
• /
• 2008

In this paper, a three-dimensional non-hydrostatic pressure model for free surface flows using a normalized vertical coordinate system is presented. To strongly couple the free surface and non-hydrostatic pressure in the momentum equations, a double predictor-corrector method is employed. This research is especially focused on implementing the dynamic boundary condition (a zero pressure condition) at the free surface. This boundary condition can be specified accurately with a small modification to existing models. Numerical results with and without this modification clearly show that a precise implementation of the dynamic boundary condition is paramountly important.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.10 no.1
• /
• pp.103-109
• /
• 2010

In this study, a three-dimensional non-hydrostatic pressure model based on a normalized vertical coordinate system for free surface flows is presented. To strongly couple the free surface and non-hydrostatic pressure with the momentum equations, a double predictor-corrector method is employed. The study is especially focused on implementing the dynamic boundary condition (a zero pressure condition) at the free surface with ignoring of the atmospheric pressure. It is shown that the boundary condition can be specified easily with a slight modification to existing models.

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

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• 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.

• International Journal of Ocean System Engineering
• /
• v.1 no.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.

• International Journal of Ocean System Engineering
• /
• v.2 no.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.

• Journal of Advanced Research in Ocean Engineering
• /
• v.2 no.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.

• Journal of Ship and Ocean Technology
• /
• v.6 no.2
• /
• pp.37-50
• /
• 2002

In the present numerical simulation of viscous free surface flow around a ship, two-fluids in-compressible Reynolds-averaged Navier-Stokes equations with the standard $\textsc{k}-\varepsilon$turbulence model are discretized on a regular grid by using a finite volume method. A local level-set method is introduced for capturing the free surface movement and the influence of the viscous layer and dynamic boundary condition of the free surface are implicitly considered. Partial differential equations in the level-set method are discretized with second order ENO scheme and explicit Euler scheme in the space and time integration, respectively. The computational results for the Series-60 model with $C_B=0.6$ show a good agreement with the experimental data, but more validation studies for commercial complicated hull forms are necessary.

• Ocean Systems Engineering
• /
• v.6 no.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.