楊 玟,王 麗 麗,2,周 海 兵,2,張 樹(shù) 道,2

(1.北 京應(yīng)用物理 與計(jì)算數(shù)學(xué) 研究所,北 京 100094; 2.北京大學(xué)應(yīng)用物理與技術(shù)研究中心,北京100871)

用浮阻力模型研究 Richtmyer-Meshkov不穩(wěn)定性誘導(dǎo)混合*

楊 玟1,王 麗 麗1,2,周 海 兵1,2,張 樹(shù) 道1,2

(1.北 京應(yīng)用物理 與計(jì)算數(shù)學(xué) 研究所,北 京 100094; 2.北京大學(xué)應(yīng)用物理與技術(shù)研究中心,北京100871)

采用浮阻力模型對(duì)激波管低壓縮和激光加載高壓縮情況下的 Richtmyer-Meshkov不穩(wěn)定性誘導(dǎo)混合現(xiàn)象進(jìn)行了研究。通過(guò)與實(shí) 驗(yàn) 和 理 論 分 析 結(jié) 果進(jìn) 行 比 較 發(fā) 現(xiàn):為 了 達(dá)到 好 的 吻 合,Richtmyer-Meshkov不穩(wěn) 定性情況下 阻力系數(shù)的 取值范圍(2.0~5.36)比 Rayleigh-Taylor不 穩(wěn) 定 性 情 況 下 的 值 (3.3~4.0)寬 得多;而在 Richtmyer-Meshkov不穩(wěn)定性情況下,高壓縮時(shí)阻力系數(shù)的不確定度(約 為3.36)明顯高于低壓縮 時(shí)的值(約為1.46),模型的進(jìn)一步完善還有待于更精確實(shí)驗(yàn)的驗(yàn)證。研究顯示:指數(shù)律經(jīng)驗(yàn)公式中指數(shù)隨工況的不同而顯著變化,目前工程設(shè)計(jì)中采用指數(shù)律經(jīng)驗(yàn)公式是粗糙的。

流體力學(xué);Richtmyer-Meshkov不穩(wěn)定性;浮阻力模型;擾動(dòng)界面;壓縮比

當(dāng)沖擊作用于不同密度或不同可壓縮性2種物質(zhì)的擾動(dòng)界面時(shí)產(chǎn)生 Richtmyer-Meshkov不穩(wěn)定性(RMI)。 這 種 不 穩(wěn) 定 性 理 論 上 由 Richtmyer發(fā) 現(xiàn) 并 描 述[1],由 Meshkov 從 實(shí) 驗(yàn) 中 證 實(shí)[2]。 該 界 面 不穩(wěn)定性問(wèn)題在 許多自 然現(xiàn)象 及科學(xué) 和工程 領(lǐng)域中 起著重 要作用[3-8],如 超新星 爆炸、磁 化等離 子 體 、磁 約束、太陽(yáng)磁化層、地下鹽礦、火山島及外殼與內(nèi)部流體混合導(dǎo)致中子收益降低的激光驅(qū)動(dòng)慣性約束聚變和沖擊波與火焰相互作用導(dǎo)致的爆燃轉(zhuǎn)爆轟等。此外,RMI也可能從受沖擊的金屬表面產(chǎn)生噴射物。

RMI的演化通常經(jīng)歷由不穩(wěn)定模式的振幅hk和波長(zhǎng)λ=2π/k描繪的若干階段。對(duì)于khk?1,擾動(dòng)隨k Ut呈正比例增長(zhǎng),U 為激波作用后的界面運(yùn)動(dòng)速度。當(dāng)khk達(dá)到某一值,非線性使增長(zhǎng)率降低,驅(qū)動(dòng)模式耦合,且增長(zhǎng)率隨著擾動(dòng)譜寬的增大而減小。然后,由于尖釘下落(重流體進(jìn)入輕流體)比氣泡上升(輕流體進(jìn)入重流體)快,界面變得不對(duì)稱。對(duì)于寬的不穩(wěn)定譜,非線性最終導(dǎo)致產(chǎn)生湍流混合層。RMI的脈沖性質(zhì)令問(wèn)題復(fù)雜,使得 RMI定性上與常見(jiàn)的 Rayleigh-Taylor不穩(wěn)定性(RTI)不同。由于沖擊的可壓縮性、復(fù)雜的物質(zhì)特性以及后期的非線性運(yùn)動(dòng)直至湍流混合,RMI演化的計(jì)算是困難的。當(dāng)然,隨著計(jì)算機(jī)技術(shù)的迅猛發(fā)展,這可以采用多維高分辨率流體力學(xué)模擬來(lái)進(jìn)行,但它們計(jì)算強(qiáng)度大,無(wú)法用于工程設(shè)計(jì)優(yōu)化研究。因此,目前實(shí)際應(yīng)用中通常采用捕捉較低分辨率時(shí)不穩(wěn)定流動(dòng)主要特征的 簡(jiǎn) 化“混 合 模 型 ”[9]。 楊 玟 等[10-11]對(duì) 此 進(jìn) 行 了 嘗 試 ,將 傳 統(tǒng) 的k-ε模 型 應(yīng) 用 于 界 面 不 穩(wěn) 定 性 引 起 的 混合,取得了令人滿意的結(jié)果。

