何惠琴，翟志剛，司 廷，羅喜勝

（中國科學(xué)技術(shù)大學(xué)近代力學(xué)系，合肥 230027）

反射激波作用下兩種重氣柱界面不穩(wěn)定性實驗研究

何惠琴，翟志剛，司 廷，羅喜勝

（中國科學(xué)技術(shù)大學(xué)近代力學(xué)系，合肥 230027）

在水平方形激波管中對兩種無膜重氣柱界面（分別是SF6和氬氣）在反射激波作用下的不穩(wěn)定性發(fā)展進行了實驗研究。氣柱界面采用射流技術(shù)形成，實驗采用連續(xù)激光片光源照射流場，乙二醇作為示蹤粒子，并用高速攝像機對流場進行拍攝，獲得了入射激波以及反射激波共同作用下，兩種不同氣柱界面的演化過程。實驗結(jié)果表明，兩種氣柱的Atwood數(shù)不同，界面演化速率不同，反射激波到達前后的界面形態(tài)不同。SF6氣柱在入射激波作用下會產(chǎn)生兩個比較明顯的反向的渦環(huán)結(jié)構(gòu)，而氬氣柱界面上由于產(chǎn)生的渦量較少，渦環(huán)結(jié)構(gòu)并不明顯。在反射激波作用下，SF6氣柱界面會出現(xiàn)明顯的次級渦對，而且次級渦對的旋轉(zhuǎn)方向與初始渦環(huán)結(jié)構(gòu)的旋轉(zhuǎn)方向相反。對于氬氣柱而言，在反射激波作用下雖然也產(chǎn)生了與初始渦環(huán)方向相反的次級渦對，但次級渦對始終未充分發(fā)展。這是因為反射激波作用時氬氣柱界面的Atwood數(shù)較小導(dǎo)致氬氣柱界面上產(chǎn)生的反向渦量較少。實驗結(jié)果充分表明了氣體Atwood數(shù)對界面不穩(wěn)定性的發(fā)展起到了較大的影響。

無膜氣柱；Richtmyer-Meshkov不穩(wěn)定性；反射激波；Atwood數(shù)；連續(xù)片光

Key words：membrane-less gas cylinder；Richtmyer-Meshkov instability；reshock；Atwood number；continuous laser sheet

0 引 言

激波作用在兩種不同密度流體界面時，壓力梯度和密度梯度的不重合會導(dǎo)致斜壓渦量的產(chǎn)生與分布，界面上的初始擾動會隨時間不斷增長，最終發(fā)展為湍流混合，這種現(xiàn)象稱為Richtmyer-Meshkov（RM）不穩(wěn)定性［1-2］。RM不穩(wěn)定性在慣性約束核聚變、超新星爆發(fā)和超燃沖壓發(fā)動機等方面有著重要的應(yīng)用價值，并且在渦動力學(xué)以及湍流形成等方面也有重要的學(xué)術(shù)意義，因而近年來受到了廣泛關(guān)注。

單次激波作用下的RM不穩(wěn)定性問題已經(jīng)得到了較多的研究。Jacobs［3-4］利用平面激光誘導(dǎo)熒光技術(shù)得到了激波與無膜重、輕氣柱相互作用的流動圖像，并獲得了流場的濃度變化；高寧等人［5］采用粒子圖像測速技術(shù)實驗測量了平面激波作用下重氣柱界面RM 不穩(wěn)定性演化的二維速度場。鄒立勇等人［6-7］也實驗研究了截面分別為圓形和橢圓形的重氣柱界面的RM不穩(wěn)定性問題。而在實際應(yīng)用中，界面通常要經(jīng)歷多次激波的作用，研究表明反射激波會加速湍流混合的進程［8］。Latini et al［9］和Hill et al［10］分別對反射激波作用下，單模和多模界面不穩(wěn)定性的演化規(guī)律開展了數(shù)值研究并與實驗結(jié)果做了對比分析。司廷等人［11］對不同時刻的反射激波條件下氣泡的演化規(guī)律進行了實驗研究。此外，關(guān)于反射激波條件下氣簾界面不穩(wěn)定性的發(fā)展也得到了大量的實驗研究［12-13］。

