A numerical study of the diffusive-thermal instability of opposed nonpremixed tubular flames

Journal Article · Thu Oct 22 00:00:00 EDT 2015 · Combustion and Flame

DOI:https://doi.org/10.1016/j.combustflame.2015.09.019· OSTI ID:1497641

Bak, Hyun Su ^[1]; Lee, Su Ryong ^[2]; Chen, Jacqueline H. ^[3]; Yoo, Chun Sang  ^[4]

---

1. Ulsan National Inst. of Science and Technology (UNIST) (Korea, Republic of). Dept. of Mechanical Engineering
2. Seoul National Univ. of Science and Technology (Korea, Republic of). Dept. of Mechanical and Automotive Engineering
3. Sandia National Lab. (SNL-CA), Livermore, CA (United States). Combustion Research Facility
4. Ulsan National Inst. of Science and Technology (UNIST) (Korea, Republic of). Dept. of Mechanical Engineering. School of Mechanical and Nuclear Engineering

The diffusive-thermal (D-T) instability of opposed nonpremixed tubular flames near extinction is investigated here using two-dimensional (2-D) direct numerical simulations together with the linear stability analysis. Two different initial conditions (IC), i.e. the perturbed IC and the C-shaped IC are adopted to elucidate the effects of small and large amplitude disturbances on the formation of flame cells, similar to conditions found in linear stability analysis and experiments, respectively. The characteristics of the D-T instability of tubular flames are identified by a critical Damköhler number, Da_C, at which the D-T instability first occurs and the corresponding number of flame cells for three different tubular flames with different flame radii. It is found that Da_C predicted through linear stability analysis shows good agreement with that obtained from the 2-D simulations performed with two different ICs. The flame cell number, N_cell, from the 2-D simulations with the perturbed IC is also found to be equal to an integer close to the maximum wavenumber, k_max, obtained from the linear stability analysis. However, N_cell from the 2-D simulations with the C-shaped IC is smaller than k_max and N_cell found from the simulations with the perturbed IC. This is primarily because the strong reaction at the edges of the horseshoe-shaped cellular flame developed from the C-shaped IC is more likely to produce larger flame cells and reduce N_cell. It is also found that for cases with the C-shaped IC, once the cellular instability occurs, the number of flame cells remains constant until global extinction occurs by incomplete reaction manifested by small Da. It is also verified through the displacement speed, S_d, analysis that the two edges of the horseshoe-shaped cellular flame are stationary and therefore do not merge due to the diffusion–reaction balance at the edges. Moreover, large negative S_d is observed at the local extinction points while small positive or negative S_d features in the movement of flame cells as they adjust their location and size towards steady state.

• Research Organization: Sandia National Lab. (SNL-CA), Livermore, CA (United States); Ulsan National Inst. of Science and Technology (UNIST) (Korea, Republic of)
• Sponsoring Organization: USDOE Office of Science (SC), Basic Energy Sciences (BES); National Research Foundation of Korea (NRF); Ministry of Science, ICT and Future Planning (MSIP) of Korea
• Grant/Contract Number: AC04-94AL85000; 2015M1A3A3A02027319; 2015R1A2A2A01007378
• OSTI ID: 1497641
• Alternate ID(s): OSTI ID: 1254880
• Report Number(s): SAND2015-10478J; 672293
• Journal Information: Combustion and Flame, Vol. 162, Issue 12; ISSN 0010-2180
• Publisher: Elsevier
• Country of Publication: United States
• Language: English

Related Subjects

42 ENGINEERING
diffusive-thermal instability
nonpremixed tubular flame
hydrogen
linear stability analysis
displacement speed