白海濤， 賴煥新（華東理工大學(xué)承壓系統(tǒng)與安全教育部重點(diǎn)實(shí)驗(yàn)室，上海 200237）

基于三種亞格子模型的空腔振蕩流動(dòng)計(jì)算

白海濤， 賴煥新
（華東理工大學(xué)承壓系統(tǒng)與安全教育部重點(diǎn)實(shí)驗(yàn)室，上海 200237）

使用三種亞格子應(yīng)力模型，對(duì)長深比（L/D）為5的三維矩形開式空腔的可壓縮流體進(jìn)行大渦模擬計(jì)算。研究得到的空腔自激振蕩頻率與Rossiter公式計(jì)算結(jié)果和實(shí)驗(yàn)結(jié)果吻合良好，結(jié)果顯示振蕩能量主要集中在較低頻率區(qū)域，壓力幅值主要出現(xiàn)在前三階模態(tài)。Dynamic Smagorinsky-Lilly（DSM）模型在空腔前后壁面附近區(qū)域的脈動(dòng)強(qiáng)度分布比Smagorinsky-Lily （SM）模型更為接近實(shí)驗(yàn)值，Wall Adapting Local Eddy Viscosity（WALE）模型的脈動(dòng)強(qiáng)度分布與實(shí)驗(yàn)值最為接近。由空腔底部監(jiān)測點(diǎn)聲壓級(jí)分布及聲壓頻譜圖可以看出：WALE模型性能最佳，DSM模型結(jié)果也與實(shí)驗(yàn)結(jié)果相符合，SM模型的預(yù)測性能略差。

開式空腔；自激振蕩；大渦模擬；亞格子應(yīng)力模型；氣動(dòng)噪聲

流體流過物體表面的空腔或缺口時(shí)，由于腔外剪切流與腔內(nèi)流動(dòng)的相互作用，會(huì)出現(xiàn)自激振蕩現(xiàn)象，同時(shí)出現(xiàn)劇烈的壓力、速度脈動(dòng)，并輻射產(chǎn)生強(qiáng)烈的噪聲，該物理現(xiàn)象稱為空腔自激振蕩?？涨蛔约ふ袷幀F(xiàn)象廣泛存在于飛行器的起落架艙、武器艙及燃燒室等部位，是典型的聲-渦干涉、非定常流和流體動(dòng)力不穩(wěn)定問題。從20世紀(jì)50年代開始，人們對(duì)空腔自激振蕩流動(dòng)特性做了大量研究。關(guān)于開式空腔自激振蕩物理機(jī)制，雖然有多種解釋，但最被人們接受的是Rossiter［1］提出的空腔流聲共振反饋模型并給出了預(yù)估振蕩頻率的半經(jīng)驗(yàn)公式，該公式在一定精度范圍內(nèi)能夠較為準(zhǔn)確地預(yù)測空腔流激振蕩的峰值頻率，成為評(píng)價(jià)數(shù)值模擬結(jié)果的重要標(biāo)準(zhǔn)。

隨著計(jì)算流體力學(xué)的發(fā)展，空腔聲學(xué)特性的數(shù)值研究得到廣泛開展［2］。湍流直接模擬（DNS）能夠分辨所有能量尺度的渦，但目前局限于低雷諾數(shù)的情況。雷諾平均（RANS）能夠較為準(zhǔn)確地預(yù)測前幾階模態(tài)頻率，但由于平均方法抹掉了湍流脈動(dòng)部分信息，對(duì)于脈動(dòng)強(qiáng)度的預(yù)測存在偏差。大渦模擬（LES）使用網(wǎng)格尺度的空間平均和濾波，可以保留大尺度渦運(yùn)動(dòng)湍流瞬態(tài)信息［3］，是目前被認(rèn)為最有前景的數(shù)值方法。目前被廣泛應(yīng)用的亞格子湍流模型有Smagorinsky-Lily（SM）模型、Dynamic Smagorinsky-Lilly（DSM）模型和Wall Adapting Local Eddy Viscosity（WALE）模型。SM模型［4］適用于各向同性的湍流流動(dòng)計(jì)算，因?yàn)樾问胶唵?，?yīng)用方便而得到廣泛使用。湍流在黏性底層的脈動(dòng)很小，相應(yīng)的雷諾應(yīng)力也很小，而壁面附近速度梯度卻很大。此時(shí)，在壁面附近SM模型產(chǎn)生過大的耗散，尤其在邊界層轉(zhuǎn)捩問題中，往往抑制了轉(zhuǎn)捩的發(fā)生或?qū)е铝鲃?dòng)重新層流化［5］。為了彌補(bǔ)SM模型的缺陷，Germano等［6］提出DSM模型，通過建立隨著時(shí)間與空間變化的亞尺度黏性系數(shù)公式，給出固體壁面上正確漸進(jìn)關(guān)系，并能預(yù)測能量逆向傳遞因而得到更廣泛關(guān)注。Moin等［7］于1991年將DSM模型應(yīng)用于可壓縮湍流。然而，DSM模型在某些流場中會(huì)得到不合理亞格子黏性系數(shù)［8］。Nicoud等［9］于1999年提出WALE模型，該模型采用速度梯度張量的平方，除了變形率張量，還考慮了旋轉(zhuǎn)張量的影響，在近壁面尺度下可以較好得到渦黏性。