在RM不穩(wěn)定性發(fā)展的過程中，斜壓渦量是導(dǎo)致界面演變的主要因素［14-15］。對于不同Atwood（A=（ρ1-ρ2）／（ρ1+ρ2），ρ1，ρ2分別是界面內(nèi)外氣體的密度）數(shù)所對應(yīng)的界面，由于密度梯度不同，單次激波作用后界面上產(chǎn)生的渦量大小不同，導(dǎo)致界面演化有較大差異。而反射激波作用在變形的界面上會再次引起渦量的產(chǎn)生與分布，從而進一步影響界面的發(fā)展。Haehn et al［16］實驗研究了SF6氣泡和氬氣泡在反射激波作用下的發(fā)展演化，討論了不同Atwood數(shù)對界面不穩(wěn)定性發(fā)展的影響。王顯圣等人［17］對反射激波作用于重氣柱的RM不穩(wěn)定性問題進行了數(shù)值研究，沙莎等人［18］以及柏勁松等人［19］也數(shù)值研究了反射激波作用下重氣柱界面不穩(wěn)定性的發(fā)展，但反射激波作用下氣柱界面演化的實驗研究相對較少。氣柱是最直觀的二維界面之一，常用于檢驗數(shù)值方法和理論模型的可靠性，因此實驗研究反射激波作用下氣柱界面的演化特征是有必要的。

本文針對兩種不同Atwood數(shù)的重氣柱（SF6和氬氣）在反射激波作用下的演化規(guī)律進行了實驗研究。實驗采用連續(xù)激光片光源照射流場，在無膜氣柱中摻入乙二醇示蹤粒子，使用高速攝像機對流場進行連續(xù)拍攝，獲得了反射激波沖擊下兩種重氣柱界面的演化過程。界面的定量變化在一定程度上揭示了Atwood數(shù)對不穩(wěn)定性發(fā)展的影響。

圖1 測試段簡圖Fig.1 Schematic of the test section

1 實驗方法

實驗所用的激波管驅(qū)動段長度為2m，被驅(qū)動段4m，截面積為95×95mm2。圖1是測試段簡圖，反射距離L定義為端壁到氣柱中心的距離，實驗中L= 90mm。入射激波馬赫數(shù)Ma=1.23±0.01，反射激波馬赫數(shù)Ma=1.22±0.01（這里考慮了反射激波前氣流的速度）。實驗前，首先將實驗氣體充入位于實驗段上方的箱體中，并用氧氣濃度檢測儀測定箱體頂部出口處氧氣的濃度，當氧氣濃度下降至0.5%以下則認為箱體中充滿了實驗氣體。然后用煙霧機對液態(tài)乙二醇進行加熱使其變成霧態(tài)（平均直徑約為0.5 μm，其在RM不穩(wěn)定性實驗研究中作為示蹤粒子的可行性已得到驗證［20］），并充入箱體中與實驗氣體混合。打開箱體底部的閥門，實驗氣體和乙二醇粒子的混合物會在重力的作用下沿著管道自上而下流入測試段，并可通過控制閥門來調(diào)節(jié)流速的大小使氣柱穩(wěn)定。這樣在實驗段中就形成了直徑約4mm的無膜重氣柱界面。采用連續(xù)激光片光結(jié)合高速攝影相機（FASTCAM SA5，Photron Limited）進行流場顯示。連續(xù)激光器（SDL-532-15000T，15W，532nm）產(chǎn)生初始光斑直徑為10mm的連續(xù)激光相繼通過一個凹柱面鏡（焦距50mm）和一個凸透鏡（焦距500mm）能夠形成厚度約1mm、寬度約80mm的連續(xù)激光片光。激光片光經(jīng)過安裝在側(cè)壁上的石英玻璃進入測試段，照射位于管口下34mm處的界面截面。高速攝像機安裝在測試段的上方，對流場進行連續(xù)拍攝。具體的實驗參數(shù)見表1。

表1 實驗用到的參數(shù)（Ma，入射激波馬赫數(shù)；Mr，反射激波馬赫數(shù)；D0，氣柱直徑；A和A′分別為入射激波前后界面的Atwood數(shù)）（p0=1 atm，T0=293 K）Table 1 Experimental parameters used（Ma，the incident shock Mach number；Mr，the reflected shock Mach number；D0，the diameter of the gas cylinder；AandA′are Atwood numbers before and after the incident shock，respectively

