• New & Renewable Energy
• /
• v.3 no.4
• /
• pp.8-15
• /
• 2007

Liquid water in flow channel is an important factor that limits the steady and transient performance of PEM fuel cells. A computational fluid dynamics study based on the volume-of-fluid [VOF] multi-phase model was conducted to understand the two-phase flow behavior of liquid water in cathode gas channels. The liquid water transport in $180^{\circ}{\Delta}$ bends was investigated, where the effects of surface characteristics (hydrophilic and hydrophobic surfaces], channel geometries (rectangular and chamfered corners], and air velocity in channel were discussed. The two-phase flow behavior of liquid water with hydrophilic channel surface and that with hydrophobic surface was found very different; liquid water preferentially flows along the corners of flow channel in hydrophilic channels while it flows in rather spherical shape in hydrophobic channels. The results showed that liquid water transport was generally enhanced when hydrophobic channel with rounded corners was used. However, the surface characteristics and channel geometries became less important when air velocity was increased over 10m/s. This study is believed to provide a useful guideline for design optimization of flow patterns or channel configurations of PEM fuel cells.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.20 no.3
• /
• pp.312-317
• /
• 2014

In this paper, the wave run-up height and depression depth around air-water interface-piercing circular cylinder have been numerically studied. The Reynolds Averaged Navier-Stokes equations (RANS) and continuity equations are solved with Reynolds Stress model (RSM) and volume of fluid (VOF) method as turbulence model and free surface modeling, respectively. A commercial Computational Fluid Dynamics (CFD) software "Star-CCM+" has been used for the current simulations. Various Froude numbers ranged from 0.2 to 1.6 are used to investigate the change of air-water interface structures around the cylinder and experimental data and theoretical values by Bernoulli are compared. The present results showed a good agreement with other studies. Kelvin waves behind the cylinder were generated and its wave lengths are longer as Froude numbers increase and they have good agreement with theoretical values. And its angles are smaller with the increase of Froude numbers.

• Journal of computational fluids engineering
• /
• v.17 no.2
• /
• pp.35-41
• /
• 2012

This paper presents incompressible Navier-Stokes solution algorithm for 2D Free-surface flow problems on the Cartesian mesh, which was implemented to run on Graphics Processing Units(GPU). The INS solver utilizes the variable arrangement on the Cartesian mesh, Finite Volume discretization along Constrained Interpolation Profile-Conservative Semi-Lagrangian(CIP-CSL). Solution procedure of incompressible Navier-Stokes equations for free-surface flow takes considerable amount of computation time and memory space even in modern multi-core computing architecture based on Central Processing Units(CPUs). By the recent development of computer architecture technology, Graphics Processing Unit(GPU)'s scientific computing performance outperforms that of CPU's. This paper focus on the utilization of GPU's high performance computing capability, and presents an efficient solution algorithm for free surface flow simulation. The performance of the GPU implementations with double precision accuracy is compared to that of the CPU code using an representative free-surface flow problem, namely. dam-break problem.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.1
• /
• pp.38-46
• /
• 2010

This study aims at investigating the effect of the baffle height on the liquid sloshing in the three-dimensional (3D) rectangular tank. In order to simulate the 3D incompressible viscous two-phase flow in the 3D tank with partially filled liquid, the present study has adopted the volume of fluid (VOF) method based on the finite-volume method which has been well verified by comparing with the results of the relevant previous researches. The ratio of the baffle height ($h_B$) to filling level (h) has been changed in the range of $0{\leq}h_B/h{\leq}1.2$ to observe the effect on the impact loads on the side wall and free surface behavior. Generally, as baffle height increases, the impact pressure on the wall decreases and the deformation of free surface becomes weaker. However it seemed that a critical ratio of the baffle height existed to reveal the lowest impact pressure on the wall. Consequently, $h_B/h=0.8$ among $h_B/hs$ considered in the study showed the lowest impact pressure.

• Transactions of the Korean hydrogen and new energy society
• /
• v.30 no.6
• /
• pp.585-592
• /
• 2019

In this study, a flashing boiling phenomenon of pressurized water jet was numerically studied and validated against an experimental data in the literatures. The volume of fluid (VOF) technique was used to consider two-phase behavior of water, while the homogeneous relaxation model (HRM) model was used to provide the velocity of phase change. During the flashing boiling through a nozzle, a mach disk was observed near nozzle exit because of pressure drop resulting from two-phase under-expansion. The flashing jet structure, local distributions of temperature/vapor volume fraction/velocity, and position of the mach disk were examined as nozzle inlet temperature changed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.10 s.253
• /
• pp.1003-1011
• /
• 2006

Droplets are ejected onto a substrate through a nozzle by pushing liquids in flow channels of drop-on-demand devices. The behavior of ejection and formation of droplets is investigated to enhance the physical understanding of the hydrodynamics involved in inkjet printing. The free surface phenomenon of a droplet is described using $CFD-ACE^{TM}$ which employs the volume-of-fluid (VOF) method with the piecewise linear interface construction (PLIC). Droplet formation characteristics are analyzed in various flow regimes with different Ohnesorge numbers. The computational results show that the droplet formations are strongly dependent on the physical properties of working fluids and the inlet flow conditions. In addition, the wetting characteristics of working fluids on a nozzle influence the volume and velocity of a droplet produced in the device. This study may provide an insight into how a liquid droplet is formed and ejected in a piezoelectric inkjet printing device.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.4
• /
• pp.401-412
• /
• 2008

