搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

热化学模型对高超声速磁流体控制数值模拟影响分析

丁明松 江涛 董维中 高铁锁 刘庆宗 傅杨奥骁

引用本文:
Citation:

热化学模型对高超声速磁流体控制数值模拟影响分析

丁明松, 江涛, 董维中, 高铁锁, 刘庆宗, 傅杨奥骁

Numerical analysis of influence of thermochemical model on hypersonic magnetohydrodynamic control

Ding Ming-Song, Jiang Tao, Dong Wei-Zhong, Gao Tie-Suo, Liu Qing-Zong, Fu Yang-Ao-Xiao
PDF
HTML
导出引用
  • 针对等离子体流场的模拟准确性问题及其对高超声速磁流体控制的影响, 通过数值求解三维非平衡Navier-Stokes流场控制方程和Maxwell电磁场控制方程, 建立了三维低磁雷诺数磁流体数值模拟方法及程序, 分析了不同空气组分化学反应模型和壁面有限催化效率等因素对高超声速磁流体控制的影响. 研究表明: 不同空气组分化学反应模型对高超声速磁流体流场结构、气动力/热特性控制的影响不容忽视; 对于本文计算条件, Park化学反应模型在组分模型一致性、等离子体模拟准确性等方面具有一定优势; 磁控热防护效果, 受壁面有限催化复合系数影响较大, 两者呈非线性关系, 不同表面区域差异较大; 磁场对磁阻力伞及其磁阻力特性影响, 受壁面催化效应的影响相对较小.
    In hypersonic flow, extremely high temperature due to shock aerodynamic heating leads plasma flow to form. By adding energy and momentum to the plas-ma flow field through the magnetic field on aircraft, the control of plasma flow field around aircraft can be realized. This has a broad prospect of applications in hypersonic aerodynamic control, aerothermal protection, and plasma distribution adjustment. Very recently, chemical reaction model and thermodynamic model were suggested to study the hypersonic magnetohydrodynamic control. However, the influence of different models and surface catalytic efficiency on hypersonic magnetohydrodynamic control are rarely analyzed in depth. In this study, a comparison of different chemical reaction models and the influence of surface catalytic efficiency are discussed. Three-dimensional (3D) nu-merical simulation method and program of extra magnetic field coupled with reentry plasma flow under the assumption of low magnetic Reynolds number are developed by solving 3D chemical non-equilibrium Navier-Stokes equations and Maxwell electromagnetic field governing equations. Based on this method, the influence of different gas component models, chemical reaction models, and surface catalytic efficiency on hypersonic magnetohydrodynamic control are analyzed. The results show that the conductivity of plasma, calculated by different gas component models and chemical reaction models, can be quite different from each other, thus can influence the accurate study on the structure of hypersonic magneto flow field as well as the aerothermal and aerodynamic characteristics. Based on the calculation conditions in this paper, the Park chemical model has advantages in the consistency and accuracy in numerical simulation. The magnetic thermal protection is greatly influenced by the surface catalytic efficiency and the correlation between the magnetic thermal protection and the surface catalytic efficiency is nonlinear and can be quite different in different region. As the surface catalytic efficiency increases, the influence of magnetic field on heat flux at stagnation point drops drastically, then increases slowly, which is a joint result of thermal conduction and chemical component diffusion. The influence of magnetic field on magnetohydrodynamic drag character is less affected by the surface catalytic efficiency. As the catalytic efficiency increases, the influence of magnetic field on magnetohydrodynamic drag character drops slowly.
      通信作者: 董维中, dms2008@qq.com
      Corresponding author: Dong Wei-Zhong, dms2008@qq.com
    [1]

    乐嘉陵 2005 再入物理(北京: 国防工业出版社)第23−43页

    Le J L 2005 Reentry Physics (Beijing: National Defence Industry Press)pp23−43 (in Chinese)

    [2]

    Maccormack R W 2005 36th AIAA Plasmadynamics and Lasers Conference Toronto, Ontario Canada, June 6-9, 2005 p4780

    [3]

    潘勇 2007博士学位论文 (南京: 南京航空航天大学)

    Pan Y 2007 Ph. D. Dissertation (Nanjing: Nanjing University of Aeronautics and Astronautics) (in Chinese)

    [4]

    董维中 1996 博士学位论文 (北京: 北京航空航天大学)

    Dong W Z 1996 Ph. D. Dissertation (Beijing: Beihang University) (in Chinese)

    [5]

    Bisek N J, Boyd I D 2010 48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition Orlando, Florida, January 4−7, 2010 p227

    [6]

    Bisek N J, Poggie J 2011 49th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition Orlando, Florida, January 4−7, 2011 p897

    [7]

    Chernyshev A, Kurakin Y, Schmidt A 2012 18th AIAA/3AF International Space Planes and Hypersonic Systems and Technologies Conference Tours, France, September 24−28, 2012 p5974

    [8]

    Andrea C, Carlo A B 2012 43rd AIAA Plasmadynamics and Lasers Conference New Orleans, Louisiana, June 25−28, 2012 p2733

    [9]

    Fujino T, Ishikawa M 2013 44th AIAA Plasmadynamics and Lasers Conference San Diego, CA, June 24−27, 2013 p3000

    [10]

    Masuda K, Shimosawa Y, Fujino T 2015 46th AIAA Plasmadynamics and Lasers Conference Dallas, TX, June 22−26, 2015 p3366

    [11]

    赫新, 陈坚强, 邓小刚 2005 空气动力学学报 23 267Google Scholar

    He X, Chen J Q, Deng X G 2005 Acta Aerodyn. Sin. 23 267Google Scholar

    [12]

    田正雨 2008 博士学位论文(长沙: 国防科学技术大学)

    Tian Z Y 2008 Ph. D. Dissertation (Changsha: National University of Defense Technology) (in Chinese)

    [13]

    陈刚, 张劲柏, 李椿萱 2008 力学学报 40 752

    Chen G, Zhang J B, Li C X 2008 Chin. J. Theoret. Appl. Mech. 40 752

    [14]

    黄富来, 黄护林 2009 航空学报 30 183

    Huang F L, Huang H L 2009 Acta Aeronaut. Astronaut. Sin. 30 183

    [15]

    黄浩, 黄护林, 张喜东, 张义宁, 刘振德 2013 推进技术 34 706

    Huang H, Huang H L, Zhang X D, Zhang Y N, Liu Z D 2013 J. Propuls. Technol. 34 706

    [16]

    何淼生, 杨文将, 郑小梅, 刘宇 2013 航空动力学报 28 365

    He M S, Yang W J, Zheng X M, Liu Y 2013 J. Aerosp. Power 28 365

    [17]

    卜少科, 薛雅心 2014 现代电子技术 37 137Google Scholar

    Bu S K, Xue Y X 2014 Modern Electron. Tech. 37 137Google Scholar

    [18]

    李开, 刘伟强 2016 物理学报 65 064701

    Li K, Liu W Q 2016 Acta Phys. Sin. 65 064701

    [19]

    李开, 柳军, 刘伟强 2017 物理学报 66 084702Google Scholar

    Li K, Liu J, Liu W Q 2017 Acta Phys. Sin. 66 084702Google Scholar

    [20]

    姚霄, 刘伟强, 谭建国 2018 物理学报 67 174702Google Scholar

    Yao X, Liu W Q, Tan J G 2018 Acta Phys. Sin. 67 174702Google Scholar

    [21]

    高铁锁, 李椿萱, 董维中, 张巧芸 2002 空气动力学学报 20 184Google Scholar

    Gao T S, Li C X, Dong W Z, Zhang Q Y 2002 Acta Aerodyn. Sin. 20 184Google Scholar

    [22]

    高铁锁, 董维中, 丁明松, 江涛 2013 空气动力学学报 31 541Google Scholar

    Gao T S, Dong W Z, Ding M S, Jiang T 2013 Acta Aerodyn. Sin. 31 541Google Scholar

    [23]

    高铁锁, 董维中, 江涛, 丁明松, 刘庆宗 2016 宇航学报 37 1193Google Scholar

    Gao T S, Dong W Z, Jiang T, Ding M S, Liu Q Z 2016 J. Astronaut. 37 1193Google Scholar

    [24]

    丁明松, 江涛, 董维中, 高铁锁, 刘庆宗 2017 航空学报 38 121030

    Ding M S, Jiang T, Dong W Z, Gao T S, Liu Q Z 2017 Acta Aeronaut. Astronaut. Sin. 38 121030

    [25]

    Dunn M G, Kang S W 1973 NASA CR-2232

    [26]

    Park C 1993 J. Thermophys. Heat Transfer 7 385Google Scholar

    [27]

    Gupta R N, Yos J M, Thompson R A 1990 NASA RP-1232

    [28]

    Gokcen T 1995 33rd AIAA Aerospace Sciences Meeting and Exhibit Reno, NV, January 9−12, 1995 p156

    [29]

    丁明松, 董维中, 高铁锁, 刘庆宗 2018 航空学报 39 12588

    Ding M S, Dong W Z, Gao T S, Liu Q Z 2018 Acta Aeronaut. Astronaut. Sin. 39 12588

    [30]

    丁明松, 董维中, 高铁锁, 江涛, 刘庆宗 2017 宇航学报 38 1361

    Ding M S, Dong W Z, Gao T S, Jiang T, Liu Q Z 2017 J. Astronaut. 38 1361

    [31]

    Yasunori N, Hirotaka O 2012 43rd AIAA Plasmadynamics and Lasers Conference New Orleans, Louisiana, June 25−28, 2012 p2734

    [32]

    Candler G V, MacCormack R W 1988 26th AIAA Aerospace Sciences Meeting Reno, Nevada, January 11−14, 1988 p511

  • 图 1  磁场配置(${B_0} = 0.2\;{\rm{T}}$)

    Fig. 1.  Magnetic field (${B_0} = 0.2\;{\rm{T}}$).

    图 2  球柱阻力特性 (a) 磁阻力伞; (b) 阻力系数

    Fig. 2.  Drag characteristics of ball-column model: (a) Magneto-resistance parachute; (b) drag coefficient.

    图 3  不同网格条件下的参数 (a) 电导率; (b) 电子数密度; (c) 热流

    Fig. 3.  Parameter of different grids: (a) Conductivity; (b) ele-ctronic number density; (c) heat flux.

    图 4  驻点线参数分布(${B_0} = 0\;{\rm{T}}$) (a) 电导率; (b) 电子数密度; (c) 温度

    Fig. 4.  Parameter along stagnation line(${B_0} = 0\;{\rm{T}}$): (a) Conductivity; (b) electronic number density; (c) temperature

    图 5  驻点线温度和电子数密度分布(${B_0} = 0.5\;{\rm{ T}}$) (a) 温度; (b) 电子数密度

    Fig. 5.  Temperature and electronic number density along stagnation line (${B_0} = 0.5\;{\rm{T}}$): (a) Temperature; (b) electronic number density.

    图 6  表面热流 (a) Dunn-Kang 模型; (b) Park模型; (c) Gupta模型

    Fig. 6.  Heat flux on wall: (a) Dunn-Kang model; (b) Park model; (c) Gupta model.

    图 7  不同模型的电子数密度 (a) 7组分; (b) 11组分

    Fig. 7.  Electron number density: (a) 7 speices; (b) 11 speices.

    图 8  流场结构云图 (a) Case 1; (b) case 2

    Fig. 8.  Flow structure map: (a) Case 1; (b) case 2.

    图 9  表面压力和热流 (a) 压力; (b) 热流

    Fig. 9.  Surface pressure and heat flux: (a) Pressure; (b) heat flux.

    图 10  不同壁面催化系数条件下驻点线参数分布 (${B_0} = 0\;{\rm{T}}$) (a) 电子数密度; (b) 电导率

    Fig. 10.  Parameter along stagnation line of different catalytic efficiency (${B_0} = 0\;{\rm{T}}$): (a) Electronic number density; (b) conductivity

    图 11  不同催化系数条件下阻力系数及磁场作用效果 (a) 阻力系数; (b) 阻力系数增大百分比

    Fig. 11.  Effect of different catalytic efficiency conditions: (a) Drag coefficient; (b) drag coefficient increment ratio.

    图 12  不同催化复合系数条件下热流受磁场作用的影响 (a) 表面热流; (b) 驻点热流; (c) 磁场效果

    Fig. 12.  Effect of different catalytic efficiency conditions: (a) Surface heat flux; (b) heat flux at stagnation point; (c) magnetic effect.

    图 13  不同催化复合系数条件下驻点线原子组分分布 (a) $\alpha = 0.001$; (b) $\alpha = 0.75$; (c) $\alpha = 0.75$

    Fig. 13.  Atomic components along stagnation line: (a) $\alpha = 0.001$; (b) $\alpha = 0.75$; (c) $\alpha = 0.75$.

    表 1  表面复合反应

    Table 1.  Surface combinative reaction.

    r化学反应式r化学反应式
    1O + O$ \Rightarrow $O25N+ + e$ \Rightarrow $N
    2N + N$ \Rightarrow $N26O2+ + e$ \Rightarrow $O2
    3NO+ + e$ \Rightarrow $NO7N2+ + e$ \Rightarrow $N2
    4O+ + e$ \Rightarrow $O
    下载: 导出CSV

    表 2  球柱阻力系数(${B_0} = 0.2\;{\rm{T}}$)

    Table 2.  Drag coefficient of ball-column model(${B_0} = 0.2\;{\rm{T}}$).

    夹角$\theta $文献结果[31]本文结果差异
    1.82021.82870.47%
    45º1.66701.67060.22%
    90º1.47551.48040.33%
    下载: 导出CSV

    表 3  不同空气化学模型阻力系数

    Table 3.  Drag coefficient under different air chemical models.

    模型${B_0} = 0\;{\rm{T}}$${B_0} = 0.5\;{\rm{T}}$增大比例/%
    7Ns Dunn-Kang0.2974720.47808161
    11Ns Dunn-Kang0.2911290.812020179
    7Ns Gupta0.3010720.51869272
    11Ns Gupta0.2587720.795648207
    7Ns Park0.2900450.48288567
    11Ns Park0.2899860.48915469
    下载: 导出CSV
  • [1]

    乐嘉陵 2005 再入物理(北京: 国防工业出版社)第23−43页

    Le J L 2005 Reentry Physics (Beijing: National Defence Industry Press)pp23−43 (in Chinese)

    [2]

    Maccormack R W 2005 36th AIAA Plasmadynamics and Lasers Conference Toronto, Ontario Canada, June 6-9, 2005 p4780

    [3]

    潘勇 2007博士学位论文 (南京: 南京航空航天大学)

    Pan Y 2007 Ph. D. Dissertation (Nanjing: Nanjing University of Aeronautics and Astronautics) (in Chinese)

    [4]

    董维中 1996 博士学位论文 (北京: 北京航空航天大学)

    Dong W Z 1996 Ph. D. Dissertation (Beijing: Beihang University) (in Chinese)

    [5]

    Bisek N J, Boyd I D 2010 48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition Orlando, Florida, January 4−7, 2010 p227

    [6]

    Bisek N J, Poggie J 2011 49th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition Orlando, Florida, January 4−7, 2011 p897

    [7]

    Chernyshev A, Kurakin Y, Schmidt A 2012 18th AIAA/3AF International Space Planes and Hypersonic Systems and Technologies Conference Tours, France, September 24−28, 2012 p5974

    [8]

    Andrea C, Carlo A B 2012 43rd AIAA Plasmadynamics and Lasers Conference New Orleans, Louisiana, June 25−28, 2012 p2733

    [9]

    Fujino T, Ishikawa M 2013 44th AIAA Plasmadynamics and Lasers Conference San Diego, CA, June 24−27, 2013 p3000

    [10]

    Masuda K, Shimosawa Y, Fujino T 2015 46th AIAA Plasmadynamics and Lasers Conference Dallas, TX, June 22−26, 2015 p3366

    [11]

    赫新, 陈坚强, 邓小刚 2005 空气动力学学报 23 267Google Scholar

    He X, Chen J Q, Deng X G 2005 Acta Aerodyn. Sin. 23 267Google Scholar

    [12]

    田正雨 2008 博士学位论文(长沙: 国防科学技术大学)

    Tian Z Y 2008 Ph. D. Dissertation (Changsha: National University of Defense Technology) (in Chinese)

    [13]

    陈刚, 张劲柏, 李椿萱 2008 力学学报 40 752

    Chen G, Zhang J B, Li C X 2008 Chin. J. Theoret. Appl. Mech. 40 752

    [14]

    黄富来, 黄护林 2009 航空学报 30 183

    Huang F L, Huang H L 2009 Acta Aeronaut. Astronaut. Sin. 30 183

    [15]

    黄浩, 黄护林, 张喜东, 张义宁, 刘振德 2013 推进技术 34 706

    Huang H, Huang H L, Zhang X D, Zhang Y N, Liu Z D 2013 J. Propuls. Technol. 34 706

    [16]

    何淼生, 杨文将, 郑小梅, 刘宇 2013 航空动力学报 28 365

    He M S, Yang W J, Zheng X M, Liu Y 2013 J. Aerosp. Power 28 365

    [17]

    卜少科, 薛雅心 2014 现代电子技术 37 137Google Scholar

    Bu S K, Xue Y X 2014 Modern Electron. Tech. 37 137Google Scholar

    [18]

    李开, 刘伟强 2016 物理学报 65 064701

    Li K, Liu W Q 2016 Acta Phys. Sin. 65 064701

    [19]

    李开, 柳军, 刘伟强 2017 物理学报 66 084702Google Scholar

    Li K, Liu J, Liu W Q 2017 Acta Phys. Sin. 66 084702Google Scholar

    [20]

    姚霄, 刘伟强, 谭建国 2018 物理学报 67 174702Google Scholar

    Yao X, Liu W Q, Tan J G 2018 Acta Phys. Sin. 67 174702Google Scholar

    [21]

    高铁锁, 李椿萱, 董维中, 张巧芸 2002 空气动力学学报 20 184Google Scholar

    Gao T S, Li C X, Dong W Z, Zhang Q Y 2002 Acta Aerodyn. Sin. 20 184Google Scholar

    [22]

    高铁锁, 董维中, 丁明松, 江涛 2013 空气动力学学报 31 541Google Scholar

    Gao T S, Dong W Z, Ding M S, Jiang T 2013 Acta Aerodyn. Sin. 31 541Google Scholar

    [23]

    高铁锁, 董维中, 江涛, 丁明松, 刘庆宗 2016 宇航学报 37 1193Google Scholar

    Gao T S, Dong W Z, Jiang T, Ding M S, Liu Q Z 2016 J. Astronaut. 37 1193Google Scholar

    [24]

    丁明松, 江涛, 董维中, 高铁锁, 刘庆宗 2017 航空学报 38 121030

    Ding M S, Jiang T, Dong W Z, Gao T S, Liu Q Z 2017 Acta Aeronaut. Astronaut. Sin. 38 121030

    [25]

    Dunn M G, Kang S W 1973 NASA CR-2232

    [26]

    Park C 1993 J. Thermophys. Heat Transfer 7 385Google Scholar

    [27]

    Gupta R N, Yos J M, Thompson R A 1990 NASA RP-1232

    [28]

    Gokcen T 1995 33rd AIAA Aerospace Sciences Meeting and Exhibit Reno, NV, January 9−12, 1995 p156

    [29]

    丁明松, 董维中, 高铁锁, 刘庆宗 2018 航空学报 39 12588

    Ding M S, Dong W Z, Gao T S, Liu Q Z 2018 Acta Aeronaut. Astronaut. Sin. 39 12588

    [30]

    丁明松, 董维中, 高铁锁, 江涛, 刘庆宗 2017 宇航学报 38 1361

    Ding M S, Dong W Z, Gao T S, Jiang T, Liu Q Z 2017 J. Astronaut. 38 1361

    [31]

    Yasunori N, Hirotaka O 2012 43rd AIAA Plasmadynamics and Lasers Conference New Orleans, Louisiana, June 25−28, 2012 p2734

    [32]

    Candler G V, MacCormack R W 1988 26th AIAA Aerospace Sciences Meeting Reno, Nevada, January 11−14, 1988 p511

  • [1] 罗仕超, 吴里银, 常雨. 高超声速湍流流动磁流体动力学控制机理. 物理学报, 2022, 71(21): 214702. doi: 10.7498/aps.71.20220941
    [2] 牛越, 包为民, 李小平, 刘彦明, 刘东林. 大功率热平衡感应耦合等离子体数值模拟及实验研究. 物理学报, 2021, 70(9): 095204. doi: 10.7498/aps.70.20201610
    [3] 丁明松, 江涛, 刘庆宗, 董维中, 高铁锁, 傅杨奥骁. 基于电流积分计算磁矢量势修正的低磁雷诺数方法. 物理学报, 2020, 69(13): 134702. doi: 10.7498/aps.69.20200091
    [4] 丁明松, 傅杨奥骁, 高铁锁, 董维中, 江涛, 刘庆宗. 高超声速磁流体力学控制霍尔效应影响. 物理学报, 2020, 69(21): 214703. doi: 10.7498/aps.69.20200630
    [5] 喻明浩. 非平衡感应耦合等离子体流场与电磁场作用机理的数值模拟. 物理学报, 2019, 68(18): 185202. doi: 10.7498/aps.68.20190865
    [6] 姜春华, 赵正予. 化学复合率对激发赤道等离子体泡影响的数值模拟. 物理学报, 2019, 68(19): 199401. doi: 10.7498/aps.68.20190173
    [7] 刘迎, 陈志华, 郑纯. 黏性各向异性磁流体Kelvin-Helmholtz不稳定性: 二维数值研究. 物理学报, 2019, 68(3): 035201. doi: 10.7498/aps.68.20181747
    [8] 成玉国, 程谋森, 王墨戈, 李小康. 磁场对螺旋波等离子体波和能量吸收影响的数值研究. 物理学报, 2014, 63(3): 035203. doi: 10.7498/aps.63.035203
    [9] 王蓬, 田修波, 汪志健, 巩春志, 杨士勤. 有限尺寸方靶等离子体离子注入动力学的三维粒子模拟研究. 物理学报, 2011, 60(8): 085206. doi: 10.7498/aps.60.085206
    [10] 杜宏亮, 何立明, 兰宇丹, 王峰. 约化场强对氮-氧混合气放电等离子体演化特性的影响. 物理学报, 2011, 60(11): 115201. doi: 10.7498/aps.60.115201
    [11] 欧阳建明, 邵福球, 邹德滨. 大气等离子体中负氧离子产生和演化过程数值模拟. 物理学报, 2011, 60(11): 110209. doi: 10.7498/aps.60.110209
    [12] 兰宇丹, 何立明, 丁伟, 王峰. 不同初始温度下H2/O2混合物等离子体的演化. 物理学报, 2010, 59(4): 2617-2621. doi: 10.7498/aps.59.2617
    [13] 杨涓, 石峰, 杨铁链, 孟志强. 电子回旋共振离子推力器放电室等离子体数值模拟. 物理学报, 2010, 59(12): 8701-8706. doi: 10.7498/aps.59.8701
    [14] 吴 翊, 荣命哲, 杨 飞, 王小华, 马 强, 王伟宗. 引入6波段P-1辐射模型的三维空气电弧等离子体数值分析. 物理学报, 2008, 57(9): 5761-5767. doi: 10.7498/aps.57.5761
    [15] 欧阳建明, 邵福球, 林明东. 含氧等离子体中臭氧形成过程数值模拟. 物理学报, 2008, 57(5): 3293-3297. doi: 10.7498/aps.57.3293
    [16] 赵国伟, 王之江, 徐跃民, 粱志伟, 徐 杰. 射频激励等离子体非线性效应的FDTD数值模拟. 物理学报, 2007, 56(9): 5304-5308. doi: 10.7498/aps.56.5304
    [17] 张 霆, 丁伯江. 原子过程对极向CXRS测量影响的数值模拟. 物理学报, 2006, 55(3): 1534-1538. doi: 10.7498/aps.55.1534
    [18] 欧阳建明, 邵福球, 王 龙, 房同珍, 刘建全. 一维大气等离子体化学过程数值模拟. 物理学报, 2006, 55(9): 4974-4979. doi: 10.7498/aps.55.4974
    [19] 黄勤超, 罗家融, 王华忠, 李 翀. EAST装置等离子体放电位形快速识别研究. 物理学报, 2006, 55(1): 281-286. doi: 10.7498/aps.55.281
    [20] 王艳辉, 王德真. 介质阻挡均匀大气压辉光放电数值模拟研究. 物理学报, 2003, 52(7): 1694-1700. doi: 10.7498/aps.52.1694
