• Journal of Mechanical Science and Technology
• /
• v.17 no.11
• /
• pp.1775-1783
• /
• 2003

The method of flame hole dynamics is demonstrated as a mean to simulate turbulent flame extinction. The core of the flame hole dynamics involves derivation of a random walk mapping for the flame holes, created by local quenching, between burning and quenched states provided that the dynamic characteristics of flame edges is known. Then, the random walk mapping is projected to a background turbulent field. The numerical simulations are carried out with further simplifications of flame string and unconditioned scalar dissipation rate. The simulation results show how the chance of partial quenching is influenced by the crossover scalar dissipation rate. Finally, a list of improvements, necessary to achieve more realistic turbulent flame quenching simulation, are discussed.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.10 no.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.

• 한국연소학회:학술대회논문집
• /
• 2004.06a
• /
• pp.102-111
• /
• 2004

Partial quenching structure of turbulent diffusion flames in a turbulent mixing layer is investigated by the method of flame hole dynamics to develope a prediction model for the turbulent lift off. The present study is specifically aimed to remedy the problem of the stiff transition of the conditioned partial burning probability across the crossover condition by adopting level-set method which describes propagating or retreating flame front with specified propagation speed. In light of the level-set simulations with two model problems for the propagation speed, the stabilizing conditions for a turbulent lifted flame are suggested. The flame hole dynamics combined with level-set method yields a temporally evolving turbulent extinction process and its partial quenching characteristics is compared with the results of the previous model employing the flame-hole random walk mapping. The probability to encounter reacting' state, conditioned with scalar dissipation rate, demonstrated that the conditional probability has a rather gradual transition across the crossover scalar dissipation rate in contrast to the stiff transition of resulted from the flame-hole random walk mapping and could be attributed to the finite response of the flame edge propagation.

• Transactions of the Korean hydrogen and new energy society
• /
• v.18 no.3
• /
• pp.265-274
• /
• 2007

A study has been conducted to clarify NO emission behavior through preferential diffusion effects of $H_2$ and H in methane-hydrogen diffusion flames. A comparison is made by employing three species diffusion models. Special concerns are focused on what is the deterministic role of the preferential diffusion effects in flame structure and NO emission. The behavior of maximum flame temperatures with three species diffusion models is not explained by scalar dissipation rate but the nature of chemical kinetics. The preferential diffusion of H into reaction zone suppresses the populations of the chain carrier radicals and then flame temperature while that of $H_2$ produces the increase of flame temperature. These preferential diffusion effects of $H_2$ and H are also discussed about NO emissions through the three species diffusion models.

• Journal of the Korean Society of Combustion
• /
• v.5 no.2
• /
• pp.19-27
• /
• 2000

A turbulent combustion model, based on edge flame dynamics, is discussed in order to predict global extinction of turbulent flames. The model is applicable to the broken flamelet regime of turbulent combustion, in which global extinction of turbulent flame is achieved by gradual expansion of flame holes. The edge flame dynamics is the key mechanism to describe the flame hole expansion or contraction. For flames with Lewis numbers near unity, there is a $Damk{\ddot{o}}hler$ number, namely the crossover $Damk{\ddot{o}}hler$ number, at which edge flame changes its direction of propagation. The parametric region between the quasi-steady extinction condition and the edge-flame crossover condition is a metastable region, in that flames without edge can stay in their burning states while flames with edge have to retract to expand quenching holes. Using the above properties of edge flame, Hartley and Dold proposed a Lagrangian hole dynamics, which allows us to simulate transient variation of quenching holes. In their model, each stoichiometric surface is subjected to a random sequence of scalar dissipation rate compatible to the equilibrium turbulence. Then, each stoichiometric surface will evolve, according to the combustion map, dependent on the scalar dissipation rate and existence of flame edge, If all the burning surfaces are annihilated, the event can be declared as a global extinction. The consequence obtained from the above model also can be used as a subgrid model to determine local extinction occurring in a calculation grid.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.11
• /
• pp.1407-1416
• /
• 2004

