• 유체기계공업학회:학술대회논문집
• /
• 1998.12a
• /
• pp.223-232
• /
• 1998

To investigate the flow inside the centrifugal impeller, computer program which can solve Three-dimensional compressible turbulent flow has been developed. The Navier-Stokes equations were chosen as the governing equation for viscous flow while Euler equations for inviscid case. Time marching method was incorporated with the Flux Difference Splitting method suggested by Roe to capture the steep gradients such as a shock. For high order of accuracy, MUSCL approach was adopted while differentiable limiter to ensure TVD property. For turbulence closure, Baldwin- Lomax model was applied due to its simplicity. To demonstrate the capabilities of present program, several validation problems have been solved and compared with experiments and other available data. From the above calculations generally good agreements were obtained. Finally, the developed code was applied to Eckardt's impeller and the performance prediction was carried out. Some important aspects on boundary condition for successful simulation were discussed and the remedy was also introduced.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.151-155
• /
• 2005

Three-dimensional structures of unsteady detonation wave propagating through a square-shaped tube is studied using computational method and parallel processing. Inviscid fluid dynamics equations coupled with variable-${\gamma}$ formulation and simplified one-step Arrhenius chemical reaction model were analysed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Results in three dimension show the two unsteady detonation wave propagating mode, the Rectangular and diagonal mode of detonation wave instabilities. Two different modes of instability showed the same cell length but different cell width and the geometric similarities in smoked-foil record.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.169-177
• /
• 2005

A high resolution scheme for solving gas-liquid two-phase flows with cavitation is described. This scheme uses the curvilinear coordinate grid and solves the density based momentum equations for mixture of gas-liquid medium with a preconditioning method to treat both compressible and incompressible flow characteristics. The present preconditioned method is based on the Runge-Kutta explicit finite-difference scheme, and is improved by using the diagonalization, the flux difference splitting and the MUSCL-TVD schemes to save computational effort and to increase stability and resolvability, especially at gas-liquid contact surfaces. A homogeneous equilibrium cavitation model is used to treat the gas-liquid two-phase medium in cavitating flow as a locally homogeneous pseudo-single-phase medium. Therefore, it is easy to solve cavitating flow, including wave propagation, large density changes and incompressible flow characteristic at low Mach number. Some numerical results obtained by the present scheme are shown.

• Journal of Mechanical Science and Technology
• /
• v.15 no.6
• /
• pp.788-795
• /
• 2001

The effects of nonequilibrium condensation on the shock boundary layer interaction over a transonic bump model were investigated experimentally and numerically. An experiment was conducted using a supersonic indraft wind tunnel. A droplet growth equation was incorporated into two-dimensional Navier-Stokes equation systems. Computations were carried out using a third-order MUSCL type TVD finite-difference scheme with a second-order fractional time step. Computation compared with the experimental results. Nonequilibirum condensation suppressed the boundary layer separation and the pressure fluctuations due to the shock boundary layer interaction. Especially the nonequilbrium condensation was helpful to suppress the high frequency components of the pressure fluctuations.

• Journal of Powder Materials
• /
• v.2 no.2
• /
• pp.120-134
• /
• 1995

The behavior of the flow about gas atomizers with a supersonic nozzle containing an under-expanded or over-expanded jet is very important with respect to performance and stability characteristics. Since detailed experiments are expensive, computational fluid mechanics have been applied recently to various relating flow field. In this study, a higher order upwind method with the 3rd order MUSCL type TVD scheme is used to solve the full Reynolds Wavier-Stokes equations. To delineate the purely exhaust jet effects, the melt flow is not considered. Comparison is made with some experimental data in terms of density fields. The influence of the exhaust-jet-to freestream pressure ratio and the effect of the protrusion length of the melt orifice are studied. The present study leads us to believe that the computational fluid mechanics should be considered as powerful tool in predicting the gas atomizer flows.

• 한국연소학회:학술대회논문집
• /
• 2005.10a
• /
• pp.15-19
• /
• 2005

Three-dimensional structures of unsteady detonation wave propagating through a square-shaped tube is studied using computational method and parallel processing. Inviscid fluid dynamics equations coupled with variable-${\gamma}$ formulation and simplified one-step Arrhenius chemical reaction model were analysed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Results in three dimension show the two unsteady detonation wave propagating mode, the Rectangular and diagonal mode of detonation wave instabilities. Two different modes of instability showed the same cell length but different cell width and the geometric similarities in smoked-foil record.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.3
• /
• pp.248-253
• /
• 2009

This paper presents a numerical study for unsteady flows in a high-pressure steam turbine with a partial admission stage. Compressible Navier-Stokes equations are solved by the high-order high-resolution finite-difference method based on the fourth-order compact MUSCL TVD scheme, Roe's approximate Riemann solver, and the LU-SGS scheme. The SST-model is also solved for evaluating the eddy-viscosity. The unsteady two-dimensional flows through whole nozzle-rotor cascade channels considering a partial admission are numerically investigated. 108 nozzle passages with two blockages and 60 rotor passages are simultaneously calculated. The influence of the flange in the nozzle box to the lift of rotors is predicted. Also the efficiency of the partial admission stage changing the number of blockages and the number of nozzles is parametrically predicted.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.184-186
• /
• 2003

In order to reduce numerical dissipation in a multi-dimensional flow, a new variable interpolation method for Roe's FDS is proposed. By introducing the Mach number weighting function w, the properties at the cell-interface are interpolated and in a non flow-aligned grid system, it can give more accurate solution because of less numerical dissipation. Various test cases including oblique contact discontinuity are simulated and compared with the results of original Roe's FDS.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.404-411
• /
• 2010

The present paper deals with the continuous work of extending multi-dimensional limiting process (MLP), which has been quite successfully proposed on two- and three-dimensional structured grids, onto the unstructured grids. The basic idea of the present limiting strategy is to control the distribution of both cell-centered and cell-vertex physical properties to mimic a multi-dimensional nature of flow physics, which can be formulated as so called the MLP condition. The MLP condition can guarantee a high-order spatial accuracy without yielding spurious oscillations. Recently, MLP slope limiter was proposed based on the MUSCL-type reconstruction in two-dimensional case and it can be readily extended to three-dimensional case. Through various numerical analyses and extensive computations, it is observed that the proposed limiters are quite effective in controlling numerical oscillations and very accurate in capturing both discontinuous and continuous multi-dimensional flow features on 3-D tetrahedral grids.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.544-549
• /
• 2000

The effects of nonequilibrium condensation on the shock boundary layer interaction over a transonic bump model were investigated experimentally and numerically. An experiment was conducted using a supersonic indraft wind tunnel. A droplet growth equation was incorporated into two-dimensional Navier-Stokes equation systems. Computations were carried out using a third-order MUSCL type TVD finite-difference scheme with a second-order fractional time step. Computations compared with the experimental results. Nonequilibirum condensation suppressed the boundary layer separation and the pressure fluctuations due to the shock boundary layer interaction. Especially the nonequilibrium condensation was helpful to suppress the high frequency components of the pressure fluctuations.