Page 132 - 《爆炸与冲击》2025年第12期
P. 132

第 45 卷               郭景琪,等: 马赫反射波系在平面重/轻界面的入射加载                                第 12 期

                                                                        Å     ã 1/2
                                                 p          2γ (      )    p
                                                                  2
                                                   −1            M −1 −      −1
                                                                  0
                                                 p 0        γ +1         p 0    
                                     tanθ = ±     Å     ã                                             (5)
                                                    p            p   γ −1       
                                               2
                                            γM −     −1             +
                                               0
                                                   p 0            p 0  γ +1
               式中:下标为      0  的参数表示波前物理量,“+”对应气体跨越激波后逆时针偏转,“−”对应气体跨越激
               波后顺时针偏转。需要注意的是,激波极曲线方法是基于平面激波假设展开的,即波系前后流场是均匀
               的。本文所有涉及的激波均为非平面激波,但假设激波在靠近三波点附近为近似平面的,则可以在三波
               点邻域内采用极曲线分析方法。本文研究问题处于                        RM  不稳定性的早期阶段(冲击压缩阶段)。在这一
               阶段,界面演化主要由冲击瞬间的压力梯度驱动,惯性力远大于黏性力,黏性的耗散作用影响十分有
               限。因此,采用无黏三激波理论仍能较好地预测界面演化前期行为。

                2    结果分析

                2.1    马赫反射波系入射加载界面后的波系演化

                   图  4  通过数值纹影图展示了马赫反射波系
                                                                                    SL 1
               入射加载界面的动态过程,定义入射激波到达界
                                                                                   SL 2
                                                                                          RS 2    RS 1
               面的时刻为初始时刻(t=0 μs)。图             4(a) 给出了                               TP 2
                                                                                   MS 2  MS 1    TP 1  IS
               t=−5 μs 时平面激波绕射刚体圆柱后形成的                2  组                                          Interface
                                                                                (a) t=−5 μs
               马赫反射波系,可以看到,该结构由平面入射激
               波  IS、 马  赫  杆  MS 和 1  MS 、 反  射  激  波  RS 和    ∇ ρ  ∇ p                         ∇ p  ∇ ρ
                                                      1
                                       2
               RS 、滑移线    SL 和 1  SL 以及连接各间断面的三              Interface                             Interface
                                   2
                 2
                                                                                           IP 2
               波点  TP 和 1  TP 组成。由于三波点        TP 和 1  TP 冲            MS 1    RW 1   RW 2
                           2
                                                      2
                                                                     IP 1
               击界面过程具有相似性且本文主要关心界面中                                         TS 1        MS 3
                                                                    (b) t= 5 μs            (c) t= 8 μs
               心位置的演化,因此将针对            TP 附近的马赫反射
                                         2
               波系(MS 、MS 、RS 、SL )与界面的作用过程展开
                                   2
                      1
                           2
                               2
               分析。
                   首先发生的是马赫杆           MS 和界面的作用,                                     RW 2         RW 1
                                          1
               t=2 μs 时,MS 开始冲击界面并发生规则折射,
                           1
                                                                                     MS 3
                                                                                              TS 1
               产生透射激波       TS 和反射稀疏波       RW ,如图    4(b)                      (d) t=11 μs
                              1
                                               1
               所示,马赫杆      MS 的加载将持续进行,直到三波
                              1
                                                                  图 4    马赫反射波系入射加载平面       SF 6 /N 2 界面
               点  TP 在 2  t=8 μs 时运动到界面上终止。同时期                              过程的数值纹影图
               发生的还有马赫杆         MS 和界面的作用,MS 冲击             Fig. 4    Numerical schlieren visualization of Mach reflection wave
                                                    2
                                  2
               界面产生透射激波         MS 和反射稀疏波         RW ,如             configuration at a planar SF 6 /N 2  gas interface
                                   3
                                                     2
               图  4(c) 所示,马赫杆    MS 也完成了对界面的加载。MS 和            1   MS 冲击界面时压力梯度与界面密度梯度不
                                    2
                                                                      2
               共线,诱导了斜压涡量的产生,图              4(b)~(c) 中红色区域展示了该过程。其中               MS 的加载诱导界面产生负
                                                                                     1
               涡量,MS 诱导界面产生正涡量。将上述                  2  次加载过程称为马赫反射波系的入射加载,经历入射加载后
                       2
               重气体和轻气体中会出现清晰的反射波系和透射波系,如图                            4(d) 所示,轻气体中的透射波系由              TS 和
                                                                                                       1
               MS 组成,重气体中的反射稀疏波系由                RW 和 1  RW 组成。
                 3
                                                           2
                2.2    马赫反射波系入射加载的理论求解
                   以上对马赫反射波系入射加载重/轻界面形成透射波系和反射波系的过程进行了定性描述,指出扰
               动激波的入射加载会在轻气体一侧形成两激波系的透射结构,在重气体一侧形成两稀疏波系的反射结
               构。本节将对于上述过程进行理论求解,以定量描述激波在重/轻界面折射的物理过程以及其对于界面
               演化的影响。
                                                         123201-5
   127   128   129   130   131   132   133   134   135   136   137