但是,由于與 RMI相關(guān)的其它物理過(guò)程非常復(fù)雜,較復(fù)雜的混合模型(如k-ε)也難以直接應(yīng)用到工程設(shè)計(jì)中。目前,很多實(shí)際應(yīng)用中對(duì) RMI誘發(fā)混合現(xiàn)象的處理都非常簡(jiǎn)單,假設(shè)混合層寬度以指數(shù)形式tθi增長(zhǎng)。而大量實(shí)驗(yàn)研究表明該比例關(guān)系仍不 確 定[3-6],因 為 考 慮 壓 縮 性 的 計(jì) 算 是 困 難 的,它 們 與 實(shí)驗(yàn)不符。即使指數(shù)律粗略滿足,但不同工況下θi的差別也很大,它顯著依賴于初始擾動(dòng)譜。由此可見(jiàn),工程設(shè)計(jì)中對(duì)RMI誘發(fā)混合現(xiàn)象的處理過(guò)于粗糙。

本文中,在簡(jiǎn)單介紹描述作用于混合層中產(chǎn)生的氣泡和尖釘?shù)母∽枇δＰ突A(chǔ)上,采用該模型對(duì)激波管低壓縮情況和激光加載高壓縮情況下的RMI誘發(fā)混合層寬度(氣泡與尖釘寬度之和)進(jìn)行計(jì)算,驗(yàn)證模型和選取參數(shù)的有效性。

1 模型介紹

目前,典型的 浮阻力 模型可 寫(xiě)為如 下形式[12]:

式 中 :下 角 標(biāo)i,j表 示 2 種 不 同 的 流 體 ,下 角 標(biāo) 為 1 時(shí) 表 示 重 流 體 (尖 釘 ),為 2 時(shí) 表 示 輕 流 體 (氣 泡 );ρi為 重 流 體/輕 流 體 的 密 度 ;vi是 尖 釘/氣 泡 的 滲 透 速 度 ,且vi=d hi/d t,hi表 示 尖 釘/氣 泡 的 瞬 時(shí) 寬 度 ;Ca是 附 加 質(zhì) 量 力 系 數(shù) ;β是 浮 力 產(chǎn) 生 的 模 型 常 數(shù) ;Cd是 阻 力 系 數(shù) ;a(t)為 激 波 脈 沖 加 速 度 ;Vi為 尖 釘/氣 泡的體積,Ai為尖釘/氣泡的截面積。方程左端第一項(xiàng)為慣性力,第二項(xiàng)為浮力,右端為阻力。關(guān)于模型的詳細(xì)論述可參考文獻(xiàn)[13-14],這里不再重復(fù)。對(duì)于 Richtmyer-Meshkov不穩(wěn)定性,通常認(rèn)為沖擊簡(jiǎn)單地給予界面上的氣泡和尖釘一個(gè)脈沖,則它們隨后的運(yùn)動(dòng)可以由慣性力和阻力相等來(lái)得到(加速度為零)。因此脈沖加速度情況是有啟發(fā)性的,可以用來(lái)研究不穩(wěn)定性的慣性特性。

本文所求解的模型方程是一組二階常微分方程,將它們簡(jiǎn)化為一階微分方程:d hi/d t=vi;d vi/d t= -fiCdvi|vi|/hi。 采 用 四 階 Runge-Kutta方 法 進(jìn) 行 求 解 。

2 結(jié)果分析與討論

采用上述模型和數(shù)值方法,對(duì)關(guān)注的激波管低壓縮情況和激光加載高壓縮情況下模型的性能進(jìn)行了考察。這2種工況下RMI產(chǎn)生的機(jī)理不同:對(duì)于弱沖擊,主要貢獻(xiàn)來(lái)自于壓力梯度和密度梯度不重合引起的旋渦沉積;對(duì)于強(qiáng)沖擊,存在激波在經(jīng)折射后產(chǎn)生了顯著的反射,這產(chǎn)生增長(zhǎng)率的振蕩,但它們最終衰減。

