柱状汇聚激波冲击单模气体界面的实验研究

在自行设计和加工的半环形汇聚激波管中,开展了柱状汇聚激波冲击单模Air/SF6气体界面的Richtmyer-Meshkov(RM)不稳定性实验研究。不同于以往的环形激波管,该激波管具有半圆形结构的实验段,使半环形管道和实验段都向外敞开,能够参考传统水平激波管的方式设置初始扰动界面和观测系统。采用线约束肥皂膜的方法形成单模初始扰动界面。利用高速纹影成像技术得到了柱形汇聚激波作用下界面演化的完整过程。为了研究初始振幅对界面演化形态的影响,实验中生成了三种不同初始振幅的单模界面,并获得了三种工况下界面位移和扰动振幅随时间的变化。结果表明,汇聚激波作用下的RM不稳定性与平面激波有很大差别,主要原因在于汇聚效应,包括结构汇聚、流动压缩以及界面反相等。     

流动肥皂膜气体界面及其在柱状汇聚激波作用下演化的实验研究

在Richtmyer-Meshkov(RM)不稳定性实验研究中,形成两种不同密度流体的初始扰动界面是前提和关健.本文提出了一种流动肥皂膜气体界面生成方法,其工作原理是由细丝构成的导流框从激波管实验段穿过,肥皂液从导流框的上端注入并在重力作用下在导流框中形成流动肥皂膜,膜的两测可以分别充入不同密度的气体从而形成稳定的气体...     

柱形和球形壳体内爆界面不稳定性数值分析

为了更好地研究柱形和球形构型下果冻界面不稳定性发展,避免内爆聚心反弹前后直角坐标网格计算导致的误差,提高对流场和界面位置的计算精度,通过应用考虑了MVFT程序的网格适应性,使其能够适用于柱形网格和球形网格下的界面不稳定性数值模拟,特别是能够保证内爆聚心反弹前后流场和界面计算的稳定性。应用改进的计算程序对两种构型下的界面不稳定性进行了数值模拟,并对二者界面演化规律进行了详细讨论和归纳。结果表明:对于内外半径相同的柱形和球形果冻,聚心反弹时前者半径较小,而后者反弹时刻早于前者,其向内聚心和向外运动的速度最大值大于前者,对内部气体的压缩强度强于前者。对于外边界带有正弦扰动情况,除遵循上述规律外,计算还给出了峰谷转换现象。该项研究结果为进一步深入进行复杂构型下界面不稳定性高精度数值模拟研究提供了一种分析工具。     

激波诱导下的气体界面不稳定性实验研究

界面不稳定性, 特别是Richtmyer–Meshkov (RM) 不稳定性, 是流体
力学中一项重要的研究内容, 无论在学术研究领域还是工程应用领域都有着
重要的研究价值和应用背景. RM 不稳定性问题自提出以来, 得到了学术界
广泛的关注, 其研究无论是在实验方法、数值模拟还是在理论分析方面都取
得了很大的进展. 在激波管中开展激波与界面相互作用的实验研究, 即研究
界面初始扰动在激波诱导下的演化规律, 是目前研究RM 不稳定性的重要手
段. RM 不稳定性实验研究包括3 个部分, 分别是激波的产生、界面的形成
以及流场的观测. 综述了RM 不稳定性的实验研究进展, 并针对目前研究的
局限性提出了RM 不稳定性今后实验研究的重点和方向: 汇聚激波作用下界
面不稳定性的发展规律; 激波冲击下多种形状及大振幅界面的演化机理; 三
维界面的RM 不稳定性发展规律; 可压缩湍流的形成与混合机理.      

激波冲击下air/SF6界面的Richtmyer-Meshkov不稳定性

实验研究了低马赫数(1.27)激波作用air/SF6界面的RM不稳定性问题.air/SF6初始正弦界面由厚度为1～2 μm的薄膜相隔得到,用阴影法测试界面演化过程.实验结果表明:由于不稳定性重流体(SF6)向轻流体(air)演化成“尖钉”结构,而轻流体演化为“气泡”结构;由于界面切向速度差的Kelvin-Helm-ho...     

