反射激波作用下重气柱界面演化的PIV研究

在水平激波管中采用PIV方法研究了反射激波作用下SF6重气柱界面的发展演化。采用射流方法形成SF6无膜气柱界面,并以乙二醇作为示踪粒子。利用连续激光片光源结合高速摄影相机对流场进行显示,得到了反射激波作用下SF6气柱界面的发展过程。结果表明,入射激波的冲击会在界面上产生反向旋转的涡环结构,而反射激波的作用会在界面上产生与初始涡环旋转方向相反的次级涡环结构。此外,对反射激波作用后的流场图像进行PIV后处理,获得了流场连续的速度场和涡量场。获得的环量与已有的理论模型进行比较,取得了较好的一致性,也验证了本文实验方法的可行性。     

激波冲击下Air/SF6斜界面不稳定性实验研究

激波在不同密度介质上的交互作用在可压缩湍流上具有重要的基础价值。激波在界面上的作用会引起Richtmyer-Meshkov不稳定性。激波不正规折射时,流场存在更多复杂的涡。研究马赫数为1.23、1.41的激波在初始倾角β=60°的Air/SF6界面上非正规折射的情况。入射激波的切向冲击和法向冲击的相互作用,在界面处产生涡,折射波在壁面发生马赫反射。利用阴影显示技术,给出了界面演化和混合的过程。     

反射激波冲击单模界面的不稳定性实验研究

激波诱导流体界面失稳问题广泛存在于惯性约束核聚变、超燃冲压发动机、武器内爆等工程应用中，相关研究具有重要意义。本文在改进的水平激波管中开展入射激波及其反射激波诱导单模气体界面失稳的实验研究，采用线约束肥皂膜技术生成较为理想的空气/六氟化硫（air/SF₆）单模气体界面，借助高速纹影技术捕捉激波冲击界面后的详细不稳定性演化过程，重点关注反射距离对不稳定性发展的影响（反射距离定义为初始界面和激波管尾端固壁的距离）。研究发现，在一定反射距离范围内，反射激波作用后的扰动增长率几乎为一个恒定值，与反射距离无关（即与反射激波作用前的扰动振幅及增长率无关），随着反射距离的进一步增大，扰动增长率降低。通过实验测量值与理论预测值的对比，发现Mikaelian模型和Charakhch’an模型（采用合适的经验系数）均能对反射激波作用界面后的扰动增长率给予有效预测，且这两种模型的经验系数都依赖于反射激波作用前的界面演化状态，如果反射距离的改变引起/（不引起）界面演化状态的变化，则需要/（不需要）改变经验系数的值。     

不同密度梯度的多层流体界面上的Richtmyer-Meshkov不稳定性研究

在气-液界面Richtmyer-Meshkov（R-M）不稳定性研究的基础上,利用不同流体间密度的差异,构造了空气-硅油-水、空气-酒精-硅油两种流体粘度方向迥异的气-液-液三相界面,实验中以氮气作为高压驱动气体,在不同激波马赫数下对这两种界面R-M不稳定性后期尖钉与气泡区发展进行了测量与统计,通过对实验数据的分析,得出了相关规律,并用这些规律与已有的单层气-液界面R-M不稳定性研究成果作比较,得出了异同点。同时,还研究了两种三相界面在R-M不稳定发展中的差异。实验结果表明：当流体的粘度梯度方向（从小到大的方向）与激波方向一致时,界面失稳更加明显,湍流混合更为显著。     

径向入射激波聚焦实验和数值模拟研究

为得到不同入射激波马赫数对凹腔内径向入射激波聚焦过程的影响,对马赫数为1.25和2.18的径向入射激波在实验段内反射聚焦的过程进行了实验和数值模拟研究。用高速CCD拍摄了凹腔中气流流场的纹影图片,并采用保持强稳定性的龙格—库塔格式、加权基本无振荡算法和块结构动态自适应网格加密对激波反射聚焦的过程进行了数值模拟,实验和数值模拟结果吻合较好。复杂的激波形状及其之间的相互作用(如激波反射,激波聚焦和不同激波的相互作用)在实验和数值模拟中都能很好地观测到。通过比较不同马赫数入射激波的反射聚焦过程,发现较高的马赫数下聚焦后能产生射流,而较低的马赫数下却没有。同时发现,入射激波马赫数能影响二次激波出现的位置和强度,进而影响聚焦过程流场的紊乱程度。     