2.1 激波管低壓縮的情況下的模擬

首先采用上述模型對(duì)4種不同激波脈沖加速度情況下氣泡和尖釘寬度進(jìn)行了計(jì)算。圖1給出了所采用的4種加速度曲線,g為重力加速度。脈沖加速 度a約 為150g,持 續(xù) 時(shí) 間t0約 為10 ms。 這 些 曲線 為 LANL 的 Dimonte 等 LEM(Linear Electric Motor)實(shí) 驗(yàn) 的 測(cè) 量 曲 線[15]。 實(shí) 驗(yàn) 中 流 體 和 脈 沖 加速 度 的 性 質(zhì) 參 數(shù) 見(jiàn) 表1,其 中R為 密 度 比,R=(1+ A)(1-A),A為 Atwood 數(shù), We為韋伯?dāng)?shù),Re為雷諾數(shù)。對(duì)于每一種情況,通過(guò)調(diào)整 阻 力 系 數(shù)Cd和 初 始 振 幅hi0來(lái) 使 隨 時(shí) 間 變 化 的解與實(shí)驗(yàn)數(shù)據(jù)相符。但是,數(shù)值實(shí)驗(yàn)發(fā)現(xiàn):在大多數(shù)情 況 下hi0對(duì) 結(jié) 果 的 影 響 遠(yuǎn) 小 于Cd的 影 響 。

圖2給出了4種加速度驅(qū)動(dòng)下氣泡和尖釘寬度隨 位 移Z的 變 化,,激 波 作 用 時(shí)Z≈Ut。

圖1 計(jì)算采用的4種不同脈沖加速度曲線Fig.1 Four kinds of impulsive accelerations used in the calculation

由圖可見(jiàn),4種加速度情況下計(jì)算的氣泡和尖釘寬度與實(shí)驗(yàn)基本吻合。計(jì)算中阻力系數(shù)Cd的取值為3.67±0.73,與文獻(xiàn)[16]中分析得到的Cd的不確定度1.2接近。從圖中還可看 出:氣泡和尖釘?shù)?不 對(duì)稱性隨著密度 比R的增大 而增大 。此外 ,本文 中還對(duì) 實(shí)驗(yàn)結(jié) 果按指 數(shù) 律hi=hi0tθi進(jìn) 行 了 擬 合,其 中hi0的 取 值 范 圍 為 0.5~1.0 cm。R=49.1 時(shí) ,θ1≈0.85,θ2≈0.33;R=23.4 時(shí) ,指 數(shù) 迅 速 下 降 ,θ1≈0.45,θ2≈0.24;R=1.57 時(shí) ,θ1≈0.28,θ2≈0.22。 由 此 可 見(jiàn) ,指 數(shù)θi隨 密 度 比 變 化 而 變 化 ,但 具 體 變 化 規(guī) 律還未從數(shù)值模擬和實(shí)驗(yàn)中最終確定,這主要是由于θi對(duì)實(shí)驗(yàn)初始條件敏感,需要計(jì)算和實(shí)驗(yàn)之間更直接的比較。

表1 實(shí)驗(yàn)中采用的流體和脈沖加速度性質(zhì)參數(shù)Table 1 Fluid combinations and characteristics for impusive accerleration experiments

圖2 氣泡和尖釘寬度隨位移的變化Fig.2 The width of bubble and spike with displacement

2.2 激光加載高壓縮情況下的模擬

為了考察模型在高壓縮情況下的性能,我們進(jìn)一步對(duì) Nova激光器上馬赫數(shù) Ma>10的實(shí)驗(yàn)進(jìn)行了模 擬。 實(shí) 驗(yàn) 采 用 一 靶 丸 裝 置 在 Nova 激 光 器 上 進(jìn) 行[17]。 流 體 1 由 厚 度 為 125μm、初 始 密 度 為1.7 g/cm3的鈹燒蝕層組成。流體2是未壓縮密度為0.12 g/cm3的泡沫。波速為46 km/s的入射沖擊與界面相互作用產(chǎn)生反射稀疏波和速度為3 km/s的透射激波。界面經(jīng)加速后速度為56 km/s,物質(zhì)被壓 縮 后 ,ρ1=2 g/cm3,ρ2=0.5 g/cm3,A=-0.6。 這 些 參 數(shù) 通 過(guò) 對(duì) 比 熱 比γ1=1.8 和γ2=1.45 的 流 體求解理想的黎曼問(wèn)題得到。