圖2 入射激波和反射激波作用下SF6氣柱界面的演化，第一幀顯示的是初始界面，入射激波從下往上運動，Δt=33.3μsFig.2 Evolution of a SF6gas cylinder accelerated by an incident shock and its reshock.The first frame gives the initial interface and the incident shock wave moves from the bottom to the top，Δt=33.3μs

2 實驗結(jié)果與分析

圖2給出了入射激波和反射激波作用下SF6氣柱演化的動態(tài)圖。第一幀是初始界面，可以看出氣柱橫截面為圓形。對于入射激波和反射激波而言，上下游界面是相反的，為了避免混淆，本文不再區(qū)分上下游界面，而統(tǒng)一以“上下界面”代替。由于透射激波在界面內(nèi)部的聚焦，在“上界面”處產(chǎn)生了射流結(jié)構(gòu)，并且隨著界面發(fā)展慢慢轉(zhuǎn)移到了“下界面”（第1～4幀）。但射流結(jié)構(gòu)始終不明顯，且僅持續(xù)200μs就消失了，這與激波-氣泡相互作用時射流結(jié)構(gòu)的發(fā)展存在較大差距［21］。這種差異可能是由于透射激波在二維結(jié)構(gòu)中聚焦效應(yīng)弱于三維結(jié)構(gòu)，產(chǎn)生的高壓峰值所致。由于斜壓渦量的正弦分布，渦量幅值最大值應(yīng)位于界面左右端點處，使得界面兩端開始形成反向的渦環(huán)結(jié)構(gòu)并逐漸發(fā)展（第5～10幀）。當反射激波經(jīng)過演化的界面時（第11幀），界面在流向上的運動受到一定的抑制。反射激波作用約60μs后，界面上衍生出次級渦環(huán)結(jié)構(gòu)，而且其旋轉(zhuǎn)方向與初始兩個渦環(huán)結(jié)構(gòu)的旋轉(zhuǎn)方向相反（第12～14幀）。次級渦環(huán)逐漸增大，愈來愈明顯（第16～27幀）。在這個過程中，“上界面”的運動速度幾乎為零，而“下界面”則發(fā)展緩慢，速度為負（這里把反射激波運動方向定義為負）。經(jīng)過充分的發(fā)展后，“上界面”獲得了正向的速度，而“下界面”獲得的仍是負向速度，這最終會導(dǎo)致“上下界面”的分離（第31～60幀）。到了后期，由于渦量大尺度的耗散，界面變得模糊不清（第61幀）。圖3給出了在入射激波和反射激波作用下，SF6氣柱界面渦量分布示意圖。在本文的反射距離下，反射激波到達界面時，初始渦環(huán)結(jié)構(gòu)的邊緣較光滑，和初始界面形狀類似，因此反射激波作用后，由于壓力梯度的反向，界面上產(chǎn)生了反向的次級渦環(huán)結(jié)構(gòu)。

圖3 入射激波（左）和反射激波（右）作用下SF6氣柱界面渦量產(chǎn)生與分布示意圖Fig.3 Schematic diagram of vorticity generation and distribution when a SF6gas cylinder is impacted by an incident planar shock（left）and reshock（right）

為了說明Atwood數(shù)對重氣柱界面不穩(wěn)定性發(fā)展的影響，本文還采用氬氣作為實驗氣體，研究了反射激波作用下氬氣柱界面的發(fā)展演化，實驗初始條件如表1所示。圖4給出了反射激波作用下氬氣柱的發(fā)展過程，這里忽略了反射激波到達前的實驗圖像。由于氬氣柱的Atwood數(shù)較小，入射激波作用之后，由于斜壓性產(chǎn)生的渦量較少，界面在反射激波到達時并沒有得到充分發(fā)展，渦環(huán)結(jié)構(gòu)尚未完全形成（第1幀）。反射激波作用之后，由于反射激波的壓縮作用，界面高度減小，這一點和SF6的結(jié)果類似，但氬氣柱隨后的演化和SF6氣柱的演化相比呈現(xiàn)出較大的區(qū)別。反射激波作用較長一段時間內(nèi)，氬氣柱界面上僅有兩個初始渦環(huán)結(jié)構(gòu)的存在（第3～26幀），并沒有像SF6界面中那樣會很快產(chǎn)生兩個次級渦環(huán)結(jié)構(gòu)。

