运动激波和气泡串相互作用的初步数值模拟.pdfVIP

物 计算 理 Vol .21.l\o. 5 目 211第 5 期 CH1NESE JOURNAL OF COMPUTATIONAL PHYSI巳 2004 帘月 单-..!!!. { 文..-liJ} I∞1-246X(2刷)0量收43-06 运动激波和气泡串相互作用的初步数值模拟 陈 菲， 张梦萍 ， 徐胜利 ( 中国剧学技术大学力学和机被工程系. 安徽合肥 23∞27) [ 摘 要] 通过对激族和li体界面相互作用诱导的大变形界面演化的数值l!Ii拟.撞证Level 制方法精确模拟多 个流体界丽的有效性. 采用2 阶迎凤TVD 求事欧位1r程得到li场解，采用5 阶 WENO 求解Level 耐方程追踪多流 体11-丽.采用 GFM 方法处理流体内界蜀，利用文 [ 仆 的计算结果佼核本文提If..在此基础上，对运动激液和气 泡峙相互作用远程进行了初步f，(伯筷似， 得到了不同时刻运动激波和圆管内的两个气泡作用后的演化图象.包 f. 压力和密度等值线分布目 计算纺』晨我叨:针对推广后的多界面Le\el 酬方穗，该方法仍可高质量地捕捉多个ifii 体界耐. [ 关键词) Level 制方法;迎风TVD 格式; GfM 方法 [申图分樊哥1 0335 [ 文献标识码) A 。 引盲 流体界面是流体力学的基本问题.在可压缩流动中，运动激波和界面相互作用问题有着重要的应用背 景.例如.可燃气体燃烧转爆轰所涉及的激波和火焰面相互作用.对水下爆炸， 当运动激波与气泡净相互作 用，激波强度会下降，降低了作用在水下航行能表面的载衔，有利于航行嚣的防护.因此，本文 目 的是研究运 动激波和气泡净的相互作用过程.该问姐包含:①流场计算和运动激波的捕捉;(2) 多个气泡界丽的捕捉. tll 目前，已有较为成熟的用于含激波流场激值模拟的高精度算法，如 TVO(2) , ENO[l [ , WENO .流体界丽模 1 拟大款有 MAC() ,VOF 刊 和近年提出的 level set 方法.本文是在 Fedki w 工作基础上开展的( 1] .Fedkiw 采用三阶 ENO- LLF 方法求解无粘流场，用 3 阶WENO 格式求解描述界丽的 level 制方程，采用 GFM(Cho附t Fluid MeÙ10d) 方法处理界面附近边界条件. 作为初步研究，基于和文[ 1)相同的思路，利用该文的计算结果验证本文计算程序，重点考察了运动激波 和两个界面相互作用后的界面演化过程.位验用 le\el set 方法描述多个流体界丽的有效性. 1 控制方程和定解条件 1. 1 徨制方程 假设爆炸流场是无粘的，可采用无粘 Euler 方程描述 豆旦 + 坦坦i+ 坦坦2=ò .h( w. r) ， (1) òt ò r 其中 ò =0， 1 分别对应方程(1 )为平丽、轴对称方程.在轴对称情况下，， Z 分别对应劣， r 坐标. w =(p , m , 11. , E) T, f ( w) = 1m ， 豆 + p ， 旦， 2{ hp)iT 飞 ρρρJ