双翼尖涡Rayleigh—Ludwieg不稳定性实验研究

采用一种结构化矩形直机翼涡发生器产生一对大小不同、方向相反的翼尖涡，调节双涡涡量的大小比例Г1／Г2及其间距b，触发两涡Rayleigh—Ludwieg不稳定性。实验采用流动显示方法定性观察双涡相互作用过程，通过二维PIV（粒子成像测速）系统定量研究双涡相互作用特征，得到双翼尖涡中主涡及次涡的运动特性、环量-时间特性。对不同实验参数下残余环量比例进行分析，发现双涡涡量大小比例Г1／Г2在1．3～1．4、b为50ram时双涡相交削弱效果良好，能够实现翼尖涡强序削弱程度扶30％～40％。     

受限空间内激波与火焰作用的三维计算

激波与火焰的相互作用常发生在超声速燃烧和燃烧转爆轰过程中,为深入了解这一现象,采用带化学反应的三维Navier-Stokes方程,对平面入射激波及其反射激波与火焰的作用进行了计算研究,其中燃烧过程使用单步反应模型描述。研究结果显示:在受限空间内模拟激波与火焰作用,能更好地符合实验结果,从而体现出受限空间的三维效应;火焰在入射激波的作用下主要经Richtmyer-Meshkov不稳定而发生变形,此时火焰变形以物理作用为主,燃烧膨胀效果相对不明显;当反射激波与变形火焰再次作用后,火焰迅速膨胀变形,放热率维持在较高水平,此时化学反应过程起主要作用;在反射激波的作用下,变形火焰复杂三维涡结构的形成能强化热量与质量的输送,提高燃烧速率。     

激波与氦气泡相互作用的实验与数值研究

提出了一种在激波管中实现运动激波与球形气一气界面相互作用的研究方案。以肥皂膜作为气体交界面,高压火花为光源,采用阴影法,在马赫数M=1．2的情况下,在矩形激波管中开展了一系列的实验研究,获得了激波作用下球形氦气泡界面发展变化的相关实验结果。还通过求解二维轴对称的Euler方程,结合网格自适应策略,对实验过程进行了数值模拟。通过与实验结果的对比分析,证实了涡量的产生和分布对于界面的变形和发展有着重要意义。     

平面激波绕30°出峰后的绕射与反射

本文针对一平面激波S绕过坡角为30°的山峰后与地面作用所产生的波系图案等问题,通过激波管光测模型实验手段进行了分析和研究。实验是在入射激波马赫数Μ_i=1.18,1.32,1.74,2.00和2.29等五种条件下进行的。实验结果表明,在入射激波S和反射激波R之间出现了一接触面C,并与Skews关于平面激波绕凸角壁面时所得到的结果做了比较。与此同时,从我们所拍摄的纹影流场照片上可以看出,平面激波S绕过30°山峰后与地面作用所产生的反射波均属于马赫反射的范围。     

基于自适应网格的超声速液体射流破碎过程研究

为研究超声速气流中液体横向射流的破碎过程,采用脉冲背景光方法和VOF方法开展了实验和数值研究。为提高液体横向射流中气液界面和气流场特征捕捉的精确性,采用自适应网格技术对于气液界面、激波出现位置进行网格细化,计算得到了较为精细的气液界面、激波特征及涡系结构。研究结果表明:在低成本仿真模拟条件下,利用自适应网格计算得到的射流轨迹和轮廓与实验吻合较好,射流轨迹的最大误差为10%;射流初始段在超声速气流条件下,仍然存在一段高度约为1.9倍喷孔直径且圆柱形态保持较好的连续光滑液柱。随着喷注压降的升高,液柱的长度逐渐增大;主流气体流经液柱发生三维绕流,在射流附近和近壁面区域形成不断演化的反转涡对,反转涡对的形成加速了液体射流一次破碎过程。     

可控肥皂膜气柱界面与激波相互作用的实验研究

