Next Article in Journal
Statistical Analysis of Bubble Parameters from a Model Bubble Column with and without Counter-Current Flow
Next Article in Special Issue
Fourier Neural Operator Networks for Solving Reaction–Diffusion Equations
Previous Article in Journal
Subgrid Turbulent Flux Models for Large Eddy Simulations of Diffusion Flames in Space Propulsion
Previous Article in Special Issue
Turbulent Channel Flow: Direct Numerical Simulation-Data-Driven Modeling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Partial Premixing and Coflow Temperature on Flame Stabilization of Lifted Jet Flames of Dimethyl Ether in a Vitiated Coflow Based on Stochastic Multiple Mapping Conditioning Approach

1
School of Engineering, Newcastle University, Claremont Road, Newcastle upon Tyne NE1 7RU, UK
2
Combustion and Energy Conversion Systems Laboratory, Department of Mechanical Engineering, Indian Institute of Technology Kanpur, Kanpur 208016, India
*
Author to whom correspondence should be addressed.
Current address: Maruti Suzuki India Limited, Gurugram 122015, India.
Fluids 2024, 9(6), 125; https://doi.org/10.3390/fluids9060125
Submission received: 13 February 2024 / Revised: 5 May 2024 / Accepted: 8 May 2024 / Published: 26 May 2024

Abstract

:
The Reynolds-averaged Navier–Stokes (RANS)-based stochastic multiple mapping conditioning (MMC) approach has been used to study partially premixed jet flames of dimethyl ether (DME) introduced into a vitiated coflowing oxidizer stream. This study investigates DME flames with varying degrees of partial premixing within a fuel jet across different coflow temperatures, delving into the underlying flame structure and stabilization mechanisms. Employing a turbulence k- ε model with a customized set of constants, the MMC technique utilizes a mixture fraction as the primary scalar, mapped to the reference variable. Solving a set of ordinary differential equations for the evolution of Lagrangian stochastic particles’ position and composition, the molecular mixing of these particles is executed using the modified Curl’s model. The lift-off height (LOH) derived from RANS-MMC simulations are juxtaposed with experimental data for different degrees of partial premixing of fuel jets and various coflow temperatures. The RANS-MMC methodology adeptly captures LOH for pure DME jets but exhibits an underestimation of flame LOH for partially premixed jet scenarios. Notably, as the degree of premixing escalates, a conspicuous underprediction in LOH becomes apparent. Conditional scatter and contour plots of OH and CH2O unveil that the propagation of partially premixed flames emerges as the dominant mechanism at high coflow temperatures, while autoignition governs flame stabilization at lower coflow temperatures in partially premixed flames. Additionally, for pure DME flames, autoignition remains the primary flame stabilization mechanism across all coflow temperature conditions. The study underscores the importance of considering the degree of premixing in partially premixed jet flames, as it significantly impacts flame stabilization mechanisms and LOH, thereby providing crucial insights into combustion dynamics for various practical applications.

1. Introduction

