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

计 算 物 理 第28卷第6期 V01．28。No．6 OF CHINESE COMPUTATIONALPHYSICS 2011年11月 JOURNAL Nov．。201l 文章编号：1001-246X(2011)06-0810-07 入射和反射激波诱导轻气泡变形和失稳的计算 朱跃进， 董 刚， 刘怡昕， 范宝春 (南京理T大学瞬态物理重点实验室，江苏南京210094) 擒要：采用NS方程对激波诱导球形轻气泡的变形失稳过程进行三维数值模拟．验证本文模型的可靠性，讨论 气泡变形、涡环形成及其三维失稳过程．结果发现：球形轻气泡在入射和反射激波先后作用下，可以变形并沿 流向形成多个旋转方向不同的涡环，而斜压效应是产生这些涡环的主要原因，涡环不仅在自诱导效应和流场作 用下运动，而且自身会受到扰动而发生方位角不稳定．当涡环形成以小尺度流向涡为主的复杂结构时，即有可 能诱发流场湍流化． 关键词：激波；斜压效应；涡量；涡环；方位角不稳定 中图分类号：0382 文献标识码：A 0 引言 不稳定)，进而可能导致非线性的湍流混合．这种现象在自然界和工程中普遍存在，比如超声速燃烧¨‘21、爆 轰一1等，因此开展相关研究非常必要．从上世纪五六十年代开始，国外开展了激波与密度界面相互作用的实 等¨-93分别采用高速摄影的手段拍摄了弱激波诱导低密度球失稳变形的过程．研究结果表明，激波掠过密度 过渡为非线性过程，最终发展到湍流混合．近十年，国内外在数值模拟方面也开展了相关工作，Bagabir和 Drikakis¨钊研究了不同马赫数的激波与圆柱形气泡的作用过程；陈坚强等…1构造了高阶的计算格式，研究了 激波诱导混合增强过程；程军波等¨纠用GhostFluid方法对激波与柱形界面相互作用的二维可压缩流场进行 了直接模拟；Layes等¨1把流场近似成二维轴对称，对激波与球形气袍的相互作用进行了模拟；本文作者¨列 曾采用二维轴对称方法，对激波诱导乙烯球形火焰变形的过程进行数值模拟． 尽管上述实验测量和数值模拟对理解激波与密度界面相互作用的机理有很大帮助，但由于密度界面失 稳变形，及其向湍流混合的过渡本质上是三维的，因此光学测量手段和二维计算(包括柱形界面和二维轴对 称的球形界面)仍不能揭示界面的三维失稳过程中的一些细节，如小尺度失稳的机制；另一方面，反射激波 与变形密度界面的相互作用可进一步加速界面失稳，而相关研究报道并不多．因此我们使用NS方程，对入 射及反射激波诱导轻质气泡的变形和失稳过程进行三维数值模拟，着重揭示涡环的形成、三维演化和小尺度 失稳的机制． 1 数学模型 对于激波与密度界面相互作用的问题，本文采用单组分理想气体模型，气体的比热比(7)为1．4．尽管不 收稿日期：2010—12-02；修回日期：201l—06-01 基金项目：国家自然科学基金资助项目 通讯作者：董刚．E·mail：dgvehicle@yahoo．corn 万方数据 第6期 朱跃进等：入射和反射激波诱导轻气泡变形和失稳的计算 8ll 同密度(组分)的气体7不同，但对于非反应流问题，单组分假设所造成的误差可以忽略¨0’14。1引，本文采用 的单组分三维NS方程如下： ’ 、一7 望+旦(p％)=0， t8)Ox二、7』71( 百+瓦‘叶)一菇+蔷， (2)‘z) 警+毒(pul叶)