圖4 反射激波作用后氬氣柱界面的演化。第一幀是反射激波到達時的界面形態(tài)，入射激波從下往上運動，Δt=28.6μsFig.4 Evolution of an argon gas cylinder accelerated by the reshock.The first frame is the initial interface for the reshock and the incident shock wave moves from the bottom to the top，Δt=28.6μs

界面的發(fā)展模式就像是僅僅經(jīng)過入射激波的沖擊一樣。然而，界面經(jīng)過了充分的發(fā)展之后，反射激波的效應(yīng)開始凸顯出來，即在初始渦環(huán)處衍生出兩個較小的反向渦環(huán)并隨時間不斷演化（第34～94幀）。相比于SF6氣柱演化中的次級渦環(huán)，氬氣條件下的次級渦環(huán)始終沒有得到充分發(fā)展。這主要是由于反射激波作用時氬氣界面的Atwood數(shù)減小較多，氬氣界面產(chǎn)生的反向渦量較少。

為了定量地說明反射激波作用下，不同Atwood數(shù)對界面不穩(wěn)定性發(fā)展的影響，圖5給出了界面位移以及高度隨時間的變化趨勢。其中，Xu是“上界面”位移，Xl是“下界面”位移，H是界面高度，并用界面的初始直徑D0進行了無量綱化，這里我們僅考慮反射激波作用后的流場。為了分析方便，令反射激波到達界面時為初始時刻并令此時“上下界面”的位移均為零。

圖5 SF6和氬氣柱結(jié)構(gòu)尺寸隨時間的變化曲線?！吧辖缑妗蔽灰芚u（a）；“下界面”位移Xl（b）；界面高度H（c）Fig.5 Changes of dimensions of the SF6and argon gas cylinders.The displacement of the upper interfaceXu（a）；the displacement of the lower interfaceXl（b）；and the height of the interfaceH（c）

反射激波經(jīng)過以后，SF6和氬氣柱的“上界面”在反射激波作用下仍向上緩慢運動，如圖5（a）所示，這是因為入射激波的作用在界面上產(chǎn)生了渦量而導(dǎo)致界面加速，經(jīng)過線性擬合得到SF6和氬氣柱的“上界面”速度分別為7.18m／s和6.94m／s。而SF6和氬氣柱“下界面”的位移則完全相反，SF6氣柱的“下界面”往下運動，速度大小約為12.43m／s，而氬氣柱的“下界面”往上運動，速度大小約為2.78m／s，如圖5（b）所示。這是由于SF6氣柱的初始渦環(huán)仍不斷發(fā)展以及后期的擴散作用致使界面向下運動，而氬氣柱的初始渦環(huán)較弱，整個界面仍是向上運動。由于SF6界面在反射激波作用下產(chǎn)生了發(fā)展較快的次級渦環(huán)，從而導(dǎo)致SF6界面的高度比氬氣界面的高度大得多，如圖5（c）所示。

3 結(jié) 論

實驗研究了兩種重氣柱界面（SF6和氬氣）在入射激波以及反射激波作用下的發(fā)展演化。采用乙二醇作為示蹤粒子，利用連續(xù)激光片光源結(jié)合高速攝像機進行流場顯示。實驗結(jié)果表明，反射激波的再次作用加速了不同流體之間的混合而且出現(xiàn)了新的流場特征。在入射激波沖擊下，SF6氣柱界面產(chǎn)生了兩個反向旋轉(zhuǎn)的初始渦環(huán)結(jié)構(gòu)。而在反射激波作用下，SF6氣柱界面會出現(xiàn)次級渦對，并成為發(fā)展后期的主要特征，而且次級渦對的旋轉(zhuǎn)方向與初始渦環(huán)的旋轉(zhuǎn)方向相反，這主要是由于反射激波在界面上產(chǎn)生了反向渦量引起的。而由于氬氣柱界面的Atwood數(shù)較小，入射激波的作用使得界面上的初始渦環(huán)結(jié)構(gòu)并沒有充分發(fā)展，而且反射激波在界面上引起的反向渦量也較少，氬氣柱界面上的次級渦環(huán)結(jié)構(gòu)并沒有像SF6氣柱界面中那樣明顯。定量結(jié)果表明，由于氣體Atwood數(shù)不同，入射激波和反射激波引起的界面演化形態(tài)有較大的區(qū)別，界面的演化速率和結(jié)構(gòu)尺寸也有較大的不同。

