• Journal of the Korean Society of Propulsion Engineers
• /
• v.9 no.1
• /
• pp.35-45
• /
• 2005

Numerical study is carried out to investigate the effects of chemistry and thermal radiation on the rocket plume flow field at various altitudes. Navier-Stokes equations for compressible flows were solved by a fully-implicit TVD code based on the finite volume method. An infinitely fast chemistry module for hydrocarbon mixture with detailed thermo-chemical properties and a thermal radiation module for optically thick media were incorporated with the fluid dynamics code. The plume flow fields of a kerosene-fueled rocket flying at Mach number zero at sea-level, 1.16 at altitude of 5.06 km and 2.90 at 17.34 km were numerically analyzed. Results showed the plume structures at different altitude conditions with the effects of chemistry and radiation. It is understood that the excess temperature by the chemical reactions in the exhaust gas may not be ignored in the view point of propulsion performance and thermal protection of the rocket base, especially at higher altitude conditions.

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

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.11-21
• /
• 1998

A vorticity-based method for the numerical solution of the two-dimensional incompressible Navier-Stokes equations is presented. The governing equations for vorticity, velocity and pressure variables are expressed in an integro-differential form. The global coupling between the vorticity and the pressure boundary conditions is fully considered in an iterative procedure when numerical schemes are employed. The finite volume method of the second order TVD scheme is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition. The velocity field is obtained by using the Biot-Savart integral. The Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well established for potential flow analysis. The present formulation is validated by comparison with data from the literature for the two-dimensional cavity flow driven by shear in a square cavity. We take two types of the cavity now: (ⅰ) driven by non-uniform shear on top lid and body forces for which the exact solution exists, and (ⅱ) driven only by uniform shear (of the classical type).