Transactions of Nanjing University of Aeronautics & Astronautics
+Advanced Search
引
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。

Effects of Chemical Reactions in the Hypersonic Reacting Flow around Blunt Bodies

  • ZHAO Faming

    Affiliation： College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, P. R. China

    ×
  • WANG Jiangfeng

    Affiliation： College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, P. R. China

    Role：Corresponding author

    Email：wangjf@nuaa.edu.cn.

    Profile：E-mail: wangjf@nuaa.edu.cn.

    ×
  • FAN Xiaofeng

    Affiliation： College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, P. R. China

    ×
  • YANG Tianpeng

    Affiliation： College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, P. R. China

    ×

College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, P. R. China

CLC: V211.3

DOI:http://dx.doi.org/10.16356/j.1005⁃1120.2019.01.006

无

中图分类号: V211.3

DOI:http://dx.doi.org/10.16356/j.1005⁃1120.2019.01.006

参考文献 1
CANDLERG, MACCORMACKR . The computation of hypersonic ionized flow in chemical and thermal non-equilibrium [C]// 26th AIAA Aerospace Sciences Meeting. Reno: AIAA, 1988.
查找原文
参考文献 2
GAOZhenxun,JIANGChongwen,ChunhianLEE .Aeroheating study of hypersonic chemical nonequilibrium flow around a reentry blunt body [C]//AIAA SPACE Conference and Exposition.San Diego:AIAA,2014.
查找原文
参考文献 3
PARKC . Review of chemical-kinetic problems of future NASA missions, I: Earth entries [J]. Journal of Thermophysics and Heat Transfer, 1993, 7(3): 385-398.
查找原文
参考文献 4
GNOFFOP, GUPTAR, SHINNJ . Conservation equations and physical models for hypersonic air flow in thermal and chemical nonequilibrium [R]. NASA-TP- 2867, 1989.
查找原文
参考文献 5
BUSSINGT, EBERHARDTS . Chemistry associated with hypersonic vehicles [J]. Journal of Thermophysics and Heat Transfer, 1989, 3(3): 245-253.
查找原文
参考文献 6
REDDYD, SINHAK . Effect of chemical reaction rates on aeroheating predictions of reentry flows [J]. Journal of Thermophysics and Heat Transfer, 2011, 25(1): 21-33.
查找原文
参考文献 7
GUPTAR, YOS J, THOMPSONR . A review of reaction rates and thermodynamic and transport properties for an 11-species air model for chemical and thermal non-equilibrium calculations to 30 000 K [R]. NASA-TM- 101528, 1989.
查找原文
参考文献 8
OUYANGShuiwu, XIEZhongqiang . Navier-Stokes computation of non-equilibrium chemically reacting flow fields [J]. Chinese Journal of Computational Physics, 1997, 14(1): 7-12. (in Chinese)
查找原文
参考文献 9
PARKC . Assessment of two-temperature kinetic model for ionizing air [J]. Journal of Thermophysics and Heat Transfer, 1989, 3(3): 233-244.
查找原文
参考文献 10
RONALDS . A century of ramjet propulsion technology evolution [J]. Journal of Propulsion and Power, 2004, 20(1): 27-58.
查找原文
参考文献 11
WANGJiangfeng, WUYizhao, JIWeidong, et al . Progress in numerical simulation techniques of hypersonic aerodynamic problems [J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(1): 159-175. (in Chinese)
查找原文
参考文献 12
CAIGuobiao, XUDajun . Hypersonic vehicle Technology[M]. Beijing: Chinese Science Press, 2012. (in Chinese)
查找原文
参考文献 13
McClintonC, RauschV, NguyenL, et al . Preliminary X-43 flight test results[J]. Acta Astronautica, 2005, 57(2): 266-276
查找原文
参考文献 14
OUYANGShuiwu, XIEZhongqiang . High temperature nonequilibrium air flow [M]. Beijing: Chinese National Defense Industry Press, 2001. (in Chinese)
查找原文
参考文献 15
ANDERSONJ . Hypersonic and high temperature gas dynamics [M]. New York: McGraw-Hill, 1989.
查找原文
参考文献 16
ZHANGXianghong, WUYizhao, WANGJiangfeng . Numerical simulation of thermo-chemical non- equilibrium hypersonic flow using HLLE+ scheme [J]. Journal of Nanjing University of Aeronautics and Astronautics, 2011, 43(2): 154-158. (in Chinese)
查找原文
参考文献 17
SHANGH, CHENY, PAULL, et al . Investigation of chemical kinetics integration algorithms for reacting flow [C]//33rd Aerospace Sciences Meeting and Exhibit. Reno: AIAA, 1995.
查找原文
参考文献 18
LEHRH . Experiments on shock-induced combustion [J]. ACTA Astronautica, 1972, 17(4): 589-597.
查找原文
参考文献 19
SOERTRISNOM, IMAYS . Simulation of the flow field of a ram accelerator [C]//27th Joint Propulsion Conference. Sacramento: AIAA, 1991.
查找原文
参考文献 20
HAYASHIK, ASO S, TANIY . Numerical study of thermal protection system by opposing jet [C]//43rd AIAA Aerospace Sciences Meeting and Exhibit. Reno: AIAA, 2005.
查找原文
目录 contents

    Abstract

    The effects of chemical reactions in the hypersonic reacting flow are investigated using an integrated algorithm considering simultaneously two different reaction mechanisms, i.e., including the high temperature air non-equilibrium chemical reactions and the H2-air combustion reactions. The program is validated by the air non-equilibrium flow at Mach number of 25.9 with the RAM C-II configuration and the shock-induced combustion flow at Mach number of 4.512 6 around a sphere, respectively. Furthermore, the mixed reacting flow with the Mach number of 10.0 with an opposing jet of hydrogen is numerically analyzed. The results show that the program is reliable, and the effects of chemical reactions engender in the decrease of peak temperature along characteristic lines, as well as on the surface. The production of water is augmented in the region with high ratio of oxygen to hydrogen and weakened in the area with low ratio of oxygen to hydrogen by the air chemical non-equilibrium effects.

    摘要

    暂无

  • 0 Introduction

    The high-temperature air surrounding the hypersonic vehicles undergoes chemical non-equilibrium processes, such as the dissociation, vibrational excitation and ionization during the reentering period[1,2,3,4]. The effects of air chemical non-equilibrium have to be considered in the numerical prediction and engineering estimation[5,6]. Advances in the numerical simulation of air non-equilibrium flow have been obtained and hence greatly enhanced the analysis of difference and influence derived from various chemical reaction models and temperature models[7,8,9]. Meanwhile, the hypersonic cruise vehicles are equipped with the propulsion system and reaction control system of rocket or scramjet engine in general. Therefore, the effects of combustion of the hydrogen or hydrocarbon fuel[10] should be considered, too. Commercial software, such as FLUENT, and self-developed program， such as UNCOMB[11], can be used to compute the combustion flow.

    Nevertheless, the numerical investigations of the air non-equilibrium flow and combustion flow are conducted independently at present. The inter-action effects of the air chemical non-equilibrium and combustion are unpredictable. For example, changing of the concentrations of oxygen and oxygen atoms near the inlet of the air-breathing vehicles because of the air chemical reactions affect the combustion reactions inside the combustion chamber. Thus, it is necessary to consider simultaneously the effects of the air chemical non-equilibrium and combustion of the advanced hypersonic vehicles, such as X-43, X-51, HTV-2, and some other air-breathing missiles[12,13].

    In this paper, a program for computing the hypersonic reacting flow is developed,considering simultaneously the air chemical reactions and combustion reactions . As no experimental or computational data can be found about the mixed reacting flow, the hypersonic air non-equilibrium flow around the RAM-C II configuration and the supersonic shock-induced combustion flow around a sphere are conducted respectively for the validation of the numerical algorithm. Furthermore, the hypersonic reacting flow with an opposing jet of hydrogen around a sphere is performed to study the interaction effects of air chemical non-equilibrium and combustion.

  • 1 Computational Method

  • 1.1 Governing equations

    The axisymmetric multi-species Navier–Stokes (N-S) equations in conservation form can be expressed in the Cartesian coordinates as[14]

    ∂U∂t+∂E∂x+∂F∂y+H=∂Ev∂x+∂Fv∂y+Hv+ω
    (1)
    U=ρ1⋮ρnsρuρvet,E=ρ1u⋮ρnsuρu2+pρuv(et+p)u,F=ρ1v⋮ρnsvρuvρv2+p(et+p)v
    (2)
    ω=ω˙1⋮ω˙ns000,Ev=ρ1D1∂c1∂x⋮ρnsDns∂cns∂xτxxτxyuτxx+vτxy+λ∂T∂xFv=ρ1D1∂c1∂y⋮ρnsDns∂cns∂yτyxτyyuτyx+vτyy+λ∂T∂y
    (3)
    H=1yρ1v⋮ρnsvρuvρv2(et+p)v,Hv=1yρ1D1∂c1∂y⋮ρnsDns∂cns∂yτyxτyy-τθθuτyx+vτyy+λ∂T∂y
    (4)
    τxx=23μ2∂u∂x-∂v∂y-uy
    (5)
    τyy=23μ2∂v∂y-∂u∂x-vy
    (6)
    τxy=μ∂u∂y+∂v∂x
    (7)
    τyy-τθθ=2μ∂v∂y-vy
    (8)

    where ρ,p,T,et,μ,λ are the density, pressure, temperature, total energy per unit volume, viscosity, and thermal conductivity of the gas mixture, respectively; ρi,ci,Di are the density, mass concentration, and coefficients of mass diffusion of species, resperctively. The pressure of the gas mixture is obtained from the Dalton law[15] and the total energy per unit volume can be expressed as

    p=∑i=1nsρiRuMiT
    (9)
    et=∑i=1nsρihi+12ρ(u2+v2)-p
    (10)

    where Mi,hi are the molecular weight and static enthalpy of species; and Ru the universal gas constant.

    An integrated algorithm for computing the mixed reacting flow with air chemical reactions and H2 -air combustion reactions is developed, using the finite volume method based on the hybrid grid techniques. The central scheme with artificial dissipation is used for calculating the inviscid flux and the mul-tistage Runge-Kutta scheme is employed for tempo-ral discretisation. The turbulence model adopted here is the Spalart-Allmaras one-equation model. The program is compiled and verified on a 128-CPUs Linux system employing the memory-distributed parallel computing techniques.

  • 1.2 Chemical kinetic model

    For NR possible chemical reactions with NS species, the chemical reactions in a general form are[7]

    ∑i=1NJαi,rXi⇄kf,rkb,r∑i=1NJβi,rXi
    (11)

    where r=1,2,⋯,NR and NJ is equal to the sum of the reacting species ( NS ) plus the number of catalytic bodies; αi,r,βi,r are the stoichiometric coefficients of the species on the rth reaction. Xi is the concentration of the chemical species and catalytic bodies in moles per unit volume, and kf,r,kb,r are the rates of forward and backward reactions, respectively. The net mass production per unit volume is

    ω˙i=Mi∑r=1NR(βi,r-αi,r)kf,r∏j=1NJ(γjρ)αj,r-kb,r∏j=1NJ(γjρ)βj,r
    (12)

    where γs is the mole mass ratio.

  • 2 Results and Discussion

  • 2.1 Hypersonic air non-equilibrium flow

    The RAM C-II configuration illustrated in Ref. [16] is used in this case, with a 129×110 structured grid shown in Fig.1. The flight altitude is 71 km and the free stream is composed by the species of N2 and O2; the associated conditions are given in Table1. The temperature of the wall is maintained at 1 500 K which is assumed to be non-catalytic. The chemical reaction model adopted is the 7-species (O2, N2, O, N, NO, NO+, e-) model with seven reaction equations proposed by Gupta[7], and the corresponding reaction coefficients are given in Ref.[2].

    Fig.1 Computational grid

    Fig.1 Computational grid

    Table 1 Free stream conditions

    Ma∞ p∞/Pa T∞/K CN2 CO2
    25.94.9219.580.791 20.208 8

    Fig.2shows the comparison of the peak electron number density as a function of the normalized axial length (x/r). A value of x/r=0 represents the nose stagnation point. The experimental data shows that the electron number density is highest at the nose and falls off rapidly around the shoulder of the body, and the computational approximation follows this trend, but the fixed-wall temperature used in the calculation may be too high that would produce much more number of electrons near the wall. Fig.3shows the comparison of the normalized surface pressure distributions ( p/p∞ ) as a function of the normalized distance from the nose (D/r). The surface pressure distributions of the test match well with the Candler’s work[1].

    Fig.2 Comparison of peak electron number density

    Fig.2 Comparison of peak electron number density

    Fig.3 Comparison of surface pressure distributions

    Fig.3 Comparison of surface pressure distributions

  • 2.2 Shock-induced combustion flow

    The radius of the sphere is 7.5 mm, and a 101×301 structured grid is generated for computation (Fig.4). The free stream is composed by the species of H2,O2, and N2, with the mole ratio n(H2 )∶n(O2 )∶n(N2 )=2∶1∶3.76; the associated conditions are given in Table2. The wall is assumed to be adiabatic and non-catalytic. The chemical reaction model adopted is the 7-species (H2 , O2 , H2O, O, H, OH, N2) model with seven reaction equations by Shang[17], and the corresponding reaction coefficients are given in Ref. [17].

    Fig.4 Computational grid

    Fig.4 Computational grid

    Table 2 Free stream conditions

    Ma∞ p∞/Pa T∞/K
    4.512 642 662250

    The comparison of the temperature distributions obtained from the test is presented in Fig.5, where the uncoupled shock-deflagration phenomenon can be observed. The temperature after the shock wave is not high enough for the ignition of H2,thus the combustion occurs after a period of induction. The shock wave and the deflagration wave are separated by the induced region clearly. The position of the shock wave matches well with the experimental results[18], and the position of the deflagration wave agrees well with that in Ref. [19].

    Fig.5 Comparison of temperature distributions

    Fig.5 Comparison of temperature distributions

    Fig.6(a) shows the comparison of the pressure and temperature distributions along the stagnation line as a function of the normalized axial length x/r. A value of x/r=－1 represents the nose stagnation point. The pressure increases rapidly behind the shock wave, but the temperature increases slightly, staying almost unchanged in the induced region, and then continues to rise rapidly at the position x/r=－1.12. The pressure and temperature distributions in this test case match well with those in Ref. [19]. The comparison of the mass fraction distributions of H2, O2 and H2O along the stagnation line is demonstrated in Fig.6(b), where we can observe clearly that the combustion occurs at the position x/r=－1.12 with the decrease of mass fraction of H2 and O2 and increase of H2O. The mass fraction distributions of the species in this test case also match well with the reference data[19].

    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F006.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F007.jpg

    Fig.6 Comparison of parameter distributions along the stagnation line

  • 2.3 Mixed reacting flow around a sphere with an opposing jet of hydrogen

    The radius of the sphere and the jet spout is 0.025 m and 0.002 m, respectively, and the computational grid is a 210×129 structured grid (Fig.7). The free stream is composed by the species of N2 and O2, and the jet flow is composed by hydrogen. The associated conditions are given in Table3. The wall is assumed to be adiabatic and non-catalytic. The chemical reactions and the rate coefficients in Arrhenius formula are given in TableTable、4, and the backward rates of reaction 8 to 14 are calculated by equilibrium constants. In addition, the premixed flow without chemical reactions and the reacting flow considering only the combustion of hydrogen are also simulated under the same conditions.

    Fig.7 Computational grid

    Fig.7 Computational grid

    Table 3 Free stream and jet conditions

    Ma∞ p∞/Pa T∞/K Maj Pjs/Pa Tjs/K
    10.01 197.03226.511.03.0E+5250

    Table 4 Chemistry model for air chemical reactions and H2‑O2 combustion

    No.ReactionForward rate coefficient/(cm3/·mol-1·s-1)Backward rate coefficient/(cm3/·mol-1·s-1)
    1O2+M⇔2O+M3.61×1018 T -1.0exp(-118 028/RT)3.01×1015 T -0.5
    2N2+M⇔2N+M1.92×1017 T -0.5exp(-224 730/RT)1.09×1016 T -0.5
    3NO+M⇔N+O+M3.97×1020 T -1.5exp(-150 217/RT)1.01×1020 T -1.5
    4N2+O⇔NO+N6.75×1013exp(-74 513/RT)1.50×1013
    5NO+O⇔O2+N3.18×109 T 1.0exp(-39 144/RT)9.63×1011 T 0.5exp(-7 153/RT)
    6N2+N⇔2N+N4.15×1022 T -1.5exp(-260 297/RT)2.32×1021 T -1.5
    7N+O⇔NO++e- 9.03×109 T 0.5exp(-64 379/RT)1.80×1019 T -1.0
    8H2+O2⇔OH+OH1.70×1013exp(-48 150/RT)-
    9H+O2⇔O+OH1.42×1014exp(-16 400/RT)-
    10OH+H2⇔H2O+H3.16×107 T 1.8exp(-3 030/RT)-
    11O+H2⇔OH+H2.07×1014exp(-13 750/RT)-
    12OH+OH⇔H2O+O5.50×1013exp(-7 000/RT)-
    13H+OH+M⇔H2O+M2.21×1022 T -2 -
    14H+H+M⇔H2+M6.53×1017 T -1 -

    Fig.8shows the main features of the flow field in terms of the Mach number. The flow field is characterized by the detached shock wave, Mach disk and recompression shock wave. And the recirculation region with separation bubble is formed due to the jet interference in the hypersonic free stream. The flow field structure obtained in this paper matches well with that in Ref. [20].

    Fig.8 Mach number distribution and representative streamlines

    Fig.8 Mach number distribution and representative streamlines

    Two characteristic lines for the comparison of the parameter distributions is illustrated in Fig.9, where the first line passes through the center of the jet with the polar angle of θ=0°, and the second line crosses over the recirculation region with the polar angle of θ=15 °.

    Fig.9 Characteristic lines with the polar angle of θ

    Fig.9 Characteristic lines with the polar angle of θ

    Fig.10(a) shows the temperature distributions along the first characteristic line. The temperature increases rapidly behind the detached shock wave, and then decreases swiftly mixing the cold jet of hydrogen. The maximum value of the temperature in the premixed flow is about 4 095.9 K, while the values of the reacting flow are lower, and the temperature peaks of the reacting flow shift slightly to the negative direction of the x-axis due to the effects of chemical reactions. Figs.10(b) , (c) show the mass fraction distributions of the major species (H2, O2, H2O) and minor species (O, H) in the reacting flow, but the peak value of the mass concentration of H2 is not given in order to display clearly the variation of the mass fraction of the species with low magnitude, which approaching to the jet condition with the value of 1.0. The value of the oxygen mass fraction decreases while the hydrogen mass fraction increases in the mixing layer, and the ratio of oxygen to hydrogen (H2/O2) is about 8.59 and 26.36 of the combustion flow and mixed reacting flow, respectively at the peak position of water. Compared with the combustion flow, the maximum values of the mass fraction of water and atomic hydrogen are greater, and the peak mass fraction of the atomic oxygen is one order higher of the mixed reacting flow. The results indicate that the dissociation of oxygen molecule is intense and the production of water is augmented by the effects of air chemical non-equilibrium.

    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F011.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F012.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F013.jpg

    Fig.10 Comparison of parameter distributions along the line with polar angle of θ=0 °

    The same parameter distributions along the second characteristic line are given in Fig.11, which the peak value of the mass fraction of H2 is also not given. The maximum value of the temperature in the premixed flow is about 3 740.7 K, which is much lower than that on the first characteristic line, and the effects of chemical reactions also result in the decrease of the peak temperature. The ratio of oxygen to hydrogen is about 0.93 and 2.16 of the combustion flow and mixed reacting flow respectively at the peak position of water. Compared with the combustion flow, the maximum values of the mass fraction of water and atomic hydrogen are less, and the peak mass fraction of the atomic oxygen is greater of the mixed reacting flow. The results indicate that the dissociation of oxygen molecule is weak and the production of water is weakened by the effects of air chemical non-equilibrium.

    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F014.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F015.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F016.jpg

    Fig.11 Comparison of parameter distributions along the line with polar angle of θ=15 °

    Fig.12shows the comparison of the pressure and temperature distributionson the surface as a function of the polar angle, where the value of θ=0° means the center of the jet. The value of the peak pressure on the surface of the premixed flow is about 82.7 kPa. The effects of chemical reactions result in the decrease of peak pressure on the surface, and the pressure peaks of the reacting flow shift slightly to the downstream of the jet. The temperature peaks on the surface locate in the region of the separation bubble, and the maximum value of the temperature of the premixed flow is about 1 376.4 K. The peak temperature decreases due to the effects of chemical reactions, especially considering only the combustion of hydrogen.

    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F017.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F018.jpg

    Fig.12 Comparison of surface parameter distributions

  • 3 Conclusions

    The effects of chemical reactions in the hypersonic reacting flow are investigated using an integrated algorithm considering simultaneously the air chemical reactions and the H2-air combustion reactions. The program is validated by the hypersonic air non-equilibrium flow of the RAM C-II configuration and the shock-induced combustion flow around a sphere, respectively. After that, the mixed reacting flow with an opposing jet of hydrogen is analyzed. The results show that the program is reliable, and the chemical reactions result in the decrease of peak temperature and peak pressure on the surface under the conditions in this paper. The production of water is augmented in the region with high ratio of oxygen to hydrogen and weakened in the area with low ratio of oxygen to hydrogen by the effects of air chemical non-equilibrium.

    The air non-equilibrium flow combined with hydrocarbon‑air combustion is under investigation, and will be used to analyze the flow field of hypersonic vehicles with the multi-species reaction control system, as well as the combustion flow field inside the scramjet engine.

  • References

    • 1

      CANDLER G, MACCORMACK R . The computation of hypersonic ionized flow in chemical and thermal non-equilibrium [C]// 26th AIAA Aerospace Sciences Meeting. Reno: AIAA, 1988.

    • 2

      GAO Zhenxun,JIANG Chongwen,LEE Chunhian .Aeroheating study of hypersonic chemical nonequilibrium flow around a reentry blunt body [C]//AIAA SPACE Conference and Exposition.San Diego:AIAA,2014.

    • 3

      PARK C . Review of chemical-kinetic problems of future NASA missions, I: Earth entries [J]. Journal of Thermophysics and Heat Transfer, 1993, 7(3): 385-398.

    • 4

      GNOFFO P, GUPTA R, SHINN J . Conservation equations and physical models for hypersonic air flow in thermal and chemical nonequilibrium [R]. NASA-TP- 2867, 1989.

    • 5

      BUSSING T, EBERHARDT S . Chemistry associated with hypersonic vehicles [J]. Journal of Thermophysics and Heat Transfer, 1989, 3(3): 245-253.

    • 6

      REDDY D, SINHA K . Effect of chemical reaction rates on aeroheating predictions of reentry flows [J]. Journal of Thermophysics and Heat Transfer, 2011, 25(1): 21-33.

    • 7

      GUPTA R, YOS J, THOMPSON R . A review of reaction rates and thermodynamic and transport properties for an 11-species air model for chemical and thermal non-equilibrium calculations to 30 000 K [R]. NASA-TM- 101528, 1989.

    • 8

      OUYANG Shuiwu, XIE Zhongqiang . Navier-Stokes computation of non-equilibrium chemically reacting flow fields [J]. Chinese Journal of Computational Physics, 1997, 14(1): 7-12. (in Chinese)

    • 9

      PARK C . Assessment of two-temperature kinetic model for ionizing air [J]. Journal of Thermophysics and Heat Transfer, 1989, 3(3): 233-244.

    • 10

      RONALD S . A century of ramjet propulsion technology evolution [J]. Journal of Propulsion and Power, 2004, 20(1): 27-58.

    • 11

      WANG Jiangfeng, WU Yizhao, JI Weidong, et al . Progress in numerical simulation techniques of hypersonic aerodynamic problems [J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(1): 159-175. (in Chinese)

    • 12

      CAI Guobiao, XU Dajun . Hypersonic vehicle Technology[M]. Beijing: Chinese Science Press, 2012. (in Chinese)

    • 13

      McClinton C, Rausch V, Nguyen L, et al . Preliminary X-43 flight test results[J]. Acta Astronautica, 2005, 57(2): 266-276

    • 14

      OUYANG Shuiwu, XIE Zhongqiang . High temperature nonequilibrium air flow [M]. Beijing: Chinese National Defense Industry Press, 2001. (in Chinese)

    • 15

      ANDERSON J . Hypersonic and high temperature gas dynamics [M]. New York: McGraw-Hill, 1989.

    • 16

      ZHANG Xianghong, WU Yizhao, WANG Jiangfeng . Numerical simulation of thermo-chemical non- equilibrium hypersonic flow using HLLE+ scheme [J]. Journal of Nanjing University of Aeronautics and Astronautics, 2011, 43(2): 154-158. (in Chinese)

    • 17

      SHANG H, CHEN Y, PAUL L, et al . Investigation of chemical kinetics integration algorithms for reacting flow [C]//33rd Aerospace Sciences Meeting and Exhibit. Reno: AIAA, 1995.

    • 18

      LEHR H . Experiments on shock-induced combustion [J]. ACTA Astronautica, 1972, 17(4): 589-597.

    • 19

      SOERTRISNO M, IMAY S . Simulation of the flow field of a ram accelerator [C]//27th Joint Propulsion Conference. Sacramento: AIAA, 1991.

    • 20

      HAYASHI K, ASO S, TANI Y . Numerical study of thermal protection system by opposing jet [C]//43rd AIAA Aerospace Sciences Meeting and Exhibit. Reno: AIAA, 2005.

    • 贡献声明和致谢

    Mr. ZHAO Faming designed the study, compiled the models, conducted the analysis, interpreted the results and wrote the manuscript. Prof. WANG Jiangfeng contributed to the discussion and background of the study. Mr. FAN Xiaofeng contributed data and model components for the combustion model. Mr. YANG Tianpeng contributed data for the analysis of reactions effects. All authors commented on the manuscript draft and approved the submission.

    This work was supported by the Fundamental Research Funds for the Central Universities (No. NZ2016101) and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

    利益冲突

    The authors declare no competing interests.

  • 基本信息

    DOI:http://dx.doi.org/10.16356/j.1005⁃1120.2019.01.006

    中图分类号: V211.3

    文献标识码: A

    引用信息

    ZHAO Faming, WANG Jiangfeng, FAN Xiaofeng, et al. Effects of Chemical Reactions in the Hypersonic Reacting Flow around Blunt Bodies[J]. Transactions of Nanjing University of Aeronautics and Astronautics, 2019,36(1):71-79.

    稿件历史

    纸刊出版 : 2019-03-30
    收稿日期 : 2017-05-09
    修回日期 : 2017-09-05
    录用日期 : 2017-09-15

    • ZHAOFaming

    Affiliation： College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, P. R. China

    Profile：Mr. ZHAO Faming is a Ph.D. candidate in the College of Aerospace Engineering，Nanjing Unlversity of Aeronautls and Astronautics. His research interests focus on the numerical simulation of hypersonic reacting flows.

    WANGJiangfeng

    Affiliation： College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, P. R. China

    角 色：通讯作者

    Role：Corresponding author

    邮 箱：wangjf@nuaa.edu.cn.

    Profile：WANG Jiangfeng received his Ph.D. degree in Computational Fluid Dynamics from Nanjing University of Aeronautics and Astronautics in 2000. His research interests focus on the layout design of hypersonic vehicles, numerical techniques of complex combustion flows, and computational techniques of aerodynamic heating.E-mail: wangjf@nuaa.edu.cn.

    FANXiaofeng

    Affiliation： College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, P. R. China

    Profile：Mr. FAN Xiaofeng is a Ph.D. candidate in the College of Aerospace Engineering，Nanjing Unlversity of Aeronautls an-d Astronautics. His research interests focus on the nume-rical simulation of combustion flows of gaseous and liquid fuel.

    YANGTianpeng

    Affiliation： College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, P. R. China

    Profile：Mr. YANG Tianpeng is a Ph.D. candidate in the College of Aerospace Engineering，Nanjing Unlversity of Aeronautls and Astronautics. His research interests focus on the numerical simulation of hypersonic aerodynamics.

    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F001.jpg
    Ma∞ p∞/Pa T∞/K CN2 CO2
    25.94.9219.580.791 20.208 8
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F002.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F003.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F004.jpg
    Ma∞ p∞/Pa T∞/K
    4.512 642 662250
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F005.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F006.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F007.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F008.jpg
    Ma∞ p∞/Pa T∞/K Maj Pjs/Pa Tjs/K
    10.01 197.03226.511.03.0E+5250
    No.ReactionForward rate coefficient/(cm3/·mol-1·s-1)Backward rate coefficient/(cm3/·mol-1·s-1)
    1O2+M⇔2O+M3.61×1018 T -1.0exp(-118 028/RT)3.01×1015 T -0.5
    2N2+M⇔2N+M1.92×1017 T -0.5exp(-224 730/RT)1.09×1016 T -0.5
    3NO+M⇔N+O+M3.97×1020 T -1.5exp(-150 217/RT)1.01×1020 T -1.5
    4N2+O⇔NO+N6.75×1013exp(-74 513/RT)1.50×1013
    5NO+O⇔O2+N3.18×109 T 1.0exp(-39 144/RT)9.63×1011 T 0.5exp(-7 153/RT)
    6N2+N⇔2N+N4.15×1022 T -1.5exp(-260 297/RT)2.32×1021 T -1.5
    7N+O⇔NO++e- 9.03×109 T 0.5exp(-64 379/RT)1.80×1019 T -1.0
    8H2+O2⇔OH+OH1.70×1013exp(-48 150/RT)-
    9H+O2⇔O+OH1.42×1014exp(-16 400/RT)-
    10OH+H2⇔H2O+H3.16×107 T 1.8exp(-3 030/RT)-
    11O+H2⇔OH+H2.07×1014exp(-13 750/RT)-
    12OH+OH⇔H2O+O5.50×1013exp(-7 000/RT)-
    13H+OH+M⇔H2O+M2.21×1022 T -2 -
    14H+H+M⇔H2+M6.53×1017 T -1 -
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F009.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F010.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F011.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F012.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F013.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F014.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F015.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F016.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F017.jpg
    html/njhkhten/201901006/alternativeImage/4471416c-513b-493b-a5bc-87f6544995b7-F018.jpg

    Fig.1 Computational grid

    Table 1 Free stream conditions

    Fig.2 Comparison of peak electron number density

    Fig.3 Comparison of surface pressure distributions

    Fig.4 Computational grid

    Table 2 Free stream conditions

    Fig.5 Comparison of temperature distributions

    Fig.6 Comparison of parameter distributions along the stagnation line

    Fig.6 Comparison of parameter distributions along the stagnation line

    Fig.7 Computational grid

    Table 3 Free stream and jet conditions

    Table 4 Chemistry model for air chemical reactions and H2‑O2 combustion

    Fig.8 Mach number distribution and representative streamlines

    Fig.9 Characteristic lines with the polar angle of θ

    Fig.10 Comparison of parameter distributions along the line with polar angle of θ=0 °

    Fig.10 Comparison of parameter distributions along the line with polar angle of θ=0 °

    Fig.10 Comparison of parameter distributions along the line with polar angle of θ=0 °

    Fig.11 Comparison of parameter distributions along the line with polar angle of θ=15 °

    Fig.11 Comparison of parameter distributions along the line with polar angle of θ=15 °

    Fig.11 Comparison of parameter distributions along the line with polar angle of θ=15 °

    Fig.12 Comparison of surface parameter distributions

    Fig.12 Comparison of surface parameter distributions

    image /

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

  • References

    • 1

      CANDLER G, MACCORMACK R . The computation of hypersonic ionized flow in chemical and thermal non-equilibrium [C]// 26th AIAA Aerospace Sciences Meeting. Reno: AIAA, 1988.

    • 2

      GAO Zhenxun,JIANG Chongwen,LEE Chunhian .Aeroheating study of hypersonic chemical nonequilibrium flow around a reentry blunt body [C]//AIAA SPACE Conference and Exposition.San Diego:AIAA,2014.

    • 3

      PARK C . Review of chemical-kinetic problems of future NASA missions, I: Earth entries [J]. Journal of Thermophysics and Heat Transfer, 1993, 7(3): 385-398.

    • 4

      GNOFFO P, GUPTA R, SHINN J . Conservation equations and physical models for hypersonic air flow in thermal and chemical nonequilibrium [R]. NASA-TP- 2867, 1989.

    • 5

      BUSSING T, EBERHARDT S . Chemistry associated with hypersonic vehicles [J]. Journal of Thermophysics and Heat Transfer, 1989, 3(3): 245-253.

    • 6

      REDDY D, SINHA K . Effect of chemical reaction rates on aeroheating predictions of reentry flows [J]. Journal of Thermophysics and Heat Transfer, 2011, 25(1): 21-33.

    • 7

      GUPTA R, YOS J, THOMPSON R . A review of reaction rates and thermodynamic and transport properties for an 11-species air model for chemical and thermal non-equilibrium calculations to 30 000 K [R]. NASA-TM- 101528, 1989.

    • 8

      OUYANG Shuiwu, XIE Zhongqiang . Navier-Stokes computation of non-equilibrium chemically reacting flow fields [J]. Chinese Journal of Computational Physics, 1997, 14(1): 7-12. (in Chinese)

    • 9

      PARK C . Assessment of two-temperature kinetic model for ionizing air [J]. Journal of Thermophysics and Heat Transfer, 1989, 3(3): 233-244.

    • 10

      RONALD S . A century of ramjet propulsion technology evolution [J]. Journal of Propulsion and Power, 2004, 20(1): 27-58.

    • 11

      WANG Jiangfeng, WU Yizhao, JI Weidong, et al . Progress in numerical simulation techniques of hypersonic aerodynamic problems [J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(1): 159-175. (in Chinese)

    • 12

      CAI Guobiao, XU Dajun . Hypersonic vehicle Technology[M]. Beijing: Chinese Science Press, 2012. (in Chinese)

    • 13

      McClinton C, Rausch V, Nguyen L, et al . Preliminary X-43 flight test results[J]. Acta Astronautica, 2005, 57(2): 266-276

    • 14

      OUYANG Shuiwu, XIE Zhongqiang . High temperature nonequilibrium air flow [M]. Beijing: Chinese National Defense Industry Press, 2001. (in Chinese)

    • 15

      ANDERSON J . Hypersonic and high temperature gas dynamics [M]. New York: McGraw-Hill, 1989.

    • 16

      ZHANG Xianghong, WU Yizhao, WANG Jiangfeng . Numerical simulation of thermo-chemical non- equilibrium hypersonic flow using HLLE+ scheme [J]. Journal of Nanjing University of Aeronautics and Astronautics, 2011, 43(2): 154-158. (in Chinese)

    • 17

      SHANG H, CHEN Y, PAUL L, et al . Investigation of chemical kinetics integration algorithms for reacting flow [C]//33rd Aerospace Sciences Meeting and Exhibit. Reno: AIAA, 1995.

    • 18

      LEHR H . Experiments on shock-induced combustion [J]. ACTA Astronautica, 1972, 17(4): 589-597.

    • 19

      SOERTRISNO M, IMAY S . Simulation of the flow field of a ram accelerator [C]//27th Joint Propulsion Conference. Sacramento: AIAA, 1991.

    • 20

      HAYASHI K, ASO S, TANI Y . Numerical study of thermal protection system by opposing jet [C]//43rd AIAA Aerospace Sciences Meeting and Exhibit. Reno: AIAA, 2005.

    • Articles

    WeChat

    Mobile website

    Website Copyright © Transactions of Nanjing University of Aeronautics & Astronautics
    Supported by:Beijing E-Tiller Technology Development Co., Ltd.