RichtmyerMeshkov不稳定性后期混合数值研讨.pdfVIP

第十一届全国激波与激波管学术会议 Ri chtmye r—ieshkov不稳定性后期混合数值研究1 周海兵 熊俊张树道伍应艳 (北京应用物理与计算数学研究所北京100088) 性增长过程、非线性增长过程，最终转捩到湍流，因此研究界面不稳定性后期混合问题对理论和工程应用都 有着重要意义。在后期混合的研究中，随着界面扰动的发展，界面形状扭曲复杂，受计算规模的限制，通常 的界面追踪算法难以获得复杂的界面形态；另一方厩，I{M不稳定性后期混合研究不可避免的要处理激波与湍 流的相互作用问题，通常的高精度数值方法(谱方法、紧致格式)难以完成强激波问题的数值模拟，而常用 的激波捕捉算法数值粘性较大，不适合湍流的数值研究。本文采用PPM方法数值求解Euler方程；采用shyue 提出的考虑压力平衡的混合网格状态方程的处理方法，完成R_M不稳定性问题后期混合的数值模拟。界面不 稳定性后期混合具有明显的三维特征，二维计算不能分辨后期混合流体团之间三维扭曲拉伸现象，因此要求 三维数值模拟。另一方面，界面不稳定性后期，通过非线性作用，小尺度运动被充分激发，必须模拟从大尺 度到小尺度的级串现象，因此数值模拟要求很高的空间分辨率，要求大规模数值计算。由此我们采用MPI、 应用区域分解方法完成程序并行化，并行程序具有较好的可扩展性。数值计算表明，本方法可以较好的模拟 有强激波作用的可压缩流动，所得结果与理论和数值实验结果一致。采用PPM方法求解Euler方程可以得到 与理论结果一致的湍流能谱，具有正确的能谱衰减指数，可以用于研究高雷诺数、强激波问题，是对N—s方 程直接数值模拟的补充。数值模拟结果说明本方法具有高精度、鲁棒性，适合数值模拟R．M不稳定性后期混 合问题。 1 引言 在自然界和人们的工程应用中，两种不相混的流体之间交界面存在不稳定性问题，其中 初始扰动经过线性、非线性增长，最终转捩到湍流，因此界面不稳定性后期混合的研究对实 际问题有着重要意义。 在后期混合的研究中，主要存在两个困难。随着界面扰动的发展，界面形状扭曲复杂， 方法)将不足以刻画界面形态。Shyue瞬】(1998)通过界面两侧压力平衡和速度一致假定导出 流体体积分数的演化方程，结合Euler方程组，得出了计算模型，该模型保证了界面附近能量 方程的一致近似并避免压力震荡，可以求解高密度比强激波的多介质可压缩流动。马东军， 压缩流动。柏劲松Iv]还将此方法推广到多项式形式状态方程。另一个困难是R．M不稳定发展 后期，不可避免的要处理激波与湍流的相互作用问题。目前激波与湍流相互作用的数值研究 主要集中在弱激波情况。而对于工程实际中的强激波问题，通常的高精度数值方法(谱方法、 紧致格式)难以完成数值模拟，而常用的激波捕捉算法数值粘性较大，不适合湍流的数值研 ParabolicMethod) 究。Porter，Woodward等人【”q【vii】[“q【调经过研究认为，PPM方法(Piecewise 求解Euler方程与N—s方程的数值结果在湍动能、能谱等方面一致，具有正确的能谱衰减指数。 采用PPM方法求解Euler方程，研究高雷诺数、强激波问题，是对N．s方程直接数值模拟的 补充。 本文采用PPM方法数值求解Euler方程，采用Shyue提出的考虑压力平衡的混合网格状 态方程的处理方法，研究R．M不稳定性后期混合问题。本文分以下几个部分。第二部分介绍 本文采用的数值方法；第三部分给出数值算例，说明本文方法适合强激波可压缩流体界面不 稳定性后期的数值模拟；第四部分是总结。 2数值方法 采用体积分数描述的一维多介质Euler方程，结合Shyue提出的多介质可压缩流体混合网 格状态方程确定方法m】，控制方程为 1中国工程物理研究院科学技术基金资助课题。 336 塑+掣：02h 掣,+_a(puZ+p)：o —O(_pE)+—O(—p_Eu-+一pu)：o f：1，…，Ⅳ。 -：by—(O+甜_3y(O：0