多介质流体力学两步欧拉方法的模型封闭性研讨.pdf

第四届全国计算爆炸力学会议论文集 多介质流体力学两步欧拉方法的模型 * 封闭性研究 梁仙红 李征 刘军 何长江 （北京应用物理与计算数学研究所，北京 100088 ） 摘 要：阐述了多介质流体力学守恒方程的模型封闭性问题，讨论了MEPH3D 程序原有封闭性方法 的优缺点，并在 Yanilkin 等人方法的基础上给出了一种根据声速特征的自适应的模型封闭性方法。 数值试验和分析均表明新方法所得到的数值结果均优于原来的方法。 关键词：多介质流体力学；模型封闭性方法，压力增量相等，压力驰豫 Closing relations of two-step Eulerian Method for Multifluid Dynamics LIANG Xian-Hong ，LI Zheng ，LIU Jun, HE Chang-Jiang (Institute of Applied Physics and Computational Mathematics, Beijing, China, 100088) Abstract: The problem that the multifluid dynamic equations are not closed is described in this paper. We discuss the disadvantages of the original closing relation in MEPH3D. And we develope an adaptive closing relation method for MEPH3D according to the method produced by Yanilkin. Numerical tests show that the new method is better than the original one. Key word: multifluid dynamics, Closing relations of two-step Eulerian method, equivalence of pressure increment, pressure relaxation 欧拉方法可以分为两大类：一类是一步欧拉方法，另一类是拉氏步加重映步的两步欧拉方法。 两步欧拉方法由于其灵活性好，因此更适合推广到多流体数值模拟的计算中。两步欧拉方法中， 拉氏阶段的解变化是由声速相关的运动引起的，而输运（重映）阶段的解变化则是由对流相关的 运动引起的[4] 。和拉氏方法相比，欧拉方法本身有着非常重要的优点，如欧拉方法网格固定，不存 在因物质界面大变形而产生网格扭曲，从而导致计算无法进行下去的现象，计算结果的网格依赖 性也较小。 尽管欧拉方法具有以上优点，但是它依然面临着一些难以解决的困难。其主要困难在于混合网 格中物理量的计算。如在拉氏步中，由于多介质流体力学守恒方程的不封闭，如何给出恰当的封闭 模型描述混合网格中各介质的运动变化就是其中的难点之一；而在重映步中，如何确定混合网格中 物质界面并计算各物质的通量（重映量）则是欧拉方法研究的另一个难点问题。 本文第一部分将描述单速度多组分流体力学守恒方程的模型不封闭问题及其研究现状；第二部 分将描述 MEPH3D 程序中原有的模型封闭性方法，并考虑拉氏计算步中解运动变化的特点，给出根 * 受国家重点基础研究专项经费（2005CB321703 ）、国家自然科学基金及计算物理国家重点实验室 基金（9140C690101070C69 ） 42 · · 多介质流体力学两步欧拉方法的模型封闭性研究 据声速特征的自适应的模型封闭性方法；第三部分则给出了数值例子及相关讨论；第四部分为结论。 1 单速度多组分流体力学拉氏守恒方程组