• Journal of Energy Engineering
• /
• v.21 no.3
• /
• pp.237-242
• /
• 2012

This study presents numerical solutions of the two-dimensional Navier-Stokes equations for supersonic unsteady flow in a convergent-divergent nozzle with a isolator. The TVD scheme in generalized coordinates is employed in order to calculate the moving shock waves caused by thermal choking. We discuss on transient characteristics, unstart phenomena, fluctuations of specific thrust caused by thermal choking and effects of isolator. The adverse pressure gradient caused by heat addition brings about separation of the wall boundary layers and formation of the oblique shock wave that proceed to upstream. The proceeding speed of the oblique shock wave to upstream direction for the convergent-divergent nozzle with isolator is lower than that for the nozzle without isolator.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.62-68
• /
• 2004

A comprehensive numerical analysis has been carried out for both non-reacting and reacting flows in a scramjet engine combustor with and without a cavity. The theoretical formulation treats the complete conservation equations of chemically reacting flows with finite-rate chemistry of hydrogen-air. Turbulence closure is achieved by means of a k-$\omega$ two-equation model. The governing equations are discretized using a MUSCL-type TVD scheme, and temporally integrated by a second-order accurate implicit scheme. Transverse injection of hydrogen is considered over a broad range of injection pressure. The corresponding equivalence ratio of the overall fuel/air mixture ranges from 0.167 to 0.50. The work features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous studies. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the .underlying physical mechanisms. Much of the flow unsteadiness is related not only to the cavity, but also to the intrinsic unsteadiness in the flow-field. The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The roles of the cavity, injection pressure, and heat release in determining the flow dynamics are examined systematically.

• Journal of computational fluids engineering
• /
• v.10 no.4 s.31
• /
• pp.59-65
• /
• 2005

A numerical analysis was made to investigate the simple silencer for high pressure blast flow fields. Reynolds-Averaged Navier-Stokes equations were solved for an axisymmetric computational domain constructed by multi block grids. A blast flow field without the silencer was also calculated to validate the present numerical method. The effect of pressure diminution for the silencer was calculated by comparing with and without silencer at the atmosphere region. It was found that the tested silencer could achieve 89.4 percent pressure diminution.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.85-94
• /
• 1999

Transient aerodynamic response of an airfoil to a moving plane-flap is numerically investigated using two-dimensional Euler equations with conservative Chimera grid method. A body moving relative to a stationary grid is treated by an overset grid bounded by a 'dynamic domain-dividing line' the concept of which is developed in this study. A conservative Chimera grid method with a dynamic domain-dividing line technique is applied and validated by solving the flowfield around circular cylinder moving supersonic speed. The unsteady and transient characteristics of the flow solver is also examined by computations of a oscillating airfoil and a ramp pitching airfoil respectively. The transient aerodynamic behavior of an airfoil with a moving plane-flap is analyzed for various flow conditions such as deflecting rate of flap and free stream Mach number.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.

• 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 for Aeronautical & Space Sciences
• /
• v.37 no.10
• /
• pp.955-966
• /
• 2009

DES method is applied to an axisymmetric base flow at supersonic mainstream. The model is based on the Spalart-Allmaras (S-A) turbulence model in the RANS mode, and is based on the subgrid scale model in the Large-eddy simulation (LES) mode. Accurate predictions of the base flowfield and base pressure are successfully achieved by using the DES methodology which is less expensive than LES. Flow properties at the edge of base, such as boundary layer thickness, momentum thickness and skin fraction are compared with Dutton et al [experimental data to proper prediction of base flowfiled. From the present results, The DES accurately resolves the physics of unsteady turbulent motions, such as shear layer rollup, large-eddy motions in the downstream region and small eddy motions inside the recirculating region. Moreover, The present results of using an empirical constant $C_{DES}$ of 1.2 shows good agreement with experimental data than conventional empirical constant $C_{DES}$ of 0.65.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.05a
• /
• pp.47-52
• /
• 2002

The supersonic flow around multi-stage rocket system is analyzed using 3-D compressible unsteady flow solver. A Chimera overset grid technique is used for the calculation of present configuration and grid around the core rocket is composed of 3 zones to represent fins in the core rocket. Flow solver is parallelized to reduce the computation time, and an efficient parallelization algorithm for Chimera grid technique is proposed. AUSMPW+ scheme is used for the spatial discretization and LU-SGS for the time integration. The flow field around multi-stage rocket was analyzed using this developed solver, and the results were compared with that of a sequential solver The speed-up ratio and the efficiency were measured in several processors. As a result, the computing speed with 12 processors was about 10 times faster than that of a sequential solver. Developed flow solver is used to predict the trajectory of booster in separation stage. From the analyses, booster collides against core rocket in free separation case. So, additional jettisoning forces and moments needed for a safe separation are examined.

• Journal of the Korean Society of Visualization
• /
• v.8 no.3
• /
• pp.56-62
• /
• 2010

The interaction patterns between shock wave and boundary layer in a rocket nozzle are mainly classified into two categories, FSS(Free Shock Separation) and RSS(Restricted Shock Separation), both of which are associated with the thrust characteristics as well as side loads of the engine. According to the previous investigations, strong side loads of the engine are produced during the period of transition from FSS to RSS or vice versa. The present work aims at investigating the unsteady behavior of the separation shock waves in a two-dimensional supersonic nozzle, using experimental method and CFD. Schlieren optical method was employed to visualize the time-mean and time-dependent shock motions in the nozzle. The unsteady, compressible N-S equations with SST K-$\omega$ turbulence closure were solved using a fully implicit finite volume scheme. The results obtained show the separation shock motions during the transition of the interaction pattern.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.19 no.3
• /
• pp.29-38
• /
• 2015

In this study, 3-dimensional URANS simulation was performed to analyze the effect of the nozzle-rotor axial gap spacing of a supersonic impulse turbine on turbine performance. The computations were conducted for four different axial gap cases corresponding to about 6%, 10%, 20% and 30% of the blade height, respectively. The results show a good agreement with previous studies and the turbine efficiency decreases drastically in certain range. It is examined that the turbine performance characteristics could change depending on the influence of leading edge shock to the nozzle outlet. It is also found that the entropy rise distributions along the span differ from each other.