韓帥斌，羅勇，張樹(shù)海

(1.中國(guó)空氣動(dòng)力研究與發(fā)展中心 空氣動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，綿陽(yáng) 621000) (2.中國(guó)空氣動(dòng)力研究與發(fā)展中心 計(jì)算空氣動(dòng)力研究所，綿陽(yáng) 621000)

0 引 言

空腔流動(dòng)是流體力學(xué)的經(jīng)典流動(dòng)問(wèn)題，包含豐富的渦結(jié)構(gòu)相關(guān)的物理現(xiàn)象，例如剪切層中渦的卷起與撞擊、渦與剪切層相互作用、渦聲的產(chǎn)生與傳播等。同時(shí)空腔流動(dòng)具有較強(qiáng)的實(shí)際應(yīng)用背景。在航空航天領(lǐng)域，剪切層內(nèi)的渦與剪切層及空腔拐角相互作用產(chǎn)生的強(qiáng)烈噪聲不僅影響身體健康，更會(huì)導(dǎo)致飛行器結(jié)構(gòu)疲勞，危及飛行安全，是目前迫切需要解決的問(wèn)題。自20世紀(jì)50年代起，空腔流動(dòng)尤其是空腔流致噪聲得到了大量研究[1-3]?？涨涣髦略肼暸c來(lái)流參數(shù)、空腔幾何形態(tài)等密切相關(guān)。J.E.Rossiter[4]通過(guò)系列實(shí)驗(yàn)總結(jié)得到Rossiter半經(jīng)驗(yàn)公式，能夠很好地預(yù)測(cè)輻射噪聲的振蕩頻率與來(lái)流馬赫數(shù)的關(guān)系。通過(guò)H.H.Heller等[5]、C.K.W.Tam等[6]的完善，公式預(yù)測(cè)范圍和精度得到進(jìn)一步提高。空腔流致噪聲的動(dòng)力學(xué)過(guò)程為：渦擾動(dòng)在剪切層中不斷增長(zhǎng)，渦撞擊空腔后緣產(chǎn)生聲波，聲波擾動(dòng)向上游傳播，聲波擾動(dòng)到達(dá)前緣誘導(dǎo)產(chǎn)生新的渦擾動(dòng)[3]。該動(dòng)力學(xué)過(guò)程的一個(gè)關(guān)鍵問(wèn)題是渦量擾動(dòng)和壓力脈動(dòng)之間如何相互轉(zhuǎn)化。M.V.Morkovin等[7]研究了壓力信號(hào)轉(zhuǎn)化為渦量擾動(dòng)的物理機(jī)制；Y.P.Tang等[8]對(duì)渦-邊緣相互作用的形式進(jìn)行了區(qū)分，對(duì)渦波轉(zhuǎn)化壓力波的物理機(jī)理作出了相應(yīng)的闡釋?zhuān)蝗f(wàn)振華[9]借助渦量及擬渦能的時(shí)空演化分析了空腔中渦結(jié)構(gòu)撞擊后拐角產(chǎn)生壓力脈動(dòng)的動(dòng)力學(xué)過(guò)程。空腔中的渦結(jié)構(gòu)作為噪聲的主要來(lái)源，其動(dòng)力學(xué)過(guò)程對(duì)空腔流動(dòng)特性及噪聲的產(chǎn)生和傳播起著重要作用，開(kāi)展空腔流動(dòng)中的渦動(dòng)力學(xué)分析能夠?yàn)橄鄳?yīng)的流動(dòng)控制和降噪提供理論依據(jù)。

識(shí)別流場(chǎng)中的渦結(jié)構(gòu)是開(kāi)展渦動(dòng)力學(xué)分析的關(guān)鍵步驟。目前存在諸多渦判據(jù)[10-11]可用來(lái)識(shí)別流場(chǎng)中的渦結(jié)構(gòu)，但仍沒(méi)有公認(rèn)的渦的嚴(yán)格定義。既往的渦動(dòng)力學(xué)研究大多采用歐拉框架下的分析工具例如閉合流線、渦量、Q判據(jù)、 準(zhǔn)則等識(shí)別流場(chǎng)中的渦，然而這些方法通常僅使用某瞬時(shí)流場(chǎng)的速度或速度梯度提取渦結(jié)構(gòu)，揭示的是流動(dòng)的瞬時(shí)狀態(tài)，缺失了動(dòng)力學(xué)系統(tǒng)中與時(shí)間相關(guān)的信息，無(wú)法準(zhǔn)確反映流動(dòng)結(jié)構(gòu)的歷史累積效應(yīng)以及相應(yīng)的動(dòng)力學(xué)特性。另外歐拉框架下的判據(jù)通常不具有客觀性，即流場(chǎng)結(jié)果會(huì)隨參考系改變而改變。

