• Journal of Advanced Marine Engineering and Technology
• /
• v.18 no.3
• /
• pp.94-101
• /
• 1994

To develope the new simulator for the analysis of fluid flow information, the influence of various convective difference schemes were evaluated. General curvilinear coordinate system with nonorthogonal grids was adopted for the successful analysis of various complex geometries. Computation results show that if we can not obtain full grid numbers within available computational environment, we need to use higher order finite difference schemes to keep the prediction accuracy.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.3
• /
• pp.87-99
• /
• 1994

Reynolds Averaged Navier-Stokes equations are solved numerically for the computation of turbulent flow around a Wigley double model. A second order finite difference method is applied for the spatial discretization on the nonstaggered grid system and 4-stage Runge-Kutta scheme for the numerical integration in time. In order to increase the time step, residual averaging scheme of Jameson is adopted. Pressure field is obtained by solving the pressure-Poisson equation with the appropriate Neumann boundary condition. For the turbulence closure, 0-equation turbulence model of Baldwin-Lomax is used. Numerical computation is carried out for the Reynolds number of 4.5 million. Comparisons of the computed results with the available experimental data show good agreements for the velocity and pressure distributions.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05b
• /
• pp.173-178
• /
• 1996

A three-dimensional thermo-hydraulic computer code is developed for simulation of incompressible flows in complex geometries. The computer code employs a body-fitted, nonorthogonal grid system in order to efficiently handle the complex geometries encountered in many engineering applications. The finite volume method is used to discretize the governing equations and the convection term is treated by higher-order bounded schemes. The cell-centered, nonstaggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is avoided by use of momentum interpolation practice. The computer code employs the SIMPLE algorithm for pressure and velocity coupling and the k-$\varepsilon$ turbulence for turbulent calculation. The computer code has been tested through application to a variety of test problems and some results are presented in this paper

• Journal of Power System Engineering
• /
• v.3 no.4
• /
• pp.16-21
• /
• 1999

Numerical simulation on two-dimensional turbine cascade flow has been performed using compressible Navier-Stokes equations. The flow equations are written in a cartesian coordinate system, then mapped into a generalized body-fitted ones. All direction of viscous terms are incoporated and turbulent effects are modeled using the extended ${\kappa}-{\epsilon}$ model. Equations are discretized using control volume SIMPLE algorithm on the nonstaggered grid sysetem. Applications are made at a VKI turbine cascade flow in atransonic wind-tunnel and compared to experimental data. Present numerical results are shown to be in good agreement with the experimental results and simulate the compressible viscous flow characteristics inside the turbine blade passage.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.160-166
• /
• 1998

In this paper, incompressible two-dimensional Navier-Stokes equations are numerically solved for the study of steady laminar flow around a body with the wall effect. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. To investigate the wall effect, numerical computations are carried out for the NACA 0012 section at the various blockage ratios. The pressure and skin friction on the foil surface, velocity pronto in its wake and drag coefficient are investigated as functions of the blockage ratio.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.175-181
• /
• 1995

Three-Dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for the various contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreements.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.8-12
• /
• 1997

The three-dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for three contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreement.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.1
• /
• pp.1-8
• /
• 2001

The computations of the turbulent flow around the ship models with the free-surface effects were carried out. Incompressible Reynolds-Averaged Navier-Stokes equations were solved by using an explicit finite-difference method with the nonstaggered grid system. The method employed second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration. For the turbulence closure, a modified Baldwin-Lomax model was exploited. The location of the free surface was determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and a free-surface conforming grid was generated at each time step so that one of the grid boundary surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition was applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method developed in the present study, the computations were carried out for beth Wigley and Series 60 $C_B=0.6$ ship model and the computational results showed good agreements with the experimental data.