高超声速三维化学非平衡流动N－S方程数值解  

NUMERICAL SOLUTION OF NAVIER-STOKES EQUATIONS FOR 3-D HYPERSONIC FLOWS IN CHEMICAL NON-EQUILIBRIUM

作　　者：欧阳水吾[1] 谢中强[1] 

机构地区：[1]北京空气动力研究所

出　　处：《力学学报》1996年第1期90-98,共9页Chinese Journal of Theoretical and Applied Mechanics

摘　　要：详细地叙述了模拟三维高超声速电离空气化学非平衡粘性流全流场（前体和底部近尾迹）的数值方法，化学反应粘性流求解是基于带有化学源项的Ｎ－Ｓ方程的守恒形式，总的连续方程由组元守恒方程组所代替．对于高温电离空气流动，存在如下７个主要组元，它们依次为Ｎ＿２，Ｏ＿２，ＮＯ，ＮＯ￣＋，Ｎ，Ｏ和ｅ￣－．研究发展了求解完全耦合的，与时间相关的偏微分方程组的数值方法．利用一类新的迎风ＴＶＤ激波捕捉有限差分格式求解守恒形式的控制方程组，对高超声速层流有攻角钝锥体绕流流场，在非催化壁条件下、对包括底部近尾迹在内的全流场各参数，如组元浓度，电子密度和温度，以及物面压力分布和热流率分布得出了计算结果．对化学反应气体和完全气体模型的计算结果进行了比较，计算的流场电子密度与飞行实验结果符合很好．We described a numerical method for the simulation of a complete 3-D hypersonic viscous flowfield (forbody and near wake) with ionization and in chemical nonequilibrium, The viscous and chemically reacting flow computations are based on a conservation form of the Navier-Stokes equations with the chemical source terms. The golbal colltinuity equation is replaced by the individual species conservation equations. For high temperature ionized air there are seven primary constituents which may be ordered as follows: N_2, O_2, NO, NO￣+, N, O and e￣-. A numerical method to solve the fully coupled time-depelldent partial differential equations has been developed. The governing equations in conservation-law form are solved using a new class of finite-difference upwind TVD shockcapturing scheme accoding to Yee, H.D. The concentration of the chemical constituents, the electron number density, and the temperature distributions are obtained over the complete flowfield assuming noncatalytic wall for a 3-D laminar hypersonic blunt body flow. The surface pressure and heat transfer rate distributiolls are also obtained. These results are compared with the perfect gas results for the 71 km case. The computed electron densities in the flowfield agree well with experimental results.

关 键 词：化学非平衡 粘性流 N-S方程 TVD 

分 类 号：V211.3[航空宇航科学与技术—航空宇航推进理论与工程]