Simulations of three-dimensional numerical wave tank are performed to investigate wave force acting on a large cylindrical structure and consequent wave deformation, which are induced by bore after breaking waves. The numerical model is based on the three-dimensional Navier-Stokes equations with a finite-difference method combined with a volume of fluid(VOF) method, which is capable of tracking the complex free surface, including wave breaking. In order to promote wave breaking of the incident wave, the approach slope was built seaward of the structure with a constant slope and a large cylindrical structure was installed on a flat bed. The incident waves were broken on the approach slope or flat bed by its wave height. In the present study, all waves acting on the large cylindrical structure were limited to breaking bore after wave breaking. The effects of the position of the structure and the incident wave height on the wave force and wave transformations were mainly investigated with the concern of wave breaking. Further, the relations between the variation of wave energy by wave propagation after wave breaking and wave force acting on the structure were discussed to give the understanding of the full-linear wave-structure interactions in three-dimensional wave fields.

• The Journal of the Acoustical Society of Korea
• /
• v.42 no.3
• /
• pp.250-261
• /
• 2023

In this study, the noise contributions by the duct, stator and rotor, which are the propulsor components, are evaluated to identify the flow noise source in cavitation and non-cavitation conditions on pumpjet propulsion and the noise levels in both conditions are compared. The unsteady incompressible Reynolds averaged Navier-Stokes (RANS) equation based on the homogeneous mixture assumption is applied on the suboff submarine hull and pumpjet propeller in the cavitation tunnel, and the Volume of Fluid (VOF) method and Schnerr-Sauer cavitation model are used to describe the two-phase flow. Based on the flow simulation results, the acoustic analogy formulated by Ffowcs Williams and Hawkings (FW-H) equation is applied to predict the underwater radiated noise. The noise contributions are evaluated by using the three types of impermeable integral surface on the duct, stator and rotor, and the two types of permeable integral surface surrounding the propulsor. As a result of noise prediction, the contribution by the stator is insignificant, but it affects the generation of flow noise source due to flow separation in the duct and rotor, and the noise is predominantly radiated into the upward and right where the flow separations are. Also, the noise is radiated into the thrust direction due to pressure fluctuation between suction and pressure sides on the rotor blades, and the it can be seen that the cavitation effect into the noise can be considered through the permeable integral surface.

• Journal of Korea Water Resources Association
• /
• v.43 no.9
• /
• pp.801-811
• /
• 2010

A three-dimensional numerical model called NEWTANK is employed to investigate solitary wave run-up with an internal wave-maker on a steep slope. The numerical model solves the spatially averaged Navier-Stokes equations for two-phase flows. The LES (large-eddy-simulation) approach is adopted to model the turbulence effect by using the Smagorinsky SGS (sub-grid scale) closure model. A two-step projection method is adopted in numerical solutions, aided by the Bi-CGSTAB (Bi-Conjugate Gradient Stabilized) method to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate VOF (volume-of-fluid) method is used to track the distorted and broken free surface. A solitary wave is first internally generated and propagated over a constant water depth in the three-dimensional domain. Numerically predicted results are compared with analytical solutions and numerical errors are analyzed in detail. The model is then applied to study solitary wave run-up on a steep slope and the obtained results are compared with available laboratory measurements.

• Wind and Structures
• /
• v.38 no.4
• /
• pp.261-275
• /
• 2024

Before performing an experimental study on the downburst-generated wave, it is necessary to examine the scale effects and corresponding corrections or compensations. Analysis of similarity is conducted to conclude the non-dimensional force ratios that account for the dynamic similarity in the interaction of downburst with wave between the prototype and the scale model, along with the corresponding scale factors. The fractional volume of fluid (VOF) method in association with the impinging jet model is employed to explore the characteristics of the downburst-generated wave numerically, and the validity of the proposed scaling method is verified. The study shows that the location of the maximum radial wind velocity in a downburst-wave field is a little higher than that identified in a downburst over the land, which might be attributed to the presence of the wave which changes the roughness of the underlying surface of the downburst. The impinging airflow would generate a concavity in the free surface of the water around the stagnation point of the downburst, with a diameter of about two times the jet diameter (D_jet). The maximum wave height appears at the location of 1.5D_jet from the stagnation point. Reynolds number has an insignificant influence on the scale effects, in accordance with the numerical investigation of the 30 scale models with the Reynolds number varying from 3.85 × 10⁴ to 7.30 × 10⁹. The ratio of the inertial force of air to the gravitational force of water, which is denoted by G, is found to be the most significant factor that would affect the interaction of downburst with wave. For the correction or compensation of the scale effects, fitting curves for the measures of the downburst-wave field (e.g., wind profile, significant wave height), along with the corresponding equations, are presented as a function of the parameter G.