• Journal of computational fluids engineering
• /
• v.18 no.1
• /
• pp.13-21
• /
• 2013

Cavitating flow is usually formed on the surface of a high speed underwater object. When a object moves near a free surface at very high speed, the cavity signature becomes one of the major factors to be overcome by sensors of military satellite. The present work was to study the free surface effect on the ventilated cavitation. The governing equations were Navier-Stokes equations based on a homogeneous mixture model. The multiphase flow solver used an implicit preconditioning method in the curvilinear coordinate system. The cavitation model used here was the one first presented by Merkle et al.(2006) and redeveloped by Park & Ha(2009). Computations considered the free surface effects were carried out with a NACA0012 hydrofoil and the corresponding results were compared with the experimental data to have a good agreement. Calculations were then performed considering the ventilated cavitation, including the effect of non-condensable gas under the free surface effects.

• 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 Mechanical Science and Technology
• /
• v.19 no.12
• /
• pp.2263-2269
• /
• 2005

The ship wave phenomena in the restricted waterway were investigated by a numerical analysis. The Euler and continuity equations were employed for the present study. The boundary fitted and moving grid system was adopted to enhance the computational efficiency. The convective terms in the governing equations and the kinematic free surface boundary condition were solved by the Constrained Interpolated Profile (CIP) algorithm in order to solve accurately wave heights in far field as well as near field. The advantage of the CIP method was verified by the comparison of the computed results by the CIP and the Maker and Cell (MAC) method. The free surface flow simulation around Wigley hull was performed and compared with the experiment for the sake of the validation of the numerical method. The present numerical scheme was applied to the free surface simulation for various canal sections in order to understand the effect of the sectional shape of waterways on the ship waves. The wave heights on the side wall and the shape of the wave patterns with their characteristics of flow are discussed.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2000.04a
• /
• pp.140-145
• /
• 2000

The spectral method which is composed of an eigenfunction expansion of free modes in the wave number domain is used to produce two dimensional unsteady inviscid wave simulation such as progressive waves in a numerical pneumatic wave tank. A spatial and time dependent free surface elevation and the potential are calculated by integrating ODE derived from fully nonlinear kinematic and dynamic free surface boundary condition at each time step. The nonlinear characteristics in the waves by this method were notable as increasing wave steepness. This method is very useful and powerful in terms of saving computational time caused by rapid convergence exponentially with increasing number of nodes, even preserving accurate numerical results. Moreover, it will given us many possibilities to apply to naval and ocean engineering fields.

• Proceedings of the Korea Water Resources Association Conference
• /
• 1987.07a
• /
• pp.223-234
• /
• 1987

A weakly nonlinear solution is presented for the two-dimensional wave kinematics forced by a generic wavemaker of variable-draft. The solution is valid for both piston and hinged wavemakers of variable-draft that may be double articulated. The second-order propagating waves generated by a planar wave board are composed of two components; viz., a Stokes second-order wave and a second-harmonic wave forced by the wavemaker which travels at a different speed. A previously neglected time-independent solution that is required to satisfy a kinematic boundary condition on the wavemaker as well as a mixed boundary condition on the free surface is included for the first time. A component of the time-independent solution is found to accurately estimate the mean return current(correct to second-order) in a closed wave flume. This mean return current is usually estimated from kinematic considerations by a conservation of mass principle.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.36 no.2
• /
• pp.70-79
• /
• 2024

In this study, we introduce a linear spectral method capable of simulating wave generation and transformation caused by a moving bottom in a 3-dimensional space. The governing equations are linear dynamic free-surface boundary conditions and linear kinematic free-surface boundary conditions, which are solved in Fourier space. Solved velocity potential and free-surface displacement should satisfy continuity equation and kinematic bottom boundary condition. For numerical analysis, a 4th order Runge-Kutta method was utilized to analyze the time integral. The results obtained in Fourier space can be converted into velocity potential and free-surface displacement in a real space using inverse Fourier transform. Regular waves generated by various types of moving bottoms were simulated with the linear spectral method. Additionally, obliquely generated regular waves using specified bottom movements were simulated. The results obtained from the spectral method were compared to analytical solutions, showing good agreement between the two.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.23 no.2
• /
• pp.1-13
• /
• 1986

The nonlinear hydrodynamic forces acting on a two-dimensional circular cylinder, oscillating with large amplitude in the free surface, are calculated by using the Semi-Lagrangian Time-Step-ping Method used by O.M. Faltinsen. In present calculation the position and the potential value of free surface are calculated using the exact kinematic and dynamic free surface boundary condition. At each time step an integral equation is solved to obtain the value of potential and normal velocity along the boundaries, consisting of both the body surface and the free surface. Some effort was devoted to the elimination of instability arising in the range of high frequency. Numerical simulations were performed up to the 3rd or 4th period which seems to be enough for the transient effect to die out. Each harmonic component and time-mean force are obtained by the Fourier transform of forces in time domain. The results are compared with others' experimental and theoretical results. Particularly, the calculation shows the tendency that the acceleration-phase 1st-harmonic component(added mass) increases as the motion amplitude increases and a reverse tendency in the velocity-phase 1st-harmonic component(damping coefficient). The Yamashita's experimental result also shows the same tendency. In general, the present result show relatively good agreement with the Yamashita's experimental result except for the time-mean force.

• 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 Ocean Engineering and Technology
• /
• v.37 no.2
• /
• pp.49-57
• /
• 2023

Since Chappelear developed a Fourier approximation method, considerable research efforts have been made. On the other hand, Fourier approximations are unsuitable for deep water waves. The purpose of this study is to provide a Fourier approximation suitable even for deep water waves and a numerical method to determine the Fourier coefficients and the wave properties. In addition, the convergence of the solution was tested in terms of its order. This paper presents a velocity potential satisfying the Laplace equation and the bottom boundary condition (BBC) with a truncated Fourier series. Two wave profiles were derived by applying the potential to the kinematic free surface boundary condition (KFSBC) and the dynamic free surface boundary condition (DFSBC). A set of nonlinear equations was represented to determine the Fourier coefficients, which were derived so that the two profiles are identical at specified phases. The set of equations was solved using Newton's method. This study proved that there is a limit to the series order, i.e., the maximum series order is N=12, and that there is a height limitation of this method which is slightly lower than the Michell theory. The reason why the other Fourier approximations are not suitable for deep water waves is discussed.

• Journal of the Society of Naval Architects of Korea
• /
• v.46 no.6
• /
• pp.650-658
• /
• 2009

Two dimensional sloshing phenomena in regularly excited liquid cargo tank are numerically simulated with finite difference method. Navier-Stokes equations and continuity equation are computed for this study. The free-surface is determined every time step satisfying kinematic boundary condition using marker-density method. And the exciting force is treated by adding the acceleration of the tank to source term. The results are compared with other existing experiment results. And the comparison results show a good agreement. The sloshing phenomena in the tank of the 138K LNG carrier in sway motion is simulated with present calculation methods in low filling level. To find the relations between impact pressure and excitation condition, the calculations are performed in various amplitudes and periods. The averaged maximum pressures are compared each other.