圖3給出了 Nova實(shí)驗(yàn)中計(jì)算的加速度曲線。由圖可見(jiàn),激光驅(qū)動(dòng)在4 ns后停止,這導(dǎo)致泡沫減壓,由于A<0而產(chǎn)生 Rayleigh-Taylor(RT)分量,因此沖擊壓縮后流動(dòng)是亞音速的,本文模型是適用的。圖4給出了混合區(qū)總寬度H隨位移Z的變化(由于實(shí)驗(yàn)不能分辨氣泡和尖釘,因此給出了總振幅H)。從圖中可看出:混合區(qū)總寬度的計(jì)算值與實(shí)驗(yàn)值吻合,而且Cd=2.0和Cd=5.36的曲線之間包括了全部的實(shí)驗(yàn)數(shù)據(jù)。但是,阻力系數(shù)Cd的不 確 定 度約為 3.36,明 顯 大 于低壓縮情況 的 值(約 為1.46)。此外,擬合得到總的混合寬度以指數(shù)為0.5的指數(shù)律增長(zhǎng),這超過(guò)了激波管低壓縮時(shí)得到的指數(shù),推測(cè)其原因可能是:(1)激光驅(qū)動(dòng)隨時(shí)間減小,使得壓力降低、界面減速,這導(dǎo)致擾動(dòng)膨脹,并引入 RT 分量(因?yàn)锳a>0)。這些影響可能顯著增加推測(cè)的指數(shù);(2)A=0.6時(shí) Nova上的初始擾動(dòng)比激波管 上 的更對(duì)稱,如果指數(shù)對(duì)初始條件敏感,這可能導(dǎo)致不同的指數(shù)。

圖3 Nova實(shí)驗(yàn)中的加速度曲線Fig.3 Acceleration history for Nova experiment

圖4 混合區(qū)寬度隨位移的變化Fig.4 Variaion of total width with displacement

3 結(jié) 論

采用浮阻力模型對(duì)激波管低壓縮和激光加載高壓縮情況下 Richtmyer-Meshkov不穩(wěn)定性誘發(fā)的物質(zhì)滲透邊界的演化過(guò)程進(jìn)行了計(jì)算,計(jì)算結(jié)果與實(shí)驗(yàn)吻合得較好。這表明本研究中模型參數(shù)的選取、方程中現(xiàn)象學(xué)比例因子的添加和模型假設(shè)是合適的。但是由于實(shí)驗(yàn)測(cè)量的局限性,模型中的一些問(wèn)題仍然是突出的,包括阻力項(xiàng)的大小和形式、壓縮的影響、“附加質(zhì)量”的描述等。為了更好地評(píng)估模型,需要一些實(shí)驗(yàn)上的完善。首先,氣泡和尖釘必須單獨(dú)分辨,因?yàn)樗鼈兊谋憩F(xiàn)相當(dāng)不同,尤其在A較大的情況。其次,實(shí)驗(yàn)持續(xù)時(shí)間應(yīng)當(dāng)延長(zhǎng)至足以揭示模型的差別為止。盡管如此,本文模型仍明顯優(yōu)于當(dāng)前實(shí)際應(yīng)用中所采用的經(jīng)驗(yàn)公式(本研究也顯示指數(shù)θi隨工況的不同而顯著變化)。

[1]Richtmyer R D.Taylor instability in shock acceleration of compressible fluids[J].Communicational Pure Applied Mathematics,1960,13(1):297-319.

[2]Meshkov E E.Instability of the interface of two gases accelerated by a shock wave[J].Soviet Fluid Dynamics, 1969,4(1):101-104.

[3]Dimonte G,Remington B.Richtmyer-Meshkov experiments on the Nova laser at high compression[J].Physical Review Letters,1993,70(12):1806-1809.

[4]Dimonte G,Frerking C E,Schneider M.Richtmyer-Meshkov instability in the turbulent regime[J].Physical Review Letters,1995,74(24):4855-4858.

[5]Vetter M,Sturtevant B.Experiments on the Richtmyer-Meshkov instability of an air/SF6interface[J].Shock Waves,1995,4(5):247-252.

[6]Jourdan G,Houas L,Haas J F,et al.Thickness and volume measurements of a Richtmyer-Meshkov instability induced mixing zone in a square shock tube[J].Journal of Fluid Mechanics,1997,349:67-94.

[7]Alon U,Hecht J,Ofer D,et al.Power laws and similarity of Rayleigh-Taylor and Richtmyer-Meshkov mixing fronts at all density ratios[J].Physical Review Letters,1995,74(4):534-537.

