平面超音速自由混合剪切流动的数值模拟.pdfVIP

平面超音速自由混合剪切流动的数值模拟 李启兵符松苏铭德 (清华大学工程力学系，北京100084) 摘要本文将简化BGK格式应用于超音速来瘴情况的平面混合剪切层的二维数值模拟．本文 详细介绍了计算结果，并与其他方法所得的结果以及实验数据进行丁比较． 关■词BGK格式；超音速混合层；对流马赫数 O引 言 平面超音速混合剪切层的发展演化过程对研究流动的物理机制、检验计算方法以及高速推进器 的设计等工程应用起着十分重要的作用．人们进行了大量的实验和数值计算研究¨、2…．由于计算条 件的限制，数值计算多为时问发展或者二维空间发展流动．无论实验还是数值计算都发现了湍流剪切 层中大尺度的复杂结构，外加扰动对混台层的发展有很大的影响，宽带谱与随机偏移相位角扰动能模 拟更“自然”的空间发展情况．对于可压缩流体的混合剪切流动，特别对于超音速流动，对数值模拟的 方法有较高的要求．近年来，由xu(1998)等人提出了一种新的计算方法——BGK方法．该方法具有丰 富的物理内涵，实际应用中也有很强的鲁棒性，因而得到越来越多的应用．本文将B(x方法应用于平 面超音速自由剪切流动的空间发展的二维数埴模拟，得到了较好的结果， 1数值方法概述 统计力学中无外加作用力的情况下，二维BGK．B01咖㈣方程为 正+识+矾=一(，一g)／f (1) 这里，为空间x、y，时间‘，微粒运动的速度“、。，内自由度}的分布函数．g为流体分子hf“weⅡ分布 函数．由于分子在碰撞过程中质量、动量和能量守恒，所以应当有相容性条件： 口=l，2，3，4 f札(，一g)／rd三=o， (2) 其中d量=dud”d}，雌=d邑…d缸，对二维问题K=3，为内自由度数．而 饥={l，M，，(Ⅱ2+口2+井+…+靠)／2}7 (3) 于是由(2)式可建立有限体积法的公式 (4) 町1一昭+士f“‘(^。r一。≯dt+古f“2(钆。一瓯川。Ⅻ：。 其中宏观守恒量和分布函数之间的关系为 蓄■■+F瞽礓 Ⅳ=(P，∥，∥，芦)‘=I觑d暑 (5) r r F=I咖d三，G=I慨d昌 (6) 具体求解时，利用前一时间步已知的守恒量构造出单元内的初始分布函数，得到初始平衡态分布 函数，根据通解求出单元边界上随时间变化的分布函数，再由该分布函数计算出边界上的守恒量通 量，最后代入有限体积法格式得到下一时间步上的守恒量．单元侧面通量的计算可参考xu(1998)．为 减少计算量，本文采用练永生(1999)等人提出的简化方法． 2边界条件 上游边界给定进口的压力、温度和密度为常数．进口的速度分布u(y)=0．5(∽+u2) 边界采用简单波的处理方法． 3计算结果及分析 x50，网格数1350×450．此外，为减少计算量，加速流动的转捩，在进口 O．001，计算域L×L：150 附加一随时间随机变化的小扰动，其幅值不大于以的2％，且限于进口附近区域，对时均量和大涡结 构不会发生大的影响． 图1给出某一时刻压力和涡量的等值线．从压力图上可以看出人口附近交界面处有比较弱的斜 以看出沿空间方向上涡的产生、发展及相互作用的情况． 40 ▲∞