［1］Richtmyer R D.Taylor instability in shock acceleration of compressible fluids［J］.Commun Pure Appl Math，1960，13：297-319.

［2］Meshkov E E.Instability of the interface of two gases accelerated by a shock wave［J］.Fluid Dyn，1969，4：151-157.

［3］Jacobs J W.Shock-induced mixing of a light-gas cylinder［J］.J Fluid Mech，1992，234：629-649.

［4］Jacobs J W.The dynamics of shock accelerated light and heavy gas cylinders［J］.Phys Fluids A，1993，5：2239-2247.

［5］高寧，張光升，林俊.PIV技術(shù)在氣體界面不穩(wěn)定性實驗研究中的應(yīng)用［J］.高壓物理學(xué)報，2011，25（2）：177-182.Gao N，Zhang G S，Lin J.An application of PIV technique in gas interfacial instability experiment［J］.Chinese Journal of High Pressure Physics，2011，25（2）：177-182.

［6］鄒立勇，劉金宏，譚多望，等.弱激波沖擊無膜重氣柱和氣簾界面的實驗研究［J］.高壓物理學(xué)報，2010，24（4）：241-247.Zou L Y，Liu J H，Tan D W，et al.Experimental study on the membraneless heavy gas cylinder and gas curtain interfaces impacted by a weak shock wave［J］.Chinese Journal of High Pressure Physics，2010，24（4）：241-247.

［7］Zou L Y，Liu C L，Tan D W，et al.On interaction of shock wave with elliptic gas cylinder［J］.J Vis，2010，13（4）：347-353.

［8］Vetter M，Sturtevant B.Experiments on the Richtmyer-Meshkov instability of an air／SF6 interface［J］.Shock Waves，1995，4：247-252.

［9］Latini M，Schilling O，Don W S.High-resolution simulations and modeling of reshocked single-mode Richtmyer-Meshkov instability：Comparison to experimental data and to amplitude growth model predictions［J］.Phys Fluids，2007，19：024104.

［10］Hill D J，Pantano C，Pullin D I.Large-eddy simulation and multiscale modeling of a Richtmyer-Meshkov instability with reshock［J］.J Fluid Mech，2006，557：29-61.

［11］Si T，Zhai Z G，Luo X S，et al.Experimental investigation of reshocked spherical gas interfaces［J］.Phys Fluids，2012，24：054101.

［12］Balakumar B J，Orlicz G C，Tomkins C D，et al.Simultaneous particle-image velocimetry-planar laser-induced fluorescence measurements of Richtmyer-Meshkov instability growth in a gas curtain with and without reshock［J］.Phys Fluids，2008，20：124103.

［13］Balakumar B J，Orlicz G C，Ristorcelli J R，et al.Turbulent mixing in a Richtmyer-Meshkov fluid layer after reshock：velocity and density statistics［J］.J Fluid Mech，2012，696：67-93.

［14］Picone J M，Boris J P.Vorticity generation by shock propagation through bubbles in a gas［J］.J Fluid Mech，1988，189：23-51.

［15］Niederhaus J H J，Greenough J A，Oakley J G，et al.A computational parameter study for the three-dimensional shock-bubble interaction［J］.J Fluid Mech，2008，594：85-124.

［16］Haehn N，Weber C，Oakley J，et al.Experimental study of the shock-bubble interaction with reshock［J］.Shock Waves，2012，22：47-56.