柱面汇聚冲击界面不稳定性的理论计算

根据流体力学界面不稳定性的线性理论,结合流体动力学零级量的计算程序,计算了柱面汇聚的界面不稳定性问题,给出的扰动发展结果与实验符合较好。     

圆柱形汇聚激波诱导 Richtmyer-Meshkov不稳定的 SPH 模拟

汇聚激波诱导不同物质界面的Richtmyer-Meshkov(RM)不稳定现象在惯性约束核聚变领域有重要的学术意义和工程背景.基于网格离散的宏观流体力学方法由于数值扩散问题往往需要高阶精度算法才能准确追踪界面演化,且对大变形和破碎合并等复杂界面追踪也极为困难.光滑粒子流体动力学(smoothed particlehydrodynamics,SPH)方法采用纯拉格朗日算法,可以有效克服上述难点.但经典SPH算法需采用人工黏性处理强间断,在激波间断处往往会出现严重的非物理振荡,对于涉及强冲击不稳定性问题,很难达到理想的模拟效果.本文采用基于HLL黎曼求解器的SPH算法,实现了对强激波和大密度比物质界面的有效分辨和追踪.一维数值校核证明了代码的可靠性、健壮性,并进一步模拟了二维圆柱形汇聚冲击波冲击四边形轻/重气界面诱导的RM不稳定性问题,与已有实验结果进行了对比,发现模拟结果与实验结果吻合.通过分析界面演化过程中的密度及压力变化,发现本文所采用的方法可准确地追踪激波与界面作用的复杂界面和波系演化规律.研究结果为进一步理解和解释汇聚冲击条件下的RM不稳定性机理奠定了基础.      

柱形汇聚几何中内爆驱动金属界面不稳定性

采用自研的高保真度爆轰与冲击动力学程序,对柱形汇聚几何中内爆驱动金属材料界面不稳定性的动力学行为,进行了数值模拟研究。结果表明,首次冲击后至约12 μs,界面发展以RM（Richtmyer-Meshkov）不稳定性为主;12 μs后至冲击波聚心反弹加载前,界面聚心运动处于加速减速状态,界面发展由RT (Rayleigh-Taylor)不稳定性主导;冲击波聚心反弹加载后,界面发展又由RM不稳定性主导。另外,还研究了初始条件（初始振幅、初始波长、钢壳初始厚度和几何构型）对柱形内爆驱动金属材料界面不稳定性的影响。结果显示：初始振幅较大时振幅增长也较大;初始波长较小（模数较大）时振幅增长较小,而且存在一个截止波长;钢壳厚度会抑制扰动增长,也存在一个截止厚度;几何汇聚效应会使扰动增长速度更快。

     

反射激波冲击重气柱的RM不稳定性数值研究

数值研究了二维气柱在入射激波以及反射激波作用下的Richtmyer-Meshkov(RM)不稳定性发展规律, 采用有限体积法结合网格自 适应技术的VAS2D程序, 精确刻画激波和界面的演化. 入射平面激波的马赫数为1.2, 气柱界面内气体为六氟化硫(SF₆), 环境气体为空气, 激波管的尾端为固壁. 通过改变气柱与尾端之间的距离调节反射激波再次作用已经变形的气柱的时间, 获得不同时刻下已经变形的气柱形态、界面尺寸以及环量演化受到反射激波的影响. 结果表明, 反射激波再次作用气柱时, 气柱所处发展阶段不同, 界面演化规律以及环量随时间的变化也不相同, 反射激波与气柱相互作用过程中的涡量产生和分布与无反射情况差异较大, 揭示了不同情况下界面演化的物理机理.     

圆柱形汇聚激波诱导Richtmyer-Meshkov不稳定的SPH模拟

