• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.215-220
• /
• 2008

A numerical method for gas-liquid two-phase flow is applied to solve shock-bubble interaction problems. The present method employs a finite-difference Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. By this method, a Riemann problem for shock tube was computed for validation. Then, shock-bubble interaction problems between cylindrical bubbles located in the liquid and incident liquid shock wave are computed.

• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.69-78
• /
• 2008

The high-order numerical method based on the adaptive mesh refinement(AMR) on the quadrilateral unstructured grids has been developed in this paper. This adaptive-grid method, originally developed with MUSCL-TVD scheme, is now extended to the WENO (weighted essentially no-oscillatory) scheme with the Runge-Kutta time integration of fifth order in spatial and temporal accuracy. The multidimensional interpolation was studied in the preliminary research, which allows us to maintain the same order of accuracy for the computation of numerical flux between two adjacent cells of different levels. Some standard benchmark tests are done to validate this method for checking the overall capacity and efficiency of the present adaptive-grid technique.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.13 no.3
• /
• pp.27-33
• /
• 2009

A computational study has been performed out to evaluate the effect of a vertical porous barrier on the pressure oscillations in a supersonic cavity. The porous barriers with different perforations were vertically installed into a rectangular cavity at Mach numbers 1.50, 1.83 and 2.50. TVD finite difference MUSCL scheme was employed to solve the two-dimensional, unsteady, compressible Navier-Stokes equations. The present vertical porous barrier considerably altered the characteristics of the time-dependent shear layers that occur at the upstream edge of cavity and remarkably reduced the pressure oscillations inside the supersonic cavity. The present results showed that the effectiveness of passive control using the present porous vertical barrier is dependent on Mach number and the perforation of the porous barrier.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.13 no.1
• /
• pp.27-33
• /
• 2009

Shock tube flow measurement has been often hampered a finite opening time of diaphragm, but there is no systematic work to investigate its effect on the shock tube flow. In the present study, both the experimental and computational works have been performed on the shock tube flows at low pressure ratios. The computational analysis has been performed using the two-dimensional, unsteady, compressible Navier-Stokes equations, based upon a TVD MUSCL finite difference scheme. It is known that the present computational results reproduce the experimental data with good accuracy and simulate successfully the process of diaphragm opening as a function of time. The concept of an imaginary center is introduced to quantify the non-centered expansion wave due to a finite opening time of diaphragm. The results obtained show that the diaphragm opening time is reduced as the initial pressure ratio of shock tube increases, leading to the effect of a finite opening time of diaphragm to be more remarkable at low pressure ratios.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.8
• /
• pp.21-29
• /
• 2002

In the present study, a passive control method, using a porous wall and cavity system, is applied to the shock wave/boundary layer interactions in transonic moist air flow. The two-dimensional, unsteady, compressible, Navier-Stokes equations, which are fully coupled with a droplet growth equation, are solved by the third-order MUSCL type TVD finite difference scheme. Baldwin-Lomax model is employed to close the governing equations. In order to investigate the effectiveness of the present control method, the total pressure loss of the flow and the time-dependent behaviour of shock motions are analyzed in detail. The computed results show that the present passive control method considerably reduces the total pressure losses due to the shock wave/boundary layer interaction in transonic moist air flow and suppresses the unsteady shock wave motions over the airfoil as well. It is also found that the location of the porous ventilation significantly affects the control effectiveness.

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

• 한국전산유체공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.128-134
• /
• 2005

In the CFD area, the numerical analysis of high Mach number flow over a blunt-body poses many issues. Various numerical schemes have been developed to cover the issues, but the traditional schemes are still used widely due to the complexities of new schemes and intricacy of modifying the established codes. In the present study, the well-known Roe's FDS based on TVD-MUSCL scheme is used for the solution of very high Mach number three-dimensional flows posing carbuncle and non-physical phenomena in numerical analysis. A parametric study was carried out to account for the effects of the entropy fixing, grid configurations and initial condition. The carbuncle phenomena could be easily overcome by the entropy fixing, and the non-physical solution could be eliminated by the use of the modified initial condition regardless of entropy fixing and grid configurations.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.177-181
• /
• 2007

Present study examines the numerical issues of cell structure simulation for various regimes of detonation phenomena ranging from weakly unstable to highly unstable detonations. Inviscid fluid dynamics equations with $variable-{\gamma}$ formulation and one-step Arrhenius reaction model are solved by a MUSCL-type TVD scheme and 4th order accurate Runge-Kutta time integration scheme. A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the detonation wave cell structure by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for the capturing cell structures of the different regimes of the detonation phenomena.

• 한국전산유체공학회:학술대회논문집
• /
• 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.

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