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

计算物理 ›› 2001, Vol. 18 ›› Issue (1): 23-26.

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

蔚喜军, 尤迎玖

北京应用物理与计算数学研究所计算物理实验室, 北京 100088

收稿日期:1999-01-28 修回日期:1999-09-01 出版日期:2001-01-25 发布日期:2001-01-25
作者简介:蔚喜军(1959-),男,内蒙古,副研,博士,从事计算流体力学方面的研究,北京8009-26信箱,100088.
基金资助:
国家自然科学基金(19771012)和中国工程物理研究院科学基金(970114)资助项目

A COMPUTATIONAL SCHEME OF DIFFERENT MEDIUM INTERFACES IN THE NUMERICAL SIMULATION OF FLOW INTERFACIAL INSTABILITY

YU Xi-jun, YOU Ying-jiu

Institute of Applied Physics and Computational Mathematics, Laboratory of Computational Physics, Beijing 100088, P R China

Received:1999-01-28 Revised:1999-09-01 Online:2001-01-25 Published:2001-01-25

摘要/Abstract

摘要： 在界面不稳定性的数值模拟中,以往的方法^[1,4~8]局限于密度不同的单种介质或绝热指数相同的两种介质问题,在计算绝热指数不同介质间的扰动时,界面上的物理量会出现非物理跳动.针对于此,在Level Set方法基础上引入了Front Tracking方法^[2]的思想,即在界面上强制加入接触间断边界条件,从而避免了界面上参数的非物理振荡.用三维程序TVD/AC对两种不同介质,在单波扰动下Rayleigh-Taylor不稳定性的发展进行了计算,得到了扰动随时间发展的图象.

关键词: 界面不稳定性, Level Set方法, Front Tracking方法, TVD格式

Abstract: In the numerical simulation of the interfacial instability, the schemes in [1,4~8] were limited probiems in one medium with different densities or two media with homologous adiabatic exponent.If the schemes are used in computing the disturbance of different medium interfaces with homologous adiabatic exponent, the physical guantities on interfaces may give rise to nonphysical oscillation.To overcome the problem, this paper introduces the front tracking scheme ^[2] into the Level Set method, name y the contact discontinuous boundary conditions are coercively added on interfaces so that the nonphysical oscillation is avoided.The evolution of Rayleigh-Taylor instability is examined by the three-dimensional TVD/AC program under the single wave.The disturbed figures with time are given.

Key words: interfacial instability, Level Set scheme, Front Tracking scheme, TVD method

中图分类号:

蔚喜军, 尤迎玖. 流体界面不稳定性数值模拟中不同介质界面的处理方法[J]. 计算物理, 2001, 18(1): 23-26.

YU Xi-jun, YOU Ying-jiu. A COMPUTATIONAL SCHEME OF DIFFERENT MEDIUM INTERFACES IN THE NUMERICAL SIMULATION OF FLOW INTERFACIAL INSTABILITY[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2001, 18(1): 23-26.

[1]	宋昱, 王飞, 郝鹏飞, 何枫. 考虑壁面接触效应的微流动数值模拟[J]. 计算物理, 2008, 25(1): 75-82.
[2]	宫翔飞, 张树道, 江松. 界面捕捉Level Set方法的(AMR)数值模拟[J]. 计算物理, 2006, 23(4): 391-395.
[3]	陈荣三, 蔚喜军. 一维多介质可压缩Euler方程的高精度RKDG有限元方法[J]. 计算物理, 2006, 23(1): 43-49.
[4]	肖志祥, 陈海昕, 李启兵, 符松. 湍流模式对转捩的初步研究[J]. 计算物理, 2006, 23(1): 61-65.
[5]	楼建锋, 于恒. CTVD格式数值计算非均质炸药爆轰问题[J]. 计算物理, 2005, 22(4): 358-364.
[6]	陈菲, 张梦萍, 徐胜利. 运动激波和气泡串相互作用的初步数值模拟[J]. 计算物理, 2004, 21(5): 443-448.
[7]	张镭, 袁礼. Ghost Fluid方法与双介质可压缩流动计算[J]. 计算物理, 2003, 20(6): 503-508.
[8]	姜洋, 赵宁, 唐维军. 界面不稳定性数值模拟中的虚拟流动方法[J]. 计算物理, 2003, 20(6): 549-555.
[9]	程军波, 唐维军, 李德元, 傅德薰. 利用Ghost Fluid方法模拟激波与柱形界面相互作用[J]. 计算物理, 2003, 20(3): 219-225.
[10]	董海涛, 张立东, 李椿萱. 基于熵条件二阶差分格式的嵌套网格分区算法[J]. 计算物理, 2003, 20(2): 102-106.
[11]	谢正桐, 李家春, 王丽丽. R-T界面不稳定性及湍流混合的大涡模拟[J]. 计算物理, 2002, 19(2): 108-114.
[12]	唐维军, 赵宁, 李晓林, 张景琳, 蔚喜军. 随机扰动下三维流体界面不稳定性的并行计算[J]. 计算物理, 2001, 18(6): 539-543.
[13]	唐维军, 张景琳, 李晓林, 赵宁. 三维流体界面不稳定性的Ghost方法[J]. 计算物理, 2001, 18(2): 163-169.
[14]	周伟江, 姜贵庆. 迎风TVD格式在粘性流计算中的应用研究与改进[J]. 计算物理, 1999, 16(4): 401-408.
[15]	张钧, 赖东显, 张维岩. 辐射烧蚀界面的单模Rayleigh-Taylor不稳定性分析[J]. 计算物理, 1998, 15(4): 483-488.