{"title":"Simulation of reaction initiation in powder compacting from the surface with composite formation in equivalent reaction cell","authors":"A. Knyazeva, N. Bukrina","doi":"10.1080/13647830.2023.2241421","DOIUrl":"https://doi.org/10.1080/13647830.2023.2241421","url":null,"abstract":"In the present work, we propose a new variant of the model of the composite synthesis under surface heating. We believe that the formation of the composition occurs at the level of reaction cell. The diffusion-controlled process of reduction of one material by another from an oxide is described in the framework of the problem with moving boundaries. It is assumed that the formation of the matrix composition is carried out by the diffusion mechanism. From the position of interfaces, we find the relative volume fractions of oxide phases and the fraction of volume occupied by the matrix. The averaging method of the analysis results at the reaction cell level makes it possible to use these data at the macroscopic level. Volume fractions of phases and average matrix composition (obtained by averaging over the area occupied by the matrix) provide values reflecting the composition of the composite at the macro level. The problem is solved numerically in dimensionless formulation. Dimensionless complexes of physical quantities are distinguished. The estimation of these parameters is performed. The conditions of correctness of the proposed approach are established from comparison of temporal and spatial scales of thermal and diffusion phenomena. A numerical algorithm for the joint solution of macro- and mesolevel problems has been developed. The proposed algorithm makes possible the investigation of the dynamics of composition changes at all points. The model is supplemented by the calculation of stresses and strains from the data on composition and temperature changes in reaction cells. Averaged values of stresses are transferred to macro level.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":"1 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2023-08-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42547632","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}
{"title":"Updated asymptotic structure of cool diffusion flames","authors":"F. Williams, V. Nayagam","doi":"10.1080/13647830.2023.2232338","DOIUrl":"https://doi.org/10.1080/13647830.2023.2232338","url":null,"abstract":"The influence of adding a seventh important elementary step to a six-step mechanism, previously employed for describing the asymptotic structure of normal-alkane droplet combustion supported by cool-flame chemistry in the negative-temperature-coefficient (NTC) range, is investigated by analytical methods. A development paralleling the classical activation-energy-asymptotic (AEA) analysis of the partial-burning regime, accompanied for the first time by an AEA analysis for a negative activation energy, to account properly for the removal of an important intermediate species, is pursued to make predictions of the combustion process, resulting in a revised asymptotic structure that agrees better with computational predictions based on detailed chemistry.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":" ","pages":""},"PeriodicalIF":1.3,"publicationDate":"2023-07-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44746145","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}
{"title":"Coflow filtration combustion waves","authors":"A. Bayliss, E. Shafirovich, V. Volpert","doi":"10.1080/13647830.2023.2226658","DOIUrl":"https://doi.org/10.1080/13647830.2023.2226658","url":null,"abstract":"Recently, it has been proposed to develop space power systems based on filtration combustion of metal powders with oxygen supplied by a chemical oxygen generator. The experiments with lithium and magnesium powders at natural infiltration of oxygen have shown propagation of both counterflow and coflow combustion waves. However, natural filtration combustion of metal powders at relatively low pressures is not sufficiently understood. In the present paper, we investigate the natural coflow combustion waves propagating through a porous medium. The porous matrix is made of metal particles that react with oxygen flowing from the open end of the sample to the reaction zone where it is consumed forming a condensed product which also has a porous structure. The gas flow is due to the pressure difference between the pressure at the open end and that in the reaction zone (the so-called natural filtration). The open end is where the sample is ignited, so that the gas flows through the reaction products, i.e. in the same direction as the combustion wave propagates (coflow filtration). Our mathematical model involves the conservation of energy equation and gas mass, solid reactant mass, and gas momentum balances, as well as an equation of state, and appropriate boundary and initial conditions. It is studied analytically under the combustion front approximation. When the reaction zone is close to the open end, there is sufficient amount of oxygen in the reaction zone and the reaction is controlled by kinetic factors (the kinetic regime of propagation). As the reaction moves away from the open end, it is gas supply that becomes a limiting factor (filtration regime). Both kinetic and filtration regimes of propagation as well the transition between them are analytically studied in this paper.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":" ","pages":""},"PeriodicalIF":1.3,"publicationDate":"2023-06-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42917635","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}
{"title":"Flamelet/transported PDF simulations of ethylene/air jet turbulent non-premixed flame using a three-equation PAH-based soot production model","authors":"F. Nmira, Antoine Bouffard, J. Consalvi","doi":"10.1080/13647830.2023.2224755","DOIUrl":"https://doi.org/10.1080/13647830.2023.2224755","url":null,"abstract":"This article reports flamelet/transported PDF (TPDF) simulations of the well-documented ethylene/air turbulent non-premixed jet flame investigated experimentally at Sandia. The transported PDF equation is solved with the Stochastic Eulerian Field method. The soot production is modelled by a validated three equation PAH-based soot model that predicts the mean soot aggregate properties at low computational time and includes a detailed description of the soot production processes. Gas and soot radiation is modelled using the rank-correlated full-spectrum k model. The turbulence/chemistry/soot production/radiation interactions are taken into account by means of the PDF method. Simulations are run by considering or not soot differential mixing. Based on recent conclusions drawn from Direct Numerical Simulations (Zhou et al., Proc. Combsut. Inst. 38 (2021) 2731–2739), soot differential mixing is modelled by neglecting soot mixing owing to sufficiently large mixing timescales. When soot differential mixing is considered, model predictions reproduce reasonably well the exhaustive set of experimental data, including flame structure, soot statistics and radiative outputs without adjusting parameters. In particular, the predictions demonstrate for the first time the capability of RANS/TPDF models to capture the soot intermittency. On the other hand, neglecting the soot differential mixing produces notable reductions in mean and fluctuating soot volume fraction and soot intermittency. Scatter plot analysis shows that the effects of soot differential mixing are more pronounced in regions of the mixture fraction space where soot surface growth and soot oxidation dominate the soot production, affecting these processes in a non-negligible manner. In an opposite way, soot nucleation and PAH condensation are much less significantly affected. Model results show also that disregarding soot differential mixing reduces the mean soot emission as well the soot emission turbulence/radiation interaction.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":"27 1","pages":"820 - 851"},"PeriodicalIF":1.3,"publicationDate":"2023-06-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43594496","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}
{"title":"High-order methods for the simulation of unsteady counterflow flames subject to stochastic forcing of large amplitude","authors":"F. Bisetti","doi":"10.1080/13647830.2023.2218621","DOIUrl":"https://doi.org/10.1080/13647830.2023.2218621","url":null,"abstract":"Unsteady counterflows are employed to understand and model the effect of turbulence on flames. We present a novel numerical approach for the simulation of one-dimensional unsteady counterflow flames with fourth order spatial discretization and up to fourth order time discretization. The approach couples a three-stage Lobatto IIIa formula for boundary value problems and variable-order, variable time step size Backward Differentiation Formulas for time integration. The framework is explained in detail, its computational performance is analysed, and its use is demonstrated for the case of stochastic forcing of premixed counterflow flames, whereby the imposed rate of strain is a multi-scale lognormal discrete random process with exponential autocorrelation. High-order spatial and temporal discretization make the approach well-suited for the accurate and computationally efficient simulation of the effect of turbulence on flames, which are characterised by large amplitude stochastic fluctuations of the local rate of strain.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":"27 1","pages":"787 - 819"},"PeriodicalIF":1.3,"publicationDate":"2023-06-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45612766","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}
{"title":"Structure and similarity properties of the laminar counterflow spray flame","authors":"Qun Hu, Lipo Wang","doi":"10.1080/13647830.2023.2212637","DOIUrl":"https://doi.org/10.1080/13647830.2023.2212637","url":null,"abstract":"The spray flame in a laminar counterflow is influenced by various setup conditions, for instance the strain rate and liquid droplet parameters, including the initial temperature, size distribution, hydrodynamics and the vaporisation history. With some reasonable simplification, the two phase governing equations can be reformulated as equations in the single gas phase with reconstructed spray-related source terms. Numerical solutions assume interesting similarity features, such as the (almost) independence of the evaporation path after mapping onto a newly defined quantity R, constructed from sensible enthalpy and mixture fraction. In this regard, the dimensionality of free parameters describing the structure of the spray flame can be hopefully reduced, which then provides a different scenario to understand the spray combustion physics. The flame regime diagram is also elaborated in the present framework.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":"27 1","pages":"768 - 786"},"PeriodicalIF":1.3,"publicationDate":"2023-05-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42628960","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}
{"title":"Conditional space evaluation of progress variable definitions for Cambridge/Sandia swirl flames","authors":"N. Sekularac, X. Fang, W. Bushe, M. Davy","doi":"10.1080/13647830.2023.2211537","DOIUrl":"https://doi.org/10.1080/13647830.2023.2211537","url":null,"abstract":"Data from all spatial locations of nine turbulent flames in the Cambridge/Sandia swirl database are combined to study how the choice of scalar variables in conditional moment closure (CMC) type approaches affect the conditional spatial fluctuations of reactive scalars. In order to investigate the influence of swirl and stratification, two additional data-sets have been constructed. Principal component analysis (PCA) is applied to help identify the number of scalar variables and the most appropriate choices to describe the composition space. Two PCA scaling methods have been adopted, namely Pareto and Auto-scaling. Regardless of the data-set investigated and the scaling method used, the results suggest that a single principal component correlated with temperature accounted for the largest variance. For the first moment hypothesis, four progress variable, c, definitions identified by PCA are selected as conditioning variables to investigate the conditional fluctuations and normalised RMS of various species and temperature from all three databases at all axial locations. The results indicate that two control variables based on mixture fraction, Z, and progress variable significantly reduce the conditional fluctuations of scalars compared to a single variable. The selection of progress variables had minimal effects on the RMS of conditional fluctuations for all tested conditions, although a slight reduction of conditional fluctuations was found for the temperature-based progress variable, which can potentially help the further extension of CMC-based models in different flame configurations. The present study also shows that using Z and c (regardless of its definition) as two conditioning scalars enables the detachment of the thermo-chemical state from space, swirl and stratification effects. This suggests that adopting a doubly conditioned source term estimation (DCSE) approach might successfully predict the considered set of flames, assuming that ensembles are divided along the axial direction.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":"27 1","pages":"736 - 767"},"PeriodicalIF":1.3,"publicationDate":"2023-05-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48205184","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}
{"title":"Quenching modes of local flame–wall interaction for turbulent premixed methane combustion in a constant volume vessel","authors":"Ye Wang, Y. Minamoto, M. Shimura, M. Tanahashi","doi":"10.1080/13647830.2023.2209047","DOIUrl":"https://doi.org/10.1080/13647830.2023.2209047","url":null,"abstract":"The quenching mode of local flame–wall interaction (FWI) is investigated for its response to different levels of turbulence intensity as well as its effect on quenching distance, wall heat flux, and near-wall reaction. For that, direct numerical simulations of turbulent premixed methane combustion in a constant volume vessel are carried with initial Karlovitz numbers (Ka) of 1.0, 10.0, and 30.0. Local flame–wall quenching positions are identified based on the local fuel consumption speed during the turbulent combustion process, and the local FWI events have been classified into four quenching modes according to the flame–wall geometric relationships of quenching positions, namely head-on quenching (HOQ), oblique-wall quenching, side-wall quenching (SWQ), and back-on quenching (BOQ). The results show that in the case with higher initial Ka, the flame surface shows a more complicated wrinkled structure due to the flame–turbulence interaction. Meanwhile, the local quenching distance defined based on the identified quenching position is strongly influenced by the near-wall flow, and the range of the local quenching mode extends further to BOQ. However, for all three cases, HOQ and near-HOQ modes account for the majority of local FWI. Wall heat flux and heat release rate (HRR) of near-wall reaction yield high values for the FWI under HOQ or BOQ and are low for SWQ. In addition, there is a discrepancy in the near-wall transportation of some species under different quenching modes, which further leads to the difference in FWI-induced near-wall reaction regarding its total and elementary HRR.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":"27 1","pages":"715 - 735"},"PeriodicalIF":1.3,"publicationDate":"2023-05-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44011529","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}
{"title":"Hydrodynamic instability of premixed tubular flame stabilized in stretched rotating flow","authors":"S. Minaev, E. Dats, K. Shtym","doi":"10.1080/13647830.2023.2204065","DOIUrl":"https://doi.org/10.1080/13647830.2023.2204065","url":null,"abstract":"The paper analyzes the hydrodynamic instability of a cylindrical flame stabilized in a converging rotating gas flow simulating the flow in a cyclone combustion chamber. The analysis was carried out within the framework of a two-dimensional model describing the gas flow with radial and tangential supply of combustible gas through the walls of a cylindrical chamber. In the stationary case, the gas flow in combustion products and the combustible mixture is described by the solution for a rotating tornado-shaped vortex. The rotation of the gas creates a centrifugal force, which has a stabilizing effect on the flame perturbations. The dependences of the perturbations growth rate on the wave vector of perturbations, the radial and tangential inlet gas velocities, and other parameters of the problem are obtained. The analytical criterion of neutral stability separating regions of stable flame with respect of any perturbations and flame instability in the plane of radial and tangential velocities is obtained. In the absence of rotation and at the large flame radii, the model converts into the Landau- Darrieus model of flame hydrodynamic instability.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":"27 1","pages":"685 - 701"},"PeriodicalIF":1.3,"publicationDate":"2023-05-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43469313","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}
{"title":"How random fluctuations can generate and suppress complex oscillatory regimes in continuous stirred tank reactors","authors":"I. Bashkirtseva, L. Ryashko","doi":"10.1080/13647830.2023.2206379","DOIUrl":"https://doi.org/10.1080/13647830.2023.2206379","url":null,"abstract":"Motivated by important chemical engineering applications, we study probabilistic mechanisms of stochastic effects in a randomly forced model of the continuous stirred tank reactor. The bifurcation analysis of the deterministic model reveals the parameter zone of bistability with coexistence of the equilibrium and oscillatory regimes. It is shown that the boundary between basins of these regimes is moving as the bifurcation parameter changes. We study noise-induced transitions across this boundary leading to the stochastic generation of spiking oscillations and backward effect of the noise-induced suppression of large-amplitude spiking. In the study of these stochastic effects, we use statistics extracted from the direct numerical simulation, and an analytical approach using confidence ellipses method. An important probabilistic phenomenon of coherence resonance is revealed and discussed.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":"27 1","pages":"702 - 713"},"PeriodicalIF":1.3,"publicationDate":"2023-04-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48590811","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}