再入飞行中DSMC与Navier-Stokes两种模型的计算与分析

王保国; 黄伟光; 钱耕; 李翔; 李学东

再入飞行中DSMC与Navier-Stokes两种模型的计算与分析

1.
北京理工大学宇航学院,北京 100081
2.
中国科学院工程热物理研究所,北京 100190

基金项目: 国家自然科学基金(50376004); 高等学校博士学科点专项基金(20030007028)

计量
- 文章访问数: 1946
- HTML浏览量: 247
- PDF量: 26
- 被引次数: 0
出版历程
- 收稿日期: 2011-02-28
- 修回日期: 2011-03-16
- 刊出日期: 2011-05-28

Computation and analysis of DSMC and Navier-Stokes model in reentry flight

1.
School of Aerospace Engineering, Beijing Institute of Technology, Beijing 100081, China
2.
Institute of Engineering Thermophysics, Chinese Academy of Sciences, Beijing 100190, China

摘要

摘要: 首次提出了再入飞行过程中“小Knudsen数特征区(即域)”的概念,并用直接模拟蒙特卡罗(DSMC)算法与求解Navier-Stokes方程两种模型详细研究了在这个区域中的流动.数值计算表明:采用王保国教授课题组自己编制的DSMC源程序和求解多组分、考虑热力学非平衡以及化学反应非平衡、弱守恒型Navier-Stokes方程组的源程序均可以完成在该区域的流动计算,而且所得流场气动参数的分布趋势大体上基本一致,对应点处相应参数的值相差较小.计算中发现:采用DSMC程序计算时,当来流工况的Knudsen数越小、越接近Kn₁时则计算所需的时间越长,使用该程序可以得到域的Kn₁值,其值近似为0.0019;采用Navier-Stokes方程的程序时,当来流工况的Knudsen数越大、越接近Kn₂时则计算越不易收敛,使用该程序能够得到Kn₂值,其值近似为0.0125.计算中还发现:在流场计算过程中,分别采用DSMC与Navier-Stokes方程的源程序,其对应的来流工况Knudsen数取值分别越接近Kn₁与Kn₂时,均会出现收敛速度明显变慢的现象.为此,文中给出了加速流场计算收敛的3种措施.正是由于小Knudsen数特征区的存在,因此整个再入过程的所有飞行工况可以通过适当选用DSMC与Navier-Stokes方程的程序完成.显然,这个结论对再入飞行过程中的数值计算来讲是非常重要的.
- 再入飞行 /
- 小Knudsen数特征区 /
- Navier-Stokes方程 /
- 直接模拟蒙特卡罗 (DSMC) /
- 高超声速气动热力学
Abstract: The concept of “low Knudsen number characteristic region”,which is ,in the process of reentry,has been presented for the first time,and the fluid in the region has been studied using direct simulation Monte Carlo (DSMC)and Navier-Stokes model respectively.It has been shown by the numerical computation that both of the DSMC code and weak conservative Navier-Stokes code with multi-species,thermodynamic non-equilibrium and chemical non-equilibrium compiled by Professor Wang Baoguo and his group can complete the computation in this region.The distribution of flowfield shows agreement with each other,and values of corresponding parameter only have very limited difference.As we can see in the computation:when DSMC is being used,the computation time increases with the decrease of Knudsen number.Value of Kn₁ can be got from this code,which is near 0.001 9.When Navier-Stokes is being used,the computation is more difficult to get convergence with a higher Knudsen number.Value of Kn₂ can be got from this code,which is near 0.012 5.It has also been found that during the computation,for Knudsen number closed to Kn₁ and Kn₂ respectively,the use of DSMC and Navier-Stokes code will lead to a decrease of the convergence rate.There are three measures to accelerate the convergence which will be introduced in this paper.Because of the low Knudsen number characteristic region,all the flight cases of the reentry process could be completed by appropriate selection between DSMC and Navier-Stokes code.Obviously,this conclusion is very important for the numerical computation in reentry flight process.
- reentry flight /
- low Knudsen number characteristic region /
- Navier-Stokes /
- direct simulation Monte Carlo (DSMC) /
- hypersonic aerothermodynamics

HTML

参考文献(59)

