Pub Date : 2024-09-30DOI: 10.1016/j.combustflame.2024.113749
Qun Hu, Lipo Wang
To overcome the modeling challenge from the coupling between liquid vaporization and chemical reaction in turbulent spray combustion, a similarity mapping approach is implemented to reduce the flamelet tabulation parameters. In the framework of Eulerian–Lagrangian multiphase large eddy simulations (LES), such a modeling idea is developed upon the conventional flamelet/progress variable model. The flamelet library is constructed from a series of quasi one-dimensional spray counterflow solutions, integrated with the multiple solution modes. Test cases, including the laminar spray counterflow flame and Sydney turbulent spray flame, indicate that this newly proposed model is in principle favorable to improve the numerical predictability with acceptable computational cost. Overall, the flame structure can be appropriately captured, showing better performance compared with some reported results.
Novelty and significance statement
A newly proposed similarity mapping spray flamelet/progress variable (SMFPV) model is implemented for turbulent spray combustion. In SMFPV, the number of entry parameters of the flamelet library is reasonably reduced and two-way coupling between flame and evaporation can be realized. Thus in principle, SMFPV is favorable to improve the numerical predictability with acceptable computational cost. Simulation results of test cases justify the modeling idea.
{"title":"Turbulent spray combustion modeling in reduced tabulation parameter space by similarity mapping","authors":"Qun Hu, Lipo Wang","doi":"10.1016/j.combustflame.2024.113749","DOIUrl":"10.1016/j.combustflame.2024.113749","url":null,"abstract":"<div><div>To overcome the modeling challenge from the coupling between liquid vaporization and chemical reaction in turbulent spray combustion, a similarity mapping approach is implemented to reduce the flamelet tabulation parameters. In the framework of Eulerian–Lagrangian multiphase large eddy simulations (LES), such a modeling idea is developed upon the conventional flamelet/progress variable model. The flamelet library is constructed from a series of quasi one-dimensional spray counterflow solutions, integrated with the multiple solution modes. Test cases, including the laminar spray counterflow flame and Sydney turbulent spray flame, indicate that this newly proposed model is in principle favorable to improve the numerical predictability with acceptable computational cost. Overall, the flame structure can be appropriately captured, showing better performance compared with some reported results.</div></div><div><h3>Novelty and significance statement</h3><div>A newly proposed similarity mapping spray flamelet/progress variable (SMFPV) model is implemented for turbulent spray combustion. In SMFPV, the number of entry parameters of the flamelet library is reasonably reduced and two-way coupling between flame and evaporation can be realized. Thus in principle, SMFPV is favorable to improve the numerical predictability with acceptable computational cost. Simulation results of test cases justify the modeling idea.</div></div>","PeriodicalId":280,"journal":{"name":"Combustion and Flame","volume":"270 ","pages":"Article 113749"},"PeriodicalIF":5.8,"publicationDate":"2024-09-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142358703","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-28DOI: 10.1016/j.combustflame.2024.113756
Ping Chen , Huichun Wang , Longxiang Qiao , Mingyan Gu , Kun Luo , Jianren Fan
Exploring the characteristics of NO formation during ammonia-coal co-combustion is crucial for achieving clean combustion. O2/CO2 atmosphere is an effective way to achieve low nitrogen combustion of pulverized coal, but the effect of O2/CO2 atmosphere on NO formation characteristics of ammonia-coal co-firing is still unclear. This work used high-temperature furnace experiments to investigate the conversion characteristics of fuel-N in ammonia-coal co-combustion under CO2 atmosphere and conducted quantum chemistry calculations to investigate the molecular pathways of fuel-N oxidation. The experimental results showed that the higher temperature and oxygen concentration results in earlier occurrence of the peak and completion of NO in ammonia-coal co-combustion. Compared to N2 atmosphere, CO2 has an inhibitory effect on NO generation in ammonia-coal co-firing and the degree of CO2 inhibition on the formation of NO per unit mass gradually decreases with increasing the temperature. The theoretical calculation results showed that CO2 reduces the energy barrier of the rate determining step for ammonia-N oxidation by about 133.81 kJ/mol, promoting the conversion of ammonia-N to nitrogen oxides. However, CO2 also has a certain degree of inhibitory effect on the coal-N oxidation, increasing the energy barrier of the rate determining step for NO generation from coal-N oxidation. The kinetic results showed that when the temperature is higher than 1673 K, CO2 has a more obvious promotion effect on ammonia-N oxidation and inhibition effect on coal-N oxidation, but a more significant inhibitory effect on coal-N oxidation, so that the CO2 atmosphere can inhibit the generation of nitrogen-containing products in ammonia-coal co-combustion to a certain extent.
探索氨煤共燃过程中 NO 的形成特征对于实现清洁燃烧至关重要。O2/CO2气氛是实现煤粉低氮燃烧的有效途径,但O2/CO2气氛对氨煤共烧NO形成特性的影响尚不清楚。本研究利用高温炉实验研究了 CO2 大气下氨煤共燃中燃料-N 的转化特性,并通过量子化学计算研究了燃料-N 氧化的分子途径。实验结果表明,温度和氧气浓度越高,氨-煤燃烧中 NO 的峰值出现和完成时间越早。与 N2 气氛相比,CO2 对氨煤共烧中 NO 的生成有抑制作用,且 CO2 对单位质量 NO 生成的抑制程度随温度升高而逐渐减小。理论计算结果表明,CO2 可使氨-N 氧化速率决定步骤的能量势垒降低约 133.81 kJ/mol,促进氨-N 向氮氧化物的转化。然而,CO2 对煤-N 氧化也有一定程度的抑制作用,增加了煤-N 氧化生成 NO 的速率决定步骤的能量势垒。动力学结果表明,当温度高于 1673 K 时,CO2 对氨-N 氧化有较明显的促进作用,对煤-N 氧化有抑制作用,但对煤-N 氧化有较明显的抑制作用,因此 CO2 大气能在一定程度上抑制氨-煤共燃中含氮产物的生成。
{"title":"Study of the NO formation characteristics and ammonia-N/coal-N transformation mechanism of ammonia-coal co-combustion in O2/CO2 atmosphere","authors":"Ping Chen , Huichun Wang , Longxiang Qiao , Mingyan Gu , Kun Luo , Jianren Fan","doi":"10.1016/j.combustflame.2024.113756","DOIUrl":"10.1016/j.combustflame.2024.113756","url":null,"abstract":"<div><div>Exploring the characteristics of NO formation during ammonia-coal co-combustion is crucial for achieving clean combustion. O<sub>2</sub>/CO<sub>2</sub> atmosphere is an effective way to achieve low nitrogen combustion of pulverized coal, but the effect of O<sub>2</sub>/CO<sub>2</sub> atmosphere on NO formation characteristics of ammonia-coal co-firing is still unclear. This work used high-temperature furnace experiments to investigate the conversion characteristics of fuel-N in ammonia-coal co-combustion under CO<sub>2</sub> atmosphere and conducted quantum chemistry calculations to investigate the molecular pathways of fuel-N oxidation. The experimental results showed that the higher temperature and oxygen concentration results in earlier occurrence of the peak and completion of NO in ammonia-coal co-combustion. Compared to N<sub>2</sub> atmosphere, CO<sub>2</sub> has an inhibitory effect on NO generation in ammonia-coal co-firing and the degree of CO<sub>2</sub> inhibition on the formation of NO per unit mass gradually decreases with increasing the temperature. The theoretical calculation results showed that CO<sub>2</sub> reduces the energy barrier of the rate determining step for ammonia-N oxidation by about 133.81 kJ/mol, promoting the conversion of ammonia-N to nitrogen oxides. However, CO<sub>2</sub> also has a certain degree of inhibitory effect on the coal-N oxidation, increasing the energy barrier of the rate determining step for NO generation from coal-N oxidation. The kinetic results showed that when the temperature is higher than 1673 K, CO<sub>2</sub> has a more obvious promotion effect on ammonia-N oxidation and inhibition effect on coal-N oxidation, but a more significant inhibitory effect on coal-N oxidation, so that the CO<sub>2</sub> atmosphere can inhibit the generation of nitrogen-containing products in ammonia-coal co-combustion to a certain extent.</div></div>","PeriodicalId":280,"journal":{"name":"Combustion and Flame","volume":"270 ","pages":"Article 113756"},"PeriodicalIF":5.8,"publicationDate":"2024-09-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142358701","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-28DOI: 10.1016/j.combustflame.2024.113753
Xiaoyuan Zhang , Maxence Lailliau , Yuyang Li , Yumeng Zhu , Zehua Feng , Wei Li , Philippe Dagaut
Ketene (CH2CO) mechanism is a building block for developing combustion kinetic models of practical fuels. To revisit the combustion chemistry related to ketene, oxidation experiments of butane-2,3‑dione (diacetyl, CH3COCOCH3), considered as an effective precursor of CH2CO, are conducted in a jet-stirred reactor (JSR) at 10 bar and temperatures ranging from 650 to 1160 K. Identification and quantification of intermediates are achieved by Fourier transform infrared spectrometry, gas chromatography, and mass spectrometry. A kinetic model of diacetyl is constructed based on recent theoretical and modeling studies on diacetyl and ketene, which has been validated against the present data and experimental data of diacetyl and CH2CO in literature. Generally, the present model can adequately predict most of them, and better predict the methyl-related intermediates under wide pyrolysis and combustion conditions than previous models. Based on modeling analyses, the unimolecular decomposition reaction of diacetyl is the dominant reaction pathway for fuel consumption under different equivalence ratio conditions, especially at high temperatures. Under lean conditions, both the H-atom abstraction reactions by methyl (i.e. CH3COCOCH3 + CH3 = CH4 + CH2CO + CH3CO, R3) and by OH (i.e. CH3COCOCH3 + OH = H2O + CH2CO + CH3CO, R5) are important for diacetyl consumption, while under rich conditions R5 becomes negligible. As the most important intermediates in diacetyl oxidation, the main consumption pathways of CH2CO and CH3 are dependent on the equivalence ratio conditions. Under lean conditions, CH2CO mainly reacts with OH to produce CH2OH and CO (i.e. CH2CO + OH = CH2OH + CO, R10), while methyl reacts with HO2 to produce CH3O and OH (i.e. CH3 + HO2 = CH3O + OH, R20). In contrast, under rich conditions, the addition-elimination reaction between CH2CO and H becomes competitive with R10, while the CH3 self-combination producing C2H6 plays a more important role than the CH3 oxidation pathway R20. Sensitivity analysis of CH2CO shows that not only the reactions of CH2CO, but also those of CH3 are sensitive to CH2CO formation. This is because CH3 related reactions influence the distribution of radical pool, which determines the oxidation reactivity of the reaction system.
{"title":"Oxidation of butane-2,3-dione at high pressure: Implications for ketene chemistry","authors":"Xiaoyuan Zhang , Maxence Lailliau , Yuyang Li , Yumeng Zhu , Zehua Feng , Wei Li , Philippe Dagaut","doi":"10.1016/j.combustflame.2024.113753","DOIUrl":"10.1016/j.combustflame.2024.113753","url":null,"abstract":"<div><div>Ketene (CH<sub>2</sub>CO) mechanism is a building block for developing combustion kinetic models of practical fuels. To revisit the combustion chemistry related to ketene, oxidation experiments of butane-2,3‑dione (diacetyl, CH<sub>3</sub>COCOCH<sub>3</sub>), considered as an effective precursor of CH<sub>2</sub>CO, are conducted in a jet-stirred reactor (JSR) at 10 bar and temperatures ranging from 650 to 1160 K. Identification and quantification of intermediates are achieved by Fourier transform infrared spectrometry, gas chromatography, and mass spectrometry. A kinetic model of diacetyl is constructed based on recent theoretical and modeling studies on diacetyl and ketene, which has been validated against the present data and experimental data of diacetyl and CH<sub>2</sub>CO in literature. Generally, the present model can adequately predict most of them, and better predict the methyl-related intermediates under wide pyrolysis and combustion conditions than previous models. Based on modeling analyses, the unimolecular decomposition reaction of diacetyl is the dominant reaction pathway for fuel consumption under different equivalence ratio conditions, especially at high temperatures. Under lean conditions, both the H-atom abstraction reactions by methyl (i.e. CH<sub>3</sub>COCOCH<sub>3</sub> + CH<sub>3</sub> = CH<sub>4</sub> + CH<sub>2</sub>CO + CH<sub>3</sub>CO, R3) and by OH (i.e. CH<sub>3</sub>COCOCH<sub>3</sub> + OH = H<sub>2</sub>O + CH<sub>2</sub>CO + CH<sub>3</sub>CO, R5) are important for diacetyl consumption, while under rich conditions R5 becomes negligible. As the most important intermediates in diacetyl oxidation, the main consumption pathways of CH<sub>2</sub>CO and CH<sub>3</sub> are dependent on the equivalence ratio conditions. Under lean conditions, CH<sub>2</sub>CO mainly reacts with OH to produce CH<sub>2</sub>OH and CO (i.e. CH<sub>2</sub>CO + OH = CH<sub>2</sub>OH + CO, R10), while methyl reacts with HO<sub>2</sub> to produce CH<sub>3</sub>O and OH (i.e. CH<sub>3</sub> + HO<sub>2</sub> = CH<sub>3</sub>O + OH, R20). In contrast, under rich conditions, the addition-elimination reaction between CH<sub>2</sub>CO and H becomes competitive with R10, while the CH<sub>3</sub> self-combination producing C<sub>2</sub>H<sub>6</sub> plays a more important role than the CH<sub>3</sub> oxidation pathway R20. Sensitivity analysis of CH<sub>2</sub>CO shows that not only the reactions of CH<sub>2</sub>CO, but also those of CH<sub>3</sub> are sensitive to CH<sub>2</sub>CO formation. This is because CH<sub>3</sub> related reactions influence the distribution of radical pool, which determines the oxidation reactivity of the reaction system.</div></div>","PeriodicalId":280,"journal":{"name":"Combustion and Flame","volume":"270 ","pages":"Article 113753"},"PeriodicalIF":5.8,"publicationDate":"2024-09-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142327801","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-28DOI: 10.1016/j.combustflame.2024.113745
N. Barléon , D.A. Lacoste , A.M. Alkhalifa , O. Vermorel , B. Cuenot
This study investigates the response of a laminar methane-air flame to Nanosecond Repetitively Pulsed (NRP) discharges in a canonical wall-stabilized burner using a combined experimental and numerical approach. The flow and flame behaviors were modeled using Direct Numerical Simulation (DNS) with an Analytically Reduced Chemistry for a precise chemical description. A phenomenological model incorporating detailed plasma kinetics and experimental observations was developed to simulate plasma effects. Zero-dimensional plasma reactor simulations were used to build up a reduced-order model describing discharge energy distribution in the specific conditions studied. Experimental measurements of electrical profiles identified two discharge regimes: a low-energy Corona discharge and a higher-energy Glow discharge, characterized by distinct spatial energy distributions. Experimental flame response analysis revealed three major phases: marginal response up to 100 pulses, a downstream shift of the flame tip, and stabilization after 400 pulses. Numerical simulations indicated that the Corona regime is crucial for explaining initial flame responses, while the Glow regime influences later stages. Adjustments in the Vibrational–Translational (VT) energy relaxation time and energy deposition ratios between fresh and burnt gases were necessary to match experimental observations. Additionally, an accurate modeling of the transient and steady-state flame responses requires integrating both the specificity of the Corona and the Glow discharge regimes. Future work should focus on measuring or theoretically calculating N2(v) relaxation times in CH4-H2O-CO2 mixtures and analyzing the spatial energy distribution of discharges interacting with flames to enhance plasma-combustion coupled models.
Novelty and significance
In this work, a phenomenological plasma-assisted combustion model has been developed, to investigate a laminar premixed stagnation plate burner, focusing on VT energy relaxation time and spatio-temporal energy distribution modeling. For the first time, not only O2, N2 and O but also the fuel and combustion intermediates and products have been considered in the VT relaxation model. It revealed their strong influence on the overall flame response in a case where the discharge crosses a flame front, highlighting the strong beneficial effect of energy deposited in vibrational form. The study questions and investigates the energy distribution from fresh to burnt gases, challenging the conventional uniform energy distribution assumption. The experimental identification of two specific plasma regimes was necessary to predict the transient flame response. Additionally, energy deposited downstream of the flame, in fresh gases, was found to more efficient than in hot gases to enhance combustion.
{"title":"Numerical investigation of lean methane flame response to NRP discharges actuation","authors":"N. Barléon , D.A. Lacoste , A.M. Alkhalifa , O. Vermorel , B. Cuenot","doi":"10.1016/j.combustflame.2024.113745","DOIUrl":"10.1016/j.combustflame.2024.113745","url":null,"abstract":"<div><div>This study investigates the response of a laminar methane-air flame to Nanosecond Repetitively Pulsed (NRP) discharges in a canonical wall-stabilized burner using a combined experimental and numerical approach. The flow and flame behaviors were modeled using Direct Numerical Simulation (DNS) with an Analytically Reduced Chemistry for a precise chemical description. A phenomenological model incorporating detailed plasma kinetics and experimental observations was developed to simulate plasma effects. Zero-dimensional plasma reactor simulations were used to build up a reduced-order model describing discharge energy distribution in the specific conditions studied. Experimental measurements of electrical profiles identified two discharge regimes: a low-energy Corona discharge and a higher-energy Glow discharge, characterized by distinct spatial energy distributions. Experimental flame response analysis revealed three major phases: marginal response up to 100 pulses, a downstream shift of the flame tip, and stabilization after 400 pulses. Numerical simulations indicated that the Corona regime is crucial for explaining initial flame responses, while the Glow regime influences later stages. Adjustments in the Vibrational–Translational (VT) energy relaxation time and energy deposition ratios between fresh and burnt gases were necessary to match experimental observations. Additionally, an accurate modeling of the transient and steady-state flame responses requires integrating both the specificity of the Corona and the Glow discharge regimes. Future work should focus on measuring or theoretically calculating N<sub>2</sub>(v) relaxation times in CH<sub>4</sub>-H<sub>2</sub>O-CO<sub>2</sub> mixtures and analyzing the spatial energy distribution of discharges interacting with flames to enhance plasma-combustion coupled models.</div><div><strong>Novelty and significance</strong></div><div>In this work, a phenomenological plasma-assisted combustion model has been developed, to investigate a laminar premixed stagnation plate burner, focusing on VT energy relaxation time and spatio-temporal energy distribution modeling. For the first time, not only O<sub>2</sub>, N<sub>2</sub> and O but also the fuel and combustion intermediates and products have been considered in the VT relaxation model. It revealed their strong influence on the overall flame response in a case where the discharge crosses a flame front, highlighting the strong beneficial effect of energy deposited in vibrational form. The study questions and investigates the energy distribution from fresh to burnt gases, challenging the conventional uniform energy distribution assumption. The experimental identification of two specific plasma regimes was necessary to predict the transient flame response. Additionally, energy deposited downstream of the flame, in fresh gases, was found to more efficient than in hot gases to enhance combustion.</div></div>","PeriodicalId":280,"journal":{"name":"Combustion and Flame","volume":"270 ","pages":"Article 113745"},"PeriodicalIF":5.8,"publicationDate":"2024-09-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142327735","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The nanosecond (NS) pulsed discharge for ammonia combustion assistance is numerically investigated with a special focus on the effects of discharge pretreatment uniformity. A method to mimic spatially non-uniform discharge pretreatment is established. Using NH3/O2/He and NH3/O2/N2 mixtures as examples, the ignition delay time (IDT), the laminar burning velocity (LBV), and the extinction strain rate (ESR) of mixtures pretreated by different discharge pretreatment methods were analyzed under a wide range of pretreatment uniformity. The results indicate that for the premixed NH3/O2/He (or N2) pretreated by NS discharge pulses, with the increase of the pretreatment non-uniformity, the IDT initially increases slightly, reaching a maximum and then decreases rapidly by 80 %, while the LBV and the ESR first decreases slowly, reaching a minimum and then increases rapidly by 10 %. The initial increment of IDT with the non-uniformity can be attributed to the deterioration of discharge-induced chemical effect but the rapid decrease of IDT as the discharge pretreatment becomes highly non-uniform is largely determined by the combustion-related chemical effect. The effects of pretreatment uniformity on LBV and ESR can be largely due to the thermal effect. When the oxidant is pretreated by the NS discharge, the minimal IDT, the maximal LBV and ESR can be observed under uniform conditions owing to the chemical effects and thus a uniform discharge pretreatment is preferred. Furthermore, a comparative analysis of the three pretreatment methods under identical energy input reveals that under uniform pretreatment, the direct pretreatment of pure NH3 is optimal, since it can generate more H2 to enhance the combustion. However, under highly non-uniform conditions, the discharge pretreatment of premixed fuel/oxidant mixtures is the most efficient for combustion enhancement.
{"title":"Numerical study on the pretreatment of ammonia by nanosecond pulsed discharge for combustion enhancement: Effects of pretreatment uniformity","authors":"Juntao Ao, Chengdong Kong, Yu Wang, Xiaojiang Wu, Zhongxiao Zhang","doi":"10.1016/j.combustflame.2024.113752","DOIUrl":"10.1016/j.combustflame.2024.113752","url":null,"abstract":"<div><div>The nanosecond (NS) pulsed discharge for ammonia combustion assistance is numerically investigated with a special focus on the effects of discharge pretreatment uniformity. A method to mimic spatially non-uniform discharge pretreatment is established. Using NH<sub>3</sub>/O<sub>2</sub>/He and NH<sub>3</sub>/O<sub>2</sub>/N<sub>2</sub> mixtures as examples, the ignition delay time (IDT), the laminar burning velocity (LBV), and the extinction strain rate (ESR) of mixtures pretreated by different discharge pretreatment methods were analyzed under a wide range of pretreatment uniformity. The results indicate that for the premixed NH<sub>3</sub>/O<sub>2</sub>/He (or N<sub>2</sub>) pretreated by NS discharge pulses, with the increase of the pretreatment non-uniformity, the IDT initially increases slightly, reaching a maximum and then decreases rapidly by 80 %, while the LBV and the ESR first decreases slowly, reaching a minimum and then increases rapidly by 10 %. The initial increment of IDT with the non-uniformity can be attributed to the deterioration of discharge-induced chemical effect but the rapid decrease of IDT as the discharge pretreatment becomes highly non-uniform is largely determined by the combustion-related chemical effect. The effects of pretreatment uniformity on LBV and ESR can be largely due to the thermal effect. When the oxidant is pretreated by the NS discharge, the minimal IDT, the maximal LBV and ESR can be observed under uniform conditions owing to the chemical effects and thus a uniform discharge pretreatment is preferred. Furthermore, a comparative analysis of the three pretreatment methods under identical energy input reveals that under uniform pretreatment, the direct pretreatment of pure NH<sub>3</sub> is optimal, since it can generate more H<sub>2</sub> to enhance the combustion. However, under highly non-uniform conditions, the discharge pretreatment of premixed fuel/oxidant mixtures is the most efficient for combustion enhancement.</div></div>","PeriodicalId":280,"journal":{"name":"Combustion and Flame","volume":"270 ","pages":"Article 113752"},"PeriodicalIF":5.8,"publicationDate":"2024-09-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142358700","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The ignition and combustion processes of boron particles are crucial to achieving high combustion efficiency in solid-fuel ramjet engines, particularly under reduced static pressure conditions in the secondary combustion chamber. This study carried out ignition and combustion experiments on amorphous boron particles with an average size of ∼3 µm, using a Hencken multi-diffusion flat flame burner under controlled pressures (0.3–1 atm) and temperatures (1900–2200 K). Optical measurements were utilized to qualify the ignition time, which increases with decreasing pressure. Then, the ignition and combustion models of small-size boron particles under reduced pressure were established, for which the Langmuir layer was introduced to calculate heat and mass transfer between particles and surrounding gases in the transition regime. Results showed that the Langmuir layer significantly lowered heat transfer rate, increasing the ignition time. A comparison of the heat fluxes of evaporation, heterogeneous reaction, and heat convection for the ∼3 µm boron particle demonstrated that the ignition process was limited by heat convection, while the combustion process is dominated by the heterogeneous reaction near the particle surface.
{"title":"Ignition and combustion characteristics of boron particles under reduced pressure","authors":"Ying Feng , Yong Tang , Dingjiang Xie , Wei Dong , Majie Zhao , Zhiwen Wu , Baolu Shi","doi":"10.1016/j.combustflame.2024.113733","DOIUrl":"10.1016/j.combustflame.2024.113733","url":null,"abstract":"<div><div>The ignition and combustion processes of boron particles are crucial to achieving high combustion efficiency in solid-fuel ramjet engines, particularly under reduced static pressure conditions in the secondary combustion chamber. This study carried out ignition and combustion experiments on amorphous boron particles with an average size of ∼3 µm, using a Hencken multi-diffusion flat flame burner under controlled pressures (0.3–1 atm) and temperatures (1900–2200 K). Optical measurements were utilized to qualify the ignition time, which increases with decreasing pressure. Then, the ignition and combustion models of small-size boron particles under reduced pressure were established, for which the Langmuir layer was introduced to calculate heat and mass transfer between particles and surrounding gases in the transition regime. Results showed that the Langmuir layer significantly lowered heat transfer rate, increasing the ignition time. A comparison of the heat fluxes of evaporation, heterogeneous reaction, and heat convection for the ∼3 µm boron particle demonstrated that the ignition process was limited by heat convection, while the combustion process is dominated by the heterogeneous reaction near the particle surface.</div></div>","PeriodicalId":280,"journal":{"name":"Combustion and Flame","volume":"270 ","pages":"Article 113733"},"PeriodicalIF":5.8,"publicationDate":"2024-09-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142358698","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-28DOI: 10.1016/j.combustflame.2024.113754
Zhen-Yang Zhou , Hong-Hao Ma , Lu-Qing Wang
Velocity deficits, which exist universally in the process of detonation propagation, were not considered in the problem of detonation Mach reflection in previous studies. Thus, an experimental study of the Mach reflections of gaseous detonations with velocity deficits is reported in this paper. Remarkable velocity deficits were achieved by a porous wall channel. The soot foil technique was utilized to monitor the cellular pattern variation and the evolution of triple-point trajectory. Two detonable mixtures, i.e., 2H2+O2+3Ar and 2H2+O2, were used to obtain stable detonation and unstable detonation. The results show that along with velocity deficits, the detonation tends to be unstable and the cell size increases. The induction length must be taken into consideration since the specific heat ratio could be different from that of a Chapman-Jouguet (CJ) detonation. Therefore, three modes of self-similarity were determined, i.e., non-self-similarity with significant velocity deficits (mode 1), global self-similarity with moderate velocity deficits (mode 2) and localized self-similarity with slight velocity deficits (mode 3). The Mach stem height of mode 1 is higher than the non-reactive three-shock theory. The triple-point trajectory of mode 2 corresponds to that of an inert shock. In the case of mode 3, the well-known frozen and equilibrium limits in detonation Mach reflections can be found. Considering the measured cell size () as the length scale, the required transition distance for the occurrence of the equilibrium limit was found to be approximately (8–10). A threshold value =0.95 was found above which the Mach reflection characteristics of a detonation agree well with those of a CJ detonation ( is the overdrive degree of the incident wave). Below the threshold value, the triple-point trajectory angle decreases linearly with the increase of , and the Mach stem overdrive degree is a constant.
{"title":"Mach reflection of detonation waves with velocity deficits","authors":"Zhen-Yang Zhou , Hong-Hao Ma , Lu-Qing Wang","doi":"10.1016/j.combustflame.2024.113754","DOIUrl":"10.1016/j.combustflame.2024.113754","url":null,"abstract":"<div><div>Velocity deficits, which exist universally in the process of detonation propagation, were not considered in the problem of detonation Mach reflection in previous studies. Thus, an experimental study of the Mach reflections of gaseous detonations with velocity deficits is reported in this paper. Remarkable velocity deficits were achieved by a porous wall channel. The soot foil technique was utilized to monitor the cellular pattern variation and the evolution of triple-point trajectory. Two detonable mixtures, i.e., 2H<sub>2</sub>+O<sub>2</sub>+3Ar and 2H<sub>2</sub>+O<sub>2</sub>, were used to obtain stable detonation and unstable detonation. The results show that along with velocity deficits, the detonation tends to be unstable and the cell size increases. The induction length must be taken into consideration since the specific heat ratio could be different from that of a Chapman-Jouguet (CJ) detonation. Therefore, three modes of self-similarity were determined, i.e., non-self-similarity with significant velocity deficits (<em>mode 1</em>), global self-similarity with moderate velocity deficits (<em>mode 2)</em> and localized self-similarity with slight velocity deficits (<em>mode 3</em>). The Mach stem height of <em>mode 1</em> is higher than the non-reactive three-shock theory. The triple-point trajectory of <em>mode 2</em> corresponds to that of an inert shock. In the case of <em>mode 3</em>, the well-known frozen and equilibrium limits in detonation Mach reflections can be found. Considering the measured cell size (<span><math><mi>λ</mi></math></span>) as the length scale, the required transition distance for the occurrence of the equilibrium limit was found to be approximately (8–10)<span><math><mi>λ</mi></math></span>. A threshold value <span><math><msub><mi>α</mi><mi>i</mi></msub></math></span>=0.95 was found above which the Mach reflection characteristics of a detonation agree well with those of a CJ detonation (<span><math><msub><mi>α</mi><mi>i</mi></msub></math></span> is the overdrive degree of the incident wave). Below the threshold value, the triple-point trajectory angle decreases linearly with the increase of <span><math><msub><mi>α</mi><mi>i</mi></msub></math></span>, and the Mach stem overdrive degree is a constant.</div></div>","PeriodicalId":280,"journal":{"name":"Combustion and Flame","volume":"270 ","pages":"Article 113754"},"PeriodicalIF":5.8,"publicationDate":"2024-09-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142358699","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-27DOI: 10.1016/j.combustflame.2024.113739
P. Prabhudeva, Srinibas Karmakar
Enhancing the energy density of hydrocarbon fuels is indeed crucial, especially in the context of addressing energy needs. Incorporating energetic particles as additives into conventional liquid fuels garners increased attention due to its potential to enhance energy density. Due to its high energy potential, boron is the most promising additive in several energetic particles. However, the study focusing on the boron-loaded slurry fuels are limited. Therefore, the present work focuses on the effect of boron loading on the spray combustion of boron-loaded slurry fuel, especially at higher loadings (10%, 20% & 30% boron + Jet A-1 fuel) using a swirl-stabilized spray combustor. The research scrutinizes the feed and exhaust particles' surface morphology, chemical composition and active boron content through diverse particle characterization techniques. Furthermore, the study investigates the influence of boron loading on the combustion characteristics of slurry fuel using methods such as colour flame imaging, spectroscopy, and BO2* Chemiluminescence imaging. The impact of boron loading on the positive thermal contribution of the slurry was examined by analyzing temperature measurements at various locations within the combustor. The BO2* chemiluminescence imaging and spectroscopy results indicate that the intensity of BO2 emission increases with a rise in particle loading. The temperature results reveal an increase in temperature compared to pure Jet A-1 for all the boron-loaded cases. In particular, JB10 and JB20 achieved nearly 90% and 85% of the theoretical temperature increase. The diffractogram and thermogravimetric results of the burnt particles collected show that the particles were burnt thoroughly for the boron-loaded cases.
{"title":"Spray combustion characteristics of boron slurry fuel in high particle loading conditions","authors":"P. Prabhudeva, Srinibas Karmakar","doi":"10.1016/j.combustflame.2024.113739","DOIUrl":"10.1016/j.combustflame.2024.113739","url":null,"abstract":"<div><div>Enhancing the energy density of hydrocarbon fuels is indeed crucial, especially in the context of addressing energy needs. Incorporating energetic particles as additives into conventional liquid fuels garners increased attention due to its potential to enhance energy density. Due to its high energy potential, boron is the most promising additive in several energetic particles. However, the study focusing on the boron-loaded slurry fuels are limited. Therefore, the present work focuses on the effect of boron loading on the spray combustion of boron-loaded slurry fuel, especially at higher loadings (10%, 20% & 30% boron + Jet A-1 fuel) using a swirl-stabilized spray combustor. The research scrutinizes the feed and exhaust particles' surface morphology, chemical composition and active boron content through diverse particle characterization techniques. Furthermore, the study investigates the influence of boron loading on the combustion characteristics of slurry fuel using methods such as colour flame imaging, spectroscopy, and BO<sub>2</sub>* Chemiluminescence imaging. The impact of boron loading on the positive thermal contribution of the slurry was examined by analyzing temperature measurements at various locations within the combustor. The BO<sub>2</sub>* chemiluminescence imaging and spectroscopy results indicate that the intensity of BO<sub>2</sub> emission increases with a rise in particle loading. The temperature results reveal an increase in temperature compared to pure Jet A-1 for all the boron-loaded cases. In particular, JB10 and JB20 achieved nearly 90% and 85% of the theoretical temperature increase. The diffractogram and thermogravimetric results of the burnt particles collected show that the particles were burnt thoroughly for the boron-loaded cases.</div></div>","PeriodicalId":280,"journal":{"name":"Combustion and Flame","volume":"270 ","pages":"Article 113739"},"PeriodicalIF":5.8,"publicationDate":"2024-09-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142327800","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The soot formation of n-heptane/iso-octane doped 2,5-dimethylfuran (DMF) with different ratios was experimentally investigated through laminar diffusion flame and numerically simulated by using the pyrolysis simulations method of reactive force field molecular dynamics (ReaxFF MD). The results showed that the laminar diffusion flame height of n-heptane/iso-octane increased with increasing DMF doping ratio. At a consistent sampling height, the primary soot particle sampling numbers collected were first decreasing and then increasing, and the size of the primary soot particle was first increasing and then decreasing accompanied by the DMF doping ratio increasing. Furthermore, the transmission electron microscope (TEM) analysis provided that DMF initially inhibited and subsequently promoted the primary soot particle growth in the n-heptane/iso-octane flame. The core-shell ratio increased and then decreased with the increase of the DMF doping ratio, which indicated that the maturity of soot decreased and then increased. In n-heptane/iso-octane doped with DMF ReaxFF MD pyrolysis simulation, it was divided into three stages 0–0.5 ns (the first stage), 0.5–3 ns (the second stage), and 3–4.5 ns (the third stage). Fuel decomposed in the first stage and reacted violently in the second stage. The soot precursors continued to react, molecule size kept to increase, and the polycyclic aromatic hydrocarbons (PAHs) continued to grow. In the third stage, the carbon number of the largest molecules were slowly increasing, and the development of soot entered mature stage, where the H/C ratio was slowly decreasing. Through ReaxFF MD simulations, it was found that the maturity of soot exhibited decreasing and then increasing trend with the increase of DMF doping ratio. The graphene-like structures were mainly concentrated around the H/C ratio of 0.297 region.
{"title":"Analysis of the nanostructure evolution of soot in n-heptane/iso-octane with 2,5-dimethylfuran addition: A combined experimental study and ReaxFF MD simulations","authors":"Wenlong Dong, Run Hong, Yuhang Yang, Dongyang Wang, Bingbing Qiu, Huaqiang Chu","doi":"10.1016/j.combustflame.2024.113751","DOIUrl":"10.1016/j.combustflame.2024.113751","url":null,"abstract":"<div><div>The soot formation of <em>n-</em>heptane/iso-octane doped 2,5-dimethylfuran (DMF) with different ratios was experimentally investigated through laminar diffusion flame and numerically simulated by using the pyrolysis simulations method of reactive force field molecular dynamics (ReaxFF MD). The results showed that the laminar diffusion flame height of <em>n-</em>heptane/iso-octane increased with increasing DMF doping ratio. At a consistent sampling height, the primary soot particle sampling numbers collected were first decreasing and then increasing, and the size of the primary soot particle was first increasing and then decreasing accompanied by the DMF doping ratio increasing. Furthermore, the transmission electron microscope (TEM) analysis provided that DMF initially inhibited and subsequently promoted the primary soot particle growth in the <em>n-</em>heptane/iso-octane flame. The core-shell ratio increased and then decreased with the increase of the DMF doping ratio, which indicated that the maturity of soot decreased and then increased. In <em>n-</em>heptane/iso-octane doped with DMF ReaxFF MD pyrolysis simulation, it was divided into three stages 0–0.5 ns (the first stage), 0.5–3 ns (the second stage), and 3–4.5 ns (the third stage). Fuel decomposed in the first stage and reacted violently in the second stage. The soot precursors continued to react, molecule size kept to increase, and the polycyclic aromatic hydrocarbons (PAHs) continued to grow. In the third stage, the carbon number of the largest molecules were slowly increasing, and the development of soot entered mature stage, where the H/C ratio was slowly decreasing. Through ReaxFF MD simulations, it was found that the maturity of soot exhibited decreasing and then increasing trend with the increase of DMF doping ratio. The graphene-like structures were mainly concentrated around the H/C ratio of 0.297 region.</div></div>","PeriodicalId":280,"journal":{"name":"Combustion and Flame","volume":"270 ","pages":"Article 113751"},"PeriodicalIF":5.8,"publicationDate":"2024-09-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142322577","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This work finds that with the assistance of filamentary plasma discharge, the ammonia combustion becomes very complicated with multiple modes, especially including a low-frequency (<2 Hz) oscillating combustion mode. This flame oscillation, which is sensitive to the discharge characteristics and the flow/fuel/boundary conditions, can be further categorized as three sub-modes denoted as Mode C1, C2 and C3. Their transitions occur depending on the flow rate (Qtot), the equivalence ratio (ϕ) and the temperature of the chamber wall. Detailed dynamic characteristics of flame kernels and flow fields during the oscillating process are further acquired to confirm the changing recirculation zone and the Karlovitz (Ka) number. A phenomenological mechanism based on recirculated energy feedback is thus proposed to explain the oscillating flame. In addition, the S-curve is used to elucidate the plasma impacts on the complicated combustion characteristics, especially the oscillating combustion modes and their transitions. It is found that the folded S-curve can be stretched with the help of plasma discharge. By considering the energy feedback effects, the S-curve can become more stretched with the Da number for extinction (i.e. ) being larger than that for ignition (i.e. ). Meanwhile, the branch with the Da number between and is unstable and the multiple oscillating combustion modes (i.e. Mode C1-C3) can be categorized on the unstable branch of the stretched S-curve.
这项研究发现,在丝状等离子体放电的帮助下,氨燃烧变得非常复杂,具有多种模式,特别是包括低频(2 赫兹)振荡燃烧模式。这种火焰振荡对放电特性和流量/燃料/边界条件非常敏感,可进一步分为三种子模式,分别称为模式 C1、C2 和 C3。它们的转换取决于流量(Qtot)、等效比(j)和腔壁温度。在振荡过程中,进一步获取了火焰内核和流场的详细动态特性,以确认不断变化的再循环区和卡尔洛维茨(Ka)数。因此,提出了一种基于再循环能量反馈的现象学机制来解释振荡火焰。此外,S 曲线还用于阐明等离子体对复杂燃烧特性的影响,特别是振荡燃烧模式及其转换。研究发现,在等离子体放电的帮助下,折叠的 S 曲线可以被拉伸。考虑到能量反馈效应,S 曲线的拉伸程度越大,熄灭的 Da 值(即 DaE′)就越大于点火的 Da 值(即 DaI′)。同时,Da 数介于 DaI′ 和 DaE′ 之间的分支是不稳定的,多重振荡燃烧模式(即模式 C1-C3)可归类于拉伸 S 曲线的不稳定分支上。
{"title":"A phenomenological understanding of the multimodal low-frequency oscillating combustion of ammonia induced by filamentary plasma discharge","authors":"Yu Wang, Chengdong Kong, Juntao Ao, Hongshen Li, Chengyi Wang, Xiaojiang Wu, Zhongxiao Zhang","doi":"10.1016/j.combustflame.2024.113748","DOIUrl":"10.1016/j.combustflame.2024.113748","url":null,"abstract":"<div><div>This work finds that with the assistance of filamentary plasma discharge, the ammonia combustion becomes very complicated with multiple modes, especially including a low-frequency (<2 Hz) oscillating combustion mode. This flame oscillation, which is sensitive to the discharge characteristics and the flow/fuel/boundary conditions, can be further categorized as three sub-modes denoted as Mode C1, C2 and C3. Their transitions occur depending on the flow rate (<em>Q</em><sub>tot</sub>), the equivalence ratio (<em>ϕ</em>) and the temperature of the chamber wall. Detailed dynamic characteristics of flame kernels and flow fields during the oscillating process are further acquired to confirm the changing recirculation zone and the Karlovitz (Ka) number. A phenomenological mechanism based on recirculated energy feedback is thus proposed to explain the oscillating flame. In addition, the <em>S</em>-curve is used to elucidate the plasma impacts on the complicated combustion characteristics, especially the oscillating combustion modes and their transitions. It is found that the folded <em>S</em>-curve can be stretched with the help of plasma discharge. By considering the energy feedback effects, the <em>S</em>-curve can become more stretched with the Da number for extinction (i.e. <span><math><mrow><mi>D</mi><msub><mi>a</mi><msup><mi>E</mi><mo>′</mo></msup></msub></mrow></math></span>) being larger than that for ignition (i.e. <span><math><mrow><mi>D</mi><msub><mi>a</mi><msup><mi>I</mi><mo>′</mo></msup></msub></mrow></math></span>). Meanwhile, the branch with the Da number between <span><math><mrow><mi>D</mi><msub><mi>a</mi><msup><mi>I</mi><mo>′</mo></msup></msub></mrow></math></span> and <span><math><mrow><mi>D</mi><msub><mi>a</mi><msup><mi>E</mi><mo>′</mo></msup></msub></mrow></math></span> is unstable and the multiple oscillating combustion modes (i.e. Mode C1-C3) can be categorized on the unstable branch of the stretched <em>S</em>-curve.</div></div>","PeriodicalId":280,"journal":{"name":"Combustion and Flame","volume":"270 ","pages":"Article 113748"},"PeriodicalIF":5.8,"publicationDate":"2024-09-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142318790","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}