• 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 the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.8-17
• /
• 2002

In this paper, the computational analysis of open-water characteristics for three model propellers(P4119, P4842 and 3 podded propeller of KRISO) is done by using a viscous-flow method based on Reynolds-Averaged Navier-Stokes equations. The results are presented for open-water performances, blade-section pressures, and circumferentially-averaged velocity profiles for the all three propeller models. Overall close agreements with available experimental data are shown. However, some discrepancies are also found in the pressure near the leading edge of the propeller blade and the open-water performance of the podded propellers.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1971-1978
• /
• 1992

A turbulent flow in a conical diffuser with total divergence angle of 8.deg. was numerically studied. The low Reynolds number k-.epsilon. model(Launder-Sharma model) was adopted to simulate the turbulence. The continuity and time averaged Navier-Stokes equations in a nonorthogonal coordinate system were solved by a finite volume method based on the fully elliptic formulation. The low Reynolds number k-.epsilon. model reasonably simulates the pressure recovery and the mean velocity components. However, there are also considerable discrepancies between predicted and measured shear stress distribution on the wall and turbulent kinetic energy distributions. It is necessary to investigate the flow structure at the entry of the diffuser, numerically as well as experimentally.

• Ocean Systems Engineering
• /
• v.4 no.3
• /
• pp.185-200
• /
• 2014

The objective of the present simulations is to evaluate the applicability of the standard $k-{\varepsilon}$ turbulence model in engineering practice in the subcritical to supercritical flow regimes. Two-dimensional numerical simulations of flow around a circular cylinder at $Re=1{\times}10^5$, $5{\times}10^5$ and $1{\times}10^6$, had been performed using Unsteady Reynolds-Averaged Navier Stokes (URANS) equations with the standard $k-{\varepsilon}$ turbulence model. Solution verification had been studied by evaluating grid and time step size convergence. For each Reynolds number, several meshes with different grid and time step size resolutions were chosen to calculate the hydrodynamic quantities such as the time-averaged drag coefficient, root-mean square value of lift coefficient, Strouhal number, the coefficient of pressure on the downstream point of the cylinder, the separation angle. By comparing the values of these quantities of adjacent grid or time step size resolutions, convergence study has been performed. Solution validation is obtained by comparing the converged results with published numerical and experimental data. The deviations of the values of present simulated quantities from those corresponding experimental data become smaller as Reynolds numbers increases from $1{\times}10^5$ to $1{\times}10^6$. This may show that the standard $k-{\varepsilon}$ model with enhanced wall treatment appears to be applicable for higher Reynolds number turbulence flow.

• Journal of Ship and Ocean Technology
• /
• v.3 no.2
• /
• pp.1-16
• /
• 1999

The incompressible Reynolds-Averaged Navier-Stokes(RANS) equations are solved using the standard $\textsc{k}-\varepsilon$ turbulence model and a finite volume method(FVM)with an O-type grid system. The computed results for its performance test are in good agreement with the published experimental data. The present method is applied to the study on the leading edge radius of a hydrofoil section Calculated results suggest that the leading edge radius has some effects on cavitation performances of a 2-D foil. A natural leading edge radius for the NACA66 section is determined from this study.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.3191-3196
• /
• 2007

The present study addresses a method to operate a fuel-cell system effectively using a recirculation ejector which recycles wasted hydrogen gas. Configuration of a recirculation ejector is changed to investigate the flow behavior through it under varying operating conditions, and how such conditions affect the fuel-cell hydrogen cycle. The numerical simulations are based on a fully implicit finite volume scheme of the axisymmetric, compressible, Reynolds-Averaged, Navier-Stokes equations for hydrogen gas, and are compared with available experimental data for validation. The results show that a hydrogen recirculation ratio is effectively controlled by a configurational alteration within the operational region in which the recirculation passage doesn't plugged by a sonic line.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.49-55
• /
• 2011

Turbulent flows with wavy wall are simulated using Residual-based Variational Multiscale Method (RB-VMS) which is proposed by Bazilves et al(2007) as new Large Eddy Simulation methodology. Incompressible Navier-Stokes equations are integrated using Isogeometric analysis which adopt the basis function as NURBS. The Reynolds number is 6760 based on the bulk velocity and averaged channel height. And the amplitude (${\alpha}/{\lambda}$) of wavy wall is 0.05. The computational domain is $2{\lambda}{\times}1.05{\lambda}{\times}{\lambda}$ in the streamwise, wall normal and spanwise direction. Mean quantities and turbulent statistics near wavy wall are compared with DNS results of Cherukat et al.(1998). The predicted results show good agreement with reference data.

• Journal of computational fluids engineering
• /
• v.12 no.1
• /
• pp.53-59
• /
• 2007

An axisymmetric supersonic screeching jet is numerically simulated to examine the length scales of screech frequency as well as screech tone generation mechanism. The axisymmetric Reynolds-averaged Navier-Stokes equations in conjuction with a modified Spalart-Allmaras turbulence model are employed. It is demonstrated that the axisymmetric jet screech tones can be simulated correctly and the numerical results are in good agreement with the experimental data. Instability waves, shock-cell structures and the phenomena of shock motion are investigated in detail to identify the screech tone generation mechanism. Shock spacings and standing wave length are analyzed to determine the dominent length scale crucial to the screech frequency formulation.

• Journal of Industrial Technology
• /
• v.17
• /
• pp.183-189
• /
• 1997

The computation of the flow around a supercritical airfoil with a divergent trailing edge(DTE) modification(DLBA 243) is compared to that of original supercritical airfoil(DLBA 186). For this computation, Reynolds-Averaged Navier-Stokes equations are solved with a linearized block implicit ADI method and a mixing length turbulence model. Results show the effects of the shock and separated flow regions on drag reduction due to DTE modification. Results also show that DTE modification accelerates the boundary layer flow near the trailing edges which has an effect similar to a chordwise extension that increases circulation and is consistent with the calculated increase in the recirculation region in the wake. Airfoil with DTE modification achieves the same lift coefficient at a lower incidence and thus at a lower drag coefficient, so that lift-to-drag ratio is increased in transonic cruise conditions compared to the original airfoil. The reduction in drag due to DTE modification is associated with weakening of shock strength and delay of shock which is greater than the increase in base drag.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.3
• /
• pp.183-188
• /
• 2017

This paper provides numerical results for the estimation of the efficiency of KRISO energy saving devices in the design stage. A finite volume method is used to solve Reynolds averaged Navier-Stokes (RANS) equations, where the SST k-$\omega$ model is selected for turbulence closure. The propeller rotating motion is determined using a rigid body motion (RBM) scheme, which is called a sliding mesh technique. The numerical analysis focuses on predicting the power reduction by the designed KRISO devices (K-DUCT) under a self-propulsion condition. The present numerical results show good agreement with the available experimental data. Finally, it is concluded that CFD can be a useful method, along with model tests, for assessing the performance of energy saving devices for propulsion efficiency improvement.