雖然前人在空腔自激振蕩方面做了很多工作［10-11］，但迄今為止，關(guān)于不同亞格子模型在預(yù)測空腔自激振蕩方面的差異鮮見報(bào)道。本文采用SM、DSM和WALE三種亞格子模型對(duì)空腔自激振蕩流動(dòng)特性和振蕩機(jī)理進(jìn)行分析，并將模擬計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比，分析這三種亞格子模型的優(yōu)劣，為進(jìn)一步研究空腔自激振蕩發(fā)聲機(jī)理提供參考和借鑒。

1　研究對(duì)象及網(wǎng)格

本文研究對(duì)象是M219三維矩形空腔［12］，長深比L/D=5。因?yàn)樵摽涨坏撞坎粫?huì)出現(xiàn)剪切層再附著情況，按照空腔靜態(tài)流動(dòng)特性［13］，該空腔為開式空腔。自由來流馬赫數(shù)為M=0.85，進(jìn)口滯止壓力p0=100 996 Pa，滯止溫度T0=309.3 K，對(duì)應(yīng)的自由流速度U∞=280.2 m/s，基于空腔長度的雷諾數(shù)ReL=6.8×106。計(jì)算區(qū)域由空腔和空腔上方主流區(qū)域兩部分組成，坐標(biāo)原點(diǎn)在空腔底板中心處，主流區(qū)長（L′）、寬（W′）、深（D′）分別平行于空腔長（L）、寬（W）、深（D）。對(duì)整個(gè)流體區(qū)域劃分結(jié)構(gòu)網(wǎng)格，網(wǎng)格采用分區(qū)劃分和局部加密技術(shù)，離壁面近處稠密，遠(yuǎn)離壁面稀疏［2］。模型尺寸及網(wǎng)格分布情況如表1和圖1所示，表1中Y＋表示第一層網(wǎng)格到壁面距離的量綱為一量，用來考察網(wǎng)格質(zhì)量的優(yōu)劣。

表1　模型尺寸及網(wǎng)格參數(shù)Table 1 Model size and mesh parameters

圖1　計(jì)算區(qū)域及網(wǎng)格分布情況Fig.1 Computational domain and distribution of the grid

在空腔底板Y=－25.4 mm的線上均勻分布P1～P10共10個(gè)監(jiān)測點(diǎn)，點(diǎn)與點(diǎn)之間的距離為50.8 mm，10個(gè)監(jiān)測點(diǎn)位置對(duì)應(yīng)英國國防科技公司QinetiQ對(duì)該空腔進(jìn)行壓力脈動(dòng)測量實(shí)驗(yàn)的監(jiān)測點(diǎn)，各個(gè)監(jiān)測點(diǎn)分布情況如圖2所示。

圖2　監(jiān)測點(diǎn)分布Fig.2 Distribution of monitoring points

2　數(shù)值計(jì)算方法

2.1控制方程與方法

可壓縮流動(dòng)的連續(xù)方程、動(dòng)量方程和能量方程的量綱為一形式如下：

其中：i，j分別表示沿x，y方向分量；ρ，t，u，p分別表示密度，時(shí)間，速度，壓力；σ，ET，q，S分別表示正應(yīng)力，總能量，熱通量和源項(xiàng)；符號(hào)頂部“－”表示普通濾波變量；“～”表示Favre濾波變量。