汇聚激波诱导不同物质界面的Richtmyer-Meshkov(RM)不稳定现象在惯性约束核聚变领域有重要的学术意义和工程背景.基于网格离散的宏观流体力学方法由于数值扩散问题往往需要高阶精度算法才能准确追踪界面演化,且对大变形和破碎合并等复杂界面追踪也极为困难.光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法采用纯拉格朗日算法,可以有效克服上述难点.但经典SPH算法需采用人工黏性处理强间断,在激波间断处往往会出现严重的非物理振荡,对于涉及强冲击不稳定性问题,很难达到理想的模拟效果.本文采用基于HLL黎曼求解器的SPH算法,实现了对强激波和大密度比物质界面的有效分辨和追踪.一维数值校核证明了代码的可靠性、健壮性,并进一步模拟了二维圆柱形汇聚冲击波冲击四边形轻/重气界面诱导的RM不稳定性问题,与已有实验结果进行了对比,发现模拟结果与实验结果吻合.通过分析界面演化过程中的密度及压力变化,发现本文所采用的方法可准确地追踪激波与界面作用的复杂界面和波系演化规律.研究结果为进一步理解和解释汇聚冲击条件下的RM不稳定性机理奠定了基础.     

Numerical investigation of Richtmyer-Meshkov instability driven by cylindrical shocks

In this paper, a numerical method with high order accuracy and high resolution was developed to simulate the Richtmyer-Meshkov(RM) instability driven by cylindrical shock waves. Compressible Euler equations in cylindrical coordinate were adopted for the cylindrical geometry and a third order accurate group control scheme was adopted to discretize the equations. Moreover, an adaptive grid technique was developed to refine the grid near the moving interface to improve the resolution of numerical solutions. The results of simulation exhibited the evolution process of RM instability, and the effect of Atwood number was studied. The larger the absolute value of Atwood number, the larger the perturbation amplitude. The nonlinear effect manifests more evidently in cylindrical geometry. The shock reflected from the pole center accelerates the interface for the second time, considerably complicating the interface evolution process, and such phenomena of reshock and secondary shock were studied. The project supported by the National Natural Science Foundation of China (10176033, 10135010 and 90205025). The English text was polished by Yunming Chen.     

激波诱导环形SF6气柱演化的机理

基于可压缩多组分Navier-Stokes控制方程,结合5阶加权本质无振荡格式以及网格自适应加密技术和level-set方法,数值模拟了平面激波（Ma=1.23）与环形SF₆气柱（内外半径分别为8和17.5 mm）界面的相互作用过程。相比于之前的实验结果,数值模拟结果揭示了入射激波在界面内4次透射过程中的复杂波系结构,观察到透射激波在内部界面传播时形成自由前导折射结构并向自由前导冯诺依曼折射结构转换的波系演变过程;另外,界面内的复杂激波结构诱导内部下游界面上的涡量发生了3次反向;在界面演化后期,内部界面形成的“射流”结构与下游界面相互作用,诱导界面形成一对主涡、一对次级涡以及一个反向“射流”结构。定量分析了环形界面长度、宽度、位移、环量以及混合率的变化情况,结果表明,内部气柱的存在减弱了前期小涡结构合并形成大涡结构过程中对界面高度与长度的影响,同时提高了重质气体与环境气体的混合率。     

Development of the Richtmyer-Meshkov instability upon interaction of a diffusion mixing layer of two gases with shock waves

Based on the previously formulated mathematical model of mechanics of a two-velocity two-temperature mixture of gases, the evolution of an initially disturbed mixing layer of two gases with different densities under the action of shock waves is considered in a two-dimensional unsteady approximation. Problems of interaction of shock waves with a sinusoidally disturbed diffuse layer are solved numerically. The predicted variation of the mixing-region width are in satisfactory agreement with experimental data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 3, pp. 3–11, May–June, 2005.     

多次激波诱导正弦扰动预混火焰界面失稳的数值研究

