可压缩多介质流体数值模拟中的Level—Set间断跟踪方法-计算物理.pdfVIP

计 算 物 理 第23卷第5期 v。1．23．No．5 乙Hl州ESE Ut乙UM上JU。lAl。lU州AL上JHY3l乙3 2006年9月 JUUK州AL sep．，2006 [文章编号】 1001．246x(2006)05．0518．07 可压缩多介质流体数值模拟中的 Level—Set间断跟踪方法 张学莹1， 赵 宁1， 王春武2 (1．南京航空航天大学航空宇航学院，江苏南京210016；2．南京航空航天大学理学院，江苏南京210016) [摘 要】 踪激波和捕捉界面，用瞰em姗问题解构造带状区域内的虚拟流体状态，对物理量的外推方法、间断附近虚拟流 体的构造、间断推进速度的计算等问题进行了研究．最后对可压缩多介质流体一维和二维守恒律方程组进行数 值模拟，数值计算采用通量重构的高精度wEN0格式，计算结果令人满意． [关键词】 多介质流体界面；kvel—set间断跟踪；虚拟流动方法；wENo格式 [中图分类号] 035；V21l [文献标识码] A O 引言 求解多介质流体的方法可以分为两大类：间断捕捉方法和间断跟踪方法，间断捕捉法在数值计算过程中 不考虑间断的存在，在整个流场采用几乎相同的计算格式自动捕获间断，期望数值解在间断处为很窄的过渡 广，然而该方法应用到多介质流计算时容易在交界面处引起虚假振荡，影响计算的准确性．间断跟踪法是把 断的移动速度．该方法在光滑区域和间断位置能同时获得高精度，但由于间断跟踪法把间断作为移动的内边 界来处理，必然会用到非正规网格，当网格变得很小时不可避免地引起cFL条件难以满足的“小网格”问题， 给数值计算造成不稳定，而且由于该方法要处理界面的拓扑结构的变化，当推广到多维问题时，交界面的重 构变得非常复杂． nuid 传统的GFM(曲ost 射造成的，波系结构依赖于介质的特性、诱导激波的强度和入射角度等因素，K锄i在文[1]中验证了这个事 方法容易推广到多维问题．传统的虚拟流动方法处理弱激波与界面的作用等问题取得了非常满意的效果，对 于界面为接触间断的情况，压强和法向速度穿过界面是连续的，虚拟流体的压强和法方向速度用该网格点上 真实流体的压强和法方向速度，密度通过等压装配技术单侧熵插值求得．对于界面为强间断的情况，上面的 处理方法显然不再适用．在这种情况下，压强、法向速度和熵在界面附近存在跳跃．为了把传统的GFM推广 到求解一般的界面问题，Fedkiw改进了虚拟流体的构造方法，针对气一液两相流问题，用气体的状态决定界 面的压力，液体的流动状态决定界面的速度b’43；针对爆轰和爆燃等强间断问题bo，根据守恒性构造虚拟流 虚拟流动状态构造虚拟流体的等压装配变量、压强和速度哺“…，Asl锄也对间断跟踪方法进行了相关的研 究㈨1引． 【收稿日期]2005一05一∞；[修回日期】2005—1l一2l [基金项目]国家自然科学基金及国防科技重点实验室基金(51479010205HK020r7)资助项目 [作者筒介】张学莹(1973一)，男，山东，博士，从事计算流体力学方面的研究，河海大学理学院210098 第5期 张学莹等：可压缩多介质流体数值模拟中的Ievel-Set间断跟踪方法 519 格式求解流体力学方程组．最后将本文的方法用于可压缩多介质流体一般间断问题的数值模拟，并对激波捕 捉方法、虚拟流动方法和LS间断跟踪方法的数值模拟结果进行分析对比．值得注意的是本文中的间断面可 以是不同介质的交界面，也可以是激波面． 1 Level．Set方程 引入一个kvd．Set函数西跟踪间断面运动变化过程，写成 西。+y·V西=O， (1) 界面的位置，则①0表示一种流体所在的区域，西0表示另一种流体所