參考進(jìn)口自由流變量進(jìn)行量綱為一化，參照長度為空腔深度。亞格子應(yīng)力張量τij有如下定義：

采用Boussinesq假定，亞格子應(yīng)力使用式（5）計(jì)算：

其中：μt為亞網(wǎng)格湍流黏性力；τkk為亞網(wǎng)格尺度各向同性的部分；δij為應(yīng)變率；S—ij為應(yīng)力張量的速率，定義為

計(jì)算時(shí)依據(jù)實(shí)驗(yàn)數(shù)據(jù)，進(jìn)口給定總溫T0和總壓p0，出口給出無反射壓力邊界條件，主流區(qū)前、后、上表面采用對(duì)稱邊界條件，其他面為無滑移固體壁面。使用Fluent軟件計(jì)算流場，采用標(biāo)準(zhǔn)k-ε模型計(jì)算得到的穩(wěn)態(tài)流場作為LES的初始場。LES模型中的亞格子湍流模型分別選取SM、DSM 和WALE模型，壓力-速度耦合采用PISO算法，動(dòng)量方程使用二階迎風(fēng)格式離散。時(shí)間推進(jìn)步長為5×10－6s，相當(dāng)于流體從流場流入到流出時(shí)間的1/1 488，計(jì)算總時(shí)長為0.3 s。

2.2三種亞格子（SGS）模型

SM模型是參照雷諾平均模式最早提出的SGS渦黏模型，以各向同性湍流為基礎(chǔ)，認(rèn)為亞網(wǎng)格湍流具有混合長度型渦黏系數(shù)：

1

DSM模型通過兩個(gè)不同的亞尺度應(yīng)變率張量，從已求解的區(qū)域信息計(jì)算得到Cs，Cs使用以下公式計(jì)算：

其中：Lij和Mij為二階張量；k和h為常數(shù)；a為濾波寬度比為格子過濾尺寸，模型中的其他參數(shù)與SM模型計(jì)算方法相同。模型系數(shù)Cs在空間上的變化會(huì)導(dǎo)致計(jì)算不穩(wěn)定，為了保證計(jì)算穩(wěn)定，通常將模型系數(shù)在空間沿均勻流動(dòng)方向進(jìn)行特定的平均處理。

在WALE模型中，渦黏模型：

其中：Cw默認(rèn)為0.325；gij為浮力產(chǎn)生項(xiàng)，其他符號(hào)與SM模型相同。WALE模型針對(duì)近壁面的流動(dòng)，優(yōu)化了亞網(wǎng)格黏性的計(jì)算方式，使預(yù)測結(jié)果更加符合真實(shí)的流動(dòng)。

3　結(jié)果分析

3.1空腔的平均特性

圖3給出了三種亞格子模型XOZ中心對(duì)稱面平均流向速度和平均縱向速度分布情況。

圖3中1.0＜Z/D＜2.0的區(qū)域表示主流區(qū)，三種模型的平均流向速度和平均縱向速度在主流區(qū)均吻合良好，表明了數(shù)值模擬的準(zhǔn)確性。由于剪切流與腔內(nèi)流動(dòng)的相互作用，空腔內(nèi)特別是空腔后壁面附近會(huì)出現(xiàn)劇烈的壓力和速度脈動(dòng)，渦黏模型的過分耗散會(huì)降低局部有效雷諾數(shù)［10］，導(dǎo)致各個(gè)亞格子模型在0＜Z/D＜1.0范圍內(nèi)的平均流向速度和平均縱向速度均有所差異。從總體上看，SM模型平均流速分布與DSM模型和WALE模型差異較大，主要是由于在空腔壁面附近存在強(qiáng)剪切作用，小尺度渦的強(qiáng)度會(huì)得到明顯加強(qiáng)，導(dǎo)致渦的黏性耗散作用比SM模型給定的要強(qiáng)，從而出現(xiàn)了上述差異［14］。在流場復(fù)雜多變的空腔后壁面處，WALE模型的平均縱向速度表現(xiàn)出與SM和DSM模型的差異，是由于上游剪切層產(chǎn)生的渦核與空腔后壁面發(fā)生碰撞，變成碎片并表現(xiàn)出高度的不規(guī)則性，造成空腔后壁面處速度分布不具備嚴(yán)格的周期性［10］，因此，即使采樣時(shí)間足夠長，不同模型的采樣時(shí)間平均值也會(huì)有一定差異。另一方面，本文LES采樣時(shí)間可能還不夠長，也是造成上述差異的原因。

