不同形状界面的RM不稳定性数值研究.pdfVIP

354 范美如等：不而形状界面的RM不稳定性数值研宠 不同形状界面的RM不稳定性数值研究 范美如，翟志刚，司廷，罗喜胜，杨基明 (近代力学系，中国科学技术大学。合肥230027。中国) 摘要：本文用VAS2D数值研究了激波与不同形状的界面相互作用下的RM不稳定性的发展。VAS2D采用自适 应非结构四边形网格的方法，可以精细刻画界面的演变。针对激波Ma=1．2，界面内气体为sF6，研究了激波 与矩形、椭圆形以及三角形界面的相互作用，详细分析了界面的波系和涡量分布，简要讨论了界面演化过程。 结果表明，不同形状的流体界面在激波作用下，流场中的波系结构以及涡量的产生与分布有所不同，各种界 面的演变过程也有较大著异，因此不同形状的流体界面引起的RM不稳定性有不同的机理。 关键词：Richtmyer-Meshkov不稳定性，界面形状，平面激波 0引言 界面的相互作用。分析了在激波作用下，不同形状 界面的伸展率，为超燃发动机的混合问题提供了丰 当激波作用在两种不同密度流体的交界面时， 富的数据。从前人研究的内容来看，不同流体界面 流场会呈现出复杂的流动现象，这种由冲击加速度 的形状主要集中于球形、椭球形、正弦以及气柱等 引起的流场不稳定性为Richtmyer-Meshkov不稳定 界面，而对于其他不同形状的二维界面在激波的加 性(简称RM不稳定性，【1，2】)。RM不稳定性的研 速下所产生的RM不稳定性的研究较少。RM不稳定 究不仅在流动稳定性、旋涡以及湍流形成机理方面 性的发展一般认为有两种主要的机理：一是斜压机 有着重要的学术意义，而且还具有很广泛的工程研 制引起的涡量；二是压力扰动机制，主要包括各种 究意义，比如：超燃冲压发动机中的燃料与空气的 波系的产生与相互作用，以及激波聚焦等，这两个 混合与燃烧、超新星运动、惯性约束核聚变(Inertial 因素显然都与界面形状密切相关。基于以上考虑， Confinement Fusion)等方面自1957年，Markstein 本文对矩形、椭圆形和三角形三种不同形状的界面 提出RM不稳定性概念以来，科学工作者己对其做 在激波冲击下的发展演化进行了数值模拟研究，以 出了大量的研究，主要表现在理论分析、实验研究 增强对这两种机理的认识。 和数值模拟三个方面。由于RM不稳定性实验研究 在界面生成及流场测量方面的困难，数值模拟就成 1数值方法 为RM不稳定性研究的一个有效的手段。 RM不稳定性的数值研究工作开始于1972年， Vectorized (2·Dimensional·Axisymmetric Meyer和Blewett[3】运用拉格朗日算法对RM不稳定 Adaptive Solve，【7，8】)。该程序采用自适应的非结构四边形 性进行了数