近年來(lái)，拉格朗日方法在復(fù)雜流動(dòng)系統(tǒng)的動(dòng)力學(xué)分析及渦識(shí)別方面取得了較大進(jìn)展[12]。拉格朗日方法在時(shí)空相空間內(nèi)跟蹤流體粒子運(yùn)動(dòng)，能夠客觀地揭示非線性動(dòng)力系統(tǒng)的動(dòng)力學(xué)特性和流動(dòng)機(jī)理?；诶窭嗜諗M序結(jié)構(gòu)(Lagrangian Coherent Structures，簡(jiǎn)稱(chēng)LCS)識(shí)別流場(chǎng)中的渦結(jié)構(gòu)并研究其時(shí)空演化規(guī)律，能夠揭示開(kāi)式空腔中的拉格朗日渦的動(dòng)力學(xué)特性。

本文采用拉格朗日擬序結(jié)構(gòu)，并結(jié)合Q判據(jù)和渦量通量分析，對(duì)開(kāi)式空腔流動(dòng)中渦的生成、脫落、對(duì)流以及撞擊破裂過(guò)程展開(kāi)研究，分析其動(dòng)力學(xué)特性，為流動(dòng)控制及降噪提供理論依據(jù)。

1 數(shù)值分析方法

1.1 高精度流場(chǎng)計(jì)算方法

通過(guò)求解非定?？蓧嚎sNavier-Stokes方程(簡(jiǎn)記為N-S方程)對(duì)空腔流動(dòng)進(jìn)行直接數(shù)值模擬，守恒形式的無(wú)量綱N-S方程為

?tρ+?i(ρvi)=0

(1)

(2)

(3)

其中，

(4)

(5)

(6)

粘性系數(shù)的Sutherland公式為

(7)

N-S方程的對(duì)流項(xiàng)和粘性項(xiàng)分別采用五階WENO格式和六階中心差分格式進(jìn)行離散。時(shí)間項(xiàng)采用三階TVD Runge-Kutta格式求解獲得高精度非定常無(wú)量綱流場(chǎng)數(shù)據(jù)。

1.2 拉格朗日擬序結(jié)構(gòu)

拉格朗日擬序結(jié)構(gòu)是基于追蹤流體粒子運(yùn)動(dòng)而提取的流場(chǎng)中主導(dǎo)流動(dòng)的吸引性、排斥性或者強(qiáng)剪切擬序結(jié)構(gòu)，是流動(dòng)的“骨架”[12]，目前已被廣泛應(yīng)用于流動(dòng)分離[13]、渦動(dòng)力學(xué)[14]、生物流體力學(xué)[15]等。諸多方法可以用來(lái)提取流場(chǎng)中的LCS，其中有限時(shí)間李雅普諾夫指數(shù)(Finite Time Lyapunov Exponent，簡(jiǎn)稱(chēng)FTLE)是被廣泛應(yīng)用的一種方法。對(duì)于在非自治系統(tǒng)

(8)

(9)

式中：λmax為矩陣的最大特征值。

在實(shí)際計(jì)算中，F(xiàn)TLE既可以前向時(shí)間積分也可以逆向時(shí)間積分，對(duì)應(yīng)的FTLE的嵴分別捕捉流場(chǎng)中的前向LCS(pLCS)和逆向LCS(nLCS)，pLCS通常是流場(chǎng)中的排斥性結(jié)構(gòu)，nLCS則通常是流場(chǎng)中的吸引性結(jié)構(gòu)。H.Yangzi等[16]采用pLCS和nLCS的邊界識(shí)別拉格朗日渦邊界并追蹤渦運(yùn)動(dòng)，本文采用該方法對(duì)空腔中的拉格朗日渦結(jié)構(gòu)進(jìn)行動(dòng)力學(xué)分析。