Lifted jet diffusion flames (LJDF) have practical applications in industrial burners and combustors, where the flame stabilizes away from the nozzle and thereby minimizes damage to the hardware [1,2]. The lift-off height (LOH) of these flames are highly sensitive to changes in coflow temperature. For example, a 1% change in the coflow temperature doubles the LOH for H2/N2 LJDF [3]. Further, at low coflow temperatures, the flame base undergoes large fluctuations, which makes the stabilization of lifted flames a complex phenomenon. Mechanisms that control flame stabilization, lift-off, and blowout have been studied extensively in the past. A comprehensive review on the stabilization of lifted jet flames by Lyons [1] presents competing theories of premixed flame propagation, autoignition, and the role of scalar dissipation as possible stabilization mechanisms of these flames. The stabilization mechanism strongly relies on turbulent transport and chemical kinetics at the flame base [4]. A detailed analysis of the flame structure and chemical processes occurring during and prior to ignition is expected to provide insights and help us to distinguish between the premixed flame propagation and autoignition stabilization mechanisms of LJDF. For autoignition-stabilized flames, HO2 and CH2O radicals form during the pre-ignition phase [5]. The subsequent decomposition of these radicals leads to the production of hydroxyl (OH) radical, which triggers high-temperature reactions. In premixed flame propagation, the formation of HO2 and CH2O radicals remains limited to the fuel-rich region, and they are quickly consumed in the high-temperature zone [6]. Due to the mixing between fuel and oxidizer streams in lifted flames at conditions not suitable for autoignition, a partially premixed mode of combustion is often observed, leading to a two-dimensional tribrachial flame (also known as triple flame) [1] consisting of a rich and a lean premixed flame and a trailing diffusion flame. The experimental setup designed by Cabra et al. [4,7] presents unique challenges for combustion models due to the delicate balance between chemical reactions and species diffusion. The detailed measurement of conditional and unconditional scalars provided by the Cabra burner series of H2/N2 [7] and CH4/air [4] flames makes it an ideal platform for developing and validating turbulent combustion models. An analysis of H2/N2 jets introduced into a vitiated coflow through PDF simulations indicates that the LOH is primarily influenced by chemical kinetics rather than mixing models [3], with minimal dependence on turbulence intensity at the inlet boundary. Numerical investigations employing PDF methods demonstrate that the LOH is strongly influenced by coflow temperature and jet velocity, with the flame being predominantly kinetically controlled, thus relying on autoignition for stabilization [8]. The experimental study of Gordon et al. [9] identifies two distinct stabilization mechanisms: turbulent premixed flame propagation and autoignition, observed at low and high coflow temperatures, respectively. Subsequent numerical analyses confirm these mechanisms through transport budget comparisons [8], while Reynolds-averaged Navier–Stokes (RANS)-based conditional moment closure (CMC) simulations further support these findings [10]. Yadav et al. [11] conducted simulations of H2/N2 lifted jet diffusion flames using a multi-environment Eulerian PDF (MEPDF) transport approach. Despite employing an advanced turbulence–chemistry interaction model, the predictions from their study did not closely correlate with the experimental data. Various research groups have reported numerical simulations of LJDF using different methodologies. These include the LES-based CMC method [12], RANS-CMC with an extinction model [13], approximated diffusion flame presumed conditional moment [14], RANS/transported PDF [4,5], LES/presumed PDF methods [15,16] with detailed and reduced chemistries, and the LES/thickened flame method [17]. Recently, an LES-based PDF simulation [18] incorporated the effects of differential diffusion and found that molecular diffusion significantly impacts the stabilization zone of lifted H2/N2 flames. The presence of a wide range of length and time scales, and highly non-linear, complex interplay between turbulence and chemistry, poses a formidable challenge towards accurate numerical modeling of LJDF. Only a few models for turbulent combustion could accurately predict the flame lift-off height (LOH), especially at low coflow temperatures. Among them, presumed PDF [15], CMC [10], MMC [19,20], and transported PDF [4,5] are a few models which could capture the variation in LOH with the RANS and large eddy simulation (LES) framework. The objectives of the present work are to investigate the complex turbulence–chemistry interaction of dimethyl ether (DME) lifted flames, which have relatively complex chemical kinetics compared to methane, and to determine the effects of partial premixing of fuel jets and coflow temperature on flame characteristics, such as LOH, structure, and stabilization mechanism. Numerical simulations are performed based on an in-house RANS-MMC solver, which was validated earlier for pilot-stabilized flames of CH4 [21] and DME [22] and lifted flames of H2/N2 [19] and CH4/air [20]. Due to complex chemical kinetics and strong turbulence–chemistry interactions, the lifted DME flame series [23] serves as an ideal candidate to evaluate the predictive capability of the RANS-MMC model. This is the first application of stochastic MMC for lifted DME flames, which display strong turbulence–chemistry interactions. In stochastic MMC, all scalars are classified into the major and the minor groups. The minor scalars, e.g., temperature and species mass fractions, are allowed to fluctuate relatively with respect to the major scalars. The mixture fraction is used as the major scalar, which is mapped with the reference space. Akin to mixture fraction-based PDF methods [10], the major scalar fluctuates with the underlying turbulent flows. The present paper is organized as follows. A brief description of the RANS-MMC approach is provided in the next section, followed by the experimental configuration and numerical details in Section 3. The numerical results from the simulations are presented in Section 4 along with an analysis of DME flames. Finally, major conclusions are drawn in Section 5.

2. Mathematical Background

