Van der Waals状态方程可压缩多流体流动的PPM方法.pdfVIP

算 物 it 理 第24 卷第3 期 CHlNESE lOURNAL OF COMPUTATIONAL PHYSICS VoI . 饵，No.3 犁渺7 年5 月 毕占3旦= {文章编每) 1∞1-刻6X(2ω7 )03-0287-08 Van der Waals 状态方程可压缩 多流体流动的 PPM 方法 郑建国， 马东军， 孙德军， 尹协远 (巾回科学技术大学工程科学学院，安徽合M! 23∞27) [111 要] 基于流体体积分数的混合烈多流体数值模型，将 pi但ewiae Parabolic MetI叫 (PPM) 方法应朋于可压缩 多流体流动的数值筷似，采用双渡近似求事多流体 van der Waale 状态方穆的Riem8nn 问题. 模拟高密度比且含有 激泼的可压缩多流体流动，典型的纯界面平移问题模拟结果表明，在接触间断的界面附近，压力和速度没有任 何的振荡且界商数值耗散都披控制在2-3 个问格之内，维和二维算例表明，该数值方法可以有效地处理接触 间断、激波和多维滑移线等物理问籍，并能够比其它多流体数值方法更精细地模拟多流体交界面. [笑.诩] 混合型多流体模lll; van der Waala 状态方程;高阶Codunov 格式; Riemann 问题 [申回分羹号J 0381; 0354 [立由标识码) A o 51 言 近年来多流体可压缩流动的数值模拟及运动界面的追踪成为计算流体力学研究的热点问题，它在天体 物理、核物理和材料科学中有着重要的应用.解决这类问题遇到很多隙碍:①不同的物质状态方程差别很大， 如气液、气固隅合问题，给模型方程的求解带来很大麻烦;②不同流体的密度比比较大.如液体和气体交界面 两侧流体的密度比一般很高，使得界面上的压力计算较为困难;@这类问题中往往有激波与界面的相互作 用，界丽的捕捉和避免压力的振荡也是很兢手的问题. 自由t对多流体可压缩流动的数值模拟主要有以下4 种 方法:①界丽追踪法(interface tracking) (ll .这是一种Lagrange 方法，界面依靠标记点来追踪;②流体体积法 (volume of fluid me由叫)(Z) • 这种方法箭要复杂的界面重构算法;③Level-Set Me曲00(3) 它用连续的Level-Set 函 数等值团来近似流体的交界面.这3 种方法各有优缺点，前两种方法精度比较高，但很难推广到三维和复杂 的拓扑结构，Level-Set 方法虽然简单但精度较低.④将混合流体的 Euler 方程组辅以体积分数的输运方程作 为控制方程组，使用高精度的数值方法进行求解，~p现在采用的方法1闻， 本文根据Saurel川和 Shyue(付关于流体界面两1Wtl 压力和速度平衡的假设，引人基于体积分数的混合型多 流体数值模型[8J 该模型形式简单且物理意义明确，保证了界面附近能量方穆的?致近似，并避免压力的振 荡.将经典的 PPM 方法应用到数值求解中，用双波近似求解 van der Waals 状态方程的Riemann 问题.一维和 二维高密度比且含强激波的数值算例显示，本方法可有效处理接触问断、激波和多维滑移线等物理问题. 1 基于体积分数形式的混合型多流体数值模型 真实流体一般采用 Mie-GrUneisen 形式的状态方程描述其本构关系{$1 . 即 p巾， e) =T(p)[ pe - pe,.r (p) ) -Pm(P)其中 r =(11ρ)(ôplôe)p 是介质的 Mie-G巾晴i盹n 系数，句，Ph!f作为一种参考的状态，只是密度的 函数.本文研究的 van der W嗣Is 状态方程是Mie-Grüneisen 方程的一种特殊形式，可以写成 p = ((y -1)/(1 - bp))(pe - p. + 印Z ) _ (P. + 叩勺 (1)