1.3 Q判據(jù)

歐拉框架下常用來(lái)識(shí)別渦的一種方法是J.C.R.Hunt等[17]提出的Q判據(jù)。速度梯度張量可以分解為對(duì)稱(chēng)的應(yīng)變率張量S和反對(duì)稱(chēng)的渦張量Ω，即

(10)

Q值的定義即為

(11)

在流場(chǎng)中Q值為正的地方，也即渦張量對(duì)流體變形的貢獻(xiàn)大于應(yīng)變率張量貢獻(xiàn)的區(qū)域被識(shí)別為渦。Q判據(jù)具有伽利略不變性，但在旋轉(zhuǎn)或加速參考系中則不能正確進(jìn)行渦識(shí)別。對(duì)于本文中的開(kāi)式空腔流動(dòng)，參考系不存在旋轉(zhuǎn)或加速，Q判據(jù)能夠給出渦的正確描述，因而被用來(lái)與LCS識(shí)別的拉格朗日渦對(duì)比驗(yàn)證，并描述相應(yīng)的渦動(dòng)力學(xué)。

2 數(shù)值計(jì)算及驗(yàn)證

2.1 計(jì)算模型及設(shè)置

本文直接數(shù)值模擬的空腔流動(dòng)計(jì)算域及空腔附近網(wǎng)格如圖1所示，其中空腔長(zhǎng)深比為2∶1，各項(xiàng)無(wú)量綱長(zhǎng)度參數(shù)為

L=2D,Ll=7D,Lr=20D,Ly=20D

(12)

腔外網(wǎng)格為810×300，腔內(nèi)網(wǎng)格為240×120，并在腔壁附近加密以捕捉邊界層。空腔流動(dòng)會(huì)產(chǎn)生強(qiáng)烈的聲波，為避免聲波在邊界的反射污染流場(chǎng)，在計(jì)算域邊界設(shè)置海綿層[18]吸收聲波，從而獲得準(zhǔn)確的高精度流場(chǎng)數(shù)據(jù)。海綿層的長(zhǎng)度參數(shù)為

Li=3D,Lo=10D,Lt=10D

(13)

本文針對(duì)亞聲速流動(dòng)開(kāi)展研究，分析所用的流場(chǎng)物理量均為無(wú)量綱量。空腔的流動(dòng)參數(shù)設(shè)置為Ma=0.8，Re=2 500，Pr=0.7，其中雷諾數(shù)基于空腔深度D。

(a) 二維空腔流動(dòng)計(jì)算區(qū)域示意圖

(b) 腔體附近網(wǎng)格分布圖1 二維空腔流動(dòng)計(jì)算域及網(wǎng)格Fig.1 Computational domain and mesh of the open cavity

2.2 計(jì)算結(jié)果及驗(yàn)證

本文直接數(shù)值模擬所用代碼已在參考文獻(xiàn)[13，19]中得到很好驗(yàn)證。針對(duì)本算例，首先進(jìn)行網(wǎng)格收斂性驗(yàn)證?？紤]一套加密網(wǎng)格，其中腔外網(wǎng)格為1 245×470，腔內(nèi)網(wǎng)格為400×200，并比較兩套網(wǎng)格下空腔底部(0.5D，0)和空腔外部(0,2D)處的壓力隨時(shí)間變化。兩套網(wǎng)格下壓力變化基本一致，如圖2所示，可以看出：幅值和相位差別均很小，說(shuō)明本計(jì)算所用網(wǎng)格密度已經(jīng)足夠。

(a) 空腔底部x=0.5D處壓力隨時(shí)間變化

(b) 空腔外x=0，y=2D處壓力隨時(shí)間變化圖2 網(wǎng)格收斂性驗(yàn)證Fig.2 Validation of grid convergency

對(duì)空腔(0.995D，D)處的速度和壓力脈動(dòng)采樣分析，結(jié)果如圖3所示，流動(dòng)呈現(xiàn)出強(qiáng)烈的周期性，且無(wú)量綱時(shí)間周期為T(mén)=3.75。對(duì)包含5個(gè)流動(dòng)周期的468個(gè)流場(chǎng)片段進(jìn)行快速傅里葉變換(FFT)，采樣點(diǎn)的頻譜如圖4所示。

