Pub Date : 2023-01-17DOI: 10.1080/13647830.2023.2165967
Z. Huo, M. Cleary, Kun Wu, A. Masri, Xuejun Fan
Due to the composition-dependent stiffness of chemistry, simulations of reactive turbulent flows may present computational load imbalance among parallel processes when spatial decomposition is used for parallelisation, causing high CPU idle time and waste of computational resources. To increase computational efficiency, a dynamic load balancing (DLB) model is proposed to redistribute computational load among computing cores. The DLB model exploits a decomposition in the mixture fraction space with two dynamic adjusting decomposition strategies to realise load redistribution. The DLB model is suitable for the integration of chemistry on stochastic particles in hybrid Eulerian/Lagrangian turbulent combustion models in which the Eulerian field is conventionally decomposed statically in physical space in a way that balances the computational load for the solution of the Navier-Stokes equation but which does not generally lead to balanced load for the computation of the composition fields. Here it is tested using an OpenFOAM-based platform, mmcFoam, which is a comprehensive object-orientated C++ library for stochastic turbulent combustion modelling. Apart from direct integration (DI) for chemistry, the DLB model is also coupled with dynamic adaptive chemistry (DAC) and in situ adaptive tabulation (ISAT), which allows for extra speedup. The performance of the coupled models is validated and assessed for two laboratory flame conditions that exhibit different levels of computational load imbalance. Overall, the DLB model effectively balances the computational load distribution and increases the effective usage of computing power, shortening the simulation wall time required. Moreover, a strong scaling test is carried out using up to 512 cores. Although all approaches have sub-ideal scalability, the scalability of each with DLB is significantly better than without DLB. While DLB-ISAT has relatively poor scalability compared to the DI- and DAC-based methods, DLB-ISAT still ranks the fastest among the algorithms in all scaling trials.
{"title":"A dynamic load balancing model coupled with DAC and ISAT for a stochastic turbulent combustion model","authors":"Z. Huo, M. Cleary, Kun Wu, A. Masri, Xuejun Fan","doi":"10.1080/13647830.2023.2165967","DOIUrl":"https://doi.org/10.1080/13647830.2023.2165967","url":null,"abstract":"Due to the composition-dependent stiffness of chemistry, simulations of reactive turbulent flows may present computational load imbalance among parallel processes when spatial decomposition is used for parallelisation, causing high CPU idle time and waste of computational resources. To increase computational efficiency, a dynamic load balancing (DLB) model is proposed to redistribute computational load among computing cores. The DLB model exploits a decomposition in the mixture fraction space with two dynamic adjusting decomposition strategies to realise load redistribution. The DLB model is suitable for the integration of chemistry on stochastic particles in hybrid Eulerian/Lagrangian turbulent combustion models in which the Eulerian field is conventionally decomposed statically in physical space in a way that balances the computational load for the solution of the Navier-Stokes equation but which does not generally lead to balanced load for the computation of the composition fields. Here it is tested using an OpenFOAM-based platform, mmcFoam, which is a comprehensive object-orientated C++ library for stochastic turbulent combustion modelling. Apart from direct integration (DI) for chemistry, the DLB model is also coupled with dynamic adaptive chemistry (DAC) and in situ adaptive tabulation (ISAT), which allows for extra speedup. The performance of the coupled models is validated and assessed for two laboratory flame conditions that exhibit different levels of computational load imbalance. Overall, the DLB model effectively balances the computational load distribution and increases the effective usage of computing power, shortening the simulation wall time required. Moreover, a strong scaling test is carried out using up to 512 cores. Although all approaches have sub-ideal scalability, the scalability of each with DLB is significantly better than without DLB. While DLB-ISAT has relatively poor scalability compared to the DI- and DAC-based methods, DLB-ISAT still ranks the fastest among the algorithms in all scaling trials.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2023-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48202425","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 : 2023-01-12DOI: 10.1080/13647830.2023.2165454
A. Tyliszczak, A. Wawrzak, K. Wawrzak
This paper presents an analysis of numerical and modelling issues based on a simulation of nitrogen-diluted hydrogen lifted flame evolving in a hot co-flow. We apply the large-eddy simulations (LES) method with the so-called ‘no combustion model’ and concentrate on the impact of chemical mechanisms and various discretisation schemes on the obtained results. The main attention is put to the latter issue in which we analyse to what extent a discretisation method of the convective terms of the scalar equations for the species and enthalpy affects the solutions. We consider eight commonly known total variation diminishing (TVD) schemes and three upwind schemes of the second order. The remaining terms in the scalar equations and all the terms in the Navier–Stokes equations are discretised applying the sixth-order compact difference method. Such an approach makes the discretisation errors of the convective terms the main factor affecting the solutions from the numerical point of view. Prior to the main analysis, the differences caused by the use of various TVD or upwind schemes are highlighted based on a single scalar transport equation with a known analytical solution. The results obtained for the flame are compared to experimental data taken from the literature. It is shown that the differences due to the application of a particular discretisation method are of similar magnitude as the differences between the simulation results and experimental data. Moreover, analysis of the impact of the chemical mechanism showed that observed differences are comparable to these originating from the use of different discretisation methods.
{"title":"Impact of a discretisation method and chemical kinetics on the accuracy of simulation of a lifted hydrogen flame","authors":"A. Tyliszczak, A. Wawrzak, K. Wawrzak","doi":"10.1080/13647830.2023.2165454","DOIUrl":"https://doi.org/10.1080/13647830.2023.2165454","url":null,"abstract":"This paper presents an analysis of numerical and modelling issues based on a simulation of nitrogen-diluted hydrogen lifted flame evolving in a hot co-flow. We apply the large-eddy simulations (LES) method with the so-called ‘no combustion model’ and concentrate on the impact of chemical mechanisms and various discretisation schemes on the obtained results. The main attention is put to the latter issue in which we analyse to what extent a discretisation method of the convective terms of the scalar equations for the species and enthalpy affects the solutions. We consider eight commonly known total variation diminishing (TVD) schemes and three upwind schemes of the second order. The remaining terms in the scalar equations and all the terms in the Navier–Stokes equations are discretised applying the sixth-order compact difference method. Such an approach makes the discretisation errors of the convective terms the main factor affecting the solutions from the numerical point of view. Prior to the main analysis, the differences caused by the use of various TVD or upwind schemes are highlighted based on a single scalar transport equation with a known analytical solution. The results obtained for the flame are compared to experimental data taken from the literature. It is shown that the differences due to the application of a particular discretisation method are of similar magnitude as the differences between the simulation results and experimental data. Moreover, analysis of the impact of the chemical mechanism showed that observed differences are comparable to these originating from the use of different discretisation methods.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2023-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48695909","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-12-28DOI: 10.1080/13647830.2022.2157753
R. Lamioni, C. Bronzoni, M. Folli, L. Tognotti, C. Galletti
Domestic condensing boilers are equipped with perforated burners ensuring short-length premixed flames issues from a series of circular holes and slits. Despite some efforts that have been devoted to understanding the effect of hole diameter, pattern, and hole-to-hole distance on the resulting flames, very little is known about the flames from a series of slits. In this work 3-dimensional numerical simulations with skeletal kinetic mechanisms are performed to determine the structure of premixed methane-air flames issuing from two-slits patterns, often recurring in practical burner designs, i.e. several equally-spaced slits and a group of four slits. The arrangement of the slits greatly influences the resulting flames. A significant change in the flame behaviour occurs depending on the inlet velocity; at low speeds, the flames issuing from the slits present a flat region, while increasing the speed they assume a conical shape along the slit length. Neighbour flames are distinct from each other at low speeds, while they interact strongly with increasing the velocity. Interestingly, a series of several slits produce a single long wedge-shaped flame, while the group of four slits generates a single conical flame.
{"title":"Effect of slit pattern on the structure of premixed flames issuing from perforated burners in domestic condensing boilers","authors":"R. Lamioni, C. Bronzoni, M. Folli, L. Tognotti, C. Galletti","doi":"10.1080/13647830.2022.2157753","DOIUrl":"https://doi.org/10.1080/13647830.2022.2157753","url":null,"abstract":"Domestic condensing boilers are equipped with perforated burners ensuring short-length premixed flames issues from a series of circular holes and slits. Despite some efforts that have been devoted to understanding the effect of hole diameter, pattern, and hole-to-hole distance on the resulting flames, very little is known about the flames from a series of slits. In this work 3-dimensional numerical simulations with skeletal kinetic mechanisms are performed to determine the structure of premixed methane-air flames issuing from two-slits patterns, often recurring in practical burner designs, i.e. several equally-spaced slits and a group of four slits. The arrangement of the slits greatly influences the resulting flames. A significant change in the flame behaviour occurs depending on the inlet velocity; at low speeds, the flames issuing from the slits present a flat region, while increasing the speed they assume a conical shape along the slit length. Neighbour flames are distinct from each other at low speeds, while they interact strongly with increasing the velocity. Interestingly, a series of several slits produce a single long wedge-shaped flame, while the group of four slits generates a single conical flame.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-12-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44021612","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-12-22DOI: 10.1080/13647830.2022.2157333
H. S. Mukunda, A. Shiva Kumar, Sachin Payannad, C. S. Bhaskar Dixit
The extensive work on wood crib fires, both experimental and model development over the last six decades is examined in some detail. The wide range of parameters of crib tests and the theories have been reviewed in the literature and a satisfactory correlation has still to emerge from these studies. From early times, the burn flux (g/m s) has been considered the most appropriate parameter to characterise the burn behaviour of the cribs. These data on a re-examination revealed a surprisingly simple behaviour of a linear variation of the mass loss rate (g/s) with the mass of the crib, particularly for smaller size sticks with the crib placed on the ground. Some insight into this behaviour is brought out and the basic idea has been pursued to reveal an alternate and a more accurate correlation for the burn rate with crib mass and the crib size as principal parameters with the crib height-to-spacing ratio providing a minor correction. The resulting correlation has been compared with over a hundred and fifty experimental data along with a modified Thomas correlation and shown to perform much better for smaller-size sticks.
{"title":"Insights into and the evolution of a novel predictive model for free burning wooden cribs","authors":"H. S. Mukunda, A. Shiva Kumar, Sachin Payannad, C. S. Bhaskar Dixit","doi":"10.1080/13647830.2022.2157333","DOIUrl":"https://doi.org/10.1080/13647830.2022.2157333","url":null,"abstract":"The extensive work on wood crib fires, both experimental and model development over the last six decades is examined in some detail. The wide range of parameters of crib tests and the theories have been reviewed in the literature and a satisfactory correlation has still to emerge from these studies. From early times, the burn flux (g/m s) has been considered the most appropriate parameter to characterise the burn behaviour of the cribs. These data on a re-examination revealed a surprisingly simple behaviour of a linear variation of the mass loss rate (g/s) with the mass of the crib, particularly for smaller size sticks with the crib placed on the ground. Some insight into this behaviour is brought out and the basic idea has been pursued to reveal an alternate and a more accurate correlation for the burn rate with crib mass and the crib size as principal parameters with the crib height-to-spacing ratio providing a minor correction. The resulting correlation has been compared with over a hundred and fifty experimental data along with a modified Thomas correlation and shown to perform much better for smaller-size sticks.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-12-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48483488","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-12-14DOI: 10.1080/13647830.2022.2156930
M. Yoon
This paper proposes a new approach for one-dimensional thermoacoustic combustor models. Our new model is a transfer function estimated from the frequency response of the linearised Euler equation to a spatially normalised and temporally impulsive input. The proposed approach can deal with combustors with varying cross-sectional areas under a non-zero mean flow, distributed heating/cooling, and outlet boundary conditions involving entropy waves, overcoming limitations of the popular network model. In addition our new approach can provide a more reliable thermoacoustic model for combustors with entropy-related boundary conditions, remedying the inaccurate entropy model of the network model. Numerical comparisons of our new model with a network model show apparent similarities between the two, validating the new model. It is also observed that, compared to our new model, the network model is more sensitive to mean flow and significantly overestimates the entropy wave effects on combustor acoustics.
{"title":"Distributed transfer function approach for one-dimensional thermoacoustic combustor models","authors":"M. Yoon","doi":"10.1080/13647830.2022.2156930","DOIUrl":"https://doi.org/10.1080/13647830.2022.2156930","url":null,"abstract":"This paper proposes a new approach for one-dimensional thermoacoustic combustor models. Our new model is a transfer function estimated from the frequency response of the linearised Euler equation to a spatially normalised and temporally impulsive input. The proposed approach can deal with combustors with varying cross-sectional areas under a non-zero mean flow, distributed heating/cooling, and outlet boundary conditions involving entropy waves, overcoming limitations of the popular network model. In addition our new approach can provide a more reliable thermoacoustic model for combustors with entropy-related boundary conditions, remedying the inaccurate entropy model of the network model. Numerical comparisons of our new model with a network model show apparent similarities between the two, validating the new model. It is also observed that, compared to our new model, the network model is more sensitive to mean flow and significantly overestimates the entropy wave effects on combustor acoustics.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-12-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46390105","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-12-12DOI: 10.1080/13647830.2022.2153739
Kun Wu, Yuting Jiang, Zhijie Huo, D. Cheng, Xuejun Fan
The integration of stiff ordinary differential equation (ODE) systems associated with detailed chemical kinetics is computationally demanding in practical combustion simulations. Despite the various approaches in expediting the computational efficiency, it is still necessary to optimise the cell-wise calculation in operator-splitting type simulations of reactive flow. In this work, we proposed an improved stiff-ODE solver framework targeting to speed up the simulation of reactive flow in OpenFOAM. This framework combines the Radau-IIA and backward differentiation formula (BDF) ODE-integration algorithms, the pyJac-based fully analytical Jacobian formulation, and dense-based LAPACK and sparse-based KLU sophisticated linear system solvers. We evaluate the performance of the efficient solver framework on various benchmark combustion problems across a wide range of chemical kinetic complexities. A comprehensive investigation of the key elements of stiff ODE solvers is conducted in the homogeneous reactor, focusing respectively on the influences of error tolerance, integration time interval, Jacobian evaluation methodology, and linear system solver on the accuracy and efficiency trade-off. More realistic simulation results are presented regarding the one-dimensional laminar flame and three-dimensional turbulent flame. The results indicate that the Radau-IIA is more preferable in both efficiency and accuracy compared with the widely used BDF and Seulex methods for large integration interval, whereas the differences between three methods diminish as the integration time interval decreases. In all cases, it is found that the full analytical Jacobian is more advantageous for small mechanisms of species number around 50–100 while the approximated formulation of Jacobian is recommended for larger ones. Furthermore, the more robust linear system solvers provide significant improvement on computational efficiency with the dense-based LAPACK solver being more suitable for small to moderate-scale mechanisms while sparse-based KLU being superior for large-scale mechanisms. The proposed efficient solver framework in its optimal configuration obtains more than 2.6 times speedup in realistic high-fidelity flame simulation with a 57 species combustion mechanism.
{"title":"An improved stiff-ODE solving framework for reacting flow simulations with detailed chemistry in OpenFOAM","authors":"Kun Wu, Yuting Jiang, Zhijie Huo, D. Cheng, Xuejun Fan","doi":"10.1080/13647830.2022.2153739","DOIUrl":"https://doi.org/10.1080/13647830.2022.2153739","url":null,"abstract":"The integration of stiff ordinary differential equation (ODE) systems associated with detailed chemical kinetics is computationally demanding in practical combustion simulations. Despite the various approaches in expediting the computational efficiency, it is still necessary to optimise the cell-wise calculation in operator-splitting type simulations of reactive flow. In this work, we proposed an improved stiff-ODE solver framework targeting to speed up the simulation of reactive flow in OpenFOAM. This framework combines the Radau-IIA and backward differentiation formula (BDF) ODE-integration algorithms, the pyJac-based fully analytical Jacobian formulation, and dense-based LAPACK and sparse-based KLU sophisticated linear system solvers. We evaluate the performance of the efficient solver framework on various benchmark combustion problems across a wide range of chemical kinetic complexities. A comprehensive investigation of the key elements of stiff ODE solvers is conducted in the homogeneous reactor, focusing respectively on the influences of error tolerance, integration time interval, Jacobian evaluation methodology, and linear system solver on the accuracy and efficiency trade-off. More realistic simulation results are presented regarding the one-dimensional laminar flame and three-dimensional turbulent flame. The results indicate that the Radau-IIA is more preferable in both efficiency and accuracy compared with the widely used BDF and Seulex methods for large integration interval, whereas the differences between three methods diminish as the integration time interval decreases. In all cases, it is found that the full analytical Jacobian is more advantageous for small mechanisms of species number around 50–100 while the approximated formulation of Jacobian is recommended for larger ones. Furthermore, the more robust linear system solvers provide significant improvement on computational efficiency with the dense-based LAPACK solver being more suitable for small to moderate-scale mechanisms while sparse-based KLU being superior for large-scale mechanisms. The proposed efficient solver framework in its optimal configuration obtains more than 2.6 times speedup in realistic high-fidelity flame simulation with a 57 species combustion mechanism.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"47566286","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-12-09DOI: 10.1080/13647830.2022.2153741
Linlin Yang, Yiqing Wang, Zheng Chen
Recently, multiple spark ignition has received great attention since it helps to increase thermal efficiency and to reduce misfire in engines. Multiple spark ignition also affects the combustion duration and thereby it can be used for knock control. However, previous studies reported opposite trends in terms of how multiple spark ignition affects engine knock. This work aims to assess and interpret the influence of flame propagation direction induced by different ignition configurations on combustion duration and end-gas autoignition/engine knock. Two simplified and idealised ignition configurations are studied theoretically and numerically. One is with infinite number of sparks at side circular wall, which induces an inwardly propagating flame (IPF); and the other is with a single central spark, which induces an outwardly propagating flame (OPF). In the asymptotic theoretical analysis, the canonical 1D formulations for IPF and OPF are reduced to 0D model. Based on the 0D model, OPF and IPF at different initial temperatures are studied and compared. Counterintuitively, it is found that the combustion duration of OPF is shorter than that of IPF when there is no end-gas autoignition. On the other hand, the combustion duration of IPF is shorter than that of OPF when end-gas autoignition occurs. Furthermore, end-gas autoignition is found to be more prone to occur in IPF than OPF. These interesting observations are interpreted through assessing the ignition delay time and different components of the absolute flame propagation speed. The theoretical results are validated by transient simulations considering detailed chemistry and transport which are conducted for IPF and OPF in an iso-octane/air mixture at different initial temperatures and pressures. Both theoretical and numerical results suggest that compared to infinite number of ignition sparks at side wall, the single central ignition has the advantages in shortening the combustion duration and reducing the tendency of end-gas autoignition.
{"title":"Comparison of combustion duration and end-gas autoignition in inwardly and outwardly propagating flames induced by different ignition configurations","authors":"Linlin Yang, Yiqing Wang, Zheng Chen","doi":"10.1080/13647830.2022.2153741","DOIUrl":"https://doi.org/10.1080/13647830.2022.2153741","url":null,"abstract":"Recently, multiple spark ignition has received great attention since it helps to increase thermal efficiency and to reduce misfire in engines. Multiple spark ignition also affects the combustion duration and thereby it can be used for knock control. However, previous studies reported opposite trends in terms of how multiple spark ignition affects engine knock. This work aims to assess and interpret the influence of flame propagation direction induced by different ignition configurations on combustion duration and end-gas autoignition/engine knock. Two simplified and idealised ignition configurations are studied theoretically and numerically. One is with infinite number of sparks at side circular wall, which induces an inwardly propagating flame (IPF); and the other is with a single central spark, which induces an outwardly propagating flame (OPF). In the asymptotic theoretical analysis, the canonical 1D formulations for IPF and OPF are reduced to 0D model. Based on the 0D model, OPF and IPF at different initial temperatures are studied and compared. Counterintuitively, it is found that the combustion duration of OPF is shorter than that of IPF when there is no end-gas autoignition. On the other hand, the combustion duration of IPF is shorter than that of OPF when end-gas autoignition occurs. Furthermore, end-gas autoignition is found to be more prone to occur in IPF than OPF. These interesting observations are interpreted through assessing the ignition delay time and different components of the absolute flame propagation speed. The theoretical results are validated by transient simulations considering detailed chemistry and transport which are conducted for IPF and OPF in an iso-octane/air mixture at different initial temperatures and pressures. Both theoretical and numerical results suggest that compared to infinite number of ignition sparks at side wall, the single central ignition has the advantages in shortening the combustion duration and reducing the tendency of end-gas autoignition.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-12-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44447563","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}
Spontaneous coal combustion is one of the most common disasters in coal mine production. In order to explore the mechanism of coal spontaneous combustion more deeply, coal samples from the Yangdong wellfield of Jizhong Energy were selected for oxidative heat energy analysis experiments. A temperature-programmed experiment was selected to study the changes in characteristic parameters during the low-temperature oxidation of coal under different air supply conditions. TG-DSC experiments were conducted to study the characteristic temperature changes and thermodynamic characteristics of coal combustion processes at different heating rates. The study results show that the coal is most easily oxidised in the low-temperature oxidation stage when the air supply is 120 ml/min. The oxygen consumption rate, CO generation rate, and maximum and minimum heat release intensity are all greater at this airflow than under other conditions. The process of spontaneous combustion of coal has six characteristic temperature points and is divided into five stages. The characteristic temperature of the coal sample increased with the increase of the heating rate, and the TG/DTG curve showed a hysteresis phenomenon. DSC temperature curve shifts toward the high temperature with the increase of the heating rate, and the exothermic region is expanded. Isokinetic analysis (F-W-O and V-W) and Coats-Redfern model for calculating thermodynamic parameters. The activation energy of the samples decreased with the increase of the heating rate in the range of 2∼20°C·min−1 and showed a decreasing trend with the increase of the conversion rate.
{"title":"Experimental study and thermodynamic analysis of coal spontaneous combustion characteristics","authors":"Minbo Zhang, Zichao Wang, Longkang Wang, Zhen Zhang, Dangyu Zhang, Chunxin Li","doi":"10.1080/13647830.2022.2153742","DOIUrl":"https://doi.org/10.1080/13647830.2022.2153742","url":null,"abstract":"Spontaneous coal combustion is one of the most common disasters in coal mine production. In order to explore the mechanism of coal spontaneous combustion more deeply, coal samples from the Yangdong wellfield of Jizhong Energy were selected for oxidative heat energy analysis experiments. A temperature-programmed experiment was selected to study the changes in characteristic parameters during the low-temperature oxidation of coal under different air supply conditions. TG-DSC experiments were conducted to study the characteristic temperature changes and thermodynamic characteristics of coal combustion processes at different heating rates. The study results show that the coal is most easily oxidised in the low-temperature oxidation stage when the air supply is 120 ml/min. The oxygen consumption rate, CO generation rate, and maximum and minimum heat release intensity are all greater at this airflow than under other conditions. The process of spontaneous combustion of coal has six characteristic temperature points and is divided into five stages. The characteristic temperature of the coal sample increased with the increase of the heating rate, and the TG/DTG curve showed a hysteresis phenomenon. DSC temperature curve shifts toward the high temperature with the increase of the heating rate, and the exothermic region is expanded. Isokinetic analysis (F-W-O and V-W) and Coats-Redfern model for calculating thermodynamic parameters. The activation energy of the samples decreased with the increase of the heating rate in the range of 2∼20°C·min−1 and showed a decreasing trend with the increase of the conversion rate.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-12-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44389286","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-12-08DOI: 10.1080/13647830.2022.2153740
Arash Mousemi, M. Jadidi, S. Dworkin, W. Bushe
The Uniform Conditional State (UCS) and the Multidimensional Flamelet Manifold (MFM) models are methods for the tabulation of chemistry in simulations of turbulent flames. The high-dimensionality of the tables these models generate and many possible combinations of the values for the input variables necessitate the allocation of a considerable size of memory during CFD calculations. This issue becomes even more problematic when adding more conditioning variables to the model. In this study, two Artificial Intelligence (AI)-based approaches referred to as Decision Tree (DT) and Artificial Neural Network (ANN) are developed and tested to provide in situ chemistry representation. The goal is to predict four parameters (outputs) accurately with low memory demand and computational cost. The trained AI models are then employed for simulation of a turbulent premixed flame. Comparison of the results from the AI-based approaches to those from the conventional UCS model shows acceptable agreement. The memory and CPU requirements from the different approaches are compared. It is found that the ANN model reduces the size of the chemistry table by around 92%. Conversely, the DT-based model reduces the size of the chemistry model by only 40%. The CPU time for using the DT model during the CFD calculations was around 10% shorter than the conventional approach while it was 8% higher for the ANN model. It was concluded that, based on the particular applications, different AI-based methods can facilitate an efficient representation of the chemistry manifold.
{"title":"Application of machine learning in low-order manifold representation of chemistry in turbulent flames","authors":"Arash Mousemi, M. Jadidi, S. Dworkin, W. Bushe","doi":"10.1080/13647830.2022.2153740","DOIUrl":"https://doi.org/10.1080/13647830.2022.2153740","url":null,"abstract":"The Uniform Conditional State (UCS) and the Multidimensional Flamelet Manifold (MFM) models are methods for the tabulation of chemistry in simulations of turbulent flames. The high-dimensionality of the tables these models generate and many possible combinations of the values for the input variables necessitate the allocation of a considerable size of memory during CFD calculations. This issue becomes even more problematic when adding more conditioning variables to the model. In this study, two Artificial Intelligence (AI)-based approaches referred to as Decision Tree (DT) and Artificial Neural Network (ANN) are developed and tested to provide in situ chemistry representation. The goal is to predict four parameters (outputs) accurately with low memory demand and computational cost. The trained AI models are then employed for simulation of a turbulent premixed flame. Comparison of the results from the AI-based approaches to those from the conventional UCS model shows acceptable agreement. The memory and CPU requirements from the different approaches are compared. It is found that the ANN model reduces the size of the chemistry table by around 92%. Conversely, the DT-based model reduces the size of the chemistry model by only 40%. The CPU time for using the DT model during the CFD calculations was around 10% shorter than the conventional approach while it was 8% higher for the ANN model. It was concluded that, based on the particular applications, different AI-based methods can facilitate an efficient representation of the chemistry manifold.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-12-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41271424","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-12-08DOI: 10.1080/13647830.2022.2153743
Aleksander D. Zakharov, R. Fursenko, S. Minaev
Optimisation approach to automate selection of global reaction mechanisms rate constants is proposed and studied. The objective of optimisation is to find reaction rate constants minimising deviation of some flame characteristics (e.g. laminar burning velocity, ignition delay time, etc.) calculated by global mechanism from their reference values known from experiments or computed by detailed mechanisms. Examples of one, two and four step mechanisms optimisation with respect to laminar burning velocity and concentration distributions in counterflow diffusion flame are given. Computer codes implementing optimisation algorithm for these examples are also afforded and can be modified and used for reaction constants selection in various applications. Uniqueness of singlecriteria and multicriteria optimisation solutions is studied numerically by computations with different initial guesses and by direct evaluation of the objective functions. Particularly, it is found that for considered global mechanisms the minimum value of objective function is reached in some subdomain of the parametric space. This means that any values of rate parameters from this subdomain results in almost the same deviations of chosen flame characteristic from its reference value.
{"title":"Optimisation method for automatic selection of rate constants of global reaction mechanisms","authors":"Aleksander D. Zakharov, R. Fursenko, S. Minaev","doi":"10.1080/13647830.2022.2153743","DOIUrl":"https://doi.org/10.1080/13647830.2022.2153743","url":null,"abstract":"Optimisation approach to automate selection of global reaction mechanisms rate constants is proposed and studied. The objective of optimisation is to find reaction rate constants minimising deviation of some flame characteristics (e.g. laminar burning velocity, ignition delay time, etc.) calculated by global mechanism from their reference values known from experiments or computed by detailed mechanisms. Examples of one, two and four step mechanisms optimisation with respect to laminar burning velocity and concentration distributions in counterflow diffusion flame are given. Computer codes implementing optimisation algorithm for these examples are also afforded and can be modified and used for reaction constants selection in various applications. Uniqueness of singlecriteria and multicriteria optimisation solutions is studied numerically by computations with different initial guesses and by direct evaluation of the objective functions. Particularly, it is found that for considered global mechanisms the minimum value of objective function is reached in some subdomain of the parametric space. This means that any values of rate parameters from this subdomain results in almost the same deviations of chosen flame characteristic from its reference value.","PeriodicalId":50665,"journal":{"name":"Combustion Theory and Modelling","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-12-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41767074","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}