计量
  • 文章访问数:  5609
  • PDF下载量:  75
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-03-16
  • 修回日期:  2019-06-30
  • 上网日期:  2019-09-01
  • 刊出日期:  2019-09-05

热化学模型对高超声速磁流体控制数值模拟影响分析

  • 中国空气动力研究与发展中心计算空气动力研究所, 绵阳 621000
  • 通信作者: 董维中, dms2008@qq.com

摘要: 针对等离子体流场的模拟准确性问题及其对高超声速磁流体控制的影响, 通过数值求解三维非平衡Navier-Stokes流场控制方程和Maxwell电磁场控制方程, 建立了三维低磁雷诺数磁流体数值模拟方法及程序, 分析了不同空气组分化学反应模型和壁面有限催化效率等因素对高超声速磁流体控制的影响. 研究表明: 不同空气组分化学反应模型对高超声速磁流体流场结构、气动力/热特性控制的影响不容忽视; 对于本文计算条件, Park化学反应模型在组分模型一致性、等离子体模拟准确性等方面具有一定优势; 磁控热防护效果, 受壁面有限催化复合系数影响较大, 两者呈非线性关系, 不同表面区域差异较大; 磁场对磁阻力伞及其磁阻力特性影响, 受壁面催化效应的影响相对较小.