In the present RANS-MMC approach, a single reference variable, ξ * , whose mapping function X ( ξ * ) represents the mixture fraction Z, is used. In this stochastic implementation, the reference variable, ξ * , is assumed to be Gaussian-distributed with zero mean and unit variance. A separate equation for the evolution of the reference variable, ξ * , is solved. The corresponding sets of equations for a stochastic particle method are as follows [24,25]:
Transport in physical space : d x * = U ( ξ * , x * , t ) d t ,
Transport in reference space : d ξ * = A ( ξ * , x * , t ) d t + b ( ξ * , x * , t ) d ω * ,
Transport in composition space : d X I * = ( W I * + S I * ) d t .
In the above equations, the asterisk is used to denote stochastic quantities. The positions x in physical space and the velocities U evolve according to a model represented by a Markov family. Here, the drift and diffusion coefficients are as follows:
A = B ξ + B ξ + 1 ρ ¯ · ( ρ ¯ U ( 1 ) ) + 2 P ˜ ξ B P ˜ ξ ξ , and
2 B = b 2 .
The chemical source terms and mixing operator are denoted by W I * and S I * , respectively. The mixing operator S I * helps to keep X I * close to its conditional mean value. This ensures that mixing is local in the reference space, which is the characteristic of the MMC mixing model. Further, micromixing models are needed for the modeling of the mixing operator S I * . The modified Curl’s [26] model was recently developed in the context of the MMC approach [19,20,21,22,27] and is utilized in the present work. In MMC-Curl, the mixing operator S I * involves an interaction between particles that are grouped in pairs (denoted by ‘p’ and ‘q’), and mixing is local in the reference space and particle pairs are mixed over a time step duration of Δ t [28]:
X I * , p ( t + Δ t ) = X I * , p ( t ) + λ ( X I * , p q ( t ) X I * , p ( t ) ) Δ t τ m i n ,
X I * , q ( t + Δ t ) = X I * , q ( t ) + λ ( X I * , p q ( t ) X I * , q ( t ) ) Δ t τ m i n ,
where λ = 1 e x p [ ( w p + w q ) W N 2 ] ( w p + w q ) W N 2 and w p , w q , and W are the weights of particles p and q and their sum, respectively. The number of particles per CFD cell is N. τ m i n is the minor mixing time scale. The particle pair p and q is not randomly selected but selected in such a manner that p and q are local to each other in the reference space. Therefore, all particles within one CFD cell are sorted by their reference value and stored in an array. Two neighboring particles within this array are then selected for mixing. This ensures localness in the composition space. The mixing operator for the selected particle pairs can now be written as
S I * = d X I * , k d t = λ ( X I * , p q X I * , k ) τ m i n .
The minor time scale, τ m i n , controls the conditional fluctuations of scalars around the conditional mean. The timescales τ m i n and τ ϕ are related by the following relation: τ m i n = C m i n · τ ϕ , where C m i n is the minor mixing time constant. Following recent nonpremixed [20] and premixed [27] RANS-MMC simulations, C m i n = 0.25 is used here.

3. Numerical Approach

The lifted DME/air jet diffusion flames studied experimentally by Macfarlane et al. [23,29] are used in the present study. The burner has been realized with certain geometrical modifications from a similar burner used by Cabra et al. [4] for autoignition studies of hydrogen and methane lifted jet flames. The setup consists of a central fuel jet having an inner diameter of 4.45 mm. An annular coflow of hot combustion products from a lean mixture of H2/air surrounds the central fuel jet. The coflow has an outer diameter of 197 mm with approximately 1800 holes of 1.6 mm diameter. Besides stabilizing the flames, these holes provide blockage to the flow and thereby promote turbulence. Fuel is issued from the central jet, which extends 70 mm downstream of the perforated plate. The temperature and bulk velocity of the fuel jet are 298 K and 50 m/s, respectively. Table 1 provides the conditions used in the present numerical simulations. In the experiments, the coflow temperature was varied by varying the equivalence ratio of the reactants while maintaining a constant coflow velocity of 4 m/s. Numerical simulations are performed for three different degrees of fuel/air premixing in the jet and at three different coflow temperatures.
The computational domain extends to 60 D and 22 D in the axial and radial directions. More details about the inflow boundary conditions are provided in Table 1. At the inlet, the Dirichlet boundary conditions are used for temperature, velocity, mixture fraction, and species mass fractions, whereas the Neumann boundary conditions are employed at the outlet for all flow variables except pressure. The Neumann boundary condition is specified at all boundaries except the outlet of the domain, where the Dirichlet boundary condition is used. The underlying turbulent flow is modeled using the k- ε model with a modified set of constants C μ = 0.09 , C ε 1 = 1.53 , and C ε 2 = 1.85 . A reduced chemical mechanism comprising 28 species and 24 chemical reactions for DME [30] derived from a full mechanism by Zhao et al. [31] is utilized. The computational domain is discretized into 120 finite-volume (FV) cells in both the axial and radial directions. At the steady state, approximately 0.6 million Lagrangian particles exist in the computational domain with an average of 50 particles per FV cell. Following our previous work on similar flames [19,20,27], the modified Curl’s model is used to describe the micromixing of the scalars with a major mixing time constant C ϕ = 3 and a minor time constant C m i n = 0.25 .

4. Results and Discussion

