• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.7
• /
• pp.1-8
• /
• 2003

Numerical characteristics of air chemistry associated with hypersonic flows are described and are compared with those of hydrogen oxygen combustion, applying the partially implicit time integration method to air chemistry. This paper reveals that the time integration of air chemistry needs a chemical Jacobian for stable calculations. However the positive real eigenvalues in air chemistry are relatively smaller than those of hydrogen combustion, and the numerical integration is less sensitive than that with combustion. lt is also found that the application of the partia1ly irnplicit method reduces the computing time without numerical instabilities.

• Journal of computational fluids engineering
• /
• v.10 no.3 s.30
• /
• pp.9-18
• /
• 2005

During a hypothetical severe accident in a nuclear power plant (NPP), hydrogen is generated by the active reaction of fuel-cladding and steam in the reactor pressure vessel and released with steam into the containment. In order to mitigate hydrogen hazards possibly occurred in the NPP containment, hydrogen mitigation system (HMS) is usually adopted. The design of the next generation NPP (APR1400) designed in Korea specifies 26 passive autocatalytic recombiners and 10 igniters installed in the containment for the hydrogen mitigation. in this study, the analysis of the hydrogen and steam behavior during a total lose of feed water (TLOFW) accident in the APR1400 containment has been conducted by using the CFD code GASFLOW. During the accident, a huge amount of hot water, steam, and hydrogen is released in the in-containment refueling water storage tank (IRWST). The current design of the APR1400 includes flap-type dampers at the IRWST vents which are operated depending on the pressure difference between inside and outside of the IRWST. it was found that the flaps strongly affects the flow structure of the steam and hydrogen in the containment. The possibilities of a flame acceleration and transition from deflagration to detonation (DDT) were evaluated by using Sigma-Lambda criteria. Numerical results indicate the DDT possibility could be heavily reduced in the IRWST compartment when the flaps are installed.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.1
• /
• pp.78-84
• /
• 2005

In the chemical laser system with a radial expansion nozzle array, the laser beam generation is achieved by mixing F atom from supersonic nozzle and $D_{2}$ molecule from holes of round-bended supply line. Based on that the fuel injection angle with main stream has a great influence of performance on supersonic combustor, the effects of $D_{2}$ injection angles with the main F flow on mixing enhancement are numerically investigated. The results are discussed by comparison with three cases of $D_{2}$ injection angles; $10^{o}$, $20^{o}$ and $40^{o}$ with the main flow direction. Major results reveal that as the $D_{2}$ injection angle increases, the maximum small signal gains and the static pressure in the laser cavity become higher. Consequently, the $D_{2}$ injection angle between $20^{o}$ and $40^{o}$ is recommended as an optimized geometric parameter in consideration of both of high gains and low cavity pressure.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.383-387
• /
• 2017

A three dimensional numerical analysis has been conducted to analyze the effects of a pipe angle, connecting a combustion chamber and a pintle nozzle, and pintle position on pintle nozzle performance. The compressibility correction of $k-{\omega}$ SST turbulent model was implemented to precisely predict the characteristics of complex flow structures inside a supersonic pintle nozzle. Due to an 3-D asymmetric pintle nozzle configuration, complex helical flow streamlines and large flow separations were observed, which resulting in significant nozzle performance losses. As the angle of connection-tube decreases, the coefficient of performance increases and Since the flow structures are evidently changed to the pintle stroke position, the performance characteristics was analyzed.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.11a
• /
• pp.299-304
• /
• 2008

Starting characteristics of the axi-symmetric second throat exhaust diffuser (STED) with zero-secondary flows are numerically investigated. Renolds-Average Navier-Stokes equations with a standard ${\kappa}-{\varepsilon}$ turbulence model incorporated with enhanced wall treatment are solved to simulate the diffusing evolutions of the nozzle plume. Minimum (optimum) starting pressure difference of 20$\sim$25% between 1-D theory and the measured data validated from previous results[5] is also applied to predict the range of an effective diffuser expansion ratio (Ad/At) in this system.

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