激波和强旋涡干扰的数值研究.pdfVIP

激波与强旋涡干扰的数值研究 张树海 (中国空气动力研究与发展中心，四川绵阳，621000) Shu Yong—TaoZhangChi—Wang Brown (DiviSionof Mathematics Applied University Island02912) Providence，Rhode 引言 激波和旋涡是空气动力学的两个基本特征，研究激波与旋涡的相互作用具有 实际价值和重要的理论意义。首先，当激波接触到旋涡的时刻，旋涡对激波产生 一种扰动，这种扰动沿着激波传播，导致激波发生变形。有趣的是，激波和旋涡 干扰引起的激波变形与激波打到固壁上产生的变形相似，都具有规则反射和马赫 反射。其次，激波碰到旋涡使旋涡产生局部压缩和膨胀，产生一种声波。研究发 现，激波和旋涡相互作用是气动噪音的一个重要源泉。 2 最早开展激波旋涡干扰实验研究的是Hol 用激波管产生平面激波，当激波扫过一个斜的机翼后，产生一个旋涡。之后，激 波碰到激波管端部墙壁后，反射回来与旋涡干扰。尽管实验中存在许多噪音波， 但他们发现了激波变形和由交替出现的压缩区和膨胀区组成的声波。此后，激波 与旋涡干扰的问题引起了广泛的研究兴趣，开展了许多实验研究“’2’”3、直接数值 模拟“‘。’””以及线性分析“5…，。这些研究主要集中于研究激波和旋涡的变形特征 以及噪音的产生机制和传播规律，给出了许多非常丰富的研究结果。 然而，大多数研究都集中在激波与弱旋涡的相互作用上，尚未发现激波与强 旋涡相互干扰的研究报告。本文研究的目的是研究激波与强旋涡干扰的流场结 构，以及噪音的产生过程。全文安排如下：第二部分介绍物理模型、汁算方法和 网格生成。第三部分为计算方法、并行程序以及网格敏感程度验证实验。第四部 分对流场结构、旋涡的变形与涡量的产生、噪音传播进行了系统分析。 2数值方法 2，1物理模型 {了 O 图2．1流动模型草图 图2．1是流动模型草图。计算区域为Ix枷≤xsXngh，，Y埘sY≤y唧】的区域。 的位置，旋涡从右向左移动。初始条件为： ％=M。rexp((1一r2)／2) “．=0 2)exp((1一r2)／2) ∞=Mu(2一r p(r)=吉【1一等M：“p(1一r2)】商 exp(1-r2)J’_ 肿郴f一等Mj 所有物理量将随r的增加而减少。但速度和涡量远比压力和密度降低的慢。 而当r=4时，速度和涡量的值已非常小，可以忽略不计。因此，初始旋涡位置选 择(Xo=4．4，Yo=O)． 2．2数值方法 采用五阶有限差分的WENO格式…1模拟二维非定常可压缩N—s方程： 等+善，舞+r／，鬻=击(文鼍+巩箐) (2．1) u=(p，P。，P。e)’为守恒型物理量，‘，，n，为坐标变换的导数值，F，G，F。，G。 为通量，可参考任何一本教科书。 对方程(2．1)的空间离散如下：非线性对流场的空间导数由五阶WENO格式 ”离散而得，该格式在光滑区以及任何极值点具有一致的五阶精度，而在间断区 域为三阶精度。对粘性项离散采用四阶中心差分格式，对时间导数的离散采用三 阶TVD型Rong—Kutta格式“”。 2．3网格生成 本问题的计算网格数在x方向为1280点，在y方向为960点。在激波附近 △y=O．05。计算网格由下面的函数生成： x：：一∞缸趣 其中o，B是控制网格间距的两个参数。