超音速平板混合层被动标量的数值模拟研究.pdf

中国工程热物理学会 燃烧学 084262 学术会议论文 编号： 超音速平板混合层被动标量 的数值模拟 1 2 薛淑艳 张会强 王希麟 （清华大学航天航空学院，北京 100084） （电话Email：xueshuyan@） 摘要： 本文首先采用空间五阶精度紧致-WENO混合格式，求解二维可压缩Euler方程，直接数值模拟 了 Mc=0.2 和 1.0 的二维空间发展超音速平板混合层被动标量混合过程，并与 clemens[3] 的实验结果进 行比较。结果表明，被动标量时均值与实验结果符合很好，而被动标量脉动均方根值与实验值相比较 小。同时还研究了可压缩性对标量混合过程的影响，随着对流马赫数的增加，被动标量时均值剖面变 化不大，而脉动均方根峰值增大。 关键字： 超音速平板混合层 直接模拟 被动标量 可压缩性 0 前言 近年来，超音速冲压发动机的发展，引起了许多学者对超音速混合层的关注[1-4] 。 在超音速冲压发动机中，燃料与空气来流的有效掺混是关键。由于燃烧室内来流的速度 非常高，而燃烧室的尺寸又非常有限，这就导致气流在燃烧室内的滞留时间非常小，要 在这非常短的时间内达到燃料和空气的充分混合，以提高燃烧效率是一个难题。因此需 要研究作为燃烧室内基本流动模型的超音速混合层的详细流动机理，以及标量混合问 题，以提出有效的增混措施。 可压缩性将影响冲压发动机中燃料和氧化剂的混合，影响发动机的效率。对流马赫 数 M =(U - U )/(a +a )，是目前人们常用的衡量混合层可压缩性的参数。许多实验研究 C 1 2 1 2 已经证明：随着可压缩性的增加，混合层增长率将减小。除此之外，还抑制脉动强度， [1,2] 。Clemens 和 Mungal[5]对不同马赫数的平板混合层进行了实验研究，结 以及雷诺应力 果表明，可压缩性抑制了标量混合。 [3] [4] Ribault 和 Sankaran 等用大涡模拟方法数值模拟了可压缩混合层的被动标量混 合过程。但其中的不同亚格子模型对模拟结果有很大的影响。因此，本文直接数值模拟 了二维空间发展超音速平板混合层被动标量混合过程，并与 Clemens 等[5] 的实验结果进 行比较。在此基础上，对不同对流马赫数超音速平板混合层被动标量混合过程进行数值 模拟，研究了可压缩性对被动标量混合的影响。 1 物理模型 本文所模拟的二维空间发展超音速平板混合层是由两股速度为 U 、U ，压力为 P 、 1 2 1 P ，密度为 ρ 、ρ ，被动标量为 f 、f 的超音速气流在尖劈后相遇发展而成。计算域 2 1 2 1 2 示意图如图 1 所示。计算域不包括尖劈，尖劈后缘构成计算域的左边界。计算域大小为： L × L 0 , 4 0 0 =× −7 5 , 7 5 。 x y [ ] [ ] y U , ρ,p f 1 1 1， 1 x U , ρ,p