入射和反射激波诱导轻气泡变形和失稳的计算.pdfVIP

算 计 物 理 Vo1. 28 ,No.6 第 28 卷第6 期 2011 年 11 月 CHINESE JOURNAL OF COMPUTATIONAL PHYSICS Nov. , 2011 文章编号:1 ∞1-246X ( 2011 ) 06.{)8 1O.{)7 入射和反射激波诱导轻气泡变形和失稳的计算 朱跃进， 董 刚， 刘怡昕， 范宝春 {南京理工大学瞬态物理重点实验室，江苏南京 21ω94) 摘 要:采用 NS 方程对激波诱导球形轻气泡的变形失稳过程进行三维数值模拟，验证本文模型的可靠性，讨论 气泡变形、涡环形成及其三维失稳过程.结果发现:球形轻气泡在入射和反射激波先后作用下，可以变形并沿 流向形成多个旋转方向不同的涡环，而斜压效应是产生这些涡环的主要原因，涡环不仅在自诱导效应和流场作 用下运动，而且自身会受到扰动而发生方位角不稳定.当涡环形成以小尺度流向涡为主的复杂结构时，即有可 能诱发流场满流化. 关键词:激波:斜压效应;涡量;涡环;方位角不稳定 中图分类号: 0382 文献标i只码 :A 。 引言 激披作用于密度界面会引起大尺度失稳(Richtmyer-Meshkov 不稳定)和小尺度失稳(如 Kelvin-Helmholtz 不稳定) ，进而可能导致非线性的揣流混合.这种现象在自然界和工程中普遍存在，比如超声速燃烧[1 -2] 、爆 轰[3] 等，因此开展相关研究非常必要.从上世纪五六十年代开始，国外开展了激披与密度界面相互作用的实 验研究. Markstein [4] 的经典实验清晰地显示了丁烧火焰在弱激波及其反射激撞作用下的失稳过程，后来为了 进一步理解激披与密度界面(如火焰表面)相互作用的失稳过程，Haas 和 Sturtevant[5] 、 Jacobs[6] 、 Layes 等[7-9) 分别采用高速摄影的手段拍摄了弱撒波诱导低密度球失稳变形的过程.研究结果表明，撒披掠过密度 界面后，可以导致界面的扰动随时间增长而增长，界面的 Richtmyer-Meshkov 不稳定性可由初期的线性过程 过渡为非线性过程，最终发展到揣流混合.近十年，国内外在数值模拟方面也开展了相关工作，8agabir 和 Drikakis [10) 研究了不同马赫数的激波与圆柱形气泡的作用过程;陈坚强等[11 J 构造了高阶的计算格式，研究了 撒波诱导混合增强过程;程军波等[12) 用 Ghost Fluid 方法对激波与柱形界面相互作用的二维可压缩流场进行 了直接模拟;Layes 等[创把流场近似成二维轴对称，对激披与球形气泡的相互作用进行了模拟;本文作者[臼] 曾采用二维轴对称方法，对激披诱导乙烯球形火焰变形的过程进行数值模拟. 尽管上述实验测量和数值模拟对理解激波与密度界面相互作用的机理有很大帮助，但由于密度界面失 稳变形，及其向揣流混合的过渡本质上是三维的，因此光学测量手段和二维计算(包括柱形界面和二维轴对 称的球形界面)仍不能揭示界面的三维失稳过程中的一些细节，如小尺度失稳的机制;另一方面，反射激波 与变形密度界面的相互作用可进一步加速界面失稳，而相关研究报道并不多.因此我们使用 NS 方程，对人 射及反射激波诱导轻质气泡的变形和失稳过程进行三维数值模拟，着重揭示涡环的形成、三维演化和小尺度 失稳的机制. 1 数学模型 对于激波与密度界面相互作用的问题，本文采用单组分理想气体模型，气体的比热比(γ) 为1. 4. 尽管不 收稿日期: 2010-12-02; 修回日期: 2011 -06 -01 基金项目:国家自然科学基金 资助项目 作者简介:朱跃进(1986 - ) ，男，江苏江都，博士生，主要从事计算流体力学和化学反应流的数值模拟， E-mail: 野jwind@yah曲，四 通讯作者:董刚， E-mail: dgvehicle@ yahω.ω田 第6 期 朱跃