圖3 速度及壓力時(shí)間演化Fig.3 Time evolution of velocity and pressure

圖4 空腔壓力及速度脈動(dòng)頻譜分布Fig.4 Frequency spectrum of velocity and pressure perturbation in the open cavity

從圖4可以看出：主頻St=0.668。

Rossiter給出的半經(jīng)驗(yàn)公式[20]

(14)

可見(jiàn)，圖4得到的主頻與式(14)中的第2模態(tài)即St2=0.686吻合較好，而且與K.Krishnamurty[21]的實(shí)驗(yàn)結(jié)果頻率St=0.656吻合良好。Krishnamurty的實(shí)驗(yàn)紋影與直接數(shù)值模擬獲得的數(shù)值紋影對(duì)比如圖5所示。

(a) Krishnamurty實(shí)驗(yàn)紋影

(b) 數(shù)值紋影圖5 實(shí)驗(yàn)紋影圖與DNS計(jì)算所得數(shù)值紋影對(duì)比Fig.5 Comparison between experimental Schlieren and Numerical Schlieren from DNS

從圖5可以看出：兩者在定性上吻合良好，波系結(jié)構(gòu)較為一致。

3 渦動(dòng)力學(xué)分析

3.1 流場(chǎng)基本結(jié)構(gòu)

一個(gè)周期內(nèi)的流線結(jié)構(gòu)如圖6所示，空腔內(nèi)部主要有四個(gè)回流區(qū)，它們構(gòu)成了流場(chǎng)的主要拓?fù)浣Y(jié)構(gòu)。在t=0.2T時(shí)刻，空腔尾緣點(diǎn)處渦量等值面撕裂，在t=0.8T時(shí)刻在前緣開(kāi)始形成閉合流線結(jié)構(gòu)，這些特征在歐拉框架下通常被看作渦破裂和渦生成。然而流線和渦量都不具有客觀性，所揭示的流場(chǎng)結(jié)果會(huì)隨參考系的變化而變化。另外瞬時(shí)流態(tài)不能揭示渦結(jié)構(gòu)的動(dòng)態(tài)演化過(guò)程，其動(dòng)力學(xué)特性也無(wú)法揭示。因此接下來(lái)本文采用拉格朗日方法對(duì)流場(chǎng)渦結(jié)構(gòu)進(jìn)行動(dòng)力學(xué)特性分析。

(a) t=0

(b) t=0.2T

(f) t=T圖6 空腔流動(dòng)的流線及渦量分布Fig.6 Streamlines and vorticity distribution in the cavity

3.2 拉格朗日渦動(dòng)力學(xué)分析

對(duì)空腔流動(dòng)的流場(chǎng)進(jìn)行FTLE的計(jì)算，積分時(shí)長(zhǎng)是影響計(jì)算結(jié)果的重要參數(shù)，通常隨著積分時(shí)長(zhǎng)的增加，F(xiàn)TLE的嵴所提取的LCS結(jié)構(gòu)更加銳利。因此在本文計(jì)算中，積分時(shí)長(zhǎng)選取為一個(gè)時(shí)間周期T=3.75，使得FTLE的嵴能夠準(zhǔn)確捕捉到流場(chǎng)中的LCS。一個(gè)周期內(nèi)不同時(shí)刻的pLCS(紅色)和nLCS(藍(lán)色)如圖7所示，pLCS為排斥性結(jié)構(gòu)，附近的流體粒子被其排斥從而流向動(dòng)力學(xué)特性不同的區(qū)域；nLCS則為吸引性結(jié)構(gòu)，在時(shí)空演化過(guò)程中吸引周?chē)鷣?lái)自不同動(dòng)力學(xué)特性區(qū)域的流體粒子向其聚集靠攏，因此pLCS和nLCS是不同動(dòng)力學(xué)特性區(qū)域的邊界。渦結(jié)構(gòu)作為流場(chǎng)中的擬序結(jié)構(gòu)，在渦邊界內(nèi)外，流動(dòng)的特性存在差異，因而渦邊界可用LCS進(jìn)行識(shí)別。將pLCS和nLCS包絡(luò)的閉合區(qū)域作為渦結(jié)構(gòu)，兩者交叉處的鞍點(diǎn)作為追蹤渦結(jié)構(gòu)運(yùn)動(dòng)的特征點(diǎn)來(lái)分析拉格朗日渦的動(dòng)力學(xué)特性。