圖3　XOZ中心對(duì)稱面平均流向速度和平均縱向速度Fig.3 Mean longitudinal velocity and mean vertical velocity of the central XOZ plane

圖4所示曲線為圖2中10個(gè)監(jiān)測點(diǎn)所在Y= －25.4 mm線上的聲壓級(jí)（SPL）分布情況，三種亞格子模型的預(yù)測結(jié)果都比實(shí)驗(yàn)值略高，是由于當(dāng)用LES計(jì)算具有強(qiáng)剪切運(yùn)動(dòng)湍流時(shí)，常用的亞格子模型都會(huì)出現(xiàn)平均值偏高的現(xiàn)象［14］。自激振蕩能量從空腔后壁面到前壁面大體上呈現(xiàn)單調(diào)遞減分布，主要是因?yàn)榭涨粌?nèi)的自激振蕩由剪切層在空腔后壁面碰撞產(chǎn)生壓力波，當(dāng)壓力波以渦旋的形式向上游運(yùn)動(dòng)時(shí)，一部分能量由于在空氣中傳播衰減而變?nèi)酰硗庖徊糠帜芰吭谕牧鬟\(yùn)動(dòng)中，由大尺度渦旋向小尺度渦旋逐級(jí)傳遞，直到有分子黏性起顯著作用的尺度，并在該尺度下耗散為熱能。因此，出現(xiàn)了下游自激振蕩能量比上游大的趨勢。

圖4　空腔底部Y=－25.4 mm處聲壓級(jí)分布Fig.4 Sound pressure level distribution of Y=－25.4 mm in the bottom of the cavity

由SM模型與實(shí)驗(yàn)結(jié)果的對(duì)比情況來看，模擬值與實(shí)驗(yàn)值差別最大的地方是在空腔前、后壁面附近區(qū)域，在這兩個(gè)區(qū)域，SM渦黏模型的過分耗散降低了局部雷諾數(shù)而使大尺度渦旋的預(yù)測結(jié)果比實(shí)驗(yàn)值偏大［15］。DSM模型由于動(dòng)態(tài)計(jì)算CS，使得空腔前、后壁面處的聲壓級(jí)得到明顯改善。WALE模型在空腔前半段與DSM模型預(yù)測結(jié)果一致，在后半段更加接近實(shí)驗(yàn)值，明顯優(yōu)于DSM模型，充分表現(xiàn)出其準(zhǔn)確預(yù)測空腔底部聲壓級(jí)分布的優(yōu)越性。

3.2空腔的瞬時(shí)特性

圖5給出了空腔自激振蕩已充分發(fā)展之后XOZ中心面一個(gè)周期（T）不同時(shí)刻的渦量等值線圖。在0.25T，空腔前緣剪切層形成了一個(gè)脫落渦（圖5 （a）），在0.50T～0.75T，脫落渦經(jīng)過空腔中部繼續(xù)向下游發(fā)展，同時(shí)渦旋強(qiáng)度不斷加強(qiáng)（圖5（b）～圖5 （c）），一個(gè)周期結(jié)束時(shí)，第1個(gè)脫落渦與空腔后壁發(fā)生碰撞，碰撞后一部分渦沿著主流區(qū)向下游運(yùn)動(dòng)，另一部分沿著空腔后壁向空腔底部運(yùn)動(dòng)，此時(shí)空腔中下游區(qū)域產(chǎn)生劇烈的自激振蕩現(xiàn)象，同時(shí)，空腔前緣剪切層又形成一個(gè)新的脫落渦（圖5（d））。因此，空腔剪切層形成的具有一定脫落頻率的渦與空腔內(nèi)流場發(fā)生相互作用，產(chǎn)生復(fù)雜的非定常特性。

圖5　空腔內(nèi)瞬時(shí)渦量發(fā)展歷程Fig.5 Instantaneous vorticity development in the cavity