The numerical results are obtained by first running the simulations for 10 characteristic flow-through times to reach a statistically stationary condition; thereafter, statistics are collected over another 10 flow-through times. The numerical results reported herein are found to be independent of any further refinement of grid and particle resolutions. The variation in flame LOH with coflow temperatures is shown in Figure 1 for three different compositions of fuel/air mixtures in the jet. During the simulation, LOH is determined based on the first axial location where OH mass fraction (YOH) reaches 100 ppm at any radial location. Thereafter, both the average and the standard deviation of LOH are estimated from the time series of LOH. The standard deviation of LOH is also plotted against the average LOH in Figure 1. A comparison with available experimental measurements [23] shows a decent agreement with predicted LOH. As seen in Figure 1, flames with low coflow temperature display a large fluctuation in the LOH for all three inlet compositions. In the experiments, the LOH was determined as the lowest distance of the main connected flame base from the nozzle exit, obtained using high-speed CH*-chemiluminescence images [23]. For the pure DME case, as shown in Figure 1, the RANS-MMC simulations almost accurately predict the average LOH for both 1400 K and 1500 K coflow temperatures, while a slight overprediction is observed for the coflow temperature of 1350 K. The temporal evolution of LOH is shown in Figure 2. For the air/DME = 1:1 case, LOH strongly depends on the coflow temperature, as a reduction of 3.5 % in the coflow temperature (i.e., from 1400 K to 1350 K) almost doubles the LOH. Further, flames with low LOH (high coflow temperature) are found to be inherently stable, while flames with high LOH (low coflow temperature) are unstable, fluctuating, and noisy [4]. For the partially premixed cases, an underprediction of LOH is observed, which continues to increase with an increase in the degree of premixing at the jet inlet. The stochastic MMC model employed herein uses the mixture fraction as the major scalar, which demonstrates a good agreement for pure DME jets (case 1). The agreement of LOH with experimental data deteriorates with an increase in the degree of premixing in the fuel jet. This may be overcome using two major scalars—one which represents the mixing of two streams and another which is the reaction progress variable. Moreover, the discrepancy in LOH prediction could also arise due to the different methodologies adopted in experiments and numerical simulations. The non-availability of excited species like CH* and OH* in the reaction mechanism does not allow us to compare LOH one-to-one with experiments. Please note that in experiments the average LOH is determined based on the chemiluminescence images predominantly from CH*. The inclusion of excited species in reaction mechanisms is possible; however, it will increase the computational cost significantly.
The temporal evolution of OH contours is shown for a air/DME = 3:1 (case 3) flame at a coflow temperature of 1400 K in Figure 3. At a low coflow temperature of 1350 K, the ignition spots form and grow in size but cannot merge with the flame base, as is evident from the OH contours shown in the last figure of this manuscript. These random spots are generated early during the start of combustion. These random spots are located within a span of 10D–30D, which causes increased fluctuations in the LOH. However, for the coflow temperature of 1400 K case, the ignition kernels form, grow, advect, and eventually mix with the flame base, suggesting that the flame stabilization is primarily due to the propagation of partially premixed flames and, to a lesser extent, to autoignition for high coflow temperature.
Figure 4 shows the instantaneous scatter plots for OH and CH2O radicals in the mixture fraction space ( η ) at several axial locations for three different degrees of fuel/air premixing in the jet and at a coflow temperature of 1500 K. The scatter plots are generated after randomly identifying 10,000 unique samples from a dataset containing instantaneous data of all stochastic particles at an axial location accumulated over ten flow-through times. In Figure 4, an overlap between OH and CH2O radicals is discernible up to z/D = 7.5 across all scenarios. Beyond this point, a clear demarcation between OH and CH2O radicals emerges in the η -space at all downstream positions. Notably, the conditional peak of CH2O is situated on the fuel-rich side initially, but it gradually transitions towards richer mixtures as the flow progresses downstream. An intriguing observation is the gradual reduction in the extent of overlap between OH and CH2O radicals as the degree of premixing in the jet increases.
In diffusion flames, autoignition stands out as the primary stabilization mechanism. During autoignition, CH2O predominantly forms within the pre-flame region. As partial premixing becomes more pronounced, the spatial separation between OH and CH2O radicals becomes more pronounced as well. OH exhibits a wide distribution within the lean region, whereas CH2O progressively confines itself to the exceedingly rich region, particularly evident at high coflow temperatures.
In the case of partially premixed flames under high coflow temperatures, turbulent premixed flame propagation emerges as the dominant flame stabilization mechanism. Consequently, the formation of CH2O becomes somewhat restricted, as depicted in Figure 4 by the diminishing peak value of CH2O. This trend is also apparent in the contour plots. Similar findings are reported in the chemiluminescence study of the same flame.
The scatter plots of OH and CH2O are shown in Figure 5 at z/D = 40 for a low coflow temperature of 1350 K. For all degrees of premixing in the jet, at a low coflow temperature, OH and CH2O overlap with each other, and a broad distribution of CH2O in η -space is found. The instantaneous contours of OH and CH2O mass fractions are shown in Figure 6 and Figure 7 for case 1 and case 3. The contours of OH and CH2O for case 2 are not shown here due to brevity. The OH radical is used here as the marker to calculate the LOH. For the pure DME case, OH is mainly present in the mixing layer between the fuel jet and the hot coflow. On the other hand, CH2O is produced in the fuel-rich region of the jet, while CH2O concentration decreases significantly in the mixing region of jet and coflow. For partially premixed flames, at high coflow temperatures, the premixing of the fuel jet plays a dominant role in flame stabilization, while at low coflow temperatures the flame is primarily controlled by autoignition [23]. A somewhat similar phenomenon is observed from the simulations, which is evident from the contours of OH and CH2O mass fractions as shown in Figure 7. At the high coflow temperature of the 1500 K case, the formation of OH kernels begins near the jet inlet. For a particular coflow temperature, the formation of OH occurs near the centerline with an increase in the degree of premixing as more O2 becomes available in the fuel jet. Consequently, CH2O concentration decreases near the centerline, as is evident in Figure 6 and Figure 7. With a decrease in coflow temperature to 1350–1400 K, the formation of OH kernels occurs at farther-downstream locations. At a low coflow temperature of 1350 K, a broad distribution of CH2O appears near the centerline due to initial preheating before the high-heat-release region is formed, and flame stabilization is expected to be controlled by autoignition only. In the absence of autoignition, CH2O would be expected to be present adjacent to the OH layer and be consumed quickly. Thus, it may be concluded that flames at low coflow temperatures are mainly stabilized by autoignition. It is observed that for the 1350 K cases, random autoignition spots appear, as is evident from the OH contours for both pure DME and partially premixed DME/air jet cases.