pLCS和nLCS包絡(luò)區(qū)域形成的渦結(jié)構(gòu)主要有三個(gè)：空腔前緣點(diǎn)剪切層卷起后形成的渦，空腔后緣點(diǎn)處即將發(fā)生撞擊和撕裂的渦以及空腔右側(cè)的主渦。這與流線所揭示的四個(gè)主要回流區(qū)有所不同，位于空腔右下拐角處的回流區(qū)Ⅱ和位于空腔左側(cè)的回流區(qū)Ⅲ未被LCS揭示出來(lái)。由于LCS揭示的是擬序結(jié)構(gòu)，能在一定生命周期內(nèi)保持動(dòng)力學(xué)及結(jié)構(gòu)特征穩(wěn)定，而回流區(qū)Ⅲ即使從流線觀察不能保持穩(wěn)定的形態(tài)且無(wú)清晰邊界，因此未被LCS識(shí)別出來(lái)；回流區(qū)Ⅱ位于空腔右側(cè)底部拐角處，速度較低，其靠近壁面的邊界由于無(wú)滑移邊界速度幾乎為零，吸引性或排斥性較弱，因而只有與主渦相交的邊界被LCS揭示出來(lái)，而完整的渦結(jié)構(gòu)未被識(shí)別出來(lái)。流場(chǎng)中Q值大于零的區(qū)域如圖7中灰黑色區(qū)域所示。Q值揭示的渦區(qū)域也主要有三個(gè)，即空腔右側(cè)的主渦和剪切層中卷起和向下游發(fā)展的兩個(gè)渦，與LCS包絡(luò)所揭示的渦結(jié)構(gòu)覆蓋的區(qū)域基本一致?？涨挥蚁鹿战翘幍慕菧u盡管被Q值所揭示，但強(qiáng)度較弱，且形態(tài)和邊界不夠清晰穩(wěn)定。流線揭示的回流區(qū)Ⅲ同樣未被Q所識(shí)別如圖6所示，表明空腔左側(cè)盡管存在瞬時(shí)封閉流線，但不能保持動(dòng)力學(xué)穩(wěn)定性及清晰邊界，該回流區(qū)不能識(shí)別為渦結(jié)構(gòu)。綜上，LCS與Q值都清晰地揭示了空腔流動(dòng)的流場(chǎng)中存在的主要渦結(jié)構(gòu)，而空腔右側(cè)的角渦則由于吸引性或排斥性較弱，未被LCS識(shí)別出來(lái)。

空腔右側(cè)的主渦在整個(gè)流動(dòng)周期中形態(tài)基本維持不變，繞渦心順時(shí)針旋轉(zhuǎn)，周期與流場(chǎng)周期一致。在運(yùn)動(dòng)過(guò)程中，空腔后緣點(diǎn)處剪切層中的渦撕裂后，一部分流體被nLCS及pLCS包絡(luò)并隨主渦運(yùn)動(dòng)進(jìn)入到主渦中；同時(shí)主渦靠近剪切層一側(cè)流體被pLCS吸引靠近剪切層進(jìn)入主流，隨主流對(duì)流至下游。

剪切層中的渦結(jié)構(gòu)則經(jīng)歷了復(fù)雜的生成、增大、脫落、對(duì)流、破裂等過(guò)程，其完整周期為流場(chǎng)中壓力和速度脈動(dòng)周期的兩倍，而且剪切層內(nèi)兩個(gè)渦結(jié)構(gòu)的完整運(yùn)動(dòng)動(dòng)力學(xué)過(guò)程和特性是相同的。在圖7中，t=0.4T時(shí)刻如圖7(c)所示，剪切層中的渦結(jié)構(gòu)在空腔前緣點(diǎn)開(kāi)始生成并逐漸增大直至t=T時(shí)刻如圖7(c)～圖7(f)所示，從前緣點(diǎn)脫落進(jìn)入到剪切層中，并隨主流一起向下游對(duì)流。在對(duì)流過(guò)程中如圖7(a)～圖7(e)所示，渦結(jié)構(gòu)發(fā)生形變扭曲，并在t=T+0.8T時(shí)刻如圖7(e)所示，開(kāi)始撞擊空腔后緣點(diǎn)，整個(gè)撞擊過(guò)程中如圖7(f)、圖7(a)、圖7(b)、圖7(c)所示，渦結(jié)構(gòu)不斷被空腔后緣點(diǎn)撕裂，最終在t=2T+0.4T如圖7(c)所示，被撕裂為兩部分，一部分隨主流在空腔外向下游流動(dòng)，另一部分則進(jìn)入主渦部分隨主渦在空腔內(nèi)部流動(dòng)。

