GhostFluid方法在一维可压缩多介质流体数值模拟中应用.pdfVIP

Ghost Fluid 方法在一维可压缩多介质流体数 值模拟中的应用 姚成宝*，浦锡锋*，田 宙*，张海波*，张柏华* *(西北核技术研究所，陕西西安，710024) 摘 要：将带有Isobaric 修正的Ghost Fluid 方法应用到多介质流体力学的数值模拟中，配合Level Set 界面识别 方法计算一维可压缩多介质无粘流体，有效地消除了在流体界面附近所产生的数值振荡和耗散。应用 WENO 格式求解Euler 方程和Level Set 方程，并对两种不同比热比的气-气交界问题进行了数值模拟，计算结果与相关 文献的结果一致，表明利用该方法来进行可压缩多介质流动的数值模拟是可行的。 关键词：Ghost Fluid，Isobaric ，Level Set ，WENO 引 言 随着科学技术的发展，可压缩多介质流动的数值模拟已经越来越受到人们的关注，尤其在 流体力学、航空航天以及核物理领域中有着重要的应用，例如激波与物质界面的相互作用、爆 炸、高速撞击、气泡的形成与运动等，都属于可压缩多介质流动的范畴。 在可压缩多介质流动问题的数值模拟中，由于现实条件的要求，流场中可能会出现多种流 体相互作用。界面两侧流体的状态方程会有很大的不同，流体的密度也会出现大梯度变化，这 给数值计算特别是交界面的数值模拟带来很大的困难。对于多介质流体，用计算单介质流体很 有效的守恒型差分格式往往不能得到理想的结果，在界面处会出现非物理振荡，不仅会降低计 算精度，而且会导致运动图像歪曲和失真，甚至使得计算难以进行。美国学者 Fedkiw 在 1999 [1] [2] 年开创性的提出了虚拟流动方法(GFM) ，并利用等压装配技术 解决了这一类间断界面的问 题，能有效的抑制界面处的非物理振荡，引起了国内外学者的广泛关注。 [2] 本文利用带Isobaric 修正的Ghost Fluid 方法 ，并和Level Set 方法相结合，用于模拟受理 想气体状态方程支配的多介质可压缩无粘流动。其中，利用Level Set方法识别流体界面 (用Level Set距离函数(x,t ) 的零等值面来表示两种流体的分界面) ，而用Ghost Fluid方法定义虚拟流体， 将具有两种介质的流场转化为两个单介质流场：流场1设定 0 为真实流体区域， 0 为 Ghost流体区域；流场2则刚好与之相反。Isobaric修正不仅将熵向Ghost 网格点延拓，而且还向 紧邻界面的真实网格点延拓，这样就可以大幅度减少所谓的“过热(overheating)”误差。 在此基础上，本文采用五阶WENO ( Weighted Essentially Non-Oscillatory) [3] 格式离散Euler [4] 方程的空间导数项，在时间维上利用具有三阶精度的TVD Runge-Kutta 方法 进行离散，并对 几个典型的具有气-气界面的多介质可压缩流场问题进行了数值模拟，得到了具有较高分辨率的 计算结果。 1 多介质流体物理模型 1.1 控制方程 二维可压缩守恒型Euler方程组为：   u   v       2    u u p   uv  (1)    0      2 