In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal e...In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.展开更多
The gas-droplet two-phase reacting flow in a model combustor with the V-gutter flame holder is studied by an Eulerian-Lagrangian large-eddy simulation (LES) approach. The k-equation subgrid-scale model is used to simu...The gas-droplet two-phase reacting flow in a model combustor with the V-gutter flame holder is studied by an Eulerian-Lagrangian large-eddy simulation (LES) approach. The k-equation subgrid-scale model is used to simulate the subgrid eddy viscosity, and the eddy-break-up (EBU) combustion subgrid-scale model is used to determine the chemical reaction rate. A two-step turbulent combustion subgrid-scale model is employed for calculating carbon monoxide CO concentration, and the NO subgrid-scale pollutant formation model for the evaluation of the rate of NO formation. The heat flux model is applied to the prediction of radiant heat transfer. The gas phase is solved with the SIMPLE algorithm and a hybrid scheme in the staggered grid system. The liquid phase equations are solved in a Lagrangian frame in reference of the particle-source-in-cell (PSIC) algorithm. From simulation results, the exchange of mass, moment and energy between gas and particle fields for the reacting flow in the afterburner with a V-gutter flame holder can be obtained. By the comparison of experimental and simulation results, profile temperature and pollutant of the outlet are quite in agreement with experimental data. Results show that the LES approach for predicting the two-phase instantaneous reacting flow and pollutant emissions in the afterburner is feasible.展开更多
Critical temperatures and pressures of nominal reacting mixture in synthesis of dimethyl carbonate （DMC） from methanol and carbon dioxide （quaternary mixture of carbon dioxide ＋ methanol ＋ water ＋ DMC） were mea...Critical temperatures and pressures of nominal reacting mixture in synthesis of dimethyl carbonate （DMC） from methanol and carbon dioxide （quaternary mixture of carbon dioxide ＋ methanol ＋ water ＋ DMC） were measured using a high-pressure view cell. The results suggested that the critical properties of the reacting mixture depended on the reaction extent as well as its initial composition （initial ratio of carbon dioxide to methanol）. Such information is essential for determining the reaction conditions when one intends to carry out the synthesis of DMC with CO2 and methanol under supercritical conditions.展开更多
In this paper, author considers a 3 x 3 system for a reacting flow model propesed by [9]. Since this model has source term, it can be considered as a relaxation approximation to 2 x 2 systems of conservation laws, whi...In this paper, author considers a 3 x 3 system for a reacting flow model propesed by [9]. Since this model has source term, it can be considered as a relaxation approximation to 2 x 2 systems of conservation laws, which include the well-known p-system. From this viewpoint, the author establishes the global existence and the nonlinear stability of travelling wave solutions by L-2 energy method.展开更多
This article focuses on the development of a discontinuous Galerkin (DG) method for simulations of multicomponent and chemically reacting flows. Compared to aerodynamic flow applications, in which DG methods have been...This article focuses on the development of a discontinuous Galerkin (DG) method for simulations of multicomponent and chemically reacting flows. Compared to aerodynamic flow applications, in which DG methods have been successfully employed, DG simulations of chemically reacting flows introduce challenges that arise from flow unsteadiness, combustion, heat release, compressibility effects, shocks, and variations in thermodynamic properties. To address these challenges, algorithms are developed, including an entropy-bounded DG method, an entropy-residual shock indicator, and a new formulation of artificial viscosity. The performance and capabilities of the resulting DG method are demonstrated in several relevant applications, including shock/bubble interaction, turbulent combustion, and detonation. It is concluded that the developed DG method shows promising performance in application to multicomponent reacting flows. The paper concludes with a discussion of further research needs to enable the application of DG methods to more complex reacting flows.展开更多
The stochastic simulation algorithm （SSA） accurately depicts spatially homogeneous wellstirred chemically reacting systems with small populations of chemical species and properly represents noise, but it is often ab...The stochastic simulation algorithm （SSA） accurately depicts spatially homogeneous wellstirred chemically reacting systems with small populations of chemical species and properly represents noise, but it is often abandoned when modeling larger systems because of its computational complexity. In this work, a twin support vector regression based stochastic simulations algorithm （TS^3A） is proposed by combining the twin support vector regression and SSA, the former is a well-known robust regression method in machine learning. Numerical results indicate that this proposed algorithm can be applied to a wide range of chemically reacting systems and obtain significant improvements on efficiency and accuracy with fewer simulating runs over the existing methods.展开更多
Presented here is an L-leap method for accelerating stochastic simulation of well-stirred chemically reacting systems, in which the number of reactions occurring in a reaction channel with the largest propensity funct...Presented here is an L-leap method for accelerating stochastic simulation of well-stirred chemically reacting systems, in which the number of reactions occurring in a reaction channel with the largest propensity function is calculated from the leap condition and the number of reactions occurring in the other reaction channels are generated by using binomial random variables during a leap. The L-leap method can better satisfy the leap condition. Numerical simulation results indicate that the L-leap method can obtain better performance than established methods.展开更多
In this paper, the supersonic chemically reacting mixing layer is simulated with the third order ENN scheme, based on the Navier-Stokes equations, containing transport equations of all species. The numerical results s...In this paper, the supersonic chemically reacting mixing layer is simulated with the third order ENN scheme, based on the Navier-Stokes equations, containing transport equations of all species. The numerical results show that the thickness of mixing layer increases gradually along the flow direction, and that the Kelvin-Helmholtz, instabilities may not exist in mixing layer flows.展开更多
The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the n...The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment （SOM） combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment （ASOM） combustion model is well validated by DNS.展开更多
We are concerned with the global existence of entropy solutions of the two-dimensional steady Euler equations for an ideal gas, which undergoes a one-step exothermic chemical reaction under the Arrhenius-type kinetics...We are concerned with the global existence of entropy solutions of the two-dimensional steady Euler equations for an ideal gas, which undergoes a one-step exothermic chemical reaction under the Arrhenius-type kinetics. The reaction rate function φ（T ） is assumed to have a positive lower bound. We first consider the Cauchy problem （the initial value problem）, that is, seek a supersonic downstream reacting flow when the incoming flow is supersonic, and establish the global existence of entropy solutions when the total variation of the initial data is sufficiently small. Then we analyze the problem of steady supersonic, exothermically reacting Euler flow past a Lipschitz wedge, generating an ad-ditional detonation wave attached to the wedge vertex, which can be then formulated as an initial-boundary value problem. We establish the global existence of entropy solutions containing the additional detonation wave （weak or strong, determined by the wedge angle at the wedge vertex） when the total variation of both the slope of the wedge boundary and the incoming flow is suitably small. The downstream asymptotic behavior of the global solutions is also obtained.展开更多
The initial boundary value problem (IBVP) for the 3×3 hyperbolic system of reacting flow with source term proposed by R.J.LeVeque and others (see [8]) is considered.It is shown, in the present paper, that if the ...The initial boundary value problem (IBVP) for the 3×3 hyperbolic system of reacting flow with source term proposed by R.J.LeVeque and others (see [8]) is considered.It is shown, in the present paper, that if the initial data are a suitable perturbation of a shiftcd shock profile which is suitably away from the boundary, then there exists a unique smooth solution in R2+ to the IBVP of the 3×3 hyperbolic system, which tends to another shifted shock profile of this system as t →∞.展开更多
In this paper,we study the zero dissipation limit with a vacuum for the reacting mixture Navier-Stokes equations.For proper smooth initial data that the initial density tends to zero as the relevant physical coefficie...In this paper,we study the zero dissipation limit with a vacuum for the reacting mixture Navier-Stokes equations.For proper smooth initial data that the initial density tends to zero as the relevant physical coefficients tend to zero,we demonstrate that the solution tends to a rarefaction wave connected to a vacuum on the left side coupled with a zero mass fraction of reactant.What is more,the uniform convergence rate is obtained.展开更多
In this paper,we study the time-asymptotically nonlinear stability of rarefaction waves for the Cauchy problem of the compressible Navier-Stokes equations for a reacting mixture with zero heat conductivity in one dime...In this paper,we study the time-asymptotically nonlinear stability of rarefaction waves for the Cauchy problem of the compressible Navier-Stokes equations for a reacting mixture with zero heat conductivity in one dimension.If the corresponding Riemann problem for the compressible Euler system admits the solutions consisting of rarefaction waves only,it is shown that its Cauchy problem has a unique global solution which tends time-asymptotically towards the rarefaction waves,while the initial perturbation and the strength of rarefaction waves are suitably small.展开更多
We consider the large time behavior of solutions of the Cauchy problem for the one-dimensional compressible Navier-Stokes equations for a reacting mixture.When the corresponding Riemann problem for the Euler system ad...We consider the large time behavior of solutions of the Cauchy problem for the one-dimensional compressible Navier-Stokes equations for a reacting mixture.When the corresponding Riemann problem for the Euler system admits a contact discontinuity wave,it is shown that the viscous contact wave which corresponds to the contact discontinuity is asymptotically stable,provided the strength of contact discontinuity and the initial perturbation are suitably small.We apply the approach introduced in Huang,Li and Matsumura(2010)[1]and the elemen tary L2-energy met hods.展开更多
This paper presents an experimental investigation of the turbulent reacting flow in a swirl combustor with staged air injection. The air injected into the combustor is composed of the primary swirling jet and the seco...This paper presents an experimental investigation of the turbulent reacting flow in a swirl combustor with staged air injection. The air injected into the combustor is composed of the primary swirling jet and the secon-dary non-swirling jet. A three dimension-laser particle dynamic analyzer (PDA) was employed to measure the in-stantaneous gas velocity. The probability density functions (PDF) for the instantaneous gas axial and tangential ve-locities at each measuring location, as well as the radial profiles of the root mean square of fluctuating gas axial and tangential velocities and the second-order moment for the fluctuating gas axial and tangential velocities are ob-tained. The measured results delineate the turbulence properties of the swirling reacting flow under the conditions of staged combustion.展开更多
This paper describes the numerical calculations of gaseous reaction flows in a model of gas turbine combustors. The profiles of hydrodynamic and thermodynamic patterns in a three-dimensional combustor model are obtain...This paper describes the numerical calculations of gaseous reaction flows in a model of gas turbine combustors. The profiles of hydrodynamic and thermodynamic patterns in a three-dimensional combustor model are obtained by solving the governing differential transport equations. The well-established numerical prediction algorithm SIMPLE, the modified k-ε turbulence model and k-ε-g turbulent diffusion flame model have been adopted in computations. The β function has been selected as probability density function. The effect of combustion process on flow patterns has been investigated. The calculated results have been verified by experiments. They are in remarkably good agreement.展开更多
A mathematical model for the analysis of a gas-solid reacting system is presented. This model is an alternative to the classical shrinking-core model. The model has a structure that can be easily transformed into a ca...A mathematical model for the analysis of a gas-solid reacting system is presented. This model is an alternative to the classical shrinking-core model. The model has a structure that can be easily transformed into a canonical control form, which is proper for controller synthesis. Analytical solution of the model to describe the open-loop behavior is expressed in terms of the Lambert function. The Lambert function is evaluated from aTaylorexpansion series. Besides, a controller is proposed to regulate the reacted layer thickness using initially the diffusion coefficient as control input. The control law is synthesized employing the feedback linearization technique. Main contributions of this work are the synthesis of the layer thickness controller, and the employment of the process temperature as substitute of the diffusion coefficient as the control input.展开更多
This paper proposes a robust and efficient oscillation-eliminating discontinuous Galerkin(OEDG)method for solving multicomponent chemically reacting flows,which is an extension and application of the recent work[M.Pen...This paper proposes a robust and efficient oscillation-eliminating discontinuous Galerkin(OEDG)method for solving multicomponent chemically reacting flows,which is an extension and application of the recent work[M.Peng,Z.Sun,and K.Wu,Math.Comput.,2024,doi.org/10.1090/mcom/3998].Following recently developed high-order bound-preserving discontinuous Galerkin method in[J.Du and Y.Yang,J.Comput.Phys.,469(2022),111548],we incorporate an OE procedure after each Runge-–Kutta time stage to suppress spurious oscillations.The OE procedure is defined by the solution operator of a damping equation,which can be analytically solved without requiring discretization,making its implementation straightforward,non-intrusive,and efficient.Through careful design of the damping coefficients,the proposed OEDG method not only achieves the essentially non-oscillatory(ENO)property without compromising accuracy but also preserves the conservative property—an indispensable aspect of the bound-preserving technique introduced in[J.Du and Y.Yang,J.Comput.Phys.,469(2022),111548].The effectiveness and robustness of the OEDG method are demonstrated through a series of one-and two-dimensional numerical tests on the compressible Euler and Navier–Stokes equations for chemically reacting flows.These results highlight the method's capability to handle complex flow dynamics while maintaining stability and high-order accuracy.展开更多
基金supported by the NSF under Grant DMS-1818467Simons Foundation under Grant 961585.
文摘In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.
文摘The gas-droplet two-phase reacting flow in a model combustor with the V-gutter flame holder is studied by an Eulerian-Lagrangian large-eddy simulation (LES) approach. The k-equation subgrid-scale model is used to simulate the subgrid eddy viscosity, and the eddy-break-up (EBU) combustion subgrid-scale model is used to determine the chemical reaction rate. A two-step turbulent combustion subgrid-scale model is employed for calculating carbon monoxide CO concentration, and the NO subgrid-scale pollutant formation model for the evaluation of the rate of NO formation. The heat flux model is applied to the prediction of radiant heat transfer. The gas phase is solved with the SIMPLE algorithm and a hybrid scheme in the staggered grid system. The liquid phase equations are solved in a Lagrangian frame in reference of the particle-source-in-cell (PSIC) algorithm. From simulation results, the exchange of mass, moment and energy between gas and particle fields for the reacting flow in the afterburner with a V-gutter flame holder can be obtained. By the comparison of experimental and simulation results, profile temperature and pollutant of the outlet are quite in agreement with experimental data. Results show that the LES approach for predicting the two-phase instantaneous reacting flow and pollutant emissions in the afterburner is feasible.
基金financial support of the State Key Fundamental Research Project(2006CB202504)the National Natural Science Foundation of China(No.20473110)Natural Science Foundation of Shanxi Province.
文摘Critical temperatures and pressures of nominal reacting mixture in synthesis of dimethyl carbonate （DMC） from methanol and carbon dioxide （quaternary mixture of carbon dioxide ＋ methanol ＋ water ＋ DMC） were measured using a high-pressure view cell. The results suggested that the critical properties of the reacting mixture depended on the reaction extent as well as its initial composition （initial ratio of carbon dioxide to methanol）. Such information is essential for determining the reaction conditions when one intends to carry out the synthesis of DMC with CO2 and methanol under supercritical conditions.
文摘In this paper, author considers a 3 x 3 system for a reacting flow model propesed by [9]. Since this model has source term, it can be considered as a relaxation approximation to 2 x 2 systems of conservation laws, which include the well-known p-system. From this viewpoint, the author establishes the global existence and the nonlinear stability of travelling wave solutions by L-2 energy method.
基金supported by an Early Career Faculty grant from NASA's Space Technology Research Grants Programprovided by the NASA High-End Computing (HEC) Program through the NASA Advanced Supercomputing (NAS) Division at Ames Research Center
文摘This article focuses on the development of a discontinuous Galerkin (DG) method for simulations of multicomponent and chemically reacting flows. Compared to aerodynamic flow applications, in which DG methods have been successfully employed, DG simulations of chemically reacting flows introduce challenges that arise from flow unsteadiness, combustion, heat release, compressibility effects, shocks, and variations in thermodynamic properties. To address these challenges, algorithms are developed, including an entropy-bounded DG method, an entropy-residual shock indicator, and a new formulation of artificial viscosity. The performance and capabilities of the resulting DG method are demonstrated in several relevant applications, including shock/bubble interaction, turbulent combustion, and detonation. It is concluded that the developed DG method shows promising performance in application to multicomponent reacting flows. The paper concludes with a discussion of further research needs to enable the application of DG methods to more complex reacting flows.
基金This work was supported by the National Natural Science Foundation of China （No.30871341）, the National High-Tech Research and Development Program of China （No.2006AA02-Z190）, the Shanghai Leading Academic Discipline Project （No.S30405）, and the Natural Science Foundation of Shanghai Normal University （No.SK200937）.
文摘The stochastic simulation algorithm （SSA） accurately depicts spatially homogeneous wellstirred chemically reacting systems with small populations of chemical species and properly represents noise, but it is often abandoned when modeling larger systems because of its computational complexity. In this work, a twin support vector regression based stochastic simulations algorithm （TS^3A） is proposed by combining the twin support vector regression and SSA, the former is a well-known robust regression method in machine learning. Numerical results indicate that this proposed algorithm can be applied to a wide range of chemically reacting systems and obtain significant improvements on efficiency and accuracy with fewer simulating runs over the existing methods.
基金Project supported by the National Natural Science Foundation of China (No.30571059)the National High-Tech Research and Development Program of China(No.2006AA02Z190).
文摘Presented here is an L-leap method for accelerating stochastic simulation of well-stirred chemically reacting systems, in which the number of reactions occurring in a reaction channel with the largest propensity function is calculated from the leap condition and the number of reactions occurring in the other reaction channels are generated by using binomial random variables during a leap. The L-leap method can better satisfy the leap condition. Numerical simulation results indicate that the L-leap method can obtain better performance than established methods.
基金The project supported by the National Natural Science Foundation of China
文摘In this paper, the supersonic chemically reacting mixing layer is simulated with the third order ENN scheme, based on the Navier-Stokes equations, containing transport equations of all species. The numerical results show that the thickness of mixing layer increases gradually along the flow direction, and that the Kelvin-Helmholtz, instabilities may not exist in mixing layer flows.
基金Supported by the National Natural Science Foundation of China （50606026, 50736006）.
文摘The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment （SOM） combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment （ASOM） combustion model is well validated by DNS.
基金Gui-Qiang CHEN was supported in part by the UK EPSRC Science and Innovation Award to the Oxford Centre for Nonlinear PDE(EP/E035027/1)the NSFC under a joint project Grant 10728101+4 种基金the Royal Society-Wolfson Research Merit Award(UK)Changguo XIAO was supported in part by the NSFC under a joint project Grant 10728101Yongqian ZHANG was supported in part by NSFC Project 11031001NSFC Project 11121101the 111 Project B08018(China)
文摘We are concerned with the global existence of entropy solutions of the two-dimensional steady Euler equations for an ideal gas, which undergoes a one-step exothermic chemical reaction under the Arrhenius-type kinetics. The reaction rate function φ（T ） is assumed to have a positive lower bound. We first consider the Cauchy problem （the initial value problem）, that is, seek a supersonic downstream reacting flow when the incoming flow is supersonic, and establish the global existence of entropy solutions when the total variation of the initial data is sufficiently small. Then we analyze the problem of steady supersonic, exothermically reacting Euler flow past a Lipschitz wedge, generating an ad-ditional detonation wave attached to the wedge vertex, which can be then formulated as an initial-boundary value problem. We establish the global existence of entropy solutions containing the additional detonation wave （weak or strong, determined by the wedge angle at the wedge vertex） when the total variation of both the slope of the wedge boundary and the incoming flow is suitably small. The downstream asymptotic behavior of the global solutions is also obtained.
文摘The initial boundary value problem (IBVP) for the 3×3 hyperbolic system of reacting flow with source term proposed by R.J.LeVeque and others (see [8]) is considered.It is shown, in the present paper, that if the initial data are a suitable perturbation of a shiftcd shock profile which is suitably away from the boundary, then there exists a unique smooth solution in R2+ to the IBVP of the 3×3 hyperbolic system, which tends to another shifted shock profile of this system as t →∞.
基金supported by the National Natural Science Foundation of China (11971193 and 12171001)。
文摘In this paper,we study the zero dissipation limit with a vacuum for the reacting mixture Navier-Stokes equations.For proper smooth initial data that the initial density tends to zero as the relevant physical coefficients tend to zero,we demonstrate that the solution tends to a rarefaction wave connected to a vacuum on the left side coupled with a zero mass fraction of reactant.What is more,the uniform convergence rate is obtained.
基金supported by the Beijing Natural Science Foundation(1182004,Z180007,1192001).
文摘In this paper,we study the time-asymptotically nonlinear stability of rarefaction waves for the Cauchy problem of the compressible Navier-Stokes equations for a reacting mixture with zero heat conductivity in one dimension.If the corresponding Riemann problem for the compressible Euler system admits the solutions consisting of rarefaction waves only,it is shown that its Cauchy problem has a unique global solution which tends time-asymptotically towards the rarefaction waves,while the initial perturbation and the strength of rarefaction waves are suitably small.
基金supported by the National Natural Science Foundation of China(11871341).
文摘We consider the large time behavior of solutions of the Cauchy problem for the one-dimensional compressible Navier-Stokes equations for a reacting mixture.When the corresponding Riemann problem for the Euler system admits a contact discontinuity wave,it is shown that the viscous contact wave which corresponds to the contact discontinuity is asymptotically stable,provided the strength of contact discontinuity and the initial perturbation are suitably small.We apply the approach introduced in Huang,Li and Matsumura(2010)[1]and the elemen tary L2-energy met hods.
基金Supported jointly by the National Natural Science Foundation of China (No.59806006) and the Laboratory Open Fund ofTsinghua University.
文摘This paper presents an experimental investigation of the turbulent reacting flow in a swirl combustor with staged air injection. The air injected into the combustor is composed of the primary swirling jet and the secon-dary non-swirling jet. A three dimension-laser particle dynamic analyzer (PDA) was employed to measure the in-stantaneous gas velocity. The probability density functions (PDF) for the instantaneous gas axial and tangential ve-locities at each measuring location, as well as the radial profiles of the root mean square of fluctuating gas axial and tangential velocities and the second-order moment for the fluctuating gas axial and tangential velocities are ob-tained. The measured results delineate the turbulence properties of the swirling reacting flow under the conditions of staged combustion.
文摘This paper describes the numerical calculations of gaseous reaction flows in a model of gas turbine combustors. The profiles of hydrodynamic and thermodynamic patterns in a three-dimensional combustor model are obtained by solving the governing differential transport equations. The well-established numerical prediction algorithm SIMPLE, the modified k-ε turbulence model and k-ε-g turbulent diffusion flame model have been adopted in computations. The β function has been selected as probability density function. The effect of combustion process on flow patterns has been investigated. The calculated results have been verified by experiments. They are in remarkably good agreement.
文摘A mathematical model for the analysis of a gas-solid reacting system is presented. This model is an alternative to the classical shrinking-core model. The model has a structure that can be easily transformed into a canonical control form, which is proper for controller synthesis. Analytical solution of the model to describe the open-loop behavior is expressed in terms of the Lambert function. The Lambert function is evaluated from aTaylorexpansion series. Besides, a controller is proposed to regulate the reacted layer thickness using initially the diffusion coefficient as control input. The control law is synthesized employing the feedback linearization technique. Main contributions of this work are the synthesis of the layer thickness controller, and the employment of the process temperature as substitute of the diffusion coefficient as the control input.
基金supported by National Key R&D Program of China(Grant No.2021YFA0719200)supported by the Shenzhen Science and Technology Program(Grant No.RCJC20221008092757098)supported by the Postdoctoral Science Foundation of China(Grant No.2024M751284)。
文摘This paper proposes a robust and efficient oscillation-eliminating discontinuous Galerkin(OEDG)method for solving multicomponent chemically reacting flows,which is an extension and application of the recent work[M.Peng,Z.Sun,and K.Wu,Math.Comput.,2024,doi.org/10.1090/mcom/3998].Following recently developed high-order bound-preserving discontinuous Galerkin method in[J.Du and Y.Yang,J.Comput.Phys.,469(2022),111548],we incorporate an OE procedure after each Runge-–Kutta time stage to suppress spurious oscillations.The OE procedure is defined by the solution operator of a damping equation,which can be analytically solved without requiring discretization,making its implementation straightforward,non-intrusive,and efficient.Through careful design of the damping coefficients,the proposed OEDG method not only achieves the essentially non-oscillatory(ENO)property without compromising accuracy but also preserves the conservative property—an indispensable aspect of the bound-preserving technique introduced in[J.Du and Y.Yang,J.Comput.Phys.,469(2022),111548].The effectiveness and robustness of the OEDG method are demonstrated through a series of one-and two-dimensional numerical tests on the compressible Euler and Navier–Stokes equations for chemically reacting flows.These results highlight the method's capability to handle complex flow dynamics while maintaining stability and high-order accuracy.
