• Journal of computational fluids engineering
• /
• v.20 no.1
• /
• pp.57-64
• /
• 2015

This paper describes a volume of fluid(VOF) method, mRHRIC for the simulation of free surface flows around ship hulls and provides its validation against benchmark test cases. The VOF method is developed on the basis of RHRIC method developed by Park et al. that uses high resolution differencing schemes to algebraically preserve both the sharpness of interface and the boundedness of volume fraction. A finite volume method is used to solve the governing equations, while the realizable ${\kappa}-{\varepsilon}$ model is used for turbulence closure. The present numerical results of the resistance performance tests for DTMB5415 and KCS hull forms show a good agreement with available experimental data and those of other free surface methods.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2558-2561
• /
• 2008

In the present study, a parallel Taylor-Galerkin/level set based two-phase flow code was developed using finite element discretization and domain decomposition method based on MPI (Message Passing Interface). The proposed method can be utilized for the analysis of a large scale free surface problem in a complex geometry due to the feature of FEM and domain decomposition method. Four-step fractional step method was used for the solution of the incompressible Navier-Stokes equations and Taylor-Galerkin method was adopted for the discretization of hyperbolic type redistancing and advection equations. A Parallel ILU(0) type preconditioner was chosen to accelerate the convergence of a conjugate gradient type iterative solvers. From the present parallel numerical experiments, it has been shown that the proposed method is applicable to the simulation of large scale free surface flows.

• Journal of Advanced Research in Ocean Engineering
• /
• v.4 no.4
• /
• pp.163-173
• /
• 2018

There are many issues in fluid structure interactions when dealing with the free surface flows in a sloshing tank. For example the problem of how yielding a highly nonlinear wave with a simple forced motion over a short duration is of concern here. Nonlinear waves are generated in a rectangular tank which is forced horizontally; its motion consists of a single cycle of oscillation. One of the objectives is to end up with a shape of the free surface yielding a wide range of critical flows by tuning few parameters. The configuration that is studied here concerns a plunging breaker accompanied with a critical jet where great kinematics are simulated. The numerical simulations are performed with a twodimensional code which solves the fully nonlinear free surface boundary conditions in Potential Theory.

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

• International Journal of Fluid Machinery and Systems
• /
• v.7 no.4
• /
• pp.174-182
• /
• 2014

The ultimate objective of this study is to develop a water turbine appropriate for low-head open channels to effectively utilize the unused hydropower energy of rivers and agricultural waterways. The application of a cross-flow runner to open channels as an undershot water turbine has been considered and, to this end, a significant simplification was attained by removing the turbine casing. However, the flow field of an undershot cross-flow water turbine possesses free surfaces, and, as a result, the water depth around the runner changes with variation in the rotational speed such that the flow field itself is significantly altered. Thus, clear understanding of the flow fields observed with free surfaces to improve the performance of this turbine is necessary. In this study, the performance of this turbine and the flow field were evaluated through experiments and numerical analysis. The particle image velocimetry technique was used for flow measurements. The experimental results reflecting the performance of this turbine and the flow field were consistent with numerical analysis. In addition, the flow fields at the inlet and outlet regions at the first and second stages of this water turbine were clarified.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.5 no.3
• /
• pp.404-413
• /
• 2013

In the present study, a towed underwater particle image velocimetry (PIV) system was validated in uniform flow and used to investigate the free-surface effects on the turbulent wake of a simple surface-piercing body. The selected test model was a cylindrical geometry formed by extruding the Wigley hull's waterplane shape in the vertical direction. Due to the constraints of the two-dimensional (2D) PIV system used for the present study, the velocity field measurements were done separately for the vertical and horizontal planes. Using the measured data at several different locations, it was possible to identify the free-surface effects on the turbulent wake in terms of the mean velocity components and turbulence quantities. In order to provide an accuracy level of the data, uncertainty assessment was done following the International Towing Tank Conference standard procedure.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.4 s.142
• /
• pp.322-329
• /
• 2005

Viscous flows behind transom stern are analyzed based on CFD simulation results. Stern wave pattern is often complicated due to the abrupt change of stern surface curvature and flow separation at transom. When a ship advances at high speed, whole transom stern is exposed out of water, resulting in the so-called 'dry transom'. However, in the moderate speed regime, stern wave development in conjunction of flow separation makes unstable wavy surface partially covering transom surface, i.e., the so-called 'wetted transom'. Transom wave formation is usually affecting the resistance characteristics of a ship, since the pressure contribution on transom surface as well as the wave-making resistance is changed. Flow modeling for 'wetted transom' is difficult, while the 'dry transom modeling' is often applied for the high-speed vessels. In the present study CFD results from the RANS equation solver using a finite volume method with level-set treatment are utilized to assess the topology of transom flow pattern for a destroyer model (DTMB5415) and a container ship (KCS). It is found that transom flow patterns are quite different for the two ships, in conformity to the shape of submerged transom. Furthermore, the existence of free surface seems to after the flow topology in case of KCS.

• Journal of Power System Engineering
• /
• v.16 no.5
• /
• pp.40-46
• /
• 2012

For free surface flow problem, a high-order spectral/boundary element method is adapted as an efficient numerical tool. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated and hydrodynamic forces also can be calculated in time domain. In this method, the velocity potential is expressed as the sum of surface potential and body potential. Then, surface potential is solved by using the high-order spectral method and body potential is solved by using the high-order boundary element method. Using the combination of these two methods, the free surface flow problems of a submerged moving body are solved in time domain. In the present study, lifting surface theory is added to the former work to include effects of lift force. Therefore, a new formulation for the basic mathematical theory is introduced to contain the lift body in calculation.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.25-30
• /
• 2011

In this study, it is performed numerical simulation on multi-phase flow by means of CIP-CSI2 scheme. It is applied In a two-phase free surface flow problem at a high density ratio equivalent to that of an air-water system, for examining the computational capability. The method that is being developed and improved is a CIP(Constrained Interpolation Profile) and CSL2(Conservative Semi-Lagrangian) based Cartesian Grid Method.

• Korea-Australia Rheology Journal
• /
• v.12 no.2
• /
• pp.107-117
• /
• 2000

Prediction of the product shape in sheet casting process is performed from the numerical simulation. A three-dimensional finite element method is used to investigate the flow behavior and to examine the effects of processing conditions on the sheet produced. Effects of inertia, gravity, surface tension and non-Newtonian viscosity on the thickness profile of the sheet are considered since the edge bead and the flow patterns in the chill roll region have great influence on the quality of the products. In the numerical simulation with free surface flows, the spine method is adopted to update the free surface, and the force-free boundary condition is imposed along the take-up plane to avoid severe singularity problems existing at the take-up plane. From the numerical results of steady isothermal flows of a generalized Newtonian fluid, it is shown that the draw ratio plays a major role in predicting the shape of the final sheet produced and the surface tension has considerable effect on the bead thickness ratio and the bead width fraction, while shear-thinning and/or tension-thickening viscosity affect the degree of neck-in.