English Abstract

    • 磁流体力学(magnetohydrodynamic, MHD)是一门研究磁场作用下流体运动规律的科学, 在天体物理、空间科学、地球科学及太阳物理学等研究领域均有非常重要的研究价值. 目前, 对磁流体力学的研究又转向了高超声速飞行领域, 即高超声速磁流体控制. 在飞行器高超声速飞行过程中, 强烈的气动加热使激波层内温度急剧上升, 流场中发生剧烈的物理化学变化, 混合气体发生离解、复合、电离和交换等化学反应, 气体分子内能模态不同程度地激发, 形成具有弱导电性等离子体流场[1]. 利用机载磁场发生装置向等离子体流场注入适当的动量和能量对流场进行控制, 可以有效地提高和改进飞行器气动特性, 在高超声速飞行器气动力控制、气动热防护和等离子体分布调整等方面具有广阔的应用前景[2,3].

      数值模拟是研究高超声速磁流体控制的主要方法之一. 要准确模拟高超声速磁流体流动, 首先依赖于等离子体流场参数的准确获得, 包括气体电离程度、组分、密度和热力学温度状态(主要是平转动温度和振动-电子温度)等参数的精确模拟. 由于等离子体是在高超声速流动中产生的, 这就存在化学反应与热力学松弛的非平衡效应[4]. 在磁流体控制数值模拟时, 必须耦合考虑高温流场中发生复杂的物理化学现象, 包括气体热化学反应、分子能量模态激活及其非平衡效应, 即高温气体非平衡效应[4], 得到详细的等离子体流场参数, 分析等离子体与电磁场耦合作用.

      在高超声速磁流体数值模拟方面, 国外有不少考虑高温气体非平衡效应的研究. 例如, 2010—2011年, Bisek等[5,6]考虑Park化学反应模型开展了非平衡磁流体磁控热防护研究; 2012年, Chernyshev等[7]将等离子体混合气体分为电子、离子和中性粒子三类来模拟磁流体电离反应非平衡流动; 同年, Andrea和Carlo [8]采用Dunn-Kang模型对钝柱体非平衡磁流体控制进行了研究; 2013年, Fujino和Ishikawa[9]采用Park化学反应模型模拟再入飞行器非平衡磁流体控制; 2015年, Masuda等[10]采用多组分化学反应模型对三维钝锥体磁流体气动力、热控制进行了研究.

      近十多年来, 国内逐渐开展了高超声速飞行器电磁流体控制方面的研究, 得到了很多有价值的研究成果. 例如, 2005年, 赫新等[11]数值模拟了理想磁流体一维MHD激波管流动和二维喷管MHD流动; 2007年, 潘勇[3]采用弱耦合方法开展了二维钝锥化学非平衡无黏磁流体数值模拟研究; 2008年, 田正雨[12]采用平衡气体模型对球头磁流体热流控制进行了数值计算分析; 陈刚等[13]对二维理想磁流体进行了数值模拟; 2009年, 黄富来和黄护林[14]采用7组分化学反应模型对钝锥高超声速磁流体控制进行了数值计算研究; 2013年, 黄浩等[15]采用完全气体模型对电子束电离的高超声速磁流体发电机进行了研究; 何淼生等[16]采用完全气体模型研究了超高超声速进气道黏性效应的磁流体控制; 2014年, 卜少科等[17]采用完全气体模型开展了高超声速磁流体数值模拟研究; 2016年, 李开、柳军等采用完全气体模型对高超声速飞行器常规螺线管磁控热防护系统进行了研究[18], 2017年, 又采用化学反应Gupta模型开展了磁控热防护霍尔电场数值模拟分析[19]. 2018年, 姚霄等[20]采完全气体模型对五种外磁场作用的磁流体增阻特性进行了分析.

      可以看出, 在高超声速飞行器外场的磁流动控制研究方面, 由于气体介质的电导率一般较低, 国外很重要的一个发展趋势是: 忽略感应磁场影响, 在低磁雷诺数假设下, 研究不同磁场配置条件对各种外形高超声速飞行器流动特性的影响规律[510]. 尽管国外高超声速非平衡磁流体控制研究得到了很大的发展, 但大多选定某个单一的化学反应模型和热力学温度模型进行研究, 尚未看到详细分析不同空气化学反应模型和壁面有限催化特性等因素对高超声速MHD控制数值模拟影响的研究. 而国内也有很多基于低磁雷诺数方法开展高超声速磁流体控制方面的研究[1320], 大多采用直接采用简化模型(如完全气体、平衡气体或无黏流体等), 也有少量考虑非平衡效应的高超声速磁流体控制研究, 但从整体上看, 仍与国外水平存在一定差距.

      作者所在研究团队对高超声速飞行器非平衡等离子体流场及磁流体控制进行了较为广泛的研究[4,2124]. 本文在此基础上, 主要针对常见等离子体空气热化学气体模型, 开展热化学非平衡气体高超声速等离子体电磁流动数值研究, 较为系统地分析了不同空气组分化学反应模型以及壁面催化效率等对高超声速磁流体控制的影响.

    • 对于高超声速飞行器绕流流场, 气体介质一般都符合低电导率的特征, 通常满足低磁雷诺数($R{e_{\rm{m}}} < < 1$)假设, 感应磁场相对于外磁场很小, 基本可以忽略, 控制方程右端出现电磁源项, 其无量纲形式为

      $\begin{split} &\frac{{\partial {{Q}}}}{{\partial t}} + \frac{{\partial {{F}}}}{{\partial x}} + \frac{{\partial {{G}}}}{{\partial y}} + \frac{{\partial {{H}}}}{{\partial z}} \\ =& \frac{1}{{Re}}\left( {\frac{{\partial {{{F}}_{\rm{V}}}}}{{\partial x}} + \frac{{\partial {{{G}}_{\rm{V}}}}}{{\partial y}} + \frac{{\partial {{{H}}_{\rm{V}}}}}{{\partial z}}} \right) \\ & +{{W}} + {{{W}}_{{\rm{MHD}}}}, \end{split} $

      守恒变量为

      ${{Q}} = {({\rho _j},\rho ,\rho u,\rho v,\rho w,\rho {E_t})^{\rm{T}}}$,

      式中${\rho _j}$为气体组分j密度; $u,v,w$为直角坐标系3个方向速度; ${E_t}$为内能; $Re$为雷诺数; ${{F}},{{G}},{{H}}$${{{F}}_{\rm{V}}},{{{G}}_{\rm{V}}},{{{H}}_{\rm{V}}}$分别为直角坐标3个方向的无黏通量与黏性通量; ${{W}}$${{{W}}_{{\rm{MHD}}}}$为非平衡源项和电磁作用源项. 电磁作用源项${{{W}}_{{\rm{MHD}}}}$

      $\begin{split} {{{W}}_{{\rm{MHD}}}} = & {Q_{\rm{m}}} \cdot [{0_j},0,{({{J}} \times {\rm B})_x},{({{J}} \times {{B}})_y},\\ & {({{J}} \times {{B}})_z},{{J}} \cdot {{E}}]^{\rm{T}}, \end{split}$

      式中${{J}}$为电流密度, ${{B}}$为磁感应强度, ${{E}}$为电场强度, ${Q_{\rm{m}}}$为磁相互作用数.

    • 为了模拟高温流场等离子体分布, 混合空气组分模型选用7组分(O2, N2, NO, O, N, NO+, e)和11组分(O2, N2, NO, O, N, NO+, e, O2+, N2+, O+, N+)的电离气体模型, 化学反应模型采用较为常见的Dunn-Kang模型[25]、Park模型[26]和Gupta模型[27]. 非平衡源项中化学反应j组分生成源项${w_j}$可写为

      $ {w_j} = {M_j}\sum\limits_{i = 1}^{{N_r}} {\left( {\gamma _{ij}^ * - {\gamma _{ij}}} \right) \,} {Q_i}, $

      式中${N_r}$为模型中化学反应个数, ${M_j}$${Q_i}$分别为第j组分的分子量和第i化学反应的生成源项, $\gamma _{ij}^ * $$\gamma _{ij}^{}$分别为第i反应中第j组分的生成物系数和反应物系数.

    • 壁面催化效应会影响表面附近等离子体电离复合, 进而影响等离子体电导率分布. 本文主要考虑完全非催化(non-catalytic wall, NCW)、完全催化(fully catalytic wall, FCW)和有限催化(partially catalytic wall, PCW). 其中NCW和FCW相对简单, 详见文献[4]. 对于PCW, 主要考虑氮原子和氧原子复合反应以及电离组分的7个复合反应[28] (见表1).

      r化学反应式r化学反应式
      1O + O$ \Rightarrow $O25N+ + e$ \Rightarrow $N
      2N + N$ \Rightarrow $N26O2+ + e$ \Rightarrow $O2
      3NO+ + e$ \Rightarrow $NO7N2+ + e$ \Rightarrow $N2
      4O+ + e$ \Rightarrow $O

      表 1  表面复合反应

      Table 1.  Surface combinative reaction.

      r反应速率常数${k_{{\rm{w,}}r}}$

      $ {k_{{\rm{w,}}r}} = \alpha \sqrt {\frac{{{R_0}{T_{\rm{w}}}}}{{2{\text{π }}{M_r}}}}, $

      式中$\alpha $是壁面催化复合系数, ${R_0}$为气体常数, ${T_{\rm{w}}}$为壁面温度, ${M_r}$为第r反应中非电子反应物的分子量.

    • 采用基于电子和离子在等离子体中碰撞和迁移的多种电离成分混合气体电导率模型, 该模型在高超声速高温非平衡等离子体磁流体控制中应用很广泛[7,9,10,19], 形式为

      $ \sigma = \frac{{{n_{\rm{e}}}{e^2}}}{{{m_{\rm{e}}}\displaystyle\sum\limits_{s \ne {\rm{e}}} {\upsilon _{{\rm{e}},s}^{\rm{m}}} }} ,$

      $\upsilon _{{\rm{e}},s}^{\rm{m}}$为电子与第s组分的有效动量传输碰撞频率. 当s为离子组分, 其表达式可写为

      $\begin{split} \upsilon _{{\rm{e}},s}^{\rm{m}} = &6{\text{π}}{\left( {\frac{{{e^2}}}{{12{\text{π}}{\varepsilon _0}{k_{\rm{B}}}{T_{}}}}} \right)^2} \\ & \times \ln \left[12{\text{π}}{\left( {\frac{{{\varepsilon _0}{k_{\rm{B}}}}}{{{e^2}}}} \right)^{1.5}}\sqrt {\frac{{{T_{}}^3}}{{{n_{\rm{e}}}}}} \right]{n_s}\sqrt {\frac{{8{k_{\rm{B}}}{T_{}}}}{{{\text{π}}{m_{\rm{e}}}}}} , \end{split}$

      s为中性组分, 其表达式为

      $ \upsilon _{{\rm{e}},s}^{\rm{m}} = \frac{4}{3}\sigma _{{\rm{e}},{\rm s}}^{\rm{m}}{n_s}\sqrt {\frac{{8{k_{\rm{B}}}{T_{}}}}{{{\text{π}}{m_{\rm{e}}}}}} , $

      这里${k_{\rm{B}}}$, ${\varepsilon _0}$分别为玻尔兹曼常数和真空介电常数; e, ${m_{\rm{e}}}$, ${n_{\rm{e}}}$, ${n_s}$分别为电子电量、电子质量、电子数密度、组分s数密度; $\sigma _{{\rm{e}},{\rm s}}^{\rm{m}}$为s电子与中性组分有效碰撞截面, 它是一般温度函数, 由试验测量拟合得到.

    • 假设等离子体微元呈电中性, 则电流密度满足连续性方程

      $ \nabla \cdot {{J}} = 0. $

      电流${{J}}$可由广义欧姆定律得到,

      $ {{J}} = \sigma ({{E}} + {{V}} \times {{B}}), $

      这里${{V}}$为混合气流速度矢量. 对于定常外磁场, 由法拉第定律$\nabla \times {{E}} = - \partial {{B}}/\partial t$ = 0, 可得电场是无旋的, 可写成电势函数形式: ${{E}} = - \nabla \phi $. 代入(8)式、(9)式可得

      $\nabla \cdot (\sigma \nabla \phi ) = \nabla \cdot [\sigma ({{V}} \times {{B}})].$

      引入虚拟时间项得到

      $ \partial \phi /\partial t + \nabla \cdot (\sigma \nabla \phi ) = \nabla \cdot [\sigma ({{V}} \times {{B}})]. $

      采用隐式数值离散(10)式, 并与流动控制方程(1)耦合求解, 得到电势函数$\phi $和感应电场${{E}}$, 再由(9)式得到感应电流${{J}}$.

    • 一般情况下, 气动力系数计算只需考虑飞行器表面的压力和黏性应力积分, 而高超声速磁流体控制气动力系数计算时, 还需考虑洛伦兹力反作用力的空间积分, 以二维磁流体阻力系数${C_{\rm{D}}}$计算为例,

      $ \begin{split} &{F_2} = - {\iiint {({{J}} \times {{B}})}_x}{\rm{d}}V\cos \beta \\ & \qquad\; - {\iiint {({{J}} \times {{B}})}_z}{\rm{d}}V\sin \beta , \\ & {{F_{\rm{D}}} = {F_1} + {F_2}}\qquad{C_{\rm{D}}} = {F_{\rm{D}}}/(0.5{q_\infty }{S_{{\rm{ref}}}}), \end{split} $

      其中${F_{\rm{D}}}$为飞行器受到的阻力, ${F_1}$为表面压力和黏性应力造成阻力, ${F_2}$为洛伦兹力反作用力造成的阻力, $\beta $为飞行攻角, ${q_\infty }$${S_{{\rm{ref}}}}$分别是动压和参考面积.

    • 在以往的工作中, 作者所在团队对数值计算方法, 从热化学非平衡等离子体流动[4,2123]、表面有限催化效应[29]、气动热环境模拟[30]以及等离子体磁流体控制[24]等多个方面进行了考察验证. 本文主要是针对高超声速磁流体控制气动力特性方面的校验.

      采用球柱外形[31], 球头半径1 m, 柱体长10 m. 计算高度为H = 69 km, 来流速度为6500 m/s. 外加磁偶极子磁场, 磁场分布方式如下:

      $ {{B}} = {B_0}\left( {\frac{{\cos (\alpha )}}{{{{(r/{r_0})}^3}}}{{r}} + \frac{{\sin (\alpha )}}{{2{{(r/{r_0})}^3}}}{{\alpha }}} \right). $

      这里(${{r}},{{\alpha }}$)为极坐标单位矢量, 偶极子中心位于坐标原点, 即头部球心; ${B_0}$为极轴上距离偶极子中心${r_0}$处磁感应强度, ${r_0} = 1.0\;{\rm{m}}$; 这里, 磁场极轴与球柱轴线的夹角$\theta $, 偶极子参考点${B_0} =0$$0.2\;{\rm{T}} $, 如图1所示.

      图  1  磁场配置(${B_0} = 0.2\;{\rm{T}}$)

      Figure 1.  Magnetic field (${B_0} = 0.2\;{\rm{T}}$).

      对于本文计算的外形, 磁场作用下, 洛伦兹力使激波外推, 形成磁阻力伞, 见图2(a). 图2(b)给出了球柱阻力系数随磁感应强度的变化分布, 其中标号${C_{\rm{D}}}$为阻力系数, ${C_{{\rm{D1}}}}$为阻力系数中压力和黏性应力分量, ${C_{{\rm{D}}2}}$为阻力系数中洛伦兹力反作用力分量. 表2给出了本文计算得到的阻力系数和文献[31]结果(由曲线取值得到)对比. 可以看出, 磁感应强度增大, 洛伦兹力反作用力带来的阻力系数分量增大, 使总的阻力系数增大; 本文计算的阻力系数, 与文献[31]几乎完全一致(图2(b)表2), 差异小于0.5%.

      图  2  球柱阻力特性 (a) 磁阻力伞; (b) 阻力系数

      Figure 2.  Drag characteristics of ball-column model: (a) Magneto-resistance parachute; (b) drag coefficient.

      夹角$\theta $文献结果[31]本文结果差异
      1.82021.82870.47%
      45º1.66701.67060.22%
      90º1.47551.48040.33%

      表 2  球柱阻力系数(${B_0} = 0.2\;{\rm{T}}$)

      Table 2.  Drag coefficient of ball-column model(${B_0} = 0.2\;{\rm{T}}$).

    • 采用RAM-C钝锥外形, 该外形有电子数密度飞行试验结果[32]. 球头半径Rn = 0.1524 m, 全长1.295 m, 半锥角为9°. 计算条件为高度71 km, 速度7650 m/s, 壁面温度Tw = 1500 K, 绝缘壁面或近似绝缘壁面条件. 外加磁偶极子磁场((12)式), 模拟实际的螺线管磁场. 这里${B_0} = 0.5\;{\rm{T}}$, ${r_0} =$$ 0.1524\;{\rm{m}} $. 磁偶极子方向指向直角坐标X轴负方向.

      为了保证计算的可靠性, 图3使用两套不同密度网格计算分析. Grid 1为稀网格结果, 网格壁面第一层距离为0.005 mm; Grid 2为密网格结果, 第一层距离为0.001 mm. 由图可以看出, 两套网格结果几乎完全重合, 计算结果受网格影响较小. 本文计算中如无特殊说明, 均采用Grid 2进行计算.

      图  3  不同网格条件下的参数 (a) 电导率; (b) 电子数密度; (c) 热流

      Figure 3.  Parameter of different grids: (a) Conductivity; (b) ele-ctronic number density; (c) heat flux.

    • 为了分析不同组分空气化学反应模型对MHD控制影响, 采用7组分、11组分的空气组分模型, 化学反应Dunn-Kang, Park, Gupta模型, 开展数值模拟分析.

      首先分析无磁场条件下(${B_0} = 0\;{\rm{T}}$)化学反应模型对高温等离子体流场数值模拟的影响. 图4给出了不同空气化学反应模型数值计算得到的驻点线电导率、电子数密度、温度分布等参数分布, R是离开壁面的法向距离, 7Ns为采用7组分模型的标识, 11Ns为采用11组分模型的标识. 由图4(a)可以看出: 化学反应模型方面, 波后很大区域内11组分的Dunn-Kang和Gupata模型计算得到的电导率远大于Park模型, 部分区域差别可达5倍以上; 组分模型方面, 7组分模型作为11组分模型的简化模型, Park模型的组分模型一致性较好(即7组分与11组分差异较小), 而Dunn-Kang和Gupata模型在这方面表现较差.

      图  4  驻点线参数分布(${B_0} = 0\;{\rm{T}}$) (a) 电导率; (b) 电子数密度; (c) 温度

      Figure 4.  Parameter along stagnation line(${B_0} = 0\;{\rm{T}}$): (a) Conductivity; (b) electronic number density; (c) temperature

      电导率受模型影响的原因, 可结合(5)—(7)式、图4(b)图4(c)进一步分析. 由图可以看出, 不同模型得到的温度分布存在差异, 但幅度较小, 其幅度变化规律与电导率存在较大不同, 这说明流场结构及温度分布差异不是电导率产生差异的主要原因; 而采用11组分Dunn-Kang和Gupata模型计算得到的电子数密度远高于Park模型, 局部可达2个数量级, 这是电导率产生较大差异的主要原因之一; Park模型电子数密度计算的组分模型一致性较好, 而Dunn-Kang和Gupta化学模型计算的电子数密度, 7组分模型结果远小于11组分模型, 这说明对于采用这两种化学反应模型, 离子组分O2+, N2+, O+, N+的相关反应较强, 对流场的电子数密度(可看作电离度的表征)影响较大, 进而影响混合气体电导率; 值得指出的是, 采用7组分模型, 化学反应Park, Dunn-Kang, Gupata模型三者流场特性(激波位置、温度、电子数密度)相差较小, 这说明这三种模型在O2/N2离解、NO电离等方面模拟具有一致性.

      在电磁流动控制中, 常用磁相互作用数${Q_{\rm{m}}}$来表征外磁场与流动耦合作用的整体效果、效率, 其表达式为${Q_{\rm{m}}} = {\sigma _{\rm{0}}}B_0^2L/({\rho _0}{u_0})$, 这里${\sigma _{\rm{0}}}$, $B_0^{}$,L, ${\rho _0}$, ${u_0}$分别为电导率、磁感应强度、磁场作用区域尺度、流体密度、流体速度的特征量. 由此可以看出电导率的差异必然影响磁流体控制效果. 但对于高超声速流场来说, 仅用${Q_{\rm{m}}}$作为控制效果的定量描述是不充分的, 这一方面是由于气体电导率受流动结构影响, 分布差异较大, 很难提取单一的电导率来综合表征磁场整体效果, 同时, 磁场作用下流动结构会发生变化, 即电导率分布与磁场作用存在耦合作用, 进一步增大了提取电导率有效表征的难度; 另一方面, 高超声速电磁流动控制效果存在“饱和”现象[24], 具体表现为气动力、热特性等随电导率、磁场强度等参量的非线性变化, 此时仅用${Q_{\rm{m}}}$不足以表征这种非线性规律. 因此, 要分析高超声速MHD控制的效果(尤其是气动热特性), 数值模拟时还需结合热化学模型进行具体分析.

      图5给出了${B_0} = 0.5\;{\rm{T}}$时不同空气化学模型数值计算得到的驻点线温度和电子数密度. 图6给出了采用不同空气化学模型模拟时磁场对表面热流的影响. 表3给出了不同空气化学模型磁场对阻力系数的影响. 可以看出, 本文计算条件下, 无磁场条件的情况下, 各种热化学模型计算得到的流场结构、表面热流、阻力系数基本一致(11组分Gupata模型除外). 磁场作用效果表现为: 激波脱体距离增大, 波后高温区和高电子数密度区增大, 峰值温度略微下降; 表面绝大部分区域热流下降, 存在局部热流反冲(X = 0—0.2 m之间, 其反冲产生的原因见文献[24]); 由于洛伦兹力形成的磁阻力伞作用(与图2类似), 总阻力系数增大. 采用不同空气化学模型, 磁场作用效果差异与图4电导率差异, 在定性上保持一致. 但在影响幅度定量方面存在差别: 流场结构方面, 采用11组分Dunn-Kang和Gupata模型磁场作用使激波脱体距离增大至5倍左右, 而采用其他模型, 磁场作用仅使激波脱体距离增大至2—3倍. 气动热环境方面, 热化学模型对磁控热流减缓幅度的影响, 不同的表面区域存在较大差异, 其变化规律较为复杂, 很难用单个数据衡量, 其中受空气化学模型影响较大区域为锥身结合处附近(X = 0 m)和锥身末端; 气动力特性方面, 11组分Dunn-Kang和Gupata模型的磁场增阻效果分别为179%和207%, 其他空气化学模型的磁场增阻效果约60%—70%.

      图  5  驻点线温度和电子数密度分布(${B_0} = 0.5\;{\rm{ T}}$) (a) 温度; (b) 电子数密度

      Figure 5.  Temperature and electronic number density along stagnation line (${B_0} = 0.5\;{\rm{T}}$): (a) Temperature; (b) electronic number density.

      图  6  表面热流 (a) Dunn-Kang 模型; (b) Park模型; (c) Gupta模型

      Figure 6.  Heat flux on wall: (a) Dunn-Kang model; (b) Park model; (c) Gupta model.

      模型${B_0} = 0\;{\rm{T}}$${B_0} = 0.5\;{\rm{T}}$增大比例/%
      7Ns Dunn-Kang0.2974720.47808161
      11Ns Dunn-Kang0.2911290.812020179
      7Ns Gupta0.3010720.51869272
      11Ns Gupta0.2587720.795648207
      7Ns Park0.2900450.48288567
      11Ns Park0.2899860.48915469

      表 3  不同空气化学模型阻力系数

      Table 3.  Drag coefficient under different air chemical models.

      由此可见, 相比于无磁场作用的高温气体流动, 空气化学模型的对高超声速MHD控制数值模拟的影响更加显著.

      为了进一步分析空气化学模型在等离子体数值模拟方面的准确性, 图7给出了不同空气化学模型计算的电子数密度与飞行试验[32]结果, 这里Rn为球头半径, YX/Rn = 8.1处远离壁面的距离(沿壁面法向方向). 可以看出, 对于7组分空气模型, 三种化学反应模型结果均与试验结果[32]符合较好, Park模型与Gupta模型表现略优; 对于11组分空气模型, Dunn-Kang模型与Gupta模型结果与试验差别较大, 而Park化学反应模型仍表现优异. 这说明, 在本文计算条件下, 相对于NO+的电离反应, 离子组分O2+, N2+, O+, N+的相关电离反应较为微弱, Park模型能较好地模拟这一现象.

      图  7  不同模型的电子数密度 (a) 7组分; (b) 11组分

      Figure 7.  Electron number density: (a) 7 speices; (b) 11 speices.

      随电导率升高, 磁扩散特征时间变长, 低磁雷诺数假设可能失效, 而本文采用11组分Dunn-Kang和Gupata化学反应模型模拟时流场中局部电导率峰值达到2000 S/m. 因此本文开展了这种条件下低磁雷诺数假设适用性分析. 图8对比给出了采用11组分Dunn-Kang模型流场结构云图, Case 1为采用低磁雷诺数假设的数值结果, case 2为不采用低磁雷诺假设(即考虑感应磁场影响)的数值模拟结果. 由图可以看出, 两种方法计算得到的流场结构, 除局部细节外, 整体上一致; 图9进一步给出了该条件下表面压力和热流的分布, 两者分布对应了飞行器的气动力、热特性, 可以看出, 这两者差异较小, 没有本质性区别. 美国斯坦福大学MacCormack针对球柱模型开展了低磁雷诺数方法适用性分析[2], 其外加磁场为偶极子磁场、来流速率为6710 m/s、气体电导率为2800—5500 S/m, 其主要结论与本文类似. 由此可见, 对于本文的计算条件来说, 基于低磁雷诺数假设, 忽略感应磁场, 对计算结果和结论的分析, 不会造成本质性影响.

      图  8  流场结构云图 (a) Case 1; (b) case 2

      Figure 8.  Flow structure map: (a) Case 1; (b) case 2.

      图  9  表面压力和热流 (a) 压力; (b) 热流

      Figure 9.  Surface pressure and heat flux: (a) Pressure; (b) heat flux.

    • 图7还可以看出, 壁面催化效应对等离子体分布存在一定影响. 为了分析壁面催化效应对等离子体及MHD控制的影响, 采用不同壁面催化系数($\alpha = $0.001—1.000), 开展磁流体控制计算对比分析.

      图10首先给出了无磁场条件下(${B_0} = 0\;{\rm{T}}$)采用不同催化复合系数计算得到的电子数密度和电导率. 可以看出, 催化复合系数在一定程度上影响了流场中电子数密度分布, 尤其是表面附近, 随催化复合系数增大($\alpha =$0.001—1.000), 电子数密度下降了约3倍; 气体电导率随催化复合系数增大呈下降趋势, 但整体上变化幅度不大(壁面附近区域除外).

      图  10  不同壁面催化系数条件下驻点线参数分布 (${B_0} = 0\;{\rm{T}}$) (a) 电子数密度; (b) 电导率

      Figure 10.  Parameter along stagnation line of different catalytic efficiency (${B_0} = 0\;{\rm{T}}$): (a) Electronic number density; (b) conductivity

      图11对比了不同催化复合系数条件下阻力系数受磁场作用的影响, 阻力系数增大百分比${D_1}$, 见(13)式, 其中${C_{{\rm{DM}}}}$为有磁场作用时的阻力系数, ${C_{{\rm{D}}0}}$为无磁场作用时的阻力系数. 可以看出, 表面催化效应对阻力系数计算有一定影响, 影响幅度较小; 随催化系数增大, 磁场增阻效果呈下降趋势, 但下降幅度较小, 这与电导率整体变化规律基本一致.

      图  11  不同催化系数条件下阻力系数及磁场作用效果 (a) 阻力系数; (b) 阻力系数增大百分比

      Figure 11.  Effect of different catalytic efficiency conditions: (a) Drag coefficient; (b) drag coefficient increment ratio.

      图12 为不同催化复合系数条件下热流受磁场作用的影响比较. 由图12(a)可以看出, 催化效应对表面热流影响明显(无论有无磁场作用), 催化系数增大, 热流升高, 这符合催化效应影响的一般规律; 磁场作用效果受表面催化复合系数影响显著, 催化系数较大($\alpha = 1.0$)时, 磁场使热流下降的幅度远小于催化系数较小时的结果, 这与“催化系数对电导率整体影响较小”的变化规律存在差别.

      图  12  不同催化复合系数条件下热流受磁场作用的影响 (a) 表面热流; (b) 驻点热流; (c) 磁场效果

      Figure 12.  Effect of different catalytic efficiency conditions: (a) Surface heat flux; (b) heat flux at stagnation point; (c) magnetic effect.

      为了分析这一现象, 图12(b)图12(c)给出了不同催化系数条件下驻点热流及其热流分量的磁控效率. 这里${D_2}$为磁场使总热流下降的百分率, 其中${Q_{{\rm{AM}}}}$为有磁场作用时的驻点总热流, ${Q_{\rm{A}}}$为无磁场作用时的驻点总热流; ${D_3}$为磁场使驻点温度传导热流下降的百分率, 其中${Q_{{\rm{1M}}}}$为有磁场作用时的温度热传导热流, ${Q_1}$为无磁场作用时的温度热传导热流; ${D_4}$为磁场使组分扩散热流下降的百分率, 其中${Q_{{\rm{2M}}}}$为有磁场作用的组分扩散热流, ${Q_2}$为无磁场作用的组分扩散热流; ${D_5}$为无磁场时组分扩散热流占总热流的百分率. 相关公式如下:

      $ \begin{split} &{D_1} = ({{{C_{{\rm{DM}}}} - {C_{{\rm{D}}0}}}})/{{{C_{{\rm{D}}0}}}} \times 100\%,\\ &{D_2} = ({{{Q_{\rm{A}}} - {Q_{{\rm{AM}}}}}})/{{{Q_A}}} \times 100\%, \\ &{D_3} = ({{{Q_1} - {Q_{1{\rm{M}}}}}})/{{{Q_1}}} \times 100\%, \\ &{D_4} = ({{{Q_2} - {Q_{2{\rm{M}}}}}})/{{{Q_2}}} \times 100\%,\\ &{D_5} = {{{Q_2}}}/{{{Q_{\rm{A}}}}} \times 100\%,\\ &{Q_{\rm{A}}} = {Q_1} + {Q_2}, \qquad {Q_{{\rm{AM}}}} = {Q_{1{\rm{M}}}} + {Q_{2{\rm{M}}}}. \end{split} $

      图12(b)图12(c)可以看出, 随催化复合系数增大, 磁场使驻点总热流下降幅度(${D_2}$)呈先剧烈下降, 然后缓慢上升; 磁场使温度传导热流的下降幅度(${D_3}$)呈下降趋势; 磁场使组分扩散热流的下降幅度(${D_4}$)呈上升趋势; 由于总热流中组分扩散热流占比不断增大(${D_5}$), 进而导致磁场使总热流下降幅度(${D_2}$)先降后升.

      值得指出的是, ${D_3}$的变化规律与图10 “随催化系数升高, 壁面附近电导率大幅下降”的规律一致, 符合磁相互作用数的理论分析, 而${D_4}$的变化规律则与其相反(在催化系数较低, 如$\alpha = 0.001$时, 磁场作用甚至使组分扩散热流增大). 为了分析这一现象的原因, 图13给出了$\alpha = 0.001$和0.75时驻点线氧原子和氮原子质量分数, 图13(c)图13(b)的壁面附近区域放大图.

      图  13  不同催化复合系数条件下驻点线原子组分分布 (a) $\alpha = 0.001$; (b) $\alpha = 0.75$; (c) $\alpha = 0.75$

      Figure 13.  Atomic components along stagnation line: (a) $\alpha = 0.001$; (b) $\alpha = 0.75$; (c) $\alpha = 0.75$.

      图4图5可知, 磁场作用使激波脱体距离增大, 高温区增大, 这意味着流场中化学高温离解反应有效时间更长. 结合图13可以看出, 磁场作用下, 氮气离解程度更大, 流场中氮原子的质量分数更高; 而氧气由于离解温度较低[4], 部分区域在无磁场作用时就离解得较为完全, 因此氧原子最大值受磁场影响较小. 由图13(a)可以看出, 在催化复合系数较低($\alpha = 0.001$)时, 由于壁面催化复合反应速率相对较慢, 壁面处氧原子、氮原子来不及复合, 有较大程度的堆积; 磁场作用下, 氮气离解度更高, 氮原子堆积的浓度也更高, 因此, 其组分扩散热流较无磁场时更大, 综合表现为磁场作用使质量扩散热流上升. 由图13(b)可见, 催化复合系数较大($\alpha = 0.75$), 壁面催化复合反应速率高, 氮原子和氧原子在壁面堆积较少, 有无磁场作用浓度差别不大; 此时在磁场作用下, 激波脱体距离大, 氮原子和氧原子扩散梯度较小, 因此组分扩散热流较无磁场时小, 符合“磁场使驻点热流下降”的一般规律. 由此可见, 随催化系数增大, 壁面附近氮原子、氧原子堆积量减少, “磁场作用使激波脱体距离增大影响氧原子和氮原子的扩散梯度”逐渐成为主导因素, 因此质量扩散热流的下降幅度${D_4}$呈上趋势.

    • 在本文的计算条件下, 有以下结论.

      1)空气化学反应模型对高超声速MHD控制影响不可忽视. 不同化学反应模型计算的等离子体电导率相差可达数倍, 严重影响了高超声速磁流体力、热控制效果. 就本文计算状态而言, Park化学反应模型在组分模型(11组分和7组分)一致性、等离子体模拟准确性等方面具有一定优势.

      2)壁面催化效应对高超声速MHD控制存在一定影响. 磁场对气动力系数的作用效果, 受表面催化效应影响幅度相对较小, 随有限催化复合系数升高, 磁场控制效率呈缓慢下降趋势; 磁场对表面热流的作用效果受催化效应影响较为显著, 催化系数$\alpha = 1.0$时, 磁场使热流下降的幅度远小于催化系数$\alpha = 0.001$时的结果, 要定量分析磁控热防护效果, 须考虑壁面材料的有限催化特性; 催化复合系数与磁场控制效果呈非线性关系, 且不同区域差异较大, 随催化复合系数升高, 驻点热流的磁场控制效率先大幅度下降, 后平缓上升, 这是温度传导热流和组分扩散热流共同作用的结果.

参考文献 (32)

目录

    /

    返回文章
    返回