Surface heterogeneous catalysis in a high-enthalpy dissociated environment leads to a remarkable enhancement of aerodynamic heating into the thermal protection surface of hypersonic aircraft.To more accurately predict...Surface heterogeneous catalysis in a high-enthalpy dissociated environment leads to a remarkable enhancement of aerodynamic heating into the thermal protection surface of hypersonic aircraft.To more accurately predict this catalytic heating,a kinetic catalytic model was constructed.This model involved four elementary reactions,the rates of which were determined on mean-field approximation and surface steady-state reaction assumption.By coupling this model into the viscous wall boundary condition of computational fluid dynamics(CFD)solver,the influences of metal material catalytic properties on heat and mass transfer into thermal protection materials were numerically investigated.Numerical results showed that atomic oxygen recombination catalyzed by surface material accounts for a major contribution to aerodynamic heating and thus variation in recombination rates from different materials leads to the significant difference in surface heat fluxes.From a comparative analysis of various materials,the catalytic activity increases from the inert platinum(Pt)to nickel(Ni)and finally to the active copper(Cu).As a result,the catalytic heating on Cu surface was more than twice of that on Pt surface.Further parametrical research revealed that the proper layout of inert material at the nose of aircraft could prevent stagnation catalytic heating from thermal damage by carrying near-wall dissociated atoms from the stagnation zone downstream.The material-relied heterogeneous catalysis mechanism in this study provides some technical support for the thermal protection system design of hypersonic aircraft.展开更多
This paper presents the extension of the global description approach of a discontinuous function, which is proposed in the previous paper, to a spectral domain decomposition method. This multi-domain spectral immersed...This paper presents the extension of the global description approach of a discontinuous function, which is proposed in the previous paper, to a spectral domain decomposition method. This multi-domain spectral immersed interlace method(IIM) divides the whole computation domain into the smooth and discontinuous parts. Fewer points on the smooth domains are used via taking advantage of the high accuracy property of the spectral method, but more points on the discontinuous domains are employed to enhance the resolution of the calculation. Two that the domain decomposition technique can placed around the discontinuity. The present reached, in spite of the enlarged computational discontinuous problems are tested to verify the present method. The results show reduce the error of the spectral IIM, especially when more collocation points are method is t:avorable for the reason that the same level of the accuracy can be domain.展开更多
The presence of Dirac delta function in differential equation can lead to a discontinuity,which may degrade the accuracy of related numerical methods.To improve the accuracy,a secondorder numerical method for elliptic...The presence of Dirac delta function in differential equation can lead to a discontinuity,which may degrade the accuracy of related numerical methods.To improve the accuracy,a secondorder numerical method for elliptic equations with singular sources is introduced by employing a local kernel flter.In this method,the discontinuous equation is convoluted with the kernel function to obtain a more regular one.Then the original equation is replaced by this fltered equation around the singular points,to obtain discrete numerical form.The unchanged equations at the other points are discretized by using a central difference scheme.1D and 2D examples are carried out to validate the correctness and accuracy of the present method.The results show that a second-order of accuracy can be obtained in the fltering framework with an appropriate integration rule.Furthermore,the present method does not need any jump condition,and also has extremely simple form that can be easily extended to high dimensional cases and complex geometry.展开更多
基金financial support of the National Key Research and Development Plan of China through the project(No.2019YFA0405202)National Natural Science Foundation of China through the project(No.12072361)。
文摘Surface heterogeneous catalysis in a high-enthalpy dissociated environment leads to a remarkable enhancement of aerodynamic heating into the thermal protection surface of hypersonic aircraft.To more accurately predict this catalytic heating,a kinetic catalytic model was constructed.This model involved four elementary reactions,the rates of which were determined on mean-field approximation and surface steady-state reaction assumption.By coupling this model into the viscous wall boundary condition of computational fluid dynamics(CFD)solver,the influences of metal material catalytic properties on heat and mass transfer into thermal protection materials were numerically investigated.Numerical results showed that atomic oxygen recombination catalyzed by surface material accounts for a major contribution to aerodynamic heating and thus variation in recombination rates from different materials leads to the significant difference in surface heat fluxes.From a comparative analysis of various materials,the catalytic activity increases from the inert platinum(Pt)to nickel(Ni)and finally to the active copper(Cu).As a result,the catalytic heating on Cu surface was more than twice of that on Pt surface.Further parametrical research revealed that the proper layout of inert material at the nose of aircraft could prevent stagnation catalytic heating from thermal damage by carrying near-wall dissociated atoms from the stagnation zone downstream.The material-relied heterogeneous catalysis mechanism in this study provides some technical support for the thermal protection system design of hypersonic aircraft.
基金National Natural Science Foundation of China(51076006)
文摘This paper presents the extension of the global description approach of a discontinuous function, which is proposed in the previous paper, to a spectral domain decomposition method. This multi-domain spectral immersed interlace method(IIM) divides the whole computation domain into the smooth and discontinuous parts. Fewer points on the smooth domains are used via taking advantage of the high accuracy property of the spectral method, but more points on the discontinuous domains are employed to enhance the resolution of the calculation. Two that the domain decomposition technique can placed around the discontinuity. The present reached, in spite of the enlarged computational discontinuous problems are tested to verify the present method. The results show reduce the error of the spectral IIM, especially when more collocation points are method is t:avorable for the reason that the same level of the accuracy can be domain.
基金supported by the National Natural Science Foundation in China(Grant Nos.51076006,11202013)BUAA SJP ‘‘111’’ Program(Grant No.B08009)+1 种基金the National Basic Research Program of China(2012CB720200)the Open Research Fund of MOE Key Lab-oratory of High-speed Railway Engineering,Southwest Jiao-tong University and the European Community’s Seventh Framework Program(FP7/2007-2013)under Grant agreement 225967‘‘NextMuSE’’
文摘The presence of Dirac delta function in differential equation can lead to a discontinuity,which may degrade the accuracy of related numerical methods.To improve the accuracy,a secondorder numerical method for elliptic equations with singular sources is introduced by employing a local kernel flter.In this method,the discontinuous equation is convoluted with the kernel function to obtain a more regular one.Then the original equation is replaced by this fltered equation around the singular points,to obtain discrete numerical form.The unchanged equations at the other points are discretized by using a central difference scheme.1D and 2D examples are carried out to validate the correctness and accuracy of the present method.The results show that a second-order of accuracy can be obtained in the fltering framework with an appropriate integration rule.Furthermore,the present method does not need any jump condition,and also has extremely simple form that can be easily extended to high dimensional cases and complex geometry.