• Journal of computational fluids engineering
• /
• v.15 no.1
• /
• pp.81-87
• /
• 2010

In this study, we conducted resistance and propulsion performance test of ship composed of the Resistance Test, Propeller Open Water Test and Self Propulsion Test using the CFD(Computational Fluid Dynamics). We used commercial RANS(Reynolds Averaged Navier Stokes equation) solver, as a calculating tool. The unstructured grids were used in a bow and stern of ship, having complex shape, for a convenience of generating grids, and the structured grids were adopted in a central hull and rest of hull having a relatively simple shape which is called hybrid grid method. In addition, The sliding mesh method was adopted to rotate a propeller directly in the Propeller Open Water and Self Propulsion Test. The Resistance Test and Self Propulsion Test were calculated using Volume of Fluid (VOF) model and considering a free surface. And all The three cases were applied realizable k-epsilon model as the turbulence model. The results of calculations were verified for the suitability of calculations by comparing MOERI's EFD results.

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

• Journal of Hydrogen and New Energy
• /
• 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.