5. Conclusions

The validation of lifted jet flames involving DME/air has been meticulously conducted through the application of the RANS-based stochastic MMC model. This comprehensive investigation seeks to elucidate the intricate interplay between coflow temperatures and the degree of premixing in the fuel jet, with a particular focus on their influence on the LOH and the underlying stabilization mechanism. Within this modeling framework, a single reference variable, namely, the mixture fraction, has been judiciously employed as the fundamental parameter on which particle mixing is conditioned. To refine the description of micromixing phenomena within the stochastic MMC approach, the modified Curl’s mixing model has been judiciously integrated. This model exhibits commendable efficacy in capturing the LOH for the pure DME jet, showcasing its capability to reproduce key flame characteristics. However, it is noteworthy that an observed underprediction of the LOH becomes apparent as the degree of premixing in the fuel jet increases. Moreover, the trend in the variation in LOH with coflow temperature is qualitatively captured. However, the predicted LOH values at low coflow temperatures differ from experiments. Further investigations are needed to enhance the capabilities of the modeling of partially premixed flames through additional reference variables, such as the reaction progress variable, in conjunction with the mixture fraction, to gain deeper insights into the phenomena under study. This nuanced behavior underscores the importance of considering varying degrees of fuel/air premixing when exploring the stabilization dynamics of lifted DME flames.
The numerical results gleaned from these simulations yield intriguing insights into the dominant factors governing flame stabilization under different operational conditions. Notably, the findings suggest that, at elevated coflow temperatures, the propagation of partially premixed flames assumes a paramount role in the overall stabilization process. This emphasizes the significance of understanding and characterizing the intricate balance between fuel and oxidizer streams, particularly in scenarios where the premixing of these components is prevalent.
Conversely, at lower coflow temperatures, the numerical results point towards autoigniting kernels emerging as the primary contributors to flame stabilization. This highlights the intricate interplay between autoignition processes and the prevailing environmental conditions, shaping the stabilization mechanisms in these complex combustion scenarios. Moreover, the investigation underscores that, in the context of pure DME jet flames, autoignition remains the dominant and prevailing mechanism governing flame stabilization.
In essence, this thorough examination employing the RANS-based stochastic MMC model provides a nuanced understanding of the multifaceted dynamics at play in lifted DME flames. As a pioneering application of this model to such flames, this study contributes valuable insights into the complex interplay between turbulent combustion, premixing, and autoignition, offering a stepping stone for further advancements in the understanding and modeling of intricate combustion processes in practical applications.

Author Contributions

Conceptualization, S.K.G. and S.D.; methodology, S.K.G. and S.D.; software, S.K.G. and S.D.; validation, R.G.; formal analysis, R.G.; investigation, R.G.; resources, S.D.; data curation, R.G.; writing—original draft preparation, S.K.G. and S.D.; writing—review and editing, S.K.G., R.G. and S.D.; visualization, S.K.G., R.G. and S.D.; supervision, S.K.G. and S.D.; project administration, S.D.; funding acquisition, S.D. All authors have read and agreed to the submitted version of the manuscript.

Funding

S.D. gratefully acknowledges the financial support from the SERB-INDIA Early Career Research Award (ECR/2016/001996). S.K.G. acknowledges the financial support of the Ministry of Education India (formerly known as MHRD) for doctoral research.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data can be obtained from the authors on reasonable request.

Acknowledgments

