• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.79-85
• /
• 2014

An air-water free surface flow simulation by using the Lattice Boltzmann Method(LBM) has not been studied a lot compared with the done by the Navier-Stoke equation. This paper shows the LBM is as one of the application tools for the free surface movement over an obstacle. The Mezo scaled application tool has been developed with two dimensional and nine discretized velocity direction using conventional lattice Bhatnagar-Gross-Krook model. Boundary conditions of a halfway-based for solid wall and a kinematic-based for interface are adopted. A validation case with a trapezoidal shape bump to make a comparison between freesurface movements from computational results and experimental ones was described with grid size dependency.

• International Journal of Ocean Engineering and Technology Speciallssue:Selected Papers
• /
• v.3 no.1
• /
• pp.14-22
• /
• 2000

The improved Green integral equation for the calculation of time-harmonic potentials in the radiation diffraction problem about a freely floating body in the presence of moderate or weak current is presented. The forward-speed Green function presented by Brard is used. The correct free surface boundary conditions on the physical free surface are employed as well as an appropriate boundary conditions on the non-physical inner free surface. The default in the existing Green integral equation as well as in the source integral equation is discussed in detail.

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

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

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.17 no.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.

• Water for future
• /
• v.29 no.2
• /
• pp.221-233
• /
• 1996

To investigate the hydrodynamic characteristics of a two-dimensional buoyant surface jet, the most important factors of the numerical analysis are the evaluation of the free surface and the turbulence transportation under the stratification. In present study, a numeriacal simulation model used with the semi-implicit method for pressure-linked equations (SIMPLE), the non-hydrostatic approximation and the algebraic stress model (ASM) is applied to investigate the vertical structure of internal flow hydrodynamically. The ASM enables to take account of anisotropy of turbulence, the damping effects of the density interface, and the free surface on the turbulence structure accurately. The ASM tested produces better agreement than the $\kappa-\varepsilon$ model with measurements by Nakatsuji (1984) on the flow development and turbulence structure. Applicability of the ASM to a two-dimensional buoyant surface jet is examined through comparison with experimental data.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.4
• /
• pp.66-74
• /
• 1992

A panel method in the spirit of Hess & Smith(1962), and also of Dawson(1977) was developed to compute the wave resistance of a submerged, or a surface piercing, body moving in the water of finite depth. As a boundary condition on the free surface what is called the Poisson equation is used, while Yasukawa(1989) chose the Dawson equation for which the double-body flow is regarded as the basic one. In order to satisfy the boundary condition on the bottom surface automatically, the sum of a Rankine source and its image with respect to the bottom surface is chosen as the Green function, and hence the singularity is distributed only on the body and on the free surface thereby decreasing the required number of panels dramatically, compared to that of Yasukawa, without the consequential loss of accuracy. Calculations were done for a submerged sphere and for the Wigley hull, and the results are compared with other existing analytical and numerical data.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.16 no.4
• /
• pp.255-262
• /
• 2013

The boltzmann equation is based on the particle distribution function while the Navire-Stokes equation based on the continuum theory. In order to simulate free surface flow, this paper used the Lattice Boltzmann Method of which is the discretized form. The detail study on the characteristics of the Lattice Boltzmann Method for the free surface simulation was investigated. The developed code was validated with the traditional dam breaking problem by tracking the front position of the water. A basic roles of density functions in the Lattice Boltzmann Method is discussed. To have an engineering applications, the simulation is also conducted the free surface behavior with an arbitrary wall geometry.

• International Journal of Ocean System Engineering
• /
• v.1 no.3
• /
• pp.136-143
• /
• 2011

In this paper, three-dimensional free-surface flows are simulated by using two different numerical methods, the constrained interpolation profile (CIP)-based and finite volume (FV)-based methods. In the CIP-based method, the governing equations are solved on stationary staggered Cartesian grids by a finite difference method, and an immersed boundary technique is applied to deal with wave-body interactions. In the FV-based method, the governing equations are solved by applying collocated finite volume discretization, and body-fitted meshes are used. A free-surface boundary is considered as the interface of the multi-phase flow with air and water, and a volumeof-fluid (VOF) approach is applied to trace the free surface. Among many variations of the VOF-type method, the tangent of hyperbola for interface capturing (THINC) and the compressive interface capturing scheme for arbitrary meshes (CICSAM) techniques are used in the CIP-based method and FV-based method, respectively. Numerical simulations have been carried out for dam-breaking and wave-body interaction problems. The computational results of the two methods are compared with experimental data and their differences are observed.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.409-413
• /
• 2005

Even though the relative importance of length scale of flow system allow us to simplify three dimensional flow problem to one or two dimensional representation, many systems still require three dimensional analysis. The objective of this study is to develop an efficient and accurate finite element model for analyzing and predicting three dimensional flow features in natural rivers and to offend to model spreading of pollutants and transport of sediments in the future. Firstly, three dimensional Reynolds averaged Navier-Stokes equations with the hydrostatic pressure assumption in generalized curvilinear coordinates were combined with the kinematic free-surface condition. Secondly. to simulate realistic high Reynolds number flow, the model employed the Streamline Upwind/Petrov-Galerkin(SU/PG) scheme as a weighting function for the finite element method in conjunction with an appropriate turbulence model(Smagorinsky scheme for the horizontal plain and Mellor-Yamada scheme for the vertical direction). Several tests is performed for the purpose of validation and verification of the developed model. A simple rectangular channel, 5-shaped and U-shaped channel are used for tests and comparisons are made with RMA-10 model. Runs for each case is converged stably without a oscillation and calculated water-surface deformation, longitudinal and transversal velocities, and velocity vector fields are in good agreement with the results of RMA-10 model.