• Journal of Korea Water Resources Association
• /
• v.57 no.8
• /
• pp.493-507
• /
• 2024

RANS-based CFD analysis is widely applied in various engineering fields, including practical hydraulic engineering, due to its high computational efficiency. However, problems of non-physical behavior in the analysis of two phase flow, such as free surfaces, have long been raised. The two-equation turbulence models used in general RANS-based analysis were developed for single phase flow and simulate unrealistically high turbulence energy at the interface where there are abrupt changes in fluid density. To solve this issue, one of the methods recently developed is the buoyancy-modified turbulence model, which has been partially validated in coastal engineering, but has not been applied to open channel flows. In this study, the applicability of the buoyancy-modified turbulence model is evaluated using the VOF method in the open-source program OpenFoam. The results of the uniform flow showed that both the buoyancy-modified k-𝜖 model and the buoyancy-modified k-ω SST model effectively simulated the reduction of turbulence energy near the free surface. Specifically, the buoyancy-modified k-ω SST model accurately simulated the vertical velocity distribution. Additionally, the model is applied to dam-break flows to examine cases with significant surface variation and cavity formation. The simulation results show that the buoyancy-modified turbulence models produce varying results depending on the VOF method and shows non-physical behavior different from experimental results. While the buoyancy-modified turbulence model is applicable in cases with stable surface shapes, it still has limitations in general application when there are rapid changes in the free surface. It is concluded that appropriate adjustments to the turbulence model are necessary for flows with rapid surface changes or cavity formation.

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

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.15 no.2
• /
• pp.111-117
• /
• 2012

Oscillating water column is the most widely used ocean energy converting systems all over the world. The operating performance is influenced by the efficiencies of the two converting stages in the OWC chamber-turbine integrated system. In order to consider the effects of the turbine, the orifice model are carried out. The VOF based Numerical Wave Tank (NWT) is utilized to simulate the water column oscillation inside the chamber and the results are compared with corresponding experimental data. This paper reviews the state of the art in interaction among wave elevation inside the chamber and air flow rate in the duct, which are considered the turbine effects. Effects of incident wave conditions and several shape parameters on the operating performance of OWC chamber are investigated numerically. The effects of the impulse turbine on the integrated system and interaction among the wave elevation, pressure and air flow velocities variations are investigated.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.3
• /
• pp.258-263
• /
• 2008

Performance of injector equiped with screw type swirler which is suitable for supersonic cruise engine combustor was investigated using theoretical, numerical, and experimental methods. Based on discharge coefficient and spray angle which represent the performance of injectors, the geometrical parameters which affect these performance parameters were defined, control group injectors were designed, and variation of performance parameters according to the geometrical parameters were examined. Within the defined range, measured value of performance of injectors was smaller than result of theoretical prediction, and prediction result from numerical simulation using VOF method agreed with the result of experiments very well. The viscous barrier was not observed, and minimum discharge coefficient and maximum spray angle, 0.05 and 104 respectively, was obtained for this type of injector.

• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.73-80
• /
• 2004

A LS(Level Set) formulation is developed for computing two-phase flows on non- orthogonal meshes. Compared with the VOF(Volume-of-Fluid) method based on a non-smooth volume-fraction function, the LS method can calculate an interfacial curvature more accurately by using a smooth distance function. Also, it is quite straightforward to implement for 3-D irregular meshes compared with the VOF method requiring much more complicated geometric calculations. The LS formulation is implemented into a general purpose program for 3-D flows and verified through several test problems.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.214-219
• /
• 2007

Both the bubble rising in a fully filled container and the droplet splash are simulated by a solution code(PowerCFD). This code employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with interface capturing method (CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The present results are compared with other numerical solutions found in the literature. It is found that the present code simulate complex free surface flows such as multi phase flows due to large density difference efficiently and accurately.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.9-13
• /
• 2006

In the present study, a numerical analysis of transient mixer flow is performed considering free surface formation. The flow patterns and free surface shape in a mixers formed by flat paddle and pitched paddle impellers are predicted. In a flat paddle mixer, two flow circulation regions are formed due to strong radial flow, whereas one large circulation is formed in a pitched paddle mixer due to axial downward flow. These differences affect the free surface evolution and shape. It is seen from the results that a flat paddle mixer gives deeper free surface at center region than a pitched paddle mixer. The free surface of 8-blades-flat-paddle mixer is also simulated to compare with the available experimental and simulation results. The present computational results agree reasonably well with the experimental data.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.95-98
• /
• 2007

A general purpose program NUFLEX has been extended for two-phase flows with topologically complex interface and cavitation flows with liquid-vapor phase change caused by large pressure drop. In analysis of two-phase flow, the phase interfaces are tracked by employing a LS(Level Set) method. Compared with the VOF(Volume-of-Fluid} method based on a non-smooth volume-fraction function, the LS method can calculate an interfacial curvature more accurately by using a smooth distance function. Also, it is quite straightforward to implement for 3-D irregular meshes compared with the VOF method requiring much more complicated geometric calculations. Also, the cavitation process is computed by including the effects of evaporation and condensation for bubble formation and collapse as well as turbulence in flows. The volume-faction and continuity equations are adapted for cavitation models with phase change. The LS and cavitation formulation are implemented into a general purpose program for 3-D flows and verified through several test problems.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.910-915
• /
• 2003

A Navier-Stokes equation solver for incompressible viscous flows with free surface is developed and tested. This is based upon a fractional time step method and a non-staggered finite volume formulation for unstructured meshes. For time advancement scheme, Adams -Bashforth method for convective term and Crank-Nicolson method for diffusive term are applied. The interface between two fluids with different fluid properties is tracked with Piecewise Linear Interface Calculation(PLIC) Volume-of-Fluid(VOF) methods. Computational results are presented for some test problems: the broken dam, the sloshing in a rectangular tank, the filling of a cylindrical tank.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.5
• /
• pp.481-486
• /
• 2011

It is important to predict the behavior of a liquid film around a rotating cylinder in the film coating process of the steel industry. When the cylinder rotates, the behavior of the liquid film on the rotating cylinder surface is influenced by the cylinder diameter, the rotation speed, the gravitational force, and the fluid properties. These parameters determine the liquid film thickness and the rise of the film on the cylinder surface. In the present study, the two-phase interfacial flow of the liquid film on the rotating cylinder were numerically investigated by using a VOF method. For various rotation speeds, cylinder diameters and fluid viscosities, the behavior of liquid film on the rotating cylinder were predicted. Thicker film around the rotating cylinder was observed with an increase in the rotation speed, cylinder diameter, and fluid viscosity. The present results for the film thickness agreed well with available experimental and analytical results.