• Journal of Mechanical Science and Technology
• /
• 제19권9호
• /
• pp.1811-1820
• /
• 2005

Thermo-acoustic instability of methane/ air flames in an industrial gas-turbine combustor is numerically investigated adopting CFD analysis. The combustor has 37 EV burners through which methane and air are mixed and then injected into the chamber. First, steady fuel! air mixing and flow characteristics established by the burner are investigated by numerical analysis with single burner. And then, based on information on the flow data, the burners are modeled numerically via equivalent swirlers, which facilitates the numerical analysis with the whole combustion system including the chamber and numerous burners. Finally, reactive flow fields within the chamber are investigated numerically by unsteady analysis and thereby, spontaneous instability is simulated. Based on the numerical results, scaling analysis is conducted to find out the instability mechanism in the combustor and the passive control method to suppress the instability is proposed and verified numerically.

• International Journal of Aeronautical and Space Sciences
• /
• 제7권2호
• /
• pp.26-32
• /
• 2006

This paper targets a direct and quantitative prediction of characteristics of unstable waves in a combustion chamber, which employs the governing equations derived in terms of amplification factors of flow variables. A freshly formulated nonlinear acoustic equation is obtained and the analysis of unsteady waves in a rocket engine is attempted. In the present formalism, perturbation method decomposes the variables into time-averaged part that can be obtained easily and accurately and time-varying part which is assumed to be harmonic. Excluding the use of conventional spatially sinusoidal eigenfunctions, a direct numerical solution of wave equation replaces the initial spatial distribution of standing waves and forms the nonlinear space-averaged terms. Amplification factor is also calculated independently by the time rate of changes of fluctuating variables, and is no longer an explicit function for compulsory representation. Employing only the numerical computation, major assumptions inevitably inherent, and in erroneous manner, in up to date analytical methods could be avoided. With two definitions of amplification factor, 1-D stable wave and 3-D unstable wave are examined, and clearly demonstrated the potentiality of a suggested theoretical-numerical method of combustion instability.

• 한국자동차공학회논문집
• /
• 제10권5호
• /
• pp.81-89
• /
• 2002

The autoignition and subsequent flame propagation of initially nonpremixed turbulent system have been numerically analyzed. The unsteady flamelet modeling based on the RIF (representative interactive flamelet) concept has been employed to account for the influences of turbulence on these essentially transient combustion processes. In this RIF approach, the partially premixed burning, diffusive combustion and formation of pollutants(NOx, soot) can be consistently modeled by utilizing the comprehensive chemical mechanism. To treat the spatially distributed inhomogeneity of scalar dissipation rate, the multiple RIFs are employed in the framework of EPFM(Eulerian particle flamelet model) approach. Computations are made for the various initial conditions of pressure, temperature, and fuel composition. The present turbulent combustion model reasonably well predicts the essential features of autoignition process in the transient gaseous fuel jets injected into high pressure and temperature environment.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2009년도 춘계학술대회 논문집
• /
• pp.79-82
• /
• 2009

초음속 엔진의 흡입구에서의 종말충격파와 연소실 화염의 상호간섭 연구를 위하여 초음속 엔진의 전영역, 즉 흡입구에서부터 연소실과 노즐까지 통합하여 비정상 연소수치해석을 수행하였다. 초음속 엔진이 상승하는 가속모드와 순항모드에서 상호간섭의 동적현상을 연구하였다. 흡입구에서의 충격파거동과 주요 위치에서 압력거동을 분석하고 초음속 엔진 전영역에서의 음향모드를 분석하여 현 시스템의 동적 거동을 파악하였다.

• 한국추진공학회지
• /
• 제4권1호
• /
• pp.74-92
• /
• 2000

A numerical analysis has been conducted to study the transient flowfield during the transition from the booster to sustainer phase in an integrated rocket ramjet (IRR) propulsion system. Emphasis is placed on the unsteady inlet aerodynamics, fuel/air mixing in an entire ramjet engine during the flow transient phase. The computational geometry consists of the entire IRR engine, including the inlet, the combustion chamber, and the exhaust nozzle. Turbulence closure is achieved using a low-Reynolds-number two-equation model. The governing equations are solved numerically by means of a finite-volume, preconditioned flux-differencing scheme over a wide range of Mach umber. Various important physical processes are investigated systemically, including terminal shock train.

• 한국연소학회:학술대회논문집
• /
• 한국연소학회 2006년도 제32회 KOSCO SYMPOSIUM 논문집
• /
• pp.205-212
• /
• 2006

Three-dimensional numerical study was carried out for the investigation of the detonation wave structures propagating in tubes. Fluid dynamics equations and conservation equation of reaction progress variable were analyzed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Chemical reaction was modeled by using a simplified one-step irreversible kinetics model. The variable gas properties between unburned and burned states were considered by using variable specific heat ratio formulation. The unsteady computational results in three-dimension show the detailed mechanisms of rectangular and diagonal mode of detonation wave instabilities resulting same cell length but different cell width in smoked-foil record. The results for the small reaction constant shows the spinning mode of three-dimensional detonation wave dynamics, which was rarely observed in the previous numerical simulation of the detonation waves.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2004년도 제22회 춘계학술대회논문집
• /
• 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.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2003년도 제20회 춘계학술대회 논문집
• /
• 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.

• 한국자동차공학회논문집
• /
• 제18권2호
• /
• pp.7-13
• /
• 2010

To satisfy the needs on fuel economy and engine performance, continuous variable valve lift systems are applying to engines. In the CVVL system, fuel economy can be improved by reducing pumping loss during the induction process, and engine performance can be also improved by controlling volumetric efficiency and the residual gas fraction. Because the residual gas fraction directly affects volumetric efficiency, engine performance, combustion efficiency and emissions in SI engines, controlling residual gas fraction is one of the important things in engine development process. This analysis investigates the residual gas fraction and volumetric efficiency with changes of intake valve lifts and intake valve timings. In this study, unsteady state solutions were solved during exhaust and induction processes. Results show variation of the residual gas fraction and volumetric efficiency by changing intake valve timing and lift. Decreasing intake valve lift leads to increase the residual gas fraction and to decrease volumetric efficiency.