流体界面不稳定性数值模拟中不同介质界面的处理方法.pdfVIP

计 　　算 　　物 　　理 第 18 卷 第 1 期 Vol . 18 ,No . 1 　 200 1 年 1 月 　CH IN ESE J OU RNAL OF COM PU TA T IONAL PH YSICS 　 J an . , 200 1 ( ) [文章编号] 　100 1246X 200 1 0 1002304 流体界面不稳定性数值模拟中 不同介质界面的处理方法 蔚喜军 , 　尤迎玖 (北京应用物理与计算数学研究所 ,计算物理实验室 ,北京 　100088) [摘 　要] 　在界面不稳定性的数值模拟中 , 以往的方法[ 1 ,4～8 ] 局限于密度不同的单种介质或绝热指数相同的两种 介质问题 ,在计算绝热指数不同介质间的扰动时 ,界面上的物理量会出现非物理跳动. 针对于此 ,在 Level Set 方法 基础上引入了 Front Tracking 方法[2 ] 的思想 ,即在界面上强制加入接触间断边界条件 ,从而避免了界面上参数的非 物理振荡. 用三维程序 TVD/ AC 对两种不同介质 ,在单波扰动下 RayleighTaylor 不稳定性的发展进行了计算 ,得到 了扰动随时间发展的图象. [ 关键词] 　界面不稳定性 ; Level Set 方法 ; Front Tracking 方法 ; TVD 格式 [ 中图分类号] 　O35 　　　[文献标识码] 　A 0 　引言 组产生任何影响. 这种方法思路清晰且容易实现 ,因 在界面不稳定性数值模拟中,所用数值格式对 而得到了广泛的应用. 随时间发展界面上的接触间断应有较高的分辨能 遗憾的是 ,当两种理想流体的绝热指数不同时 , 力. 传统的二阶 TVD 格式较好地解决了激波传播 用 TVD/ AC 格式附以L evel Set 方法计算所得界面 问题 ,但由于其格式耗散较大 ,对接触间断的求解往 上物理量出现了非物理振荡[ 1 ] . 对于这个 问题 , 目 往难以得到较高的精度. 因此用这种格式计算界面 前的看法和认识有分歧 ,有关的研究也较少 ,但这个 不稳定性问题 ,几个时间步长后初始极薄的扰动界 问题的解决对界面不稳定性的研究有着十分重要的 面就会弥散到几个网格上. 加人工压缩项的 TVD 意义. 本文认为 , 由于 L evel 函数的引入并不影响实 格式 ,即 TVD/ AC 格式[ 1 ] 能较好地解决这个问题 , 际流场中物理量的求解 , 因而 L evel Set 方法并不是 计算中数值耗散并不随时间步长的增加而增大. 界面上物理量产生非物理振荡的原因 ,其真正原因 在过去几年 中 , 对运动界面 的研 究使 Front 在于根据差分格式计算界面一侧点时 ,用到界面另 Tracking 方法[2 ,3 ] 和 L evel Set [4 ,5 ] 方法得到了充分 一侧满足不同状态方程点上的物理量. 鉴于此 , 本文尝试在处理界面时 , 引入 Front 的发展. Front Tracking 方法是一种 L agrange 方法 , 界面的运动是通过界面上质点的运动来表征 ,而这 Tracking 的思想 ,在界面上强制给出接触