• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.1
• /
• pp.1-13
• /
• 2008

In the present work, wave transformation by vertical columns and its wave forces acting on them are discussed using a direct 3-D numerical model based on the VOF (Volume Of Fluid) method. The numerical results for wave transformations and wave forces are critically compared to an advanced experimental data, and provide the verification of the numerical model used in the present study. Overall model-data comparisons are good. After verification of the numerical model, it is used to simulate wave fields around dual vertical columns with arbitrary cross section, and the characteristics of nonlinear wave forces and wave transformations according to the variations of different cross section types of vertical columns, an interval of vertical columns and incident wave angle are discussed.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.21 no.1
• /
• pp.72-83
• /
• 2017

A thermo-fluid dynamic analysis was performed through the numerical simulation of a missile canister. Calculation was made in a fixed analytical volume and fully evaporated water was used as a coolant. To analyze the interaction among the hot gas, coolant, and mixture flow, Realizable $k-{\varepsilon}$ turbulence and VOF(Volume Of Fluid) model were chosen and parametric study was performed with the change of coolant flow rate. It could be found that the pressure on the canister top nonlinearly increased with the increase of coolant flow rate. Temperature and coolant distribution were closely related to the flow behavior in canister. Temperature on the canister bottom indicated a decrease being proportional to coolant flow rate in early times but after a specific time, the temperature increased with the tendency being reversed. In addition, the early part of temperature showed a fluctuating phenomenon because of the overall circulatory flow of mixture gas.

• Nuclear Engineering and Technology
• /
• v.53 no.4
• /
• pp.1119-1126
• /
• 2021

Spacer grid with mixing vane had been widely used in nuclear reactor core. One of the main feather of spacer grid with mixing vane was that strong swirl flow was formed after the spacer grid. The swirl flow not only changed the bubble generation in the near wall field, but also affected the bubble behaviors in the center region of the subchannel. The interaction between bubble and the swirl flow was one of the basic phenomena for the two phase flow modeling in fuel assembly. To obatin better understanding on the bubble behaviors in swirl flow, full three dimension numerical simulations were conducted in the present paper. The swirl flow was assumed in the cylindral calculation domain. The bubble interface was captured by Volume Of Fluid (VOF) method. The properties of saturated water and steam at different pressure were applied in the simulation. The bubble trajectory, motion, shape and force were obtained based on the bubble parameters captured by VOF. The simulation cases in the present study included single bubble with different size, at different angular velocity conditions and at different pressure conditions. The results indicated that bubble migrated to the center in swirl flow with spiral motion type. The lateral migration was mainly related to shear stress magnitude and bubble size. The bubble moved toward the center with high velocity when the swirl magnitude was high. The largest bubble had the highest lateral migration velocity in the present study range. The effect of pressure was small when bubble size was the same. The prelimenery simulation result would be beneficial for better understanding complex two phase flow phenomena in fuel assembly with spacer grid.

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

• Journal of Korea Water Resources Association
• /
• v.50 no.2
• /
• pp.89-97
• /
• 2017

Turbulent flow structure in the high amplitude meandering channel is complex due to secondary recirculation with helicoidal motions and shear layers formed by flow separation from the curved sidewall. In this work, the secondary flow and the superelevation of the water surface produced in the high-amplitude Kinoshita channel are reproduced by the unsteady Reynolds-averaged Navier-Stokes (RANS) computations using the VOF technique for resolving the variation of water surface elevation and three statistical turbulence models ($k-{\varepsilon}$, RNG $k-{\varepsilon}$, $k-{\omega}$ SST). The numerical results computed by a second-order accurate finite volume method are compared with an existing experimental measurement. Among applied turbulence models, $k-{\omega}$ SST model relatively well predicts overall distribution of the secondary recirculation in the Kinoshita channel, while all three models yield similar prediction of water superelevation transverse slope. The secondary recirculation driven by the radial acceleration in the upstream bend affects the flow structure in the downstream bend, which yields a pair of counter-rotating vortices at the bend apex. This complex flow pattern is reasonably well reproduced by the $k-{\omega}$ SST model. Both $k-{\varepsilon}$ based models fail to predict the clockwise-rotating vortex between a pair of counter-rotating vortices which was observed in the experiment. Regardless of applied turbulence models, the present computations using the VOF method appear to well reproduce the superelevation of water surface through the meandering channel.

• Proceedings of the KSME Conference
• /
• 2005.11a
• /
• pp.84.1-84.1
• /
• 2005

• Wind and Structures
• /
• v.14 no.5
• /
• pp.435-447
• /
• 2011

A Computational Fluid Dynamics model is presented in this study for the simulation of the complex fluid flows with free surfaces inside the Tuned Liquid Column Dampers in horizontal motion. The characteristics of the fluid model of the TLCD in horizontal motion include the free surface of the multiphase flow and the horizontal moving frame. In this study, the time depend unsteady Standard ${\kappa}-{\varepsilon}$ turbulent model based on Navier-Stokes equations is chosen. The volume of fluid (VOF) method and sliding mesh technique are adopted to track the free surface of water inside the vertical columns of TLCD and treat the moving boundary of the walls of TLCD in horizontal motion. Several model solution parameters comprising different time steps, mesh sizes, convergence criteria and discretization schemes are examined to establish model parametric independency results. The simulation results are compared with the experimental data in the dimensionless amplitude of the water column in four different configured groups of TLCDs with four different orifice areas. The predicted natural frequencies and the head loss coefficient of TLCDs from CFD model are also compared with the experimental data. The predicted numerical results agree well with the available experimental data.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.215-222
• /
• 2009

In the present study, a numerical analysis on the sloshing in a tank with the harmonic motion was investigated. A VOF method was used for two-phase flows inside the sloshing tank and a source term of the momentum equation was applied for the harmonic motion. This numerical method was verified by comparing its results with the available experimental data. The sloshing in a tank causes the instability of the fluid flows and the fluctuation of the impact pressure on the tank. By these phenomena of the tank sloshing, the sloshing problems such as the failure and the noise of system can be generated. For the reduction of these sloshing problems, the various baffles such as the horizontal/vertical plate baffles and the porous baffles inside the tank are installed. With the installations of these baffles, the characteristics of the liquid behavior in the sloshing tank, the impact pressure on the wall, the amplitude of the free surface near the wall and the sloshing noise were numerically analyzed.

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