• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.19-26
• /
• 1999

This paper describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows with free-surface. The Navier-Stokes equations governing the flows have been discretized by means of finite-difference approximations, and the resulting equations have been solved via the SIMPLE-C algorithm. The free-surface is defined by the motion of a set of marker particles and the interface behaviour was investigated by means 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. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.13 no.1
• /
• pp.1-11
• /
• 2021

Performance of a NACA 63₄-021 hydrofoil in motion under and in close proximity of a free surface for a large range of angles of attack is studied. Lift and drag coefficients of the hydrofoil at different submergence depths are investigated both numerically and experimentally, for 0° ≤ AoA ≤ 30° at a Reynolds number of 10⁵. The results of the numerical study are in good agreement with the experimental results. The agreement confirms the new finding that for a submerged hydrofoil operating at high angles of attack close to a free surface, the interaction between the hydrofoil-motion induced waves on the free surface and the hydrofoil results in mitigation of the flow separation characteristics on the suction side of the foil and delay in stall, and improvement in hydrofoil performance. In comparing with a baseline case, results suggest a 55% increase in maximum lift coefficient and 90% average improvement in performance for, based on the lift-to-drag ratio, but it is also observed significant decrease of lift-to-drag ratio at lower angles of attack. Flow details obtained from combined finite volume and volume of fluid numerical methods provide insight into the underlying enhancement mechanism, involving interaction between the hydrofoil and the free surface.

• Journal of the Society of Naval Architects of Korea
• /
• v.40 no.4
• /
• pp.8-15
• /
• 2003

This paper presents the method for developing an optimum hull form with minimum wave resistance using SQP(sequential quadratic programming) as an optimization technique. The wave resistance is evaluated by a Rankine source panel method with non-linear free surface conditions and the ITTC 1957 friction line is used to predict the frictional resistance coefficient. The geometry of the hull surface is represented and modified using B-spline surface patches. The optimization method is applied to Series 60 hull and KCS(KRISO 3600 TEU Container Ship). The obtained results prove that the method is appropriate for preliminary hull form design.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.9
• /
• pp.1230-1243
• /
• 1997

A numerical technique for simulating incompressible viscous flow with free surface is presented. The flow field is obtained by penalty finite element formulation. In this work, a modified volume of fluid (VOF) method which is compatible with 4-node element is proposed to track the moving free surface. This scheme can be applied to irregular mesh system, and can be easily extended to three dimensional geometries. Numerical analyses were done for two benchmark examples, namely the broken dam problem and the solitary wave propagation problem. The numerical results were in close agreement with the existing data. Illustrative examples were studied to show the effectiveness of the proposed numerical scheme.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.1
• /
• pp.15-23
• /
• 1998

Characteristics of turbulence beneath the free surface around a blunt bow are numerically-investigated. Three dimensional wavier-Stokes and continuity equations are solved for the simulations.. The Large Eddy Simulation(LES) with the external disturbance is performed to simulate the turbulent free surface flow called sub-breaking wave. The result shows that the free surface fluctuates beyond a certain critical condition and the characteristics of the fluctuation are similar to the turbulent boundary layer flow around a solid body.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.17 no.4
• /
• pp.13-20
• /
• 1980

The integral equation method to solve the boundary-value problem of a 2-dimensional body oscillating in the presence of a free surface generally breaks down at and near irregular frequencies due to the hypothetical flow inside the body. In this paper singularity distributions were extended to an inner free surface to remove the irregular frequency as Ohmatsu's work in 1978, and the solution for the above problem was found by using stream function. For various bodies including Lewis form cylinders, the hydrodynamic forces were calculated numerically at various wave numbers. From the results we concluded that the irregular frequencies can be removed even for the Lewis form cylinder as Ohmatsu done for circular cylinders, and calculated hydrodynamic forces by the present method are little higher than those of Ohmatsu's when the singularities are put on the inner free surface of the body. We specially point out that the solution for heaving motion converges in an oscillatory manner but not for swaying and rolling motions.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.19 no.2
• /
• pp.33-39
• /
• 1982