[8]Aglitskiy Y,Velikovich A L,Karasik M,et al.Basic hydrodynamics of Richtmyer-Meshkov-type growth and os-cillations in the inertial confinement fusion-relevant conditions[J].Philosophical Transactions of Royal Society A, 2010,368(1916):1739-1768.

[9]Llor A.Statistical hydrodynamic models for developed mixing instability flows[S].Springer,2005.

[10]楊 玟,王 麗麗 ,張樹(shù) 道.Rayleigh-Taylor不穩(wěn)定性誘 導(dǎo)湍流混合 的數(shù)值模擬[J].工程力 學(xué),2011,28(6):236-241. Yang Min,Wang Li-li,Zhang Shu-dao.Numerical simulation of turbulent mixing induced by Rayleigh-Taylor instability[J].Engineering Mechanics,2011,28(6):236-241.

[11]楊玟,王麗麗,張樹(shù)道,等.用湍流模型 研 究 Richtmyer-Meshkov不穩(wěn)定性 誘 導(dǎo)的湍流 混 合[J].空氣動(dòng) 力 學(xué)學(xué)報(bào), 2010,28(1):119-123. Yang Min,Wang Li-li,Zhang Shu-dao,et al.The study of turbulent mixing induced by Richtmyer-Meshkov instability using turbulence model[J].Acta Aerodynamica Sinica,2010,28(1):119-123.

[12]Cheng B.Modeling chaotic mixing[J].Nuclear Weapons Journal,2010,1:8-17.

[13]Layzer D.On the gravitational instability of two superposed fluids in a gravitational field[J].Astrophysics Journal,1955,122(1):1-12.

[14]楊玟,王麗 麗,周 海兵 ,張樹(shù) 道.用浮阻 力模型研究 Rayleigh-Taylor不穩(wěn)定性誘 發(fā)混合現(xiàn)象[J].工程力 學(xué),2013,30 (4):385-391. Yang Min,Wang Li-li,Zhang Shu-dao.Study on mixing induced by Rayleigh-Taylor instability using buoyancydrag model[J].Engineering Mechanics,2013,30(4):385-391.

[15]Dimonte G,Schneider M.Density ratio dependence of Rayleigh-Taylor mixing for sustained and impulsive acceleration histories[J].Physics of Fluids,2000,12(2):304-321.

[16]Dimonte G,Schneider M.Turbulent Richtmyer-Meshkov instability experiments with strong radiatively driven shocks[J].Physics of Plasmas,1997,4(12):4347-4357.

[17]Dimonte G.Nonlinear evolution of the Rayleigh-Taylor and Richtmyer-Meshkov instabilities[J].Physics of Plasmas,1999,6(5):2009-2015.

Study on mixing induced by Richtmyer-Meshkov instability by using buoyancy-drag model

Yang Min1,Wang Li-li1,2,Zhou Hai-bing1,2,Zhang Shu-dao1,2
(1.Insitute of Applied Physics and Computational Mathematics,Beijing 100094,China; 2.Center for Applid Physics and Tchnology,Peking University,Beijing 100871,China)

The mixing induced by Richtmyer-Meshkov instability under lower and high compression ratios is studied by using buoyancy-drag model.It is found by the comparison between the experimental and theoretical analytic results that in order to achieve good agreement,the range of the drag coefficient value chosen in Richtmyer-Meshkov instability is much wider than that in Rayleigh-Taylor instability;while in Richtmyer-Meshkov instability the uncertainty of the drag coefficient under high compression is larger than that under lower compression.It is pointed that the further improvement of the model needs to be validated by more accurate experiment.Moreover the study reveals that the exponent in empirical expression varies with calculation conditions remarkably and the use of the empirical expression in current engineering design is rough.

fluid mechanics;Richtmyer-Meshkov instability;buoyancy-drag model;disturbed interface;compression ratio

O358國(guó)標(biāo)學(xué)科代碼:13025

:A

10.11883/1001-1455(2015)03-0423-05

(責(zé)任編輯 王小飛)

2013-11-28;

2014-06-28

國(guó)家 自然科學(xué)基 金項(xiàng)目(11072040);裝 備預(yù)研重點(diǎn) 實(shí)驗(yàn)室基金 項(xiàng)目(9140C690204120C69259);中國(guó) 工程物理研 究院科學(xué)技 術(shù)發(fā)展基金 項(xiàng)目(2012B0201028、2012B0201030)

楊 玟(1970- ),女,博士,副研究員,yang_min@iapcm.ac.cn。