Analysis on capabilities of density-based solvers within OpenFOAM to distinguish aerothermal variables in difusion boundary layer

Open source feld operation and manipulation（OpenFOAM）is one of the most prevalent open source computational fluid dynamics（CFD）software.It is very convenient for researchers to develop their own codes based on the class library toolbox within OpenFOAM.In recent years,several density-based solvers within OpenFOAM for supersonic/hypersonic compressible flow are coming up.Although the capabilities of these solvers to capture shock wave have already been verifed by some researchers,these solvers still need to be validated comprehensively as commercial CFD software.In boundary layer where diffusion is the dominant transportation manner,the convective discrete schemes＇capability to capture aerothermal variables,such as temperature and heat flux,is different from each other due to their own numerical dissipative characteristics and from viewpoint of this capability,these compressible solvers within OpenFOAM can be validated further.In this paper,frstly,the organizational architecture of density-based solvers within OpenFOAM is analyzed.Then,from the viewpoint of the capability to capture aerothermal variables,the numerical results of several typical geometrical felds predicted by these solvers are compared with both the outcome obtained from the commercial software Fastran and the experimental data.During the computing process,the Roe,AUSM＋（Advection Upstream Splitting Method）,and HLLC（Harten-Lax-van Leer-Contact）convective discrete schemes of which the spatial accuracy is 1st and 2nd order are utilized,respectively.The compared results show that the aerothermal variables are in agreement with results generated by Fastran and the experimental data even if the1st order spatial precision is implemented.Overall,the accuracy of these density-based solvers can meet the requirement of engineering and scientifc problems to capture aerothermal variables in diffusion boundary layer.

A new surface catalytic model for silica-based thermal protection material for hypersonic vehicles

Silica-based materials are widely employed in the thermal protection system for hypersonic vehicles, and the investigation of their catalytic characteristics is crucially important for accurate aerothermal heating prediction. By analyzing the disadvantages of Norman＇s high and low temperature models, this paper combines the two models and proposes an eight-reaction combined surface catalytic model to describe the catalysis between oxygen and silica surface. Given proper evaluation of the parameters according to many references, the recombination coefficient obtained shows good agreement with experimental data. The catalytic mechanisms between oxygen and silica surface are then analyzed. Results show that with the increase of the wall temperature, the dominant reaction contributing to catalytic coefficient varies from Langmuir Hinshelwood （LH） recombination （Tw 〈 620 K） to Eley Rideal （ER） replacement （620 K 〈 Tw 〈 1350 K）, and then to 02 desorption （Tw 〉 1350 K）. The surface coverage of chemisorption areas varies evidently with the dominant reactions in the high temperature （HT） range, while the surface coverage of physisorption areas varies within quite low temperature （LT） range （Tw 〈 250 K）. Recommended evaluation of partial parameters is also given.