• 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 Korea Water Resources Association
• /
• v.44 no.6
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

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

• Journal of Ocean Engineering and Technology
• /
• v.26 no.5
• /
• pp.55-62
• /
• 2012

In this paper, the rigid lid boundary condition is applied to simulate the influence of floating structures such as ships or pontoons, and the pressure term in both the momentum equations and continuity equation are modified. The pressure of a floating structure under the free surface is dependent on the draft of the structure, generally called a ship. If the free surface is covered by a floating structure, the free surface cannot move freely. The water level should be fixed, using a rigid lid boundary condition. This boundary condition is implemented by reducing the storage area of the grid cell with a factor between zero and one. The numerical model developed by Hong (2009) is verified through a comparison with experimental results, and the influence of the reduction factor is investigated using the verified numerical model.

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

• Journal of the Society of Naval Architects of Korea
• /
• v.40 no.2
• /
• pp.11-20
• /
• 2003

In this paper, the high-order spectral/boundary-element method is developed to calculate the 3-dimensional water waves generated by a submerged body. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated Tn time-domain. Three-dimensional free-surface flows generated by a submerged sphere which is moving under the free-surface are calculated. Through example calculations, nonlinear effects on free-surface profiles and hydrodynamic forces are shown. Comparisons with others' results show good agreements.

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

• Geophysics and Geophysical Exploration
• /
• v.3 no.2
• /
• pp.61-69
• /
• 2000

The free-surface boundary conditions are persistent problem in elastic wave equation modeling by finite-difference method, which can be summarized with the degradation of the accuracy of the solution and limited stability range in Poisson's ratio. In this paper, we propose the mid-point averaging scheme as an alternative way of implementing the free-surface boundary conditions, and present the solution to Lamb's problem to verify our approach.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.1
• /
• pp.233-243
• /
• 2019

In this paper, numerical and experimental studies are performed to investigate the liquid impact on a free falling wedge. In the numerical simulation, the structure is assumed to be rigid and the elastic response is ignored. The fully nonlinear coupling between wedge and water is considered by an auxiliary function method based on the Boundary Element Method (BEM). At the intersection of the wedge surface and liquid surface, two coincident nodes are used to decouple the boundary conditions. The Eulerian free surface conditions in the local coordinate system are adopted to update the deformed free surface. In the experiments, five pressure sensors are fixed on each side of the wedge which is released from an experimental installation. Steel and aluminum wedges that have different structural elasticity are used in the experiments to investigate the influence of structural elasticity on the impact force. Numerical results are compared with experimental data and they agree very well. The influence of fluid gravity, body mass, initial entry speed and deadrise angle on the impact pressure are further investigated.

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