It is one of the basic problems of naval architecture and ocean engineering how to describe the wave kinematics normally under the assumption of an ideal fluid. At present, there are many wave theories available for design purposes. These can be classified into two groups: One is the analytic theory and the other is the numerical theory. This paper briefly introduces the stream function method of R.G. Dean which belongs to the latter group and shows its numerical evaluations exemplified for two cases: One is applied to observed waves and the other is for design waves. In the former case, the wave profiles are calculated by the stream function method and compared with those of the observed waves and also with the results of R.G. Dean. They show good agreement. In the latter case, the wave kinematics and wave loads on a column of diameter 1m are calculated by the stream function method and these are compared with those resulted from the 5th-order gravity wave theory. As a result of comparison the values by the stream function method are slightly larger than those by the 5th-order gravity wave theory but the difference are negligible. From this it is concluded that the stream function method is very useful. And as characteristics of the numerical theories, the stream function method of R.G. Dean can be easily extended to the higher order terms and can include easily the current velocity and the pressure distribution on the free surface. In addition, when the data of observed wave profile are given, this method can reproduced the observed wave profile as closely as possible so that this method seems to describe the ocean wave more realistically. And from standpoint of a mathematical principle the stream function method exactly satisfies the kinematic free-surface boundary condition.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.14 no.2
• /
• pp.171-181
• /
• 2002

Among various numerical methods of wave transformations including wave breaking by structure, models using VOF(Volume Of Fluid) method to trace free surface are getting into the spotlight recently. In order to analyze wave transformations and velocity of the wave fields due to the permeable submerged breakwater(PSB), This study applied VOF method to the two-dimensional wave channel installed line-source to generate waves and added dissipation zone to offer a non-reflective boundary. Hydraulic experiments was performed to obtain the application of two-dimensional numerical wave channel. The results of numerical experiments using the two-dimensional wave channel agree well with the experimental data. It was shown that vortices are formed behind the PSB, and in case of the 2-rowed PSB they also are occurred in between PSBs, strongly non-linear waves are developed on the crown of the PSB, and the direction of velocities in porous media is determined by the shape of free surface.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.12
• /
• pp.1555-1569
• /
• 2000

Numerical simulation of fluid flow with moving free surface has been carried out. For the free surface flow, a VOF(Volume of Fluid)-based algorithm utilizing a fixed grid system has been investigated. In order to reduce numerical smearing at the free surface represented on a fixed grid system, a new free surface tracking algorithm based on the donor-acceptor scheme has been presented. Novel features of the proposed algorithm are characterized as two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed algorithm can be easily implemented in any irregular non-uniform grid systems that are usual in finite element method (FEM). Moreover, the proposed algorithm can be extended and applied to the 3-D free surface flow problem without additional efforts. For computation of unsteady incompressible flow, a finite element approximation based on the explicit fractional step method has been adopted. In addition, the SUPG(streamline upwind/Petrov-Galerkin) method has been implemented to deal with convection dominated flows. Combination of the proposed free surface tracking scheme and explicit fractional step formulation resulted in an efficient solution algorithm. Validity of the present solution algorithm was demonstrated from its application to the broken dam and the solitary wave propagation problems.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.4
• /
• pp.32-42
• /
• 1990

The method using Rankine Soure distribution over the hull surface and undisturbed free surface was applied to calculate the free surface flow around a ship. The ship hull as well as a local portion of the undisturbed free surface arc geometrically represented by quadrilateral panels and the source density is determined so as to satisfy the linearized free surface condition based on the double model flow. The pressure distribution, wave resistance, wave profile and hydrodynamic sinkage force and trim moment for the Wigley hull and the Series 60 hull with $C_B=0.60$ were calculated in the fixed condition. The calculated results were compared with the measured values. The dependance of the solution on the panel arrangement, particularly on the free suraface, was also studied through 11 numerical test cases for the Wigley hull.