流体界面不稳定性数值模拟中不同介质界面的处理方法-计算物理.PDF

计 算 物 理 第 卷 第 期 ， 18 1 VoI . 18 No . 1 年 月 CHINESE JOURNAL OF COMPUTATIONAL PHYSICS ， 2001 1 Jan . 2001 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ［文章编号］ （ ） 1001-246X- 2001 01-0023-04 流体界面不稳定性数值模拟中 不同介质界面的处理方法 蔚喜军， 尤迎玖 （北京应用物理与计算数学研究所，计算物理实验室，北京 100088） ［， ］ ［摘 要］ 在界面不稳定性的数值模拟中，以往的方法 1 4 ~ 8 局限于密度不同的单种介质或绝热指数相同的两种 介质问题，在计算绝热指数不同介质间的扰动时，界面上的物理量会出现非物理跳动 针对于此，在 方法 . LeveI Set ［］ 2 基础上引入了Front Tracking 方法 的思想，即在界面上强制加入接触间断边界条件，从而避免了界面上参数的非 物理振荡. 用三维程序TVD/ AC 对两种不同介质，在单波扰动下RayIeigh-TayIor 不稳定性的发展进行了计算，得到了 扰动随时间发展的图象. ［关键词］ 界面不稳定性； 方法； 方法； 格式 LeveI Set Front Tracking TVD ［中图分类号］ O35 ［文献标识码］ A 响 这种方法思路清晰且容易实现，因而得到了广泛 ! 引言 . 在界面不稳定性数值模拟中，所用数值格式对 的应用. 随时间发展界面上的接触间断应有较高的分辨能 遗憾的是，当两种理想流体的绝热指数不同时， 力. 传统的二阶TVD 格式较好地解决了激波传播问 用TVD/ AC 格式附以LeveI Set 方法计算所得界面上 ［］ 1 物理量出现了非物理振荡 对于这个问题，目前的 题，但由于其格式耗散较大，对接触间断的求解往往 . 难以得到较高的精度. 因此用这种格式计算界面不 看法和认识有分歧，有关的研究也较少，但这个问题 稳定性问题，几个时间步长后初始极薄的扰动界面