采用Q準(zhǔn)則可以對(duì)流場的渦結(jié)構(gòu)進(jìn)行識(shí)別和顯示，定義為

其中，Dij和Ωij分別為速度梯度的對(duì)稱和反對(duì)稱量，因此Q是速度矢量梯度的第二Halmiton不變量，Q等值面包絡(luò)區(qū)定義了渦核分布情況。圖6給出了使用Q指標(biāo)等值面定義的渦結(jié)構(gòu)?？梢钥闯鲈诳涨磺胺郊罢瓜騾^(qū)域都沒有渦核的存在，由于Kelvin-Helmholtz作用，準(zhǔn)周期性渦核在空腔前緣剪切層處產(chǎn)生并快速向空腔中下游發(fā)展，渦核不斷變大，直至與空腔后壁面發(fā)生碰撞變成碎片并表現(xiàn)出高度的不規(guī)則性。

由圖5中渦量等值線的變化歷程及圖6的渦結(jié)構(gòu)情況，對(duì)空腔自激振蕩機(jī)理分析：剪切層脫落渦與空腔后壁發(fā)生碰撞產(chǎn)生一次壓力波，一次壓力波向空腔前緣傳播，擾動(dòng)剪切層并激發(fā)剪切層更大的不穩(wěn)定性，產(chǎn)生新的脫落渦，形成空腔自激振蕩，完成一次反饋過程。

圖6　Q指標(biāo)定義的瞬時(shí)三維渦結(jié)構(gòu)Fig.6 Instantaneous three-dimensional vortexstructure defined by Q criterion

在本文計(jì)算中，使用的工作站為3.0 GHz主頻、20核心的DellTMPower EdgeTMT7610。計(jì)算表明在同樣網(wǎng)格數(shù)目和相同邊界條件下計(jì)算300個(gè)時(shí)間推進(jìn)步長，SM模型需要40 min，DSM模型和WALE模型均需要41 min。因此可以看出這些模型在計(jì)算效率上沒有明顯差異。

經(jīng)過三種亞格子模型計(jì)算，各個(gè)監(jiān)測點(diǎn)的聲壓史被記錄。由于瞬態(tài)計(jì)算的流場是一個(gè)由不穩(wěn)定狀態(tài)逐步發(fā)展到準(zhǔn)穩(wěn)態(tài)的過程，最終選用經(jīng)過0.1 ～0.3 s達(dá)到準(zhǔn)穩(wěn)態(tài)的215個(gè)數(shù)據(jù)進(jìn)行快速傅里葉變換（FFT）。實(shí)驗(yàn)采樣時(shí)間為本文FFT變換采樣時(shí)間的20倍，因此將實(shí)驗(yàn)數(shù)據(jù)均分為20段，圖7中的實(shí)驗(yàn)值的壓力頻譜是20段數(shù)據(jù)的整體平均。圖7將對(duì)空腔內(nèi)有代表性的P1、P4、P7和P10進(jìn)行重點(diǎn)分析。為清晰表明三種亞格子應(yīng)力模型在預(yù)測空腔自激振蕩方面的差異，在P10點(diǎn)分別列出三種亞格子模型結(jié)果與實(shí)驗(yàn)結(jié)果的對(duì)比，P1，P4和P7點(diǎn)處為綜合對(duì)比。

圖7　聲壓頻譜圖Fig.7 Sound pressure spectrum

三種亞格子模型的模態(tài)頻率和聲壓級(jí)與實(shí)驗(yàn)值吻合良好，都在合理預(yù)測范圍之內(nèi)?？涨蛔约ふ袷幠芰恐饕性?～1 000 Hz的低頻區(qū)域，前三階模態(tài)頻率對(duì)應(yīng)的聲壓級(jí)較大，是空腔自激振蕩的主要振蕩模式。在預(yù)測精準(zhǔn)度上，三種亞格子模型仍有一定差異。表2、表3給出了P10點(diǎn)處三種亞格子模型前四階模態(tài)頻率及對(duì)應(yīng)聲壓級(jí)與Rossiter公式和實(shí)驗(yàn)值的對(duì)比情況，括號(hào)內(nèi)數(shù)值為計(jì)算值與實(shí)驗(yàn)值的相對(duì)誤差。Rossiter公式如下：

