国际期刊投稿平台
登录 | 注册
当前位置: 首页 > 亚太科研论坛 > 马赫数对激波—气泡相互作用界面演化理论分析
亚太科研论坛

亚太科研论坛

Asia-Pacific Research Forum

  • 主办单位: 
    未來中國國際出版集團有限公司
  • ISSN: 
    3079-3645(P)
  • ISSN: 
    3079-9945(O)
  • 期刊分类: 
    人文社科
  • 出版周期: 
    月刊
  • 投稿量: 
    4
  • 浏览量: 
    587

相关文章

暂无数据

马赫数对激波—气泡相互作用界面演化理论分析

Theoretical Analysis of the Interface Evolution in Shock-Wave Bubble Interaction by Mach Number

发布时间:2026-04-15
作者: 王尚 :北华航天工业学院 河北廊坊;
摘要: 激波—气泡相互作用(Shock–Bubble Interaction,SBI)是指激波与具有密度差异的气泡相互作用所产生的一类复杂流动现象。在这一过程中,会伴随激波反射与折射、界面变形以及多相流体混合等一系列物理过程。该现象广泛存在于自然界及多种工程应用中,例如超新星爆发、惯性约束核聚变、高超声速推进以及超声速燃烧等。因此,深入研究SBI的演化规律并探索其调控机制,对于揭示复杂流动的物理机理以及推动相关工程应用的发展具有重要意义。SBI属于典型的非线性、非平衡、多尺度复杂流动问题,涉及激波—界面相互作用、激波—激波相互作用以及多种波系的耦合演化过程。其主要特征包括激波的反射、透射、绕射与汇聚,气泡界面的涡量沉积,以及界面结构的剧烈变形与破碎等复杂动力学行为。本文重点论述激波马赫数对流动演化及界面动力学的影响,并进一步分析其对后续流动结构演化及混合过程的重要作用。
Abstract: Shock–bubble interaction (SBI) refers to a complex flow phenomenon that occurs when a shock wave interacts with a bubble that has a density difference relative to the surrounding medium. During this process, a series of physical phenomena arise, including shock reflection and refraction, interface deformation, and multiphase fluid mixing. This phenomenon is widely observed in nature and in various engineering applications, such as supernova explosions, inertial confinement fusion, hypersonic propulsion, and supersonic combustion. Therefore, investigating the evolution characteristics of SBI and exploring effective control strategies are of great significance for understanding the physical mechanisms of complex flows and promoting related engineering applications.SBI is a typical nonlinear, nonequilibrium, and multiscale flow problem involving shock–interface interaction, shock–shock interaction, and the coupled evolution of multiple wave systems. Its main features include shock reflection, transmission, diffraction, and focusing processes, vorticity deposition on the bubble interface, and severe deformation and fragmentation of the interface, leading to complex flow dynamics. In this study, particular attention is paid to the influence of the shock Mach number on flow evolution and interface dynamics, and its important effects on the subsequent flow structure development and mixing process are further analyzed.
关键词: 激波—气泡相互作用;马赫数效应;界面演化;流动混合
Keywords: shock–bubble interaction; Mach number; interface; flow mixing

引言

SBI是可压缩流动研究中的重要问题之一,其理论研究最早可追溯到20世纪60年代。1960年,Rudinger等人首次提出了气泡对加速度响应的解析模型,为描述冲击波作用下气泡界面的演化提供了初步的理论框架。随后,研究者逐渐从界面动力学和涡量生成机制等方面对SBI展开深入研究。1988年,Picone等人提出了新的理论模型,阐明了激波与气泡相互作用过程中气泡边界涡度的生成机制。研究表明,当激波通过气泡时,气泡界面处的密度梯度与压力梯度发生错位,从而产生条形涡度(baroclinic vorticity),最终形成与气泡几何形状相对应的涡结构,例如圆柱气泡产生涡丝对,而球形气泡则形成涡环结构。随着实验技术的发展,研究者开始通过实验手段深入探究SBI过程中的流动结构演化。2003年,Layes等人在方形截面激波管中利用多重曝光阴影照相与高速摄影技术,对不同强度激波作用下氦气(He)、氮气(N₂)和氪气(Kr)气泡的演化过程进行了实验研究。实验结果表明,轻气泡在激波作用下会发生明显的界面反转并形成涡环结构,而重气泡则表现为收缩并逐渐发展为伞状结构,近密度气泡仅发生轻微压缩。随后,Ranjan等人系统开展了强激波作用下不同密度球形气泡的实验研究。除了理论分析和实验研究之外,数值模拟也成为研究SBI的重要手段。1988年,Picone首次对激波与高密度和低密度气泡相互作用进行了数值模拟,其结果与Haas等人的实验结果基本一致,表明数值方法能够较好地捕捉实验中观察到的涡结构,并提出了描述气泡和气柱环量变化的模型。基于上述研究进展,本文将从基本流体力学公式与理论模型出发,探讨激波马赫数对气泡形态演化的调控作用。重点围绕斜压涡量生成、Rankine–Hugoniot关系、气泡加速度及界面环量沉积等关键公式,分析马赫数变化如何影响密度梯度与压力梯度的耦合强度,进而决定气泡界面的变形模式、涡结构生成及整体形态发展。通过上述理论分析,将从机制层面阐明不同马赫数下轻气泡、重气泡呈现不同演化特征的原因,为相关工程流动问题提供理论参考。

