• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.1
• /
• pp.58-66
• /
• 2014

Numerical simulations for the flows containing free surface remain difficult problems because the drastic differences of physical properties of water and air, The difference of densities makes the solution instable in particular. For the stabilities of the solutions, the most typical methods to simulate free surface flows, such as Volume Of Fluid(VOF) and Level-Set(LS) methods, impose transient zones where the physical prosperities are continuously distributed. The thickness of the transient zone is the source of the numerical errors. The other side, marker-density method does not use such a transient zone. In the traditional marker-density method, however, the air velocities of free surface cells are extrapolated from the water velocity, and the pressures on the free surface are extrapolated from the air pressures for the stability of the solution. In this study, the marker-density method is modified for the decrease of such numerical errors. That is, the pressure on the free surface is determined to coincide with the pressure gradient terms of the governing equations, and the velocity of free surface cells are calculated with the governing equations. Two-dimensional steady spilling breakers behind of a submersed hydrofoil and three-dimensional spilling breaker near a wedge shaped ship model are simulated using INHAWAVE-II including the modified marker-density(MMD) method. The results are compared with the results of Fluent V6.3 including VOF method and several published research results.

• Journal of Ocean Engineering and Technology
• /
• v.26 no.6
• /
• pp.7-13
• /
• 2012

This paper presents a numerical study of the wave loads acting on offshore structures using a Cartesian-grid-based flow simulation method. Finite volume discretization with a volume-of-fluid (VOF) method is adopted to solve two-phase Navier-Stokes equations. Among the many variations of the VOF method, the CICSAM scheme is applied. The body boundary conditions are satisfied using a porosity function, and wave generation is carried out by using transient (wave or damping) zone approaches. In order to validate the present numerical method, three different basic offshore structures, including a sphere, Pinkster barge, and Wigley model, are numerically investigated. First, diffraction and radiation problems are solved using the present numerical method. The wave exciting and drift forces from the diffraction problems are compared with potential-based solutions. The added mass and wave damping forces from the radiation problems are also compared with the potential results. Next, the wave-induced motion responses of the structures are calculated and compared with the existing experimental data. The comparison results are fairly good, showing the validity of the present numerical method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.11
• /
• pp.1169-1176
• /
• 2011

Liquid film flows are classified into waveless laminar, wavy laminar, and turbulent flows depending on the Reynolds number or the flow stability. Since the wavy motions of the film flows are so intricate and nonlinear, studies on them have largely been experimental. Most numerical approaches have been limited to the waveless flow regime. The various free surface-tracking schemes adopted for this problem were used to more accurately estimate the average film thickness, rather than to capture the unsteady wavy motion. In this study, the wavy motions in laminar wavy liquid film flows with Reynolds numbers of 200-1000 were simulated with various numerical schemes based on the volume of fluid (VOF) method for interface tracking. The results from each numerical scheme were compared with the experimental results in terms of the average film thickness, the wave velocity, and the wave amplitude.

• Journal of computational fluids engineering
• /
• v.20 no.4
• /
• pp.93-101
• /
• 2015

This paper provides numerical results of the simulation for the flow around the hull and the propeller of KCS model ship advancing in shallow water conditions. A finite volume method is used to solve the unsteady Reynolds averaged Navier-Stokes(RANS) equations, where the wave-making problem is solved by using a volume-of-fluid(VOF) method. The wave formed near the hull surface in shallow water conditions shows a deep trough dominant pattern that causes the loss of buoyancy followed by hull squat. The flow past the hull increases as the depth of water decreases. However, the axial flow velocity around the stern shows a reduction in magnitude by the effect of shallow water accompanied by the hull-propeller interaction. As a results, the thrust and torque coefficient increase about 8.3% and 6.2%, respectively for a depth of h/T=3.0 corresponding to a depth Froude number of $F_h=0.693$. The resistance coefficient increases about 11.6% at this Froude number condition.

• Journal of computational fluids engineering
• /
• v.21 no.1
• /
• pp.110-116
• /
• 2016

The present study numerically investigated the deformation of the interface of two-phase fluids flow in a blast furnace. To simulate three-dimensional(3D) incompressible viscous two-phase flow in the furnace filled with the air and molten iron, the volume of fluid(VOF) method based on the finite volume method has been utilized. In addition, the porous medium with the porosity has been considered as the bed of the particles such as cokes and char etc. For the comparison, the single phase flow and the two-phase flow without the porosity have been simulated. The two-phase flow without porosity condition revealed the smooth parabolic profile of the free surface near the outlet. However, the free surface under the porosity condition formed the viscous finger when the free surface was close to the outlet. This viscous finger accelerated the velocity of the free surface falling and the outflow velocity of the fluids near the outlet.