其中：fn為空腔振蕩的模態(tài)頻率，Hz；L為空腔長度，m；U∞為自由流速度；M∞為自由流馬赫數(shù)；n為模態(tài)數(shù)；κ為渦速度與自由流速度比相關(guān)的常數(shù)，κ=0.57；γ為渦通過與產(chǎn)生聲壓之間的時(shí)間延遲因子，Lionel Larcheveque認(rèn)為γ是與空腔長寬比相關(guān)的一個(gè)常數(shù)，當(dāng)L/D=5時(shí)，γ=0.29［16］。

由表2和表3可以看出，若以3%相對(duì)誤差為界限作為評(píng)價(jià)誤差大小依據(jù)，SM模型的第一、二和第四階模態(tài)頻率預(yù)測偏差較大，而WALE模型的第二階模態(tài)頻率預(yù)測準(zhǔn)確度較低。同時(shí)，SM模型和DSM模型的第一、四階模態(tài)對(duì)應(yīng)的聲壓級(jí)與實(shí)驗(yàn)值差別較大。上述差異也體現(xiàn)在P1、P4和P7這三個(gè)點(diǎn)上。因此，可以看出，由于SM模型采用單一的Cs常數(shù)，不能準(zhǔn)確模擬有剪切流的固體壁面附近流動(dòng)情況；DSM模型動(dòng)態(tài)計(jì)算Cs常數(shù)的方法可以有效改善這一情況；WALE模型使用了正確的固體壁面流動(dòng)的漸進(jìn)行為，各監(jiān)測點(diǎn)的聲壓頻譜預(yù)測與實(shí)驗(yàn)值最為接近，表明WALE模型和DSM模型在預(yù)測監(jiān)測點(diǎn)頻率聲壓級(jí)分布方面要優(yōu)于SM模型。

表2　P10點(diǎn)前四階模態(tài)頻率Table 2 Frequency of the first four modes at P10

表3　P10點(diǎn)前四階模態(tài)聲壓級(jí)Table 3 SPL of the first four modes at P10

4　結(jié) 論

本文使用SM、DSM和WALE三種亞格子模型對(duì)三維矩形空腔自激振蕩流動(dòng)進(jìn)行大渦數(shù)值模擬，并與實(shí)驗(yàn)結(jié)果對(duì)比分析，得出以下結(jié)論：

（1）三種亞格子模型計(jì)算的前四階模態(tài)頻率與Rossiter公式計(jì)算結(jié)果和實(shí)驗(yàn)數(shù)據(jù)吻合良好，空腔自激振蕩能量主要集中在0～1 000 Hz的低頻區(qū)域，前三階模態(tài)是主要的振蕩模式，空腔自激振蕩并沒有嚴(yán)格的周期性。

（2）從計(jì)算的時(shí)間效率上來看，三種亞格子模型并無顯著差異。由空腔底部監(jiān)測點(diǎn)聲壓級(jí)分布及聲壓頻譜圖可以看出：WALE模型性能最佳，DSM模型也能給出與實(shí)驗(yàn)相符合的結(jié)果，SM模型的預(yù)測性能略差。

（3）由于WALE模型在近壁區(qū)域計(jì)算的優(yōu)越性，相比SM和DSM模型，能更準(zhǔn)確地預(yù)測空腔自激振蕩的模態(tài)頻率和壓力幅值。該結(jié)果也間接說明亞格子模型對(duì)數(shù)值模擬結(jié)果的重要性。

［1］ ROSSITER J E.Wind tunnel experiments of the flow over rectangular cavities at subsonic andtransonic speeds［R］. London：Ministry of Aviation，1967：3438.

［2］ LARCHEVêQUE L，SAGAUT P，LE T H，et al.Large-eddy simulation of a compressible flow in a three-dimensional open cavity at high Reynolds number［J］.Journal of Fluid Mechanics，2004，516：265-301.

［3］ GALPERIN A，ORSZAG S A.Large Eddy Simulation of Complex Engineering and Geophysical Flows［M］.USA：Cambridge University Press，1993.

［4］ SMAGORINSKY J.General circulation experiments with the primitive equations：I.The basic experiment［J］.Monthly Weather Review，1963，91（3）：99-164.