1论述内容

1.1 激波基本关系与马赫数的影响

激波的强度由马赫数Ma直接决定,Ma越大,激波对周围流体的压缩作用越强,激波前后的压力、密度跃变也越显著,而这种跃变是气泡界面受力、产生变形的根本驱动力。对于可压缩理想气体,激波前后的压强比、密度比满足Rankine–Hugoniot基本关系,这是分析马赫数影响的基础公式。

设激波前流体的压强为()、密度为(),激波后流体的压强为()、密度为(),气体比热比为(),空气(.4),氦气(.67),氪气(.66),则激波前后的密度比与压强比分别为:

(1)
(2)

对上述公式进行分析可知,马赫数Ma对激波前后的参数跃变具有显著调控作用:当(Ma<1)时,流体处于亚声速状态,无激波产生,气泡仅受微弱扰动,形态基本保持不变;当(Ma>1)时,形成激波,且随着(Ma)的增大,压强比与密度比均单调递增,激波的压缩效应显著增强。

具体而言,低马赫数条件下,激波前后的压强差较小,对气泡的冲击作用较弱,气泡界面仅发生轻微压缩或拉伸,形态无明显畸变;中高马赫数条件下,压强差急剧增大,激波对气泡的冲击力显著增强,气泡界面会发生剧烈变形,甚至出现界面反转、涡结构生成等现象,这也是不同马赫数下气泡形态差异的核心原因。

此外,结合气泡内外气体的密度差异,重气泡,马赫数的影响会进一步放大:对于轻气泡,较低马赫数即可引发明显的界面变形,而重气泡则需要更高的马赫数才能产生显著的收缩或畸变,这一规律可通过后续气泡加速度公式进一步验证。

1.2斜压涡量生成公式与马赫数对气泡涡结构的影响

气泡形态的演化与涡量的生成、沉积密切相关,而涡量的产生主要源于斜压效应,其核心机制可通过涡量输运方程描述。对于可压缩黏性流体,涡量(=▽u),其输运方程为:

=(▽)u++v(3)

式(3)侧共三项,其中第一项为涡量拉伸项,主要影响涡量的大小与方向;第二项为斜压项是SBI中涡量生成的主要来源;第三项为黏性耗散项,对于高马赫数激波作用下的短时间演化过程,黏性耗散效应可忽略不计,因此涡量的生成主要由斜压项决定。当激波穿过气泡界面时,气泡内外的密度差异导致界面处存在显著的密度梯度沿界面法向分布,而激波传播产生的压力梯度沿激波传播方向(与界面法向不共线),二者叉乘不为零,从而产生斜压涡量。且马赫数Ma越大,激波前后的压强差越大,压力梯度的强度也越大,进而导致斜压项的幅值增大,涡量生成量显著增加。

为进一步量化马赫数对涡量生成的影响,可采用界面环量生成率公式。环量是表征涡结构强度的重要参数,其生成率可表示为:具体表现为:低马赫数下,涡量生成量少,气泡界面仅发生轻微变形,无明显涡结构形成,气泡形态基本保持初始形状(球形、柱形等);中高马赫数下,涡量生成量显著增加,轻气泡会因涡量的作用发生界面反转,形成清晰的涡环结构,而重气泡则会在涡量与压力的共同作用下发生剧烈收缩,逐渐发展为伞状结构。

1.3 气泡加速度公式与马赫数对气泡运动及形态的调控

气泡在激波作用下的运动与变形,本质上是气泡在压力梯度作用下的加速过程,其加速度特征可通过 Rudinger等人(1960)提出的经典模型描述,该模型直接关联了激波压力、气泡密度与气泡加速度,是分析马赫数对气泡形态影响的关键公式。

设气泡外流体密度为,气泡内气体密度为,气泡界面的加速度为,则气泡加速度公式为:

=(4)

由式(4)知,气泡加速度的大小与压力梯度成正比,与气泡内外密度比负相关。故马赫数Ma越大,压力梯度越大,气泡加速度的幅值也越大,气泡的运动与变形越剧烈。