[1]	卞荫贵,徐立功.气动热力学[M].合肥:中国科学技术大学出版社,1997.
[2]	Anderson J D,Jr.Hypersonic and high temperature gas dynamics[M].New York:Mc Graw-Hill,1989.
[3]	王保国,黄伟光.高超声速气动热力学[M].北京:科学出版社,2011.
[4]	Bird G A.Molecular gas dynamics and the direct simulation of gas flow[M].Oxford:Clarendon Press,1994.
[5]	沈青.稀薄气体动力学[M].北京:国防工业出版社,2003.
[6]	王保国,刘淑艳.稀薄气体动力学计算[M].北京:北京航空航天大学出版社,2011.
[7]	Cercignani C.The Boltzmann equation and its appdlcations[M].New York:Springer-Verlag,1988.
[8]	Bird G A.Application of the DSMC method to the full shuttle geometry .AIAA Paper 90-1692,1990.
[9]	Bertin J J.Hypersonic aerothermodynamics[M].Washington D C:American Institute of Aeronautics and Astronautics,1994.
[10]	王保国,李翔.多工况下高超声速飞行器再入时流场的计算[J].西安交通大学学报,2010,44(1):71-76. WANG Baoguo,LI Xiang.Prediction of hypersonic vehicle reentry flowfields in multiple cases[J].Journal of Xi'an Jiaotong University,2010,44(1):71-76.(in Chinese)
[11]	王保国,李翔,黄伟光.激波后高温高速流场中的传热特性研究[J].航空动力学报,2010,25(5):963-980. WANG Baoguo,LI Xiang,HUANG Weiguang.Analytical study of heat transfer of high-velocity flow field behind shock waves[J].Journal of Aerospace Power,2010,25(5):963-980.(in Chinese)
[12]	Druguet M,Candler G V,Nompelis I.Effect of numerics on Navier-Stokes computations of hypersonic double-cone flows[J].AIAA Journal,2005,43(3):616-623.
[13]	Gnoffo P A.Computational aerothermodynamics in aeroassist applications[J].Journal of Spacecraft and Rockets,2003,40(3):305-312.
[14]	WANG Baoguo.On general form of Navier-Stokes equations and implicit factored scheme[J].Applied Mathematics and Mechanics,1988,9(2):179-188.
[15]	王保国.叶栅流基本方程组特征分析及矢通量分裂[J].中国科学院研究生院学报,1987,4(2):54-65. WANG Baoguo.Characteristic analysis and flux vector splitting for the basic govering equations of cascade flow[J].Journal of the Graduate School of Chinese Academy of Sciences,1987,4(2):54-65.(in Chinese)
[16]	王保国,卞荫贵.求解三维欧拉流的隐-显式格式及改进的三维LU算法[J].计算物理,1992,9(4):423-425. WANG Baoguo,BIAN Yingui.An efficient implicit-explicit scheme and an improved LU scheme and for 3D Euler solver[J].Chinese Journal of Computation Physics,1992,9(4):423-425.(in Chinese)
[17]	王保国,卞荫贵.关于三维Navier-Stokes方程的黏性项计算[J].空气动力学报,1994,12(4):375-382. WANG Baoguo,BIAN Yingui.On numerical computation of viscous terms for 3-D Navier-Stokes equations[J].Acta Aerodynamica Sinica,1994,12(4):375-382.(in Chinese)
[18]	王保国.新的解跨音速Euler方程的隐式杂交方法[J].航空学报,1989,10(7):309-315. WANG Baoguo.A new implicit hybrid schemes for the Euler equation of transonic flow[J].Acta Aeronautica et Astronautica Sinica,1989,10(7):309-315.(in Chinese)
[19]	王保国,沈孟育.高速粘性内流的高分辨率高精度迎风型杂交格式[J].空气动力学报,1995,13(4):365-373. WANG Baoguo,SHEN Mengyu.High resolution and high order accurate upwind TVD algorithm of high speed viscous internal flows[J].Acta Aerodynamica Sinica,1995,13(4):365-373.(in Chinese)
[20]	王保国,刘秋生,卞荫贵.三维Euler方程的两种高效高分辨率算法及其在高速进气道中的应用[J].应用基础与工程科学学报,1994,2(4):360-370. WANG Baoguo,LIU Qiusheng BIAN Yingui.Two efficient and high-resolution hybrid schemes for solving 3-d Euler equations and its application to high speed inlet flows[J].Journal of Basic Science and Engineering,1994,2(4):360-370.(in Chinese)
[21]	王保国,郭延虎,沈孟育.恢复函数的三点迎风紧致格式构造方法及应用[J].计算物理,1997,14(4):666-668. WANG Baoguo,GUO Yanhu,SHEN Mengyu.Third order accurate reconstruction with upwind compact difference and its application[J].Chinese Journal of Computation Physics,1997,14(4):666-668.(in Chinese)
[22]	王保国,郭延虎,沈孟育.贴体曲线坐标系中恢复函数的一种构造方法及应用[J].空气动力学报,1999,17(1):117-122. WANG Baoguo,GUO Yanhu,SHEN Mengyu.Third order accurate reconstruction with upwind compact difference in curvilinear coordinate and its application[J].Acta Aerodynamica Sinica,1999,17(1):117-122.(in Chinese)
[23]	王保国,孙成海,李荣先.快速生成三维非结构网格的一种方法[J].清华大学学报,1999,39(8):91-95. WANG Baoguo,SUN Chenghai,LI Rongxian.An efficient method for generation of three dimensional unstructured grid[J].Journal of Tsinghua University(Science and Technology),1999,39(8):91-95.(in Chinese)
[24]	王保国,刘淑艳,闫为革.高精度强紧致三点格式的构造及边界条件的处理[J].北京理工大学学报,2003,23(1):13-18. WANG Baoguo,LIU Shuyan,YAN Weige.High-order strongly compact schemes and the treatment of boundary conditions[J].Journal of Beijing Institute of Technology,2003,23(1):13-18.(in Chinese)
[25]	王保国,刘淑艳,潘美霞.强紧致六阶格式的构造及其应用[J].工程热物理学报,2003,24(5):761-763. WANG Baoguo,LIU Shuyan,PAN Meixia.Sixth order accurate strongly compact scheme and its application[J].Journal of Engineering Thermophysics,2003,24(5):761-763.(in Chinese)
[26]	王保国,刘淑艳,张雅.非结构网格下非定常流场双时间步长加权ENO-强紧致杂交高分辨格式[J].工程热物理学报,2005,26(6):941-943. WANG Baoguo,LIU Shuyan,ZHANG Ya.A new weighted ENO-strongly compact difference schemes and dual-time method for the solution of the unsteady flows on unstructured grids[J].Journal of Engineering Thermophysics,2005,26(6):941-943.(in Chinese)
[27]	王保国,刘淑艳.三点强紧致加权高分辨率高阶格式及其应用[J].工程热物理学报,2006,27(6):929-932. WANG Baoguo,LIU Shuyan.Three point strongly compact weighting high resolution high order scheme and its application[J].Journal of Engineering Thermophysics,2006,27(6):929-932.(in Chinese)
[28]	王保国,刘淑艳,李翔.基于Nash-Pareto策略的两种改进算法及其应用[J].航空动力学报,2008,23(2):374-382. WANG Baoguo,LIU Shuyan,LI Xiang.Two improved algorithms based on Nash-Pareto strategy and their applications[J].Journal of Aerospace Power,2008,23(2):374-382.(in Chinese)
[29]	王保国,刘淑艳,姜国义.高阶格式及其在内外流场计算中的应用[J].航空动力学报,2008,23(1):55-63. WANG Baoguo,LIU Shuyan,JIANG Guoyi.High-order scheme and its application to the computation of the internal and external flow field[J].Journal of Aerospace Power,2008,23(1):55-63.(in Chinese)
[30]	王保国,刘淑艳,钱耕.一种小波神经网络与遗传算法结合的优化方法[J].航空动力学报,2008,23(11):1953-1960. WANG Baoguo,LIU Shuyan,QIAN Geng.Optimization method by combination of wavelet neural networks and genetic algorithm[J].Journal of Aerospace Power,2008,23(11):1953-1960.(in Chinese)
[31]	王保国,吴俊宏,刘淑艳.流场特性预测的两类高效方法[J].航空动力学报,2010,25(8):1763-1767. WANG Baoguo,WU Junhong,LIU Shuyan.Two effective methods for flow-field performance forecast[J].Journal of Aerospace Power,2010,25(8):1763-1767.(in Chinese)
[32]	王保国,吴俊宏,朱俊强.基于小波奇异分析的流场计算方法及应用[J].航空动力学报,2010,25(12):2728-2747. WANG Baoguo,WU Junhong,ZHU Junqiang.Method based on wavelet sigularity analysis for complicated flow and its application[J].Journal of Aerospace Power,2010,25(12):2728-2747.(in Chinese)
[33]	王保国,卞荫贵.超声速和高超声速进气道的数值模拟[J].力学进展,1992,22(3):318-323. WANG Baoguo,BIAN Yingui.On numerical simulations of supersonic and hypersonic inlet[J].Advances in Mechanics,1992,22(3):318-323.(in Chinese)
[34]	WANG Baoguo,LIU Qiusheng,BIAN Yingui.A high-resolution hybrid scheme for solving three dimensional Euler equations of high speed inlet flows[J].Journal of Thermal Science,1996,5(3):164-167.
[35]	WANG Baoguo,LIU Qiusheng,BIAN Yingui.Numerical simulation of high speed inlets using 3D Reynolds averaged Navier-Stokes equations and SIP finite volume scheme //First Asian Computational Fluid Dynamics Conference.Hong Kong:Hong Kong University Press,1995,3:1061-1066.
[36]	WANG Baoguo,GUO Yanhu.High order accurate and high resolation implicit upwind finite volume scheme for solving Euler/Reynolds averaged Navier-Stokes equations[J].Tsinghua Science and Technology,2000,5(1):47-53.
[37]	WANG Baoguo,BIAN Yingui.An LU-TVD finite volume scheme for solving 3D Reynolds averaged Navier-Stokes equations of high speed inlet flows //First Asian Computational Fluid Dynamics Conference.Hong Kong:Hong Kong University Press,1995,3:1055-1060.
[38]	王保国,李耀华,钱耕.四种飞行器绕流的三维DSMC计算与传热分析[J].航空动力学报,2011,26(1):1-20. WANG Baoguo,LI Yaohua,QIAN Geng.Three-dimensional DSMC calculation and heat transfer analysis of four capsules for hypersonic rarefied conditions[J].Journal of Aerospace Power,2011,26(1):1-20.(in Chinese)
[39]	李翔,王保国.高超声速飞行器气动热力学环境的研究 //龙升照,Dhillon B S.第10届人-机-环境系统工程大会论文集.美国,纽约:美国科研出版社,2010:294-298. LI Xiang,WANG Baoguo.Research of aerothermodynamics environments of hypersonic vehicle //Long Shengzhao,Dhillon B S.Proceedings of the 10th Conference on Man-Machine-Environment System Engineering.New York,USA:Scientific Research Publishing,USA,2010:294-298.(in Chinese)
[40]	LI Xiang,WANG Baoguo.CFD prediction of hypersonic blunt-cone configurations of multiple cases .AIAA Paper 2009-7385,2009.
[41]	王保国,孙业萍,钱耕.两类典型高速飞行器壁面热流的工程算法 //龙升照,Dhillon B S.第10届人-机-环境系统工程大会论文集.美国,纽约:美国科研出版社,2010:299-305. WANG Baoguo,SUN Yeping,QIAN Geng.An engineering algorithm for wall heat flux of two typical hypersonic vehicles //Long Shengzhao,Dhillon B S.Proceedings of the 10th Conference on Man-Machine-Environment System Engineering.New York,USA:Scientific Research Publishing,USA,2010:299-305.(in Chinese)
[42]	QIAN Geng,WANG Baoguo.A comparative study of Navier-Stokes and DSMC simulation of hypersonic flowfields .AIAA Paper 2011-765,2011.
[43]	钱耕,王保国.70°圆锥体行星探测器绕流热化学非平衡流场的Navier-Stokes方程与热环境的分析 //龙升照,Dhillon B S.第10届人-机-环境系统工程大会论文集.美国,纽约:美国科研出版社,2010:306-310. QIAN Geng,WANG Baoguo.Navier-Stokes solution of thermo-chemical non-equilibrium flow of 70° cone section of planetary and the analysis of thermal environment //Long Shengzhao,Dhillon B S.Proceedings of the 10th Conference on Man-Machine-Environment System Engineering.New York,USA:Scientific Research Publishing,USA,2010:306-310.(in Chinese)
[44]	李翔,王保国.轨道再入飞行器气动热力学环境研究[J].科技导报,2010,28(9):73-75. LI Xiang,WANG Baoguo.Aerothermodynamics environments of orbital return vehicle[J].Science & Technology Review,2010,28(9):73-75.(in Chinese)
[45]	钱耕,王保国.高超声速双锥体绕流的数值计算与流场分析[J].科技导报,2010,28(14):49-55. QIAN Geng,WANG Baoguo.Numerical computation and flowfield analysis of hypersonic flow over a double-cone body[J].Science & Technology Review,2010,28(14):49-55.(in Chinese)
[46]	王保国,李学东,刘淑艳.高温高速稀薄流的DSMC算法与流场传热分析[J].航空动力学报,2010,25(6):1203-1220. WANG Baoguo,LI Xuedong,LIU Shuyan.DSMC algorithm and heat transfer analysis of high temperature and high velocity rarefied gas flow[J].Journal of Aerospace Power,2010,25(6):1203-1220.(in Chinese)
[47]	李学东,王保国.稀薄气体高超声速流动的热环境研究 //龙升照,Dhillon B S.第10届人-机-环境系统工程大会论文集.美国,纽约:美国科研出版社,2010:203-206. LI Xuedong,WANG Baoguo.Research of thermal environments of hypersonic rarefied gas flow //Long Shengzhao,Dhillon B S.Proceedings of the 10th Conference on Man-Machine-Environment System Engineering.New York,USA:Scientific Research Publishing,USA,2010:203-206.(in Chinese)
[48]	孙成海,王保国,沈孟育.格子Boltzmann方法的质量扩散模型[J].计算物理,1997,14(4):671-673. SUN Chenghai,WANG Baoguo,SHEN Mengyu.Lattice Boltzmann models for mass diffusion[J].Chinese Journal of Computation Physics,1997,14(4):671-673.(in Chinese)
[49]	SUN Chenghai,WANG Baoguo,SHEN Mengyu.Adaptive lattice Boltzmann model for compressible flows[J].Tsinghua Science and Technology,2000,5(1):43-46.
[50]	SUN Chenghai,WANG Baoguo,SHEN Mengyu.Lattice Boltzmann models for heat transfer[J].Communications in Nonlinear Science and Numerical Simulation,1997,2(4):212-216.
[51]	孙成海,王保国,沈孟育.完全气体格子Boltzmann热模型[J].清华大学学报,2000,40(4):51-54. SUN Chenghai,WANG Baoguo,SHEN Mengyu.Thermal lattice-Boltzmann model for perfect gas[J].Journal of Tsinghua University(Science and Technology),2000,40(4):51-54.(in Chinese)
[52]	王保国,刘淑艳,姜国义,等.高超声速化学非平衡流动的数值计算[J].气体物理:理论与应用,2007,2(2):150-153. WANG Baoguo,LIU Shuyan,JIANG Guoyi,et al.Numerical computation of hypersonic chemical nonequilibrium flows[J].Physical of Gases:Theory and Application,2007,2(2):150-153.(in Chinese)
[53]	李学东,王保国.稀薄气体高超声速流动的非结构DSMC并行化计算[J].科学导报,2010,28(4):64-67. LI Xuedong,WANG Baoguo.Parallelization of DSMC method on unstructured grids for hypersonic rarefied gas[J].Science & Technology Review,2010,28(4):64-67.(in Chinese)
[54]	王保国,刘淑艳,黄伟光.气体动力学[M].北京:北京理工大学、北京航空航天大学、西北工业大学、哈尔滨工业大学、哈尔滨工程大学出版社,2005.
[55]	Akey N D,Cross A E.Radio blackout alleviation and plasma diagnostic results from a 25000 foot per second blunt-body reentry .NASA TN D-5615,1970.
[56]	Grantham W L.Flight results of 25000 foot per second reentry experiment using microwave reflectometers to measure plasma electron density and standoff distance .NASA TN D-6062,1970.
[57]	Inouye Y.OREX flingt-quick report and lessons learned //Proceedings of the 2nd European Symposium on Aerothermodynamics of Space Vehicles,ESTEC.Paris:ESA,1994.
[58]	王保国,黄虹宾.叶轮机械跨声速及亚声速流场的计算方法[M].北京:国防工业出版社,2000.
[59]	王保国.近20年AMME Lab在人机环境工程中的研究与进展 //龙升照,Dhillon B S.第11届人-机-环境系统工程大会论文集.美国,纽约:美国科研出版社,2011. WANG Baoguo.Research and progress of AMME Lab on man-machine-environment engineering in the past 20 years //Long Shengzhao,Dhillon B S.Proceedings of the 11th Conference on Man-Machine-Environment System Engineering.New York,USA:Scientific Research Publishing,USA,2011.(in Chinese)