• Journal of computational fluids engineering
• /
• v.17 no.4
• /
• pp.75-80
• /
• 2012

Due to the effect of surface tension, liquid film around a curved edge of solid surface moves from the corner to the flat surface. During this behavior of liquid film, film sagging phenomenon is easily occurred at the solid surface. Behavior of liquid film is determined with the effects of the properties of liquid film and the geometric factors of solid surface. In the present study, 2-D transient CFD simulations were conducted on the behavior of liquid film around a curved edge. The two-phase interfacial flow of liquid film was numerically investigated by using a VOF method in order to predict the film sagging around a curved edge. In the steady state of behavior of liquid film, the liquid film thickness of numerical result showed a good agreement with experimental data. After verifying the numerical results, the characteristics of behavior of liquid film were numerically analyzed with various properties of liquid film such as surface tension coefficient and viscosity. The effects of geometric factors on film sagging were also investigated to reduce the film sagging around a curved edge.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.7
• /
• pp.641-647
• /
• 2016

A drop-on-demand inkjet is used widely for various applications. Therefore, it is important to understand the jetting behavior of the drop from the piezo inkjet. In this study, to predict the jetting behavior, VOF (Volume-of-Fluid) simulation techniques were used and compared with the experimental results. The experimentally measured meniscus movement was used as the input data for the simulation. To verify the simulation, the measured jetting behavior of the mixture fluids of ethylene glycol and IPA (isopropyl alcohol), which has a mixing ratio of 50:50, was used. The numerical simulation of the drop formation using various mixture ratios and its comparison with the measured drop formation confirmed that the proposed method can predict the actual jetting. On the other hand, the satellite drop behavior showed slight differences because the small sized droplet is subject to a more aerodynamic effect during flight because the kinetic energy of the satellite droplet is far smaller than that of the main droplet.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.9
• /
• pp.571-577
• /
• 2017

In this numerical study, the separation efficiency of three-phase separator in an oil-sand plant was studied with various inlet head configurations. The free water knockout (FWKO) vessel was designed with a flow rate of $15.89m^3/day$ (100 bbl/day) and the SOR(stream-to-oil ratio)=3.5 was derived using Stokes' law. For modularization, optimization of the design of the inlet head configuration was performed with parallel-connected dual FWKO vessels. The feed condition of bitumen emulsion was API=17, $T_{in}=150^{\circ}C$ and $P_{in}=50bar$. A mean residence time was determined the time when 95% of the oil and water in FWKO vessel was separated. The combination between the volume of fluid (VOF) and the discrete phase model (DPM) was used to simulate the phase separation phenomenon in a multi-phase separator. Furthermore, in order to calculate multi-phase flow the pseudo-transient method was adopted.

• Journal of Korea Water Resources Association
• /
• v.36 no.5
• /
• pp.741-749
• /
• 2003

This study presents a combined experimental and numerical effort to investigate experimentally and numerically the Bragg reflection of sinusoidal waves due to trapezoidal submerged porous breakwaters. Numerical predictions of the study are verified by comparing to laboratory measurements. In the numerical model, the flow in porous structures is described by the spatially averaged Navier-Stokes equations and the volume of fluid method is employed to track the free surface displacements. Numerical solutions are agree well with laboratory measurements. The reflection coefficients of porous structures are smaller than those of non-porous structures and become stronger in proportion to the increase of number of submerged breakwaters.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.13 no.12
• /
• pp.5661-5668
• /
• 2012

In the era of water scarcity, saving toilet water is one of the most effective ways to save water. In this study, two-phase flow for the development of highly-efficient toilet has been analysed with VOF(Volume Of Fluid) method. Since the whole model requires extensive computing time, part of the whole model has also been adopted to reduce the analysis time. Four different rim models were considered for flow distribution analysis and sizes of rim holes were found to become more important rather than locations or numbers of rim holes in achieving flow distribution effectively. In addition, velocity and pressure of two-phase flow due to siphon phenomenon have been studied through the analysis of whole model. Therefore, this study provides a variety of fundamental data for the development of highly-efficient toilet.