Pub Date : 2022-05-20DOI: 10.1080/13647830.2022.2072237
Seyed Mehdi Ashrafizadeh, C. Devaud
The modelling of soot formation is investigated for two turbulent flames, at atmospheric and 3 atm pressure conditions. For the first time, a semi-empirical soot formulation that accounts for soot inception, coagulation, surface growth, and oxidation processes is coupled with the turbulent combustion model, Conditional Source-term Estimation (CSE) using Reynolds-Averaged Navier–Stokes equations. Detailed chemistry is included and an optically thin radiation model is considered. Non-adiabatic chemistry tabulations are created. Good agreement with the experiments is found for turbulent mixing and temperature fields in both flames, with some discrepancies believed to be due to the turbulence modelling approach. At 1 atm, the soot volume fractions are in reasonable agreement with the experiments, but typically smaller than the measurements with the centerline peak locating closer to the fuel exit. At 3 atm, good agreement between the numerical predictions and experimental data is achieved for the soot volume fraction within the experimental error. The centerline peak location is observed slightly farther downstream. Possible sources of discrepancies are examined and comparison with previously published numerical results is also undertaken. Differential diffusion and modified soot chemistry constants may bring further improvement. Without any particular tuning of soot chemistry, soot modelling within CSE is shown to be a promising approach.
{"title":"Investigation of conditional source-term estimation coupled with a semi-empirical model for soot predictions in two turbulent flames","authors":"Seyed Mehdi Ashrafizadeh, C. Devaud","doi":"10.1080/13647830.2022.2072237","DOIUrl":"https://doi.org/10.1080/13647830.2022.2072237","url":null,"abstract":"The modelling of soot formation is investigated for two turbulent flames, at atmospheric and 3 atm pressure conditions. For the first time, a semi-empirical soot formulation that accounts for soot inception, coagulation, surface growth, and oxidation processes is coupled with the turbulent combustion model, Conditional Source-term Estimation (CSE) using Reynolds-Averaged Navier–Stokes equations. Detailed chemistry is included and an optically thin radiation model is considered. Non-adiabatic chemistry tabulations are created. Good agreement with the experiments is found for turbulent mixing and temperature fields in both flames, with some discrepancies believed to be due to the turbulence modelling approach. At 1 atm, the soot volume fractions are in reasonable agreement with the experiments, but typically smaller than the measurements with the centerline peak locating closer to the fuel exit. At 3 atm, good agreement between the numerical predictions and experimental data is achieved for the soot volume fraction within the experimental error. The centerline peak location is observed slightly farther downstream. Possible sources of discrepancies are examined and comparison with previously published numerical results is also undertaken. Differential diffusion and modified soot chemistry constants may bring further improvement. Without any particular tuning of soot chemistry, soot modelling within CSE is shown to be a promising approach.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-05-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43759005","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2022-05-19DOI: 10.1080/13647830.2022.2070549
F. Chitgarha, F. Ommi, M. Farshchi
The reaction progress variable is a crucial concept in the advanced flamelet combustion models. As a controlling variable, a well-defined progress variable must consider the essential features of the combustion process. It is usually a heuristically defined linear combination of some major chemical species mass fractions. However, such a simple definition could lead to inaccurate results for the fuel-rich reactive mixtures or complicated fuels, due to the vast number of chemical species in the combustion process. In this paper, a new method for generating a reaction progress variable is proposed through solving a constrained optimisation problem. The proposed method uses a genetic algorithm with new constraints. The major new constraint is the minimisation of the inverse of a progress variable-based Damköhler number in addition to the minimisation of the gradients of a collection of chemical species concentrations, as used in the previous methods. Hence, this scheme increases the Damköhler number defined based on the progress variable. The applicability and performance of the current optimised progress variable are evaluated for ethanol–air partially premixed flames in an axisymmetric two-dimensional counterflow burner and a two-dimensional plugged flow triple-flame burner. The effects of the number of chemical species included in the progress variable and the flow field strain rate on a partially premixed ethanol–air flame prediction are investigated. Results indicate that including the progress variable Damköhler number in the determination of the progress variable has a considerable effect on the accuracy of Flamelet Generated Manifold (FGM) model prediction of fuel-rich and lean reactive mixtures, especially at higher strain rates. Also, it is shown that the inclusion of the critical chemical species for ignition and fuel decomposition processes, such as CH3O2, CH3CHO, sC2H4OH, HO2, H and H2O2, in the definition of progress variable has a significant effect on the accuracy of the ethanol–air flame predictions.
{"title":"Assessment of optimal reaction progress variable characteristics for partially premixed flames","authors":"F. Chitgarha, F. Ommi, M. Farshchi","doi":"10.1080/13647830.2022.2070549","DOIUrl":"https://doi.org/10.1080/13647830.2022.2070549","url":null,"abstract":"The reaction progress variable is a crucial concept in the advanced flamelet combustion models. As a controlling variable, a well-defined progress variable must consider the essential features of the combustion process. It is usually a heuristically defined linear combination of some major chemical species mass fractions. However, such a simple definition could lead to inaccurate results for the fuel-rich reactive mixtures or complicated fuels, due to the vast number of chemical species in the combustion process. In this paper, a new method for generating a reaction progress variable is proposed through solving a constrained optimisation problem. The proposed method uses a genetic algorithm with new constraints. The major new constraint is the minimisation of the inverse of a progress variable-based Damköhler number in addition to the minimisation of the gradients of a collection of chemical species concentrations, as used in the previous methods. Hence, this scheme increases the Damköhler number defined based on the progress variable. The applicability and performance of the current optimised progress variable are evaluated for ethanol–air partially premixed flames in an axisymmetric two-dimensional counterflow burner and a two-dimensional plugged flow triple-flame burner. The effects of the number of chemical species included in the progress variable and the flow field strain rate on a partially premixed ethanol–air flame prediction are investigated. Results indicate that including the progress variable Damköhler number in the determination of the progress variable has a considerable effect on the accuracy of Flamelet Generated Manifold (FGM) model prediction of fuel-rich and lean reactive mixtures, especially at higher strain rates. Also, it is shown that the inclusion of the critical chemical species for ignition and fuel decomposition processes, such as CH3O2, CH3CHO, sC2H4OH, HO2, H and H2O2, in the definition of progress variable has a significant effect on the accuracy of the ethanol–air flame predictions.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-05-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45543332","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2022-05-10DOI: 10.1080/13647830.2022.2069601
T. Shen, Huahua Xiao
This paper studies premixed flame dynamics in half-open tubes by solving the two-dimensional, fully compressible, reactive Navier-Stokes equations on a dynamically adapting mesh using a high-order algorithm. A simplified chemical-diffusive model was used to describe the chemical reaction and diffusive transports in a stoichiometric hydrogen-air mixture. The influence of the length scale was examined by considering four tube heights at a fixed aspect ratio α = 7. The numerical simulations show that the flame evolves into a tulip flame (TF) in all the tubes shortly after being ignited at the open end. Variation in tube size leads to differences in the evolution of TF and generation of expansion waves. In a sufficiently large tube (d > 0.5 cm), the TF further develops into a series of more unstable distorted tulip flames (DTFs). By contrast, in a small tube (d < 0.5 cm), the TF shape remains until the end of the combustion. In addition, both the flame and pressure oscillate significantly almost in the entire process of flame propagation in the large tubes, while the oscillating behaviour in flame or pressure is negligible in the small tube after TF forms. It was found that the TF formation mechanism is length-scale dependent even for the same type of geometry and condition. A detailed examination of the interactions between flame, boundary layer, and pressure waves showed two mechanisms of TF formation: (1) boundary layer effect for the larger tubes (d ≥ 0.5 cm), and (2) Rayleigh–Taylor instability driven by compression waves for the smallest tube (d = 0.25 cm). The DTF formation in the half-open tubes is closely related to the expansion waves generated by the collapse of the TF cusp. The expansion waves cause a reverse flow in the boundary layer ahead of the flame front and consequently initiate the DTF.
{"title":"Numerical study of the stability of premixed flames propagating in half-open tubes","authors":"T. Shen, Huahua Xiao","doi":"10.1080/13647830.2022.2069601","DOIUrl":"https://doi.org/10.1080/13647830.2022.2069601","url":null,"abstract":"This paper studies premixed flame dynamics in half-open tubes by solving the two-dimensional, fully compressible, reactive Navier-Stokes equations on a dynamically adapting mesh using a high-order algorithm. A simplified chemical-diffusive model was used to describe the chemical reaction and diffusive transports in a stoichiometric hydrogen-air mixture. The influence of the length scale was examined by considering four tube heights at a fixed aspect ratio α = 7. The numerical simulations show that the flame evolves into a tulip flame (TF) in all the tubes shortly after being ignited at the open end. Variation in tube size leads to differences in the evolution of TF and generation of expansion waves. In a sufficiently large tube (d > 0.5 cm), the TF further develops into a series of more unstable distorted tulip flames (DTFs). By contrast, in a small tube (d < 0.5 cm), the TF shape remains until the end of the combustion. In addition, both the flame and pressure oscillate significantly almost in the entire process of flame propagation in the large tubes, while the oscillating behaviour in flame or pressure is negligible in the small tube after TF forms. It was found that the TF formation mechanism is length-scale dependent even for the same type of geometry and condition. A detailed examination of the interactions between flame, boundary layer, and pressure waves showed two mechanisms of TF formation: (1) boundary layer effect for the larger tubes (d ≥ 0.5 cm), and (2) Rayleigh–Taylor instability driven by compression waves for the smallest tube (d = 0.25 cm). The DTF formation in the half-open tubes is closely related to the expansion waves generated by the collapse of the TF cusp. The expansion waves cause a reverse flow in the boundary layer ahead of the flame front and consequently initiate the DTF.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-05-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43777883","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2022-05-10DOI: 10.1080/13647830.2022.2071170
W. Jayasuriya, T. C. Mulky, Kyle E. Niemeyer
Smouldering combustion plays a key role in wildfires in forests, grasslands, and peatlands due to its common occurrence in porous fuels like peat and duff. As a consequence, understanding smouldering behaviour in these fuels is crucial. Such fuels are generally composed of cellulose, hemicellulose, and lignin. Here we present an updated computational model for simulating smouldering combustion in cellulose and hemicellulose mixtures. We used this model to examine changes in smouldering propagation speed and peak temperatures with varying fuel composition and density. For a given fuel composition, increases in density decrease the propagation speed and increase mean peak temperature; for a given density, increases in hemicellulose content increase both propagation speed and peak temperature. We also examined the role of natural fuel expansion with the addition of water. Without expansion, addition of moisture content reduces the propagation speed primarily due to increasing (wet) fuel density. However, with fuel expansion similar to that observed in peat, the propagation speed increases due to the overall drop in fuel density. Finally, we studied the influence of fuel composition on critical moisture content of ignition and extinction: mixtures dominated by hemicellulose have 10% higher critical moisture content due to the increase in peak temperature.
{"title":"Smouldering combustion in cellulose and hemicellulose mixtures: Examining the roles of density, fuel composition, oxygen concentration, and moisture content","authors":"W. Jayasuriya, T. C. Mulky, Kyle E. Niemeyer","doi":"10.1080/13647830.2022.2071170","DOIUrl":"https://doi.org/10.1080/13647830.2022.2071170","url":null,"abstract":"Smouldering combustion plays a key role in wildfires in forests, grasslands, and peatlands due to its common occurrence in porous fuels like peat and duff. As a consequence, understanding smouldering behaviour in these fuels is crucial. Such fuels are generally composed of cellulose, hemicellulose, and lignin. Here we present an updated computational model for simulating smouldering combustion in cellulose and hemicellulose mixtures. We used this model to examine changes in smouldering propagation speed and peak temperatures with varying fuel composition and density. For a given fuel composition, increases in density decrease the propagation speed and increase mean peak temperature; for a given density, increases in hemicellulose content increase both propagation speed and peak temperature. We also examined the role of natural fuel expansion with the addition of water. Without expansion, addition of moisture content reduces the propagation speed primarily due to increasing (wet) fuel density. However, with fuel expansion similar to that observed in peat, the propagation speed increases due to the overall drop in fuel density. Finally, we studied the influence of fuel composition on critical moisture content of ignition and extinction: mixtures dominated by hemicellulose have 10% higher critical moisture content due to the increase in peak temperature.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-05-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42177999","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2022-04-26DOI: 10.1080/13647830.2022.2066024
A. Bayliss, E. Shafirovich, V. Volpert
Combustion of a porous solid fuel is considered. An exothermic reaction takes place between the fuel and a gaseous oxidiser which is delivered to the reaction zone by filtration through the pores in the sample from an open end toward which the combustion wave propagates (counterflow filtration). The gas reacts with the solid fuel to form a solid product. The gas filtration is due to the pressure difference between the ambient pressure at the open end and the pressure in the reaction zone where the gas is being consumed (referred to as natural filtration). A 1D mathematical model based on equations describing conservation of energy, gas mass, solid reactant mass, and gas momentum, as well as an equation of state, and appropriate boundary and initial conditions is formulated and analytically studied taking advantage of the separation of length scales in the process. When the reaction zone is sufficiently far from the open end, the combustion wave propagates at a constant speed and has a time-independent structure, while when the reaction is close to the open end (closer than the filtration length), the structure of the combustion wave and its speed become time dependent. Both cases are discussed in the paper though the main emphasis is on short samples, in which the combustion wave is affected by the gas flow from the open end during the entire propagation process. A specific example of interest involves magnesium as the solid fuel and oxygen as the gaseous oxidiser.
{"title":"Counterflow combustion waves in short samples of metal powders at natural filtration of oxygen","authors":"A. Bayliss, E. Shafirovich, V. Volpert","doi":"10.1080/13647830.2022.2066024","DOIUrl":"https://doi.org/10.1080/13647830.2022.2066024","url":null,"abstract":"Combustion of a porous solid fuel is considered. An exothermic reaction takes place between the fuel and a gaseous oxidiser which is delivered to the reaction zone by filtration through the pores in the sample from an open end toward which the combustion wave propagates (counterflow filtration). The gas reacts with the solid fuel to form a solid product. The gas filtration is due to the pressure difference between the ambient pressure at the open end and the pressure in the reaction zone where the gas is being consumed (referred to as natural filtration). A 1D mathematical model based on equations describing conservation of energy, gas mass, solid reactant mass, and gas momentum, as well as an equation of state, and appropriate boundary and initial conditions is formulated and analytically studied taking advantage of the separation of length scales in the process. When the reaction zone is sufficiently far from the open end, the combustion wave propagates at a constant speed and has a time-independent structure, while when the reaction is close to the open end (closer than the filtration length), the structure of the combustion wave and its speed become time dependent. Both cases are discussed in the paper though the main emphasis is on short samples, in which the combustion wave is affected by the gas flow from the open end during the entire propagation process. A specific example of interest involves magnesium as the solid fuel and oxygen as the gaseous oxidiser.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-04-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45611452","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2022-04-16DOI: 10.1080/13647830.2022.2029947
P. Breda, Chunkan Yu, U. Maas, M. Pfitzner
Detailed chemistry simulations of turbulent reacting flows involving combustion of hydrocarbons can easily exceed the available computational resources, depending on the dimensions of the chemical system. Previous work of the authors showed how the combination of the Eulerian Stochastic Fields (ESF) model with tabulated chemistry based on 2-dimensional Reaction-Diffusion Manifolds (REDIM) provided a significant computational speed-up, compared to the finite rate ESF solver. In this work, the behaviour for flame F, featuring a strong degree of extinction, is further investigated. A comparison is performed for 2D and 3D databases, both using simplified and detailed transport, where the scalar dissipation rate is included as the third table parameter. The results show that the upstream sections are well captured by the REDIM built for detailed transport, while the downstream sections are better captured by the simplified transport database. While a 3D-REDIM based on simplified transport seems to be necessary to capture the extinction events, a 2D-REDIM with differential diffusion already provides satisfactory results. Overall, the use of a 3D-REDIM with differential diffusion better describes the global behaviour of flame F.
{"title":"Reaction-Diffusion Manifolds including differential diffusion applied to methane/air combustion in strong extinction regimes","authors":"P. Breda, Chunkan Yu, U. Maas, M. Pfitzner","doi":"10.1080/13647830.2022.2029947","DOIUrl":"https://doi.org/10.1080/13647830.2022.2029947","url":null,"abstract":"Detailed chemistry simulations of turbulent reacting flows involving combustion of hydrocarbons can easily exceed the available computational resources, depending on the dimensions of the chemical system. Previous work of the authors showed how the combination of the Eulerian Stochastic Fields (ESF) model with tabulated chemistry based on 2-dimensional Reaction-Diffusion Manifolds (REDIM) provided a significant computational speed-up, compared to the finite rate ESF solver. In this work, the behaviour for flame F, featuring a strong degree of extinction, is further investigated. A comparison is performed for 2D and 3D databases, both using simplified and detailed transport, where the scalar dissipation rate is included as the third table parameter. The results show that the upstream sections are well captured by the REDIM built for detailed transport, while the downstream sections are better captured by the simplified transport database. While a 3D-REDIM based on simplified transport seems to be necessary to capture the extinction events, a 2D-REDIM with differential diffusion already provides satisfactory results. Overall, the use of a 3D-REDIM with differential diffusion better describes the global behaviour of flame F.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-04-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42748964","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2022-03-19DOI: 10.1080/13647830.2022.2049881
A. Varma, U. Ahmed, N. Chakraborty
The effects of body force on the statistical behaviour of turbulent scalar flux and its closure in the context of Reynolds Averaged Navier–Stokes simulations have been studied using Direct Numerical Simulations (DNS) of statistically planar turbulent premixed flames under different turbulence intensities and Froude numbers. An increase in body force magnitude in the case of unstable density stratification has been found to augment flame wrinkling, burning rate and gradient transport in comparison to a case without body force but with statistically similar unburned gas turbulence. By contrast, an increase in body force magnitude in the case of stable stratification reduces the flame wrinkling, burning rate and gradient transport in comparison to the flame without body force subjected to statistically similar unburned gas turbulence. Based on a-priori DNS analysis, an algebraic closure for turbulent scalar flux has been identified where the Froude number effects are explicitly accounted for. The body force has been found to have significant influence on the statistical behaviours and magnitudes of the terms of the scalar flux transport equation and this effect is particularly strong for the mean pressure gradient term in the scalar flux transport equation. Based on a detailed a priori DNS analysis, suitable model expressions have been identified for the turbulent transport, pressure gradient, dissipation and reaction rate-velocity correlation terms of the scalar flux transport equation by incorporating the effects of body force (e.g. Froude number effects) for improved model performance.
{"title":"Effects of buoyancy on turbulent scalar flux closure for turbulent premixed flames in the context of Reynolds Averaged Navier–Stokes simulations","authors":"A. Varma, U. Ahmed, N. Chakraborty","doi":"10.1080/13647830.2022.2049881","DOIUrl":"https://doi.org/10.1080/13647830.2022.2049881","url":null,"abstract":"The effects of body force on the statistical behaviour of turbulent scalar flux and its closure in the context of Reynolds Averaged Navier–Stokes simulations have been studied using Direct Numerical Simulations (DNS) of statistically planar turbulent premixed flames under different turbulence intensities and Froude numbers. An increase in body force magnitude in the case of unstable density stratification has been found to augment flame wrinkling, burning rate and gradient transport in comparison to a case without body force but with statistically similar unburned gas turbulence. By contrast, an increase in body force magnitude in the case of stable stratification reduces the flame wrinkling, burning rate and gradient transport in comparison to the flame without body force subjected to statistically similar unburned gas turbulence. Based on a-priori DNS analysis, an algebraic closure for turbulent scalar flux has been identified where the Froude number effects are explicitly accounted for. The body force has been found to have significant influence on the statistical behaviours and magnitudes of the terms of the scalar flux transport equation and this effect is particularly strong for the mean pressure gradient term in the scalar flux transport equation. Based on a detailed a priori DNS analysis, suitable model expressions have been identified for the turbulent transport, pressure gradient, dissipation and reaction rate-velocity correlation terms of the scalar flux transport equation by incorporating the effects of body force (e.g. Froude number effects) for improved model performance.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45148606","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2022-03-11DOI: 10.1080/13647830.2022.2103452
W. Sirignano
A new rotational flamelet model with inward swirling flow through a stretched vortex tube is developed for sub-grid modelling to be coupled with the resolved flow for turbulent combustion. The model has critical new features compared to existing models. (i) Non-premixed flames, premixed flames, or multi-branched flame structures are determined rather than prescribed. (ii) The effects of vorticity and the related centrifugal acceleration are determined. (iii) The strain rates and vorticity applied at the sub-grid level can be directly determined from the resolved-scale strain rates and vorticity without a contrived progress variable. (iv) The flamelet model is three-dimensional. (v) The effect of variable density is addressed. (vi) The inward swirl is created by vorticity combined with two compressive normal strain components; this feature distinguishes the model from counterflow flamelet models. Solutions to the multicomponent Navier–Stokes equations governing the flamelet model are obtained. By coordinate transformation, a similar solution is found for the model, through a system of ordinary differential equations. Vorticity creates a centrifugal force on the sub-grid counterflow that modifies the molecular transport rates, burning rates, and flammability limits. Sample computations of the inward swirling rotational flamelet model without coupling to the resolved flow are presented to demonstrate the importance of the new features. Premixed, nonpremixed, and multi-branched flame structures are examined. Parameter surveys are made with rate of normal strain, vorticity, Damköhler number, and Prandtl number. The centrifugal effect has interesting consequences when combined with the variable-density field. Flow direction can reverse; burning rates can be modified; flammability limits can be extended.
{"title":"Inward swirling flamelet model","authors":"W. Sirignano","doi":"10.1080/13647830.2022.2103452","DOIUrl":"https://doi.org/10.1080/13647830.2022.2103452","url":null,"abstract":"A new rotational flamelet model with inward swirling flow through a stretched vortex tube is developed for sub-grid modelling to be coupled with the resolved flow for turbulent combustion. The model has critical new features compared to existing models. (i) Non-premixed flames, premixed flames, or multi-branched flame structures are determined rather than prescribed. (ii) The effects of vorticity and the related centrifugal acceleration are determined. (iii) The strain rates and vorticity applied at the sub-grid level can be directly determined from the resolved-scale strain rates and vorticity without a contrived progress variable. (iv) The flamelet model is three-dimensional. (v) The effect of variable density is addressed. (vi) The inward swirl is created by vorticity combined with two compressive normal strain components; this feature distinguishes the model from counterflow flamelet models. Solutions to the multicomponent Navier–Stokes equations governing the flamelet model are obtained. By coordinate transformation, a similar solution is found for the model, through a system of ordinary differential equations. Vorticity creates a centrifugal force on the sub-grid counterflow that modifies the molecular transport rates, burning rates, and flammability limits. Sample computations of the inward swirling rotational flamelet model without coupling to the resolved flow are presented to demonstrate the importance of the new features. Premixed, nonpremixed, and multi-branched flame structures are examined. Parameter surveys are made with rate of normal strain, vorticity, Damköhler number, and Prandtl number. The centrifugal effect has interesting consequences when combined with the variable-density field. Flow direction can reverse; burning rates can be modified; flammability limits can be extended.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42260821","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2022-03-10DOI: 10.1080/13647830.2022.2049370
V. Kurdyumov, C. Jiménez
We present an investigation of the stabilisation of premixed laminar flames by means of an isolated highly conductive bluff-body, a circular cylinder, placed in a uniform flow of a combustible mixture. It is shown that the problem has non-unique steady-state solutions for certain values of the parameters. Moreover, we solve the time-dependent equations to check the stability of the solutions and demonstrate the possibility of controlling the convergence to a certain steady-state solution.
{"title":"Flame stabilisation by a highly conductive body: multiple steady-state solutions and time-dependent dynamics","authors":"V. Kurdyumov, C. Jiménez","doi":"10.1080/13647830.2022.2049370","DOIUrl":"https://doi.org/10.1080/13647830.2022.2049370","url":null,"abstract":"We present an investigation of the stabilisation of premixed laminar flames by means of an isolated highly conductive bluff-body, a circular cylinder, placed in a uniform flow of a combustible mixture. It is shown that the problem has non-unique steady-state solutions for certain values of the parameters. Moreover, we solve the time-dependent equations to check the stability of the solutions and demonstrate the possibility of controlling the convergence to a certain steady-state solution.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-03-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43593837","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2022-02-18DOI: 10.1080/13647830.2022.2036373
Li Ma, F. Nmira, J. Consalvi
The main objective of this article is to investigate the capability of the flamelet progress variable (FPV) model to capture the extinction processes observed in under-ventilated fire scenarios. To this end, large eddy simulation (LES) of the methane line fire plumes in oxygen-reduced environments down to global extinction, investigated experimentally at the University of Maryland (UMD), is performed. Two experimental burner configurations, that differ by the presence (anchored) or not (non-anchored) of an oxygen anchor to stabilise the flame base, are considered leading to two different extinction modes. Both the FPV and the steady laminar flamelet (SLF) model coupled with a presumed filtered density function (FDF) are considered. The Rank Correlated Full Spectrum k-distribution (RCFSK) model is used as a gas radiative property model. In both non-anchored and anchored scenarios, the FPV model reproduces with fidelity the evolution of the fire plume structure, radiative loss, and combustion efficiency with decreasing down to global extinction, without introducing any adjustable constant. The extinction in the non-anchored scenario occurs owing to flame-based detachment coupled to the generation of a buoyancy-driven vortex and is found to be very sensitive to the grid resolution in the near burner region. The present results suggest that these processes can be adequately resolved with a spatial resolution of 2.5 mm in this region. The SLF model, for its part, provides reliable predictions comparable to the FPV as long as no local extinction/re-ignition process occurs.
{"title":"Modelling extinction/re-ignition processes in fire plumes under oxygen-diluted conditions using flamelet tabulation approaches","authors":"Li Ma, F. Nmira, J. Consalvi","doi":"10.1080/13647830.2022.2036373","DOIUrl":"https://doi.org/10.1080/13647830.2022.2036373","url":null,"abstract":"The main objective of this article is to investigate the capability of the flamelet progress variable (FPV) model to capture the extinction processes observed in under-ventilated fire scenarios. To this end, large eddy simulation (LES) of the methane line fire plumes in oxygen-reduced environments down to global extinction, investigated experimentally at the University of Maryland (UMD), is performed. Two experimental burner configurations, that differ by the presence (anchored) or not (non-anchored) of an oxygen anchor to stabilise the flame base, are considered leading to two different extinction modes. Both the FPV and the steady laminar flamelet (SLF) model coupled with a presumed filtered density function (FDF) are considered. The Rank Correlated Full Spectrum k-distribution (RCFSK) model is used as a gas radiative property model. In both non-anchored and anchored scenarios, the FPV model reproduces with fidelity the evolution of the fire plume structure, radiative loss, and combustion efficiency with decreasing down to global extinction, without introducing any adjustable constant. The extinction in the non-anchored scenario occurs owing to flame-based detachment coupled to the generation of a buoyancy-driven vortex and is found to be very sensitive to the grid resolution in the near burner region. The present results suggest that these processes can be adequately resolved with a spatial resolution of 2.5 mm in this region. The SLF model, for its part, provides reliable predictions comparable to the FPV as long as no local extinction/re-ignition process occurs.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-02-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46377333","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}