激波诱导火焰失稳是实际中常见的现象,为深入研究火焰失稳特性,采用三维单步化学反应的Navier-Stokes方程和9阶weighted essentially non-oscillatory (WENO)的高精度格式,对不同马赫数的入射激波及其反射激波多次诱导正弦型预混火焰界面失稳的现象进行了三维数值模拟,并对计算结果的可靠性进行了验证。研究结果显示,在激波的多次作用下,火焰界面的演变过程主要受Richtmyer-Meshkov (RM)不稳定特性和化学反应特性的双重影响,且随着入射激波强度的增强,上述2种特性均得到进一步强化。为构造体现反应性RM不稳定特性的参数,根据火焰界面混合区平均涡量和化学反应速率,提出了表征界面受不稳定性和化学反应影响的量纲一参数。通过分析发现,在同一入射激波强度下,该参数的对数形式随入射激波和反射激波的多次作用呈基本相同的线性增长趋势;对不同马赫数的入射激波,该参数对数形式的线性增长率也基本一致。这样的变化表明该量纲一参数能够反映反应性RM不稳定过程中火焰界面发展的内在规律。     

反射激波作用下三维凹气柱界面演化的数值研究

激波与气柱相互作用是Richtmyer-Meshkov不稳定性研究的经典案例. 单次激波与二维气柱相互作用已得到广泛关注, 但是反射激波再次冲击气柱 (尤其是三维气柱) 的研究较少, 相关演化规律和机理尚不清楚. 反射激波再次冲击演化中的气柱界面会产生新的斜压涡量, 影响涡量的输运和分布, 从而影响界面的演化. 本文采用自主开发的HOWD (high-order WENO and double-flux methods) 程序, 研究了马赫数为1.29的平面激波冲击N激波与气柱相互作用是Richtmyer-Meshkov不稳定性研究的经典案例. 单次激波与二维气柱相互作用已得到广泛关注, 但是反射激波再次冲击气柱 (尤其是三维气柱) 的研究较少, 相关演化规律和机理尚不清楚. 反射激波再次冲击演化中的气柱界面会产生新的斜压涡量, 影响涡量的输运和分布, 从而影响界面的演化. 本文采用自主开发的HOWD (high-order WENO and double-flux methods) 程序, 研究了马赫数为1.29的平面激波冲击N激波与气柱相互作用是Richtmyer-Meshkov不稳定性研究的经典案例.单次激波与二维气柱相互作用已得到广泛关注,但是反射激波再次冲击气柱(尤其是三维气柱)的研究较少,相关演化规律和机理尚不清楚.反射激波再次冲击演化中的气柱界面会产生新的斜压涡量,影响涡量的输运和分布,从而影响界面的演化.本文采用自主开发的HOW...     

流场非均匀性对非平面激波诱导的Richtmyer-Meshkov不稳定性影响的数值研究

基于Navier-Stokes方程组,采用可压缩多介质黏性流动和湍流大涡模拟程序MVFT (multi-viscousflow and turbulence),模拟了均匀流场与初始密度呈现高斯函数分布的非均匀流场中马赫数为1.25的非平面激波加载初始扰动air/SF6界面的Richtmyer-Meshkov (RM)不稳定性现象。数值模拟结果表明,初始流场非均匀性将会影响非平面激波诱导的RM不稳定性演化过程。反射激波加载前,非平面激波导致的界面扰动振幅随着流场非均匀性增强而增大;反射激波加载后,非均匀流场与均匀流场条件下的界面扰动振幅差异有所减小。进一步,定量分析流场中环量分布及脉动速度统计量揭示了前述规律的原因。此外,还与平面激波诱导的RM不稳定性进行了简单对比,发现由于非平面激波波阵面区域的涡量与激波冲击界面时产生的涡量的共同作用,使得非平面激波与平面激波诱导的界面失稳过程存在差异。     

收敛柱激波在粉尘气体中的传播

取稀颗粒群气固两相耦合的双流体简化模型,采用高解数值方法研究了收敛柱激波在粉尘气体中的传播和波后流场特性。通过与纯气体情况比较,揭示了固体颗粒出现对收敛柱激波传播特性的影响。