(a) t=0

(b) t=0.2T

(c) t=0.4T

(d) t=0.6T

(e) t=0.8T

(f) t=T圖7 LCS及Q值分布Fig.7 LCS and Q value distribution in the cavity

前緣點(diǎn)和后緣點(diǎn)附近區(qū)域的渦通量為

(15)

式中：ω為渦量。

剪切層中渦的生成、脫落、撞擊和破裂時(shí)刻的量化分析通過(guò)式(15)的時(shí)間演化特性分析給出，分析結(jié)果與LCS分析給出的特征時(shí)刻相互對(duì)照驗(yàn)證。選取的前后緣點(diǎn)附近區(qū)域分別如圖7中A、B所示，兩個(gè)時(shí)間周期內(nèi)的渦通量時(shí)間演化如圖8所示，注意在兩個(gè)區(qū)域內(nèi)渦通量都是負(fù)值。由A區(qū)域內(nèi)的渦通量時(shí)間演化可知：在t=0.35T時(shí)刻，渦通量開(kāi)始逐漸積累增加，意味著渦結(jié)構(gòu)開(kāi)始形成；在t=0.87T時(shí)刻渦通量達(dá)到極大值并開(kāi)始減小，意味著此時(shí)渦結(jié)構(gòu)開(kāi)始從前緣點(diǎn)脫落。由B區(qū)域內(nèi)的渦通量時(shí)間演化可知：在t=0.76T時(shí)刻，渦通量開(kāi)始迅速增加并在t=1.12T時(shí)刻達(dá)到極大值，意味著從上游而來(lái)的渦結(jié)構(gòu)開(kāi)始對(duì)流至后緣點(diǎn)處并撞擊后緣點(diǎn)，隨后渦通量開(kāi)始減小，直至t=1.43T時(shí)刻渦通量減小至極小值，此時(shí)渦結(jié)構(gòu)撕裂成為兩部分，一部分隨主流對(duì)流而下，另一部分則進(jìn)入空腔中的主渦內(nèi)。以上基于渦通量的渦結(jié)構(gòu)動(dòng)力學(xué)特性分析與LCS分析結(jié)果基本一致，印證了基于LCS分析渦動(dòng)力學(xué)的準(zhǔn)確性。

圖8 渦通量時(shí)間演化Fig.8 Time evolution of vorticity flux

4 結(jié) 論

(1) 吸引性和排斥性拉格朗日擬序結(jié)構(gòu)邊界的包絡(luò)作為空腔流動(dòng)中的渦邊界，可有效識(shí)別出空腔中形態(tài)基本維持不變的主渦和剪切層內(nèi)發(fā)生復(fù)雜動(dòng)力學(xué)過(guò)程的渦。空腔右側(cè)的角渦則由于吸引性和排斥性較弱，未被識(shí)別出來(lái)。

(2)Q判據(jù)以及渦量通量的時(shí)間演化分析表明，在流動(dòng)周期內(nèi)，拉格朗日擬序結(jié)構(gòu)在空腔前緣點(diǎn)的卷起與脫離伴隨著前緣渦的生成與脫落，并向下游發(fā)展；拉格朗日擬序結(jié)構(gòu)在后緣點(diǎn)的撞擊破裂伴隨著Q與渦的破裂。

(3) 剪切層中的渦結(jié)構(gòu)經(jīng)歷了復(fù)雜的生成、增大、脫落、對(duì)流、破裂等過(guò)程，空腔內(nèi)部的主渦則在整個(gè)流動(dòng)周期內(nèi)都維持相對(duì)穩(wěn)定形態(tài)。