利用Ghost+Fluid方法模拟激波与柱形界面相互作用.pdfVIP

[文章编号]l001．246x(2003)03—0219一07 利用GhostFluid方法模拟激波与 柱形界面相互作用 程军波1， 唐维军1， 李德元1， 傅德薰2 (】．北京应用物理与计算数学研究所，计算物理实验室，北京 100088；2．LNM，中国科学院力学研究所．北京100080) [摘 要] 利用Ghostnuid方法(后面简称Ghost方法)和7．model方法，在同样的时空离散精度条件下，对激 波与柱形界面相互作用的二维可压缩流场进行了直接模拟，并与实验结果相比较．从模拟结果看，在短时间内， Ghost方法和*model方法模拟的结果与实验结果基本相同，两种方法均正确地模拟出界面的位置、激波的强度和 速度．但随着时间的发展，具有较大数值耗散的卜model方法的计算结果与实验差别越来越大；而数值耗散较小 的Ghost方法能较为正确地模拟界面的运动． Level [关键词] set方法；Ghostnuid方法；7．model方法；R．M失稳 [中图分类号】 035 [文献标识码] A 于界面之间存在接触间断，二者的压力和法向速度 O 引言 相等，而切向速度和密度不等，为维持接触间断的性 针对激波的捕捉，有很多行之有效的格式，如 质，虚拟流场的压力和法向速度为该点真实的压力 TvD，EN0，wENO¨’23以及带群速度控制的紧致差分 和法向速度；而切向速度和密度则最好是由界面另 一边的真实流场延拓过来．这样处理使得构造的流 格式b’41等．但对于接触间断的模拟，即使在这些格 式的基础上加上人工压缩法或亚网格子法，由于始 场既保留了界面的基本特征(压力和法向速度)，又 终存在数值耗散而导致接触间断被抹平．而且，如果 保证了界面处各种物理量的连续，这使得高精度捕 接触间断两边是两种不同的流体，用Euler方法进行捉界面成为可能． 数值模拟时往往需要进行特别的处理．Ghost方法的 将物理量由真实流场延拓到虚拟流场的延拓公 出现能较好地解决这个问题． 式为 Ghost方法是Fedkiw¨1在kvelSet方法№o的基 ，，土n·V，=0， (1) 础上引入的．kvelSet方法的主要思想是利用kvel，表示需延拓的物理量；n表示单位法向矢量，n= set函数9的零等值面来表示两种流体的分界面，品；±为延拓的方向，+表示从口o的流体延 Set函数值为该点到界面的 界面周围其他点的kvel 有符号的垂直距离，不同符号的kvelset函数表示 拓到90的流体，而一正好相反． set 为了使kvelset函数缈始终为有符号的距离 不同的流体．Ghost方法的主要思想是利用kvel 函数将含有两种流体的流场分成各含有一种流体的 函数，必须在每个时间步后对其进行重新初始化，使 其满足方程 两个流场，每个流场由真实流场和虚拟流场(Ghost nuid)组成，每个流场可单独进行计算．这样可成功 V9l=1． (2) 地避开界面两边流体状态方程不同的问题，这是 方程(2)的边值问题有两种求解方法，其一是简单、 Ghost方法的最大优点． 计算时间长的迭代方法1；其二是复杂、计算时间短 09 al Ghost方法的关键是构造虚拟流场的物理量．由的Fast Metho(1[8]．后来，su