强正激波后非平衡流动的计算方法研究.pdfVIP

强正激波后非平衡流动的计算方法 强正激波后非平衡流动的计算方法 张若凌乐嘉陵 中国空气动力研究与发展中心吸气式高超声速技术实验室，四川，绵阳，621000 摘要：对于求解正激波后非平衡流动，本文基于Euler方程，给出一种新颖的、接近解析的计算方法。本文 利用这个计算方法来说明空气非平衡流动的重要特点，分析球模型驻点线上和正激波波后流动参数的关系， 及非平衡流动里温度场和浓度场的相似性。 1．引言 求解冻结激波和平衡激波流动，已有常规的解析方法¨’21。对于求解正激波后非平衡流动， 本文基于Euler方程，采用一种新颖的、接近解析的计算方法∞1。经实践证明，这种方法是有 效的。如果采用常规的一维CFD技术，必须求解一个非定常问题。由于刚性问题的存在，对于 高效的计算方法要求较高。 本文利用这个计算方法来说明非平衡流动的重要特点，分析球模型驻点线上和正激波波后 流动参数的关系，及非平衡流动里温度场和浓度场的相似性。通过计算可以看出，这个简单实 用的求解方法是高效的，对于理解高超声速流动的机理，对于分析实验结果，都会有所帮助。 2．强正激波后非平衡流动的计算方法 当气体的耗散性质被忽略时，冻结和平衡正激波关系，均有解析的求解方法，这在普通的 气体动力学教材里，已经给出了详细的说明。具体地说，波前和波后物理量的关系 (Rankine—Hugoniot关系)为 2 p2甜2 P1甜l (1) P2+P2材2’=Pl+Plul‘ (2) ㈥ o、 己 o、 ”去+扣钏·咯+新2 其中，下标“1”和“2”分别表示激波间断前后的物理量。当(卜3)式表示冻结激波关系时， 在波后的组分质量百分比、振动温度与电子温度和波前相同。当(卜3)式表示平衡激波关系 时，波后气体处于平衡状态。需要求解的三个未知数为p、￡，和乃这三个方程已经构成封闭 的数学问题。此时波前和波后都是均匀流动。 对于非平衡情形，一维定常Euler方程是常微分方程组，即 dR 二==Q (4) dr 以三温度模型为例，F和Q为 p1“ pnⅣ E= ，Q=．q；．‰o (5) puu+P phu 肛，U pee纠 o．q．皱 把各组分方程叠加可以得到质量守恒方程，即 P“=Cl (6) 根据动量方程，可以得到 2 (7) UU+P P C2 根据总能量方程，有 38l 第十一届全国激波与激波管学术会议 一一————————————————