The authors would like to acknowledge the IITK computer center (https://www.iitk.ac.in/cc/) for providing the computing facility.

Conflicts of Interest

Rajat Gupta is currently employed at the Maruti Suzuki India Limited. The present study was conducted at IIT Kanpur. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Lyons, K.M. Toward an understanding of the stabilization mechanisms of lifted turbulent jet flames: Experiments. Prog. Energy Combust. Sci. 2007, 33, 211–231. [Google Scholar] [CrossRef]
  2. Ghai, S.K.; De, S. A Review on Autoignition in Laminar and Turbulent Nonpremixed Flames. In Combustion for Power Generation and Transportation: Technology, Challenges and Prospects; Agarwal, A.K., De, S., Pandey, A., Singh, A.P., Eds.; Springer: Singapore, 2017; pp. 11–37. [Google Scholar]
  3. Cao, R.R.; Pope, S.B.; Masri, A.R. Turbulent lifted flames in a vitiated coflow investigated using joint PDF calculations. Combust. Flame 2005, 142, 438–453. [Google Scholar] [CrossRef]
  4. Cabra, R.; Cabra, R.; Chen, J.; Dibble, R.W.; Karpetis, A.N.; Karpetis, A.N.; Barlow, R.S. Lifted methane–air jet flames in a vitiated coflow. Combust. Flame 2005, 143, 491–506. [Google Scholar] [CrossRef]
  5. Gkagkas, K.; Lindstedt, R.P. Transported PDF modelling with detailed chemistry of pre- and auto-ignition in CH4/air mixtures. Proc. Combust. Inst. 2007, 31, 1559–1566. [Google Scholar] [CrossRef]
  6. Jangi, M.; Zhao, X.; Haworth, D.C.; Bai, X.S. Stabilization and liftoff length of a non-premixed methane/air jet flame discharging into a high-temperature environment: An accelerated transported PDF method. Combust. Flame 2015, 162, 408–419. [Google Scholar] [CrossRef]
  7. Cabra, R.; Myhrvold, T.; Chen, J.; Dibble, R.; Karpetis, A.; Barlow, R. Simultaneous laser raman-rayleigh-lif measurements and numerical modeling results of a lifted turbulent H2/N2 jet flame in a vitiated coflow. Proc. Combust. Inst. 2002, 29, 1881–1888. [Google Scholar] [CrossRef]
  8. Gordon, R.L.; Masri, A.R.; Pope, S.B.; Goldin, G.M. Transport budgets in turbulent lifted flames of methane autoigniting in a vitiated co-flow. Combust. Flame 2007, 151, 495–511. [Google Scholar] [CrossRef]
  9. Gordon, R.; Stårner, S.; Masri, A.; Bilger, R. Further characterisation of lifted hydrogen and methane flames issuing into a vitiated coflow. In Proceedings of the 5th Asia-Pacific Conference on Combustion, Adelaide, Australia, 17–20 July 2005; pp. 333–336. [Google Scholar]
  10. Patwardhan, S.S.; De, S.; Lakshmisha, K.N.; Raghunandan, B.N. CMC simulations of lifted turbulent jet flame in a vitiated coflow. Proc. Combust. Inst. 2009, 32, 1705–1712. [Google Scholar] [CrossRef]
  11. Yadav, R.; Kushari, A.; De, A. Modeling of turbulent lifted flames in vitiated co-flow using multi environment Eulerian PDF transport approach. Int. J. Heat Mass Transf. 2014, 77, 230–246. [Google Scholar] [CrossRef]
  12. Navarro-Martinez, S.; Kronenburg, A. LES–CMC simulations of a lifted methane flame. Proc. Combust. Inst. 2009, 32, 1509–1516. [Google Scholar] [CrossRef]
  13. Roy, R.N.; Kumar, S.; Sreedhara, S. A new approach to model turbulent lifted CH4/air flame issuing in a vitiated coflow using conditional moment closure coupled with an extinction model. Combust. Flame 2014, 161, 197–209. [Google Scholar] [CrossRef]
  14. Michel, J.B.; Colin, O.; Angelberger, C.; Veynante, D. Using the tabulated diffusion flamelet model ADF-PCM to simulate a lifted methane–air jet flame. Combust. Flame 2009, 156, 1318–1331. [Google Scholar] [CrossRef]
  15. Domingo, P.; Vervisch, L.; Veynante, D. Large-eddy simulation of a lifted methane jet flame in a vitiated coflow. Combust. Flame 2008, 152, 415–432. [Google Scholar] [CrossRef]
  16. Enjalbert, N.; Domingo, P.; Vervisch, L. Mixing time-history effects in Large Eddy Simulation of non-premixed turbulent flames: Flow-Controlled Chemistry Tabulation. Combust. Flame 2012, 159, 336–352. [Google Scholar] [CrossRef]
  17. Schulz, O.; Jaravel, T.; Poinsot, T.; Cuenot, B.; Noiray, N. A criterion to distinguish autoignition and propagation applied to a lifted methane–air jet flame. Proc. Combust. Inst. 2017, 36, 1637–1644. [Google Scholar] [CrossRef]
  18. Han, W.; Raman, V.; Chen, Z. LES/PDF modeling of autoignition in a lifted turbulent flame: Analysis of flame sensitivity to differential diffusion and scalar mixing time-scale. Combust. Flame 2016, 171, 69–86. [Google Scholar] [CrossRef]
  19. Ghai, S.K.; De, S.; Kronenburg, A. Numerical simulations of turbulent lifted jet diffusion flames in a vitiated coflow using the stochastic multiple mapping conditioning approach. Proc. Combust. Inst. 2019, 37, 2199–2206. [Google Scholar] [CrossRef]
  20. Ghai, S.K.; De, S. Numerical investigation of auto-igniting turbulent lifted CH4/air jet diffusion flames in a vitiated co-flow using a RANS based stochastic multiple mapping conditioning approach. Combust. Flame 2019, 203, 362–374. [Google Scholar] [CrossRef]
  21. Straub, C.; De, S.; Kronenburg, A.; Vogiatzaki, K. The effect of timescale variation in multiple mapping conditioning mixing of PDF calculations for Sandia Flame series D–F. Combust. Theory Model. 2016, 20, 894–912. [Google Scholar] [CrossRef]
  22. Ghai, S.K.; De, S. Numerical investigation of flow and scalar fields of piloted, partially-premixed dimethyl ether/air jet flames using stochastic multiple mapping conditioning approach. Combust. Flame 2019, 208, 480–491. [Google Scholar] [CrossRef]
  23. Macfarlane, A.W.; Dunn, M.J.; Juddoo, M.; Masri, A.R. Stabilisation of turbulent auto-igniting dimethyl ether jet flames issuing into a hot vitiated coflow. Proc. Combust. Inst. 2017, 36, 1661–1668. [Google Scholar] [CrossRef]
  24. Ghai, S.K.; De, S.; Vogiatzaki, K.; Cleary, M.J. Theory and Application of Multiple Mapping Conditioning for Turbulent Reactive Flows. In Modeling and Simulation of Turbulent Combustion; Springer: Singapore, 2018; pp. 447–474. [Google Scholar]
  25. Klimenko, A.Y.; Pope, S.B. The modeling of turbulent reactive flows based on multiple mapping conditioning. Phys. Fluids 2003, 15, 1907–1925. [Google Scholar] [CrossRef]
  26. Janicka, J.; Kolbe, W.; Kollmann, W. Closure of the Transport Equation for the Probability Density Funcfion of Turbulent Scalar Fields. J. Non-Equilib. Thermody 1979, 4, 47–66. [Google Scholar] [CrossRef]
  27. Ghai, S.K.; De, S. Numerical modeling of turbulent premixed combustion using RANS based stochastic multiple mapping conditioning approach. Proc. Combust. Inst. 2019, 37, 2519–2526. [Google Scholar] [CrossRef]
  28. Nooren, P.A.; Wouters, H.A.; Peeters, T.; Roekaerts, D.J.; Maas, U.; Schmidt, D. Monte Carlo PDF modelling of a turbulent natural-gas diffusion flame. Combust. Theory Model. 1997, 1, 79–96. [Google Scholar] [CrossRef]
  29. Macfarlane, A.R.; Dunn, M.; Juddoo, M.; Masri, A. The evolution of autoignition kernels in turbulent flames of dimethyl ether. Combust. Flame 2018, 197, 182–196. [Google Scholar] [CrossRef]
  30. Chin, G.T.; Chen, J.Y.; Rapp, V.H.; Dibble, R.W. Development and Validation of a Reduced DME Mechanism Applicable to Various Combustion Modes in Internal Combustion Engines. J. Combust. 2011, 2011, 630580. [Google Scholar] [CrossRef]
  31. Zhao, Z.; Chaos, M.; Kazakov, A.F.; Dryer, F.L. Thermal decomposition reaction and a comprehensive kinetic model of dimethyl ether. Int. J. Chem. Kinet. 2008, 40, 1–18. [Google Scholar] [CrossRef]
Figure 1. Variation in average LOH vs. coflow temperature: (left) and standard deviation plotted against average LOH (right) for case 1 (top row), case 2 (middle row), and case 3 (bottom row). Asterisk symbols represent experimental mean, and square symbols represent standard deviation from measurements [23]. Blue dashed lines represents results from RANS-MMC.
Figure 1. Variation in average LOH vs. coflow temperature: (left) and standard deviation plotted against average LOH (right) for case 1 (top row), case 2 (middle row), and case 3 (bottom row). Asterisk symbols represent experimental mean, and square symbols represent standard deviation from measurements [23]. Blue dashed lines represents results from RANS-MMC.
Fluids 09 00125 g001
Figure 2. Temporal variation in LOH for different coflow temperatures: 1500 K (top), 1400 K (middle), 1350 K (bottom). Solid lines represent instantaneous LOH, and dash-dot lines represent average LOH.
Figure 2. Temporal variation in LOH for different coflow temperatures: 1500 K (top), 1400 K (middle), 1350 K (bottom). Solid lines represent instantaneous LOH, and dash-dot lines represent average LOH.
Fluids 09 00125 g002
Figure 3. Temporal evolution of OH contours of air/DME = 3:1 (case 3) flame at coflow temperature of 1400 K.
Figure 3. Temporal evolution of OH contours of air/DME = 3:1 (case 3) flame at coflow temperature of 1400 K.
Fluids 09 00125 g003
Figure 4. Instantaneous scatter plots of OH and CH2O for case 1 (left), case 2 (middle), and case 3 (right) at coflow temperature of 1500 K. Stoichiometric mixture fractions and ξ s t values for case 1, case 2, and case 3 are 0.06, 0.11, and 0.2.
Figure 4. Instantaneous scatter plots of OH and CH2O for case 1 (left), case 2 (middle), and case 3 (right) at coflow temperature of 1500 K. Stoichiometric mixture fractions and ξ s t values for case 1, case 2, and case 3 are 0.06, 0.11, and 0.2.
Fluids 09 00125 g004
Figure 5. Instantaneous scatter plots of OH and CH2O for coflow temperature of 1350 K.
Figure 5. Instantaneous scatter plots of OH and CH2O for coflow temperature of 1350 K.
Fluids 09 00125 g005
Figure 6. Instantaneous contours of OH (left) and CH2O (right) mass fraction for pure DME jet (case 1) at different coflow temperatures: 1500 K (left), 1400 K (middle), and 1350 K (right).
Figure 6. Instantaneous contours of OH (left) and CH2O (right) mass fraction for pure DME jet (case 1) at different coflow temperatures: 1500 K (left), 1400 K (middle), and 1350 K (right).
Fluids 09 00125 g006
Figure 7. Instantaneous contours of OH (left) and CH2O (right) mass fraction for 3:1 air/DME jet (case 3) at different coflow temperatures: 1500 K (left), 1400 K (middle), and 1350 K (right).
Figure 7. Instantaneous contours of OH (left) and CH2O (right) mass fraction for 3:1 air/DME jet (case 3) at different coflow temperatures: 1500 K (left), 1400 K (middle), and 1350 K (right).
Fluids 09 00125 g007
Table 1. Inlet boundary conditions are derived from experimental data from Macfarlane et al. [23]. Here, Y represents species mass fraction.
Table 1. Inlet boundary conditions are derived from experimental data from Macfarlane et al. [23]. Here, Y represents species mass fraction.
ParameterFuel JetCoflow
Case 1Case 2Case 3
Jet composition (by volume)Pure DMEAir/DME (1:1)Air/DME (3:1)
Temperature (K)298135014001500
Velocity (m/s)504
Reynolds number (Re)38,61022,99016,650- - -
Diameter (mm)4.45197
Equivalence ratio- - -0.360.380.43
Stoichiometric mixture fraction, ξ s t 0.060.110.20- - -
Y O 2 0.00.09410.15780.14590.1410.1307
Y N 2 0.00.29190.49480.75890.75850.7575
Y H 2 O 0.00.00.00.09520.10050.1118
Y CH 3 OCH 3 1.00.6140.34740.00.00.0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ghai, S.K.; Gupta, R.; De, S. Effects of Partial Premixing and Coflow Temperature on Flame Stabilization of Lifted Jet Flames of Dimethyl Ether in a Vitiated Coflow Based on Stochastic Multiple Mapping Conditioning Approach. Fluids 2024, 9, 125. https://doi.org/10.3390/fluids9060125

AMA Style

Ghai SK, Gupta R, De S. Effects of Partial Premixing and Coflow Temperature on Flame Stabilization of Lifted Jet Flames of Dimethyl Ether in a Vitiated Coflow Based on Stochastic Multiple Mapping Conditioning Approach. Fluids. 2024; 9(6):125. https://doi.org/10.3390/fluids9060125

Chicago/Turabian Style

Ghai, Sanjeev Kumar, Rajat Gupta, and Santanu De. 2024. "Effects of Partial Premixing and Coflow Temperature on Flame Stabilization of Lifted Jet Flames of Dimethyl Ether in a Vitiated Coflow Based on Stochastic Multiple Mapping Conditioning Approach" Fluids 9, no. 6: 125. https://doi.org/10.3390/fluids9060125

APA Style

Ghai, S. K., Gupta, R., & De, S. (2024). Effects of Partial Premixing and Coflow Temperature on Flame Stabilization of Lifted Jet Flames of Dimethyl Ether in a Vitiated Coflow Based on Stochastic Multiple Mapping Conditioning Approach. Fluids, 9(6), 125. https://doi.org/10.3390/fluids9060125

Article Metrics

Back to TopTop