A two-dimensional direct numerical simulation was performed to investigate the flame structure of CH$_4$/$N_2$-Air counterflow nonpremixed flame interacting with a single vortex. The detailed transport properties and a modified 16-step augmented reduced mechanism based on Miller and Bowman's detailed chemistry were adopted in this computation. The results showed that an initially flat stagnation plane, on which an axial velocity was zero, was deformed into a complex-shaped plane, and an initial stagnation point was moved far away from a vortex head when the counterflow field was perturbed by the vortex. It was noted that the movement of stagnation point could alter the species transport mechanism to the flame surface. It was also identified that the altered species transport mechanism affected the distributions of the mixture fraction and the scalar dissipation rate.

• 한국연소학회:학술대회논문집
• /
• 2007.05a
• /
• pp.118-125
• /
• 2007

The firtst-order conditional moment closure (CMC) model is applied to CH4/air swirl diffusion flame in a gas turbine model combustor. The flow and mixing fields are calculated by fast chemistry assumption with SLFM library and a beta function pdf for mixture fraction. RNG k-e model is used to consider the swirl flame in a confined wall. Reacting scalar fields are calculated by elliptic CMC formulation with chemical kinetic mechanism, GRI Mech 3.0. Validation is done against measurement data for mean flow and scalar fields in the model combustor [1]. Results show reasonable agreement with the mean mixture fraction and its variance, while temperature is overpredicted as the level of local extinction increases. The second-order CMC model is needed to consider local extinction with considerable conditional fluctuations near the nozzle.

• 한국연소학회:학술대회논문집
• /
• 2002.11a
• /
• pp.11-17
• /
• 2002

The first-order conditional moment closure (CMC) model is applied to $CH_4$/Air turbulent jet diffusion flames(Sandia Flame D, E and F). The flow and mixing fields are calculated by fast chemistry assumption and a beta function pdf for mixture fraction. Reacting scalar fields are calculated by elliptic CMC formulation. The results for Flame D show reasonable agreement with the measured conditional mean temperature and mass fractions of major species, although with discrepancy on the fuel rich side. The discrepancy tends to increase as the level of local extinction increases. Second-order CMC may be needed for better prediction of these near-extinction flames.

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