针对不同类型的气泡,马赫数的影响存在显著差异:

对于轻气泡(如氦气泡,()),式(5)中,()气泡加速度可近似为:

(5)

此时马赫数的增大对加速度的影响最为显著。低马赫数下,轻气泡加速度较小,界面变形平缓;中高马赫数下,加速度急剧增大,气泡界面会快速加速、反转,形成涡环结构,最终呈现出明显的“膨胀—反转—涡环”形态。对于重气泡,气泡加速度相对较小,因此需要更高的马赫数才能产生显著的变形。低马赫数下,重气泡仅发生轻微收缩,形态基本保持稳定;高马赫数下,压力梯度足够大,气泡加速度增大,界面会发生剧烈收缩,涡量沉积导致界面出现不规则畸变,最终发展为伞状结构。对于低密度气泡,气泡内外密度梯度较小,斜压涡量生成量少,且气泡加速度适中,因此无论马赫数处于何种范围,气泡形态仅发生轻微压缩或拉伸,无明显畸变与涡结构生成。

1.4 马赫数对气泡形态的影响规律

结合上述Rankine–Hugoniot关系、斜压涡量生成公式和气泡加速度公式的分析,可总结出马赫数对气泡形态的影响规律如下:

1.4.1低马赫数(1

激波强度较弱,压力梯度与密度梯度跃变较小,涡量生成量少,气泡加速度小。此时,轻气泡仅发生轻微变形,重气泡基本保持稳定,近密度气泡无明显形态变化,所有类型气泡均未出现显著涡结构。

1.4.2 中马赫数(2

激波强度中等,压力梯度显著增大,涡量生成量增加,气泡加速度明显提升。此时,轻气泡发生明显界面反转并形成涡环结构,重气泡开始出现收缩变形,近密度气泡出现轻微压缩,气泡形态差异逐渐显现。

1.4.3 高马赫数(Ma>4)

激波强度极强,压力梯度与涡量生成量急剧增大,气泡加速度达到最大值。此时,轻气泡的涡环结构更加清晰,界面反转更为剧烈;重气泡发生剧烈收缩,形成典型的伞状结构;近密度气泡出现明显压缩变形,整体形态发生显著畸变。

通过对SBI综合分析可以看出,马赫数是决定气泡受力强度、界面变形程度与流场结构演化的关键参数。在低马赫数条件下,激波对气泡的压力扰动较弱,气泡整体变形平缓,形态变化具有可逆性,涡量生成有限,流场结构基本保持稳定。随着马赫数的提高,激波前后的压力跃变显著增强,气泡所受的压力梯度与加速度明显增大,界面处斜压涡量快速生成并累积,促使气泡从轻微压缩逐步发展为明显变形、界面翻转乃至结构破碎。

不同密度特征的气泡对马赫数的响应存在明显差异。轻气泡在较低马赫数下即可出现显著变形与涡环结构,而重气泡由于自身密度较大,需要更高的马赫数才能产生明显的形态改变。近密度气泡则在较宽马赫数范围内保持相对稳定,仅在高马赫数条件下才表现出可观测的畸变。这一规律表明,气泡内外密度比与马赫数具有耦合作用,共同决定了气泡的最终形态。

通过对激波—气泡相互作用过程的理论分析,可以发现该类流动现象本质上是激波压缩效应、界面密度梯度以及流场动力学结构共同作用的结果。激波在传播过程中携带大量动量和能量,当其与具有明显密度差异的气泡相遇时,会在界面处触发一系列复杂的流动响应。这些响应不仅体现在气泡几何形态的变化上,还包括局部流场结构的重组、涡结构的生成与演化以及多相流体之间的混合增强。通过对相关理论关系的综合分析,可以更加清晰地认识到马赫数在整个相互作用过程中的核心调控作用。

首先,从宏观流动结构角度来看,马赫数决定了激波对周围流体的压缩强度,从而影响激波传播后的流场结构。当马赫数较低时,激波的能量输入有限,流体状态变化相对平缓,气泡界面受到的扰动主要表现为整体压缩与局部形变,其发展过程具有较强的稳定性。在这一阶段,流动结构仍保持相对简单的特征,界面附近仅出现弱小的涡旋结构,流体之间的混合过程较为缓慢。随着马赫数的逐渐增大,激波所携带的动量与压力跃变迅速增强,流场中形成更强的压力梯度与速度梯度,气泡界面由稳定变形逐渐进入不稳定发展阶段,界面结构开始呈现出明显的非线性演化特征。

其次,从界面动力学角度来看,气泡界面的变形过程体现了外部冲击与内部流体惯性之间的动态平衡。当激波到达气泡界面时,气泡前缘首先受到压缩并发生局部形变,而界面两侧流体密度差异的存在使得这一形变具有明显的不对称性。在后续演化过程中,界面附近的流体被不断加速并发生卷曲,从而逐渐形成具有一定尺度的涡旋结构。这些涡结构不仅改变了界面形态,还对周围流体产生卷吸作用,使得气泡内部与外部流体之间的物质交换逐渐增强。随着时间的推进,初始较为规则的界面逐渐发展为复杂的多尺度结构,表现出典型的非线性不稳定特征。

从流动混合角度来看,SBI为不同流体之间的快速混合提供了有效机制。由于激波能够在极短时间内产生强烈的速度差异与压力变化,界面附近会出现明显的剪切流动,从而促进界面扰动的增长。在这一过程中,涡结构起到了关键作用:一方面,涡旋可以将气泡内部流体卷入外部流场,使得原本清晰的界面逐渐模糊;另一方面,涡旋之间的相互作用会形成更加复杂的流动结构,例如涡对合并、涡环拉伸以及局部破碎等现象。这些过程共同加速了流体混合的发展,使得流动逐渐从初始的界面主导阶段转变为湍流混合阶段。

此外,气泡内部流体性质的差异也会对相互作用过程产生重要影响。不同密度气体在激波作用下所表现出的惯性特征不同,这会直接影响界面的响应速度与变形方式。例如,低密度气泡更容易受到激波驱动而发生快速形变,而高密度气泡则由于惯性较大,其界面变形过程相对滞后。这种差异使得不同类型气泡在同一马赫数条件下呈现出不同的演化模式,并进一步影响整体流场结构的发展。因此,在研究SBI问题时,不仅需要考虑激波强度,还需要综合分析气泡密度、气体性质以及界面初始形态等多种因素。

从工程应用角度来看,深入理解马赫数对气泡界面演化的影响,对于多个领域具有重要意义。在惯性约束核聚变研究中,界面扰动的增长会直接影响燃料压缩的均匀性,从而决定聚变效率;在高超声速飞行与超声速燃烧技术中,激波诱导的混合过程能够显著提高燃料与氧化剂之间的混合效率;在爆轰物理与冲击波工程领域,激波与不均匀介质的相互作用也会对能量传递和冲击载荷分布产生重要影响。因此,对SBI过程的系统研究不仅具有基础科学价值,同时也能够为相关工程技术的发展提供理论指导。

随着计算流体力学和实验诊断技术的不断发展,对激波—气泡相互作用的研究正逐渐向更加精细化和多尺度方向推进。未来研究可以从以下几个方面进一步展开:首先,在理论模型方面,需要建立更加完善的界面动力学描述方法,以更准确地刻画复杂界面形态的演化规律;其次,在数值模拟方面,应发展更高精度的界面捕捉算法,以提高对涡结构和界面破碎过程的模拟能力;再次,在实验研究方面,可以结合高速成像、粒子图像测速等先进测量技术,对流场结构进行更加精确的观测与分析。此外,还可以通过改变气泡形状、密度比以及多气泡排列方式等条件,进一步探索不同参数组合下的流动响应规律。

2结语

综合上述分析可得出结论:马赫数通过控制压力梯度、涡量生成与气泡加速度,确实从根本上决定了SBI的强度与模式;密度比则决定了气泡对马赫数变化的敏感程度。在实际流动控制、界面扰动与多相流工程问题中,可通过合理调控马赫数,实现对气泡变形、破碎与界面结构的有效控制,从而为相关数值模拟、实验设计及工程应用提供直接的理论依据。

参考文献:

  1. [1] Rudinger G, Somers L M. Behaviour of small regions of different gases carried in accelerated gas flows[J]. Journal of Fluid Mechanics,1960,7(02):161-176.
  2. [2] Picone J M, Boris J P. Vorticity generation by shock propagation through bubbles in a gas[J]. Journal of Fluid Mechanics,1988(189):235-251.
  3. [3] Layes G, Jourdan G, Houas L. Distortion of a spherical gaseous interface accelerated by a plane shock wave[J]. Physical Review Letters,2003,91(17):174502.
  4. [4] Layes G, Jourdan G, Houas L. Experimental investigation of the shock wave interaction with a spherical gas inhomogeneity[J]. Physics of Fluids,2005,17(02):028103.
  5. [5] Ranjan D, Niederhaus J H J, Oakley J G, et al. Experimental investigation of shock-induced distortion of a light spherical gas inhomogeneity[C]//Shock Waves. Proceedings of the18th International Shock Waves Symposium,2009.

联系我们
人工客服,稿件咨询
投稿
扫码添加微信
客服
置顶