基于肥皂膜技术提出了一种生成气柱界面的简单方法。形成的界面不需要支撑网格,因此初始条件可以很好地控制。在水平激波管中结合高速纹影开展了平面激波与二维或三维气柱界面作用的实验。纹影照片显示实验结果杂波更少,演变中的界面也比文献的结果更对称。此外,还特别研究了由极小曲面特征导致的三维效应,发现三维效应使界面的演变和发展变慢了。本研究将有助于更多了解三维性对Richtmyer-Meshkov不稳定性发展的影响。     

一种用于RM不稳定性研究的竖直环形激波管的设计与验证

设计并加工了一套竖直环形同轴无膜激波管,可用于环形汇聚激波诱导下的Richtmyer-Meshkov不稳定性实验研究。与前人工作相比,本文在流体界面的形成以及流场的观测方法上做了较大的改进。通过实验和数值方法,对该竖直激波管产生的环形柱状汇聚激波的参数进行测量和分析,验证了同轴激波管形成柱状汇聚激波方法的可行性和可靠性。在界面形成方面,采用细丝约束肥皂膜技术形成正八边形气体界面,并利用数值方法考察了细丝对界面发展的影响。结果表明在界面发展的前期,细丝的影响几乎可以忽略。利用连续激光片光结合高速摄影相机对流场进行观测,获得了正八边形air/SF6气体界面在环形汇聚激波及其反射激波冲击下的演化过程,并与数值结果进行了对比,获得了较好的一致性,进一步验证了汇聚激波的对称性以及细丝约束肥皂膜技术用于形成多边形气体界面的可靠性。     

垂直交叉双圆柱绕流强迫振动数值模拟

采用虚拟体法数值模拟了平面单圆柱绕流强迫振动,以及空间垂直交叉双圆柱绕流下游圆柱强迫振动流场。验证了单圆柱强迫振动中的锁定状态以及相位突变现象,从而证实了该数值方法模拟振动流场的可靠性。研究了在雷诺数Re=150、间距为5倍圆柱直径、下游圆柱按正弦曲线振动时,对尾流场的影响。下游圆柱两端尾流场在振动的作用下,涡的横向间距增大,而中心尾流场由于受上游圆柱尾流的影响而保持原先的状态。下游圆柱在锁定区的振动使尾流场变得稳定,表现在流向二次涡结构的减少,以及尾流场的速度场频谱趋向单一化。     

Shock Wave Interaction with the Magnetopause

The magnetopause is in continuous motion and shock waves and impulsive acceleration events can occur. As an example, we show that the interaction of an interplanetary shock with the bow shock can generate a shock wave that after passing through the magnetosheath can interact with the magnetopause. In fluid dynamics, when a shock wave encounters a fluid discontinuity, the interface may become unstable and form bubbles and spikes. We consider this Richtmyer-Meshkov instability in magnetohydrodynamics. At the dayside magnetopause, the instability tends to be stabilized by the magnetic field. However, the shock wave interaction can initiate magnetic field reconnection for the southward IMF, which may be important in strong interplanetary shock events. This revised version was published online in August 2006 with corrections to the Cover Date.     

自由剪切层中的三维不稳定性

本文是在文[1]的基础上研究自由剪切层中由Kelvin-Helmholtz不稳定波发展而形成的展向大涡结构的三维不稳定性。以大涡结构为基本流动,将稳定性分析归结为二维特征值问题,用pseudo spectral(伪谱)方法数值求解。研究发现:在没有亚谐波存在的情况下,大涡结构的最不稳定的扰动波是流向波长和其相同。展向波数较高,有对流特性的三维扰动波。它在剪切层中的发展与展向涡量的分布有关,大涡结构的涡核不稳定性和辫子不稳定性是流向涡形成的主要力学机制。本文还给出了不同雷诺数下三维扰动波增长率与展向波数的关系,这些结论与实验及数值模拟结果基本一致。     

小振辐振荡圆柱诱导的二次旋涡流动

本文给出充满粘性不可压缩流体时同心圆柱间内柱作小振幅振荡所诱导的二次旋涡流动的解析解,对典型的柱径比给出容器外边界对二次旋涡流动流线、涡量分布和动能分布的影响。给出了等涡量线和等动能线分布的鞍点。采用数值求解非定常Stokes流动考察了二次旋涡流动的发展过程。二次定常旋涡流动的显示实验表明,解析解、数值计算和实验显示所示流谱基本一致。