［17］王顯圣，司廷，羅喜勝，等.反射激波沖擊重氣柱的RM不穩(wěn)定性數(shù)值研究［J］.力學(xué)學(xué)報，2012，44（4）：664-672.Wang X S，Si T，Luo X S，et al.Numerical study on the RM instability of a heavy-gas cylinder interacted with reshock［J］.Chinese Journal of Theoretical and Applied Mechanics，2012，44（4）：664-672.

［18］沙莎，陳志華，薛大文.激波沖擊R22重氣柱所導(dǎo)致的射流與混合研究［J］.物理學(xué)報，2013，62（14）：144701：1-10.Sha S，Chen Z H，Xue D W.The generation of jet and mixing induced by the interaction of shock wave with R22 cylinder［J］.Acta Phys Sin，2013，62（14）：144701：1-10.

［19］柏勁松，王濤，李平，等.激波反射對擾動界面湍流混合區(qū)影響的數(shù)值分析［J］.中國科學(xué)，2009，39（11）：1646-1653.Bo J S，Wang T，Li P，et al.Numerical analyse of the turbulence mixing of re-shocks with a perturbation interface［J］.Science in China：G，2009，39（11）：1646-1653.

［20］Prestridge K，Rightley P，Vorobieff P，et al.Simultaneous density-field visualization and PIV of a shock-accelerated gas curtain［J］.Exp Fluids，2000，29：339-346.

［21］Zhai Z G，Si T，Luo X S，et al.On the evolution of spherical gas interface accelerated by a planar shock wave［J］.Phys Fluids，2011，23：084104.

Experimental study on thereshocked RM instability of two kinds of heavy gas cylinder

He Huiqin，Zhai Zhigang，Si Ting，Luo Xisheng
（Department of Modern Mechanics，University of Science and Technology of China，Hefei 230027，China）

The evolutions of the Richtmyer-Meshkov instability of two kinds of membraneless heavy gas cylinder（SF6and argon），under reshock condition with a specific reflected distance，are experimentally studied in a horizontal square shock tube.Based on the jet technique，SF6and argon gas cylinders are generated respectively to study the effect of Atwood number on the development of the Richtmyer-Meshkov instability.For the visualization of the flow，the tested gases are mixed with glycol droplets which are generated by the fog generator.Illuminated by a continuous laser sheet with a width of 80mm and a thickness of 1 mm，the interface morphologies after incident shock and reshock impact are captured in a single test run with the help of the high speed camera.The results show that different evolving rates and interface morphologies before and after reshock are observed for SF6and argon gas cylinders due to the different Atwood numbers.For the same reflected end wall，the visible two reversed vortex rings are generated for the SF6gas cylinder after the incident shock passage，which are not apparently observed for the argon gas cylinder because of less vorticity deposition on the interface.Moreover，after the reshock impact，secondary vortex rings which have reversed rotating directions to the original vortex rings are quickly generated in the SF6gas cylinder and dominate the flow field at the later stage，becoming the primary feature of the flow.However，in the argon gas cylinder，the secondary vortex rings，though generated with opposite directions of rotation to the original ones，are not fully developed all the time due to the smaller Atwood number that results in less production of reversed baroclinic vorticity on the argon interface.The present results illustrate the pronounced influence of the Atwood number on the development of the Richtmyer-Meshkov instability.

O354.5

：A

1672-9897（2014）06-0056-05doi：10.11729／syltlx20130077

（編輯：張巧蕓）

2013-09-03；

：2013-10-25

國家自然科學(xué)基金（11272308，11302219）；中央高?；究蒲袠I(yè)務(wù)費專項資金（WK2090050020）和中國博士后科學(xué)基金（BH2090050031）

HeHQ，ZhaiZG，SiT，etal.ExperimentalstudyonthereshockedRMinstabilityoftwokindsofheavygascylinder.JournalofExperimentsinFluidMechanics，2014，28（6）：56-60.何惠琴，翟志剛，司 廷，等.反射激波作用下兩種重氣柱界面不穩(wěn)定性實驗研究.實驗流體力學(xué)，2014，28（6）：56-60.

何惠琴（1991-），女，安徽池州人，碩士研究生。研究方向：實驗流體力學(xué)。通信地址：安徽省合肥市黃山路中國科學(xué)技術(shù)大學(xué)近代力學(xué)系（230027）。E-mail：hhq＠m(xù)ail.ustc.edu.cn