고체나 액체 추진로켓에 비하여 하이브리드 추진 시스템은 작동조건의 안정성과 안전함등의 많은 장점을 가지고 있다. HTPB와 같은 고체연료는 제작 및 저장, 운송 그리고 장착상의 안정성을 가지고 있으며 하이브리드 로켓의 고체연료로의 산화제의 유입을 제어하면서 추력의 변화와 엔진내부의 연소중단과 재 점화를 용이하게 할 수 있다. 이러한 이유로 인하여 하이브리드 엔진은 좀 더 경제적인 장치로 기대를 모으고 있다. 그러나, 기존의 하이브리드 로켓 엔진은 고체 추진 로켓에 비하여 낮은 연료 regression 율과 연소효율을 가지는 단점이 있다. 이러한 단점을 해결하고 요구되어지는 추력값과 연료유량을 증가시키기 위하여 고체연료의 표면적을 증가시킬 필요가 있다. 기존의 하이브리드 엔진에서는 연료 그레인에 다수의 연소포트를 만들어 표면적을 증가시켰으나 이는 비 활용 공간의 증가와 추진제의 질량 및 체적분율의 상당한 감소를 초래한다. 지난 수십년간에 걸쳐 하이브리드 엔진에서 연료의 regression 특성 및 엔진 성능 향상을 위한 연구가 계속되어 왔으며 최근에 엔진의 체적 규제를 경감시키고 연료의 regression율을 향상시키기 위하여 선회유동을 이용하는 하이브리드 로켓 엔진들이 제안되고 있다. 이러한 선회유동을 가지는 하이브리드 로켓은 고체연료 그레인에 대하여 평행하게 유입되는 기존의 하이브리드 로켓에 비하여 고체연료 벽면에서의 대류열전달이 현저하게 증가하게 되어 아주 높은 고체연료의 regression율을 얻을 수 있는 이점이 있다. 선회유동 하이브리드 로켓의 연소과정은 고체 연료의 열분해과정, 대류 열전달, 난류 혼합, 난류와 화학반응의 상호작용, soot의 생성 및 산화과정, soot 입자 및 연소가스에 의한 복사 열전달, 연소장과 음향장의 상호작용 등의 복잡한 물리적 과정을 포함하고 있다. 이러한 물리적 과정 중 난류연소, 고체연료 벽면 근방에서의 대류 열전달 및 연소과정에서 생성되는 soot 입자로부터의 복사 열전달, 그리고 고체연료 열 분해시 표면반응들은 고체연료의 regression율에 큰 영향을 미친다. 특히 고체연료의 난류화염면의 위치와 폭, 그리고 비 예혼합 난류화염장에서 생성되는 soot의 체적분율의 예측은 난류연소모델, 열전달 모델, 그리고 regression율 모델에 의해 크게 영향을 받기 때문에 수치모델의 예측 능력 향상시키기 위하여 이러한 물리적 과정을 정확히 모델링해야 할 필요가 있다. 특히 vortex hybrid rocket내의 난류연소과정은 아래와 같은 Laminar Flamelet Model에 의해 모델링 하였다. 상세 화학반응 과정을 고려한 혼합분율 공간에서의 화염편의 화학종 및 에너지 보존 방정식은 다음과 같다. 화염편 방정식과 혼합분률과 scalar dissipation rate의 관계식을 이용하여 혼합분률과 scalar dissipation rate에 따른 모든 reactive scalar들을 구하게 된다. 이러한 화염편 방정식들을 mixture fraction space에서 이산화시켜서 얻은 비선형 대수방정식은 TWOPNT(Grcar, 1992)로 계산돼 flamelet Library에 저장되게 된다. 저장된 laminar flamelet library를 이용하여 난류화염장의 열역학 상태량 평균치는 presumed PDF approach에 의해 구해진다. 본 연구에서는 강한 선회유동을 가지는 Hybrid Rocket 연소장내의 난류와 화학반응의 상호작용을 분석하기 위하여 Laminar Flamelet Model, 화학평형모델, 그리고 Eddy Dissipation Model을 이용한 수치해석결과를 체계적으로 비교하였다. 또한 Laminar Flamelet Model과 state-of-art 물리모델들을 이용하여 선회 유동을 갖는 하이브리드 로켓 엔진의 연소 및 Soot 생성 및 산화과정을 살펴보았으며 복사 열전달이 고체 연료 표면의 regression율에 미치는 영향도 살펴보았다. 특히 swirl강도, 산화제의 유입위치 그리고 선회유동의 형성방식이 하이브리드 로켓의 연소특성 및 regression rate에 미치는 영향을 상세히 해석하였다.

• Journal of the Korean Society of Combustion
• /
• v.4 no.2
• /
• pp.51-62
• /
• 1999

NOx formation in turbulent flames is strongly coupled with temperature, superequilibrium concentration of O radical, and residence time. This implies that in order to accurately predict NO level, it is necessary to develop sophisticated models able to account for the complex turbulent combustion processes including turbulence/chemistry interaction and radiative heat transfer. The present study numerically investigates the turbulent nonpremixed hydrogen jet flames using the laminar flamelet model. Flamelet library is constructed by solving the modified Peters equations and the turbulent combustion model is extended to nonadiabatic flame by introducing the enthalpy defect. The effects of turbulent fluctuation are taken into account by the presumed joint PDFs for mixture fraction, scalar dissipation rate, and enthalpy defect. The predictive capability of the present model has been validated against the detailed experimental data. Effects of nonequilibrium chemistry and radiative heat loss on the thermal NO formation are discussed in detail.