［5］ 鄧小兵.不可壓縮湍流大渦模擬研究［D］.四川綿陽：中國空氣動(dòng)力研究與發(fā)展中心，2008.

［6］ GERMANO M，PIOMELLI U，MOIN P，et al.A dynamic subgrid-scale eddy viscosity model［J］.Physics of Fluids A：Fluid Dynamics，1991，3（7）：1760-1765.

［7］ MOIN P，SQUIRES K，CABOT W，et al.A dynamic subgridscale model for compressible turbulence and scalar transport ［J］.Physics of Fluids A：Fluid Dynamics，1991，3（11）：2746-2757.

［8］ 周磊，解茂昭，賈明，等.不同亞網(wǎng)格尺度應(yīng)力模型在燃油噴霧大渦模擬中的應(yīng)用［J］.內(nèi)燃機(jī)學(xué)報(bào)，2011，29（1）：29-35.

［9］ NICOUD F，DUCROS F.Subgrid-scale stress modelling based on the square of the velocity gradient tensor［J］.Flow，Turbulence and Combustion，1999，62（3）：183-200.

［10］ LAI Huanxin，LUO Kai.A three-dimensional hybrid LES-acoustic analogy method for predicting open-cavity noise［J］. Flow，Turbulence and Combustion，2007，79（1）：55-82.

［11］ XU Lan，CUI Guixiang，WANG Zhishi，et al.High accurate finite volume method for large eddy simulation of complex turbulent flows［J］.International Journal of Turbo and Jet Engines，2006，23（3）：191-210.

［12］ HENSHAW M J.M219 cavity case：Verification and validation data for computational unsteady aerodynamics［R］. ［s.l.］：British Aerospace（Operations）Ltd.，2000：453-472.

［13］ PLENTOVICH E B，STALLINGS R L，TRACY M B. Experimental cavity pressure measurements at subsonic and transonic speeds［J］.NASA Technical Paper，1993，3358：1-128.

［14］ 肖紅林，羅紀(jì)生.大渦模擬中亞格子模型的改進(jìn)及其在槽道湍流中的應(yīng)用［J］.航空動(dòng)力學(xué)報(bào)，2007，22（4）：583-587.

［15］ 賴煥新，周邵萍，羅開紅.空腔的非定?？蓧嚎s過流及相關(guān)氣動(dòng)聲學(xué)問題［J］.工程熱物理學(xué)報(bào)，2007，28（5）：755-758.

［16］ LARCHEVêQUE L，SAGAUT P，MARY I，et al.Largeeddy simulation of a compressible flow past a deep cavity［J］. Physics of Fluids，2003，15（1）：193-210.

Calculation of Oscillation Flow in Open Cavity Using Three Sub-grid Scale Models

BAI Hai-tao， LAI Huan-xin
（Key Laboratory of Pressurized Systems and Safety，Ministry of Education，East China University of Science and Technology，Shanghai 200237，China）

Three sub-grid scale models are employed to calculate the compressible fluid in a three dimensional rectangular open cavity，which has a length-to-depth ratio of 5.The frequencies of sustainedoscillation in the cavity are in agreement with experimental results and the Rossiter formula.The oscillating energy mainly concentrates on the low frequency areas.The first three modes are the basic frequencies of the flow induced oscillation.The main discrepancies of sound pressure level（SPL）occur in the areas near the fore and aft walls.The results of dynamic Smagorinsky-Lilly（DSM）model are more closer to experimental data than that of Smagorinsky-Lily（SM）model.The results of Wall Adapting Local Eddy Viscosity（WALE）model are the closest to the experimental data.The results of sound pressure level and sound pressure spectrum at the bottom of the cavity demonstrate that WALE model is the best one to predict the self-sustained oscillation of the cavity，while DSM model is fairly good，SM model is slightly inaccurate.

open cavity；self-sustained oscillation；large eddy simulation；sub-grid scale model；aerodynamic noise

O353.4

A

1006-3080（2016）01-0125-07 DOI：10.14135/j.cnki.1006-3080.2016.01.020

2015-05-08

國家自然科學(xué)基金（51576067）

白海濤（1989-），男，河南南陽人，碩士生，主要從事氣動(dòng)聲學(xué)的研究。E-mail：bht119@126.com

賴煥新，E-mail：hlai@ecust.edu.cn