高速飛行器空腔脈動壓力主動控制與非線性數(shù)值模擬

王一丁， 郭亮， 童明波， 張杰

1. 南京航空航天大學(xué) 航空宇航學(xué)院， 南京 210016

高速飛行器空腔脈動壓力主動控制與非線性數(shù)值模擬

王一丁1,2， 郭亮1,2， 童明波1，*， 張杰3

1. 南京航空航天大學(xué) 航空宇航學(xué)院， 南京 210016

2. 成都飛機設(shè)計研究所， 成都 610091

3. 北京銀景科技有限公司， 北京 100107

空腔脈動壓力(空腔噪聲)預(yù)測是高速飛行器內(nèi)埋彈艙的關(guān)鍵技術(shù)之一。非線性噪聲求解方法是近年來新提出的一種噪聲求解方法，為研究該方法對空腔噪聲的預(yù)測性能，將雷諾平均Navier-Stokes(RANS)方程與之相結(jié)合。首先，通過RANS求解空腔周圍流場，得到初始湍流統(tǒng)計平均解，其中包含平均流場基本特征及強制設(shè)定的湍流脈動的統(tǒng)計描述。然后，采用非線性噪聲求解方法重構(gòu)噪聲源并高精度模擬壓力脈動的傳播，計算了馬赫數(shù)Ma=1.5和Ma=5條件下的空腔噪聲。結(jié)果表明，噪聲特性計算值與試驗結(jié)果基本吻合，說明非線性噪聲求解方法對于高速空腔流動噪聲具有較好的預(yù)測能力。在此基礎(chǔ)上，研究了馬赫數(shù)Ma=1.5和Ma=5條件下在空腔前緣加入氣簾噴流主動控制措施對噪聲的抑制作用，并得出在超聲速和高超聲速條件下，氣簾噴流對于空腔脈動壓力都有較好的抑制作用。

脈動； 非線性； 空腔； 湍流； 主動控制

空腔結(jié)構(gòu)在航空航天領(lǐng)域有著廣泛的應(yīng)用。作為飛行器上最常用的空腔結(jié)構(gòu)，內(nèi)埋彈艙有效地提高了飛行器的低可探測性和戰(zhàn)場生存率，同時避免了武器外掛所產(chǎn)生的阻力，優(yōu)化了飛行性能，是新一代戰(zhàn)斗機的最基本特征[1]。飛行器內(nèi)埋彈艙艙門在高速條件下打開后，艙內(nèi)及周圍將產(chǎn)生嚴(yán)重的壓力脈動，包含由湍流剪切層決定的寬帶噪聲以及由艙體幾何尺寸、外流條件決定的離散噪聲。強烈的脈動壓力會激發(fā)腔體及艙內(nèi)儲藏物產(chǎn)生振動，甚至導(dǎo)致材料失效。不穩(wěn)定的脈動流場同樣會嚴(yán)重影響武器的安裝投放[2]。隨著新一代航空航天器的發(fā)展，空腔脈動壓力(空腔噪聲)控制顯得尤為重要。目前，國內(nèi)外對于從低馬赫數(shù)到超聲速的空腔流動已經(jīng)做了大量研究，但是對于大馬赫數(shù)特別是高超聲速下的空腔壓力脈動及其控制研究還較少[3-4]。

流動控制分為被動控制和主動控制兩大類。其中，被動控制沒有外部能量輸入，一般靠改變幾何外形來控制流場，如加入擋板、擾流板、斜坡等；在主動控制中，外部能量可輸入到可調(diào)驅(qū)動器中，將工質(zhì)所具有的相應(yīng)動量直接作用在局部流場中，從而改變局部流場附面層的相互作用，并由此改變整個流場形態(tài)來達到控制流動的目的[5-7]。傳統(tǒng)內(nèi)埋武器彈艙一般采用被動控制措施來改變彈艙剪切層以降低噪聲水平，并確保武器安全分離。但在超聲速條件下，被動控制的實際效果并不理想。2001年，美國波音公司啟動了高頻激勵主動流動控制項目(HIFEX)用以發(fā)展其全球打擊飛行器全包線武器投放系統(tǒng)[8]，通過在武器艙內(nèi)安裝微型噴氣發(fā)動機陣列來實現(xiàn)對流場的控制。Bower等[9]在此基礎(chǔ)上開展了一系列主動控制方法研究。

計算氣動聲學(xué)(CAA)可在風(fēng)洞試驗前對設(shè)計修改進行分析，從而大大節(jié)省研制費用，縮短研制周期[10]。常用的數(shù)值方法包括直接數(shù)值模擬(DNS)、大渦模擬(LES)以及離散渦模擬(DES)等[11-13]，但對于許多工程實際問題來說，這些方法的計算代價大且存在不確定性。例如，DES方法經(jīng)常需要與有限傳輸算法結(jié)合使用，大大增加了亞格子尺度模型的耗散，將會導(dǎo)致有效黏度過大，且雷諾平均Navier-Stokes(RANS)方程和LES方法之間的統(tǒng)計學(xué)湍流能量傳輸也存在困難，從而限制了DES方法的應(yīng)用。2002年，Batten等[14-15]提出了一種非線性噪聲求解方法，通過對湍流物理量進行重構(gòu)可以兼顧亞格子尺度聲源的影響，在保持計算精度的同時降低了網(wǎng)格需求。目前，此方法已應(yīng)用于汽車、噴流等領(lǐng)域的噪聲預(yù)測[16-17]。

本文首次將Batten提出的非線性噪聲求解器應(yīng)用于內(nèi)埋彈艙的噪聲求解程序中，對超聲速條件下的空腔噪聲進行數(shù)值模擬，并與試驗結(jié)果進行比對以驗證方法的準(zhǔn)確性。在此基礎(chǔ)上，對超聲速和高超聲速工況下氣簾噴流對空腔噪聲的抑制效果開展研究，為內(nèi)埋彈艙的設(shè)計及優(yōu)化提供參考。

1 非線性噪聲求解方法

由初始統(tǒng)計平均湍流造成的噪聲可通過RANS求解得到，而非線性噪聲求解器主要用于模擬該噪聲的產(chǎn)生和傳播。所得的統(tǒng)計平均解不僅能夠提供平均流場的基本特征，也可以給出強制設(shè)定的湍流脈動的統(tǒng)計描述。非線性噪聲求解器在此統(tǒng)計結(jié)果上重構(gòu)噪聲源并高精度地模擬壓力脈動的傳播。非線性求解器可以使用原來的RANS計算網(wǎng)格來模擬壓力脈動的傳播，也可以在一個獨立的噪聲計算網(wǎng)格上插值噪聲進行數(shù)據(jù)求解。

(1)

(2)

(3)

(4)

(5)

(6)

(7)

式中：i、j、k的取值均為 1、2、3(其中，1、2、3分別表示x軸、y軸和z軸方向)；ρ為來流密度；ui(uj、uk)為擾動沿x軸(y軸、z軸)方向的速度；p為壓強；e為單位體積能；δij為克羅內(nèi)克函數(shù)；τij為剪切應(yīng)力項；θ為熱傳導(dǎo)項。忽略密度脈動項，對式(2)取時間平均可得

(8)

其中，

(9)

(10)

其中：

(11)

2 空腔驗證計算

本節(jié)將馬赫數(shù)Ma=1.5下空腔噪聲的計算結(jié)果與試驗值進行對比驗證。

為了捕捉流動在近壁面位置的細(xì)微結(jié)構(gòu)，對空腔進行了局部加密，使近壁面的第1層網(wǎng)格寬度為0.05 mm。彈艙壁面均采用無滑移壁面條件，黏性與溫度滿足Sutherland定律。

RANS計算時，入口條件為固定超聲速來流；遠場為特征線條件；出口邊界不指定，由內(nèi)層網(wǎng)格物理量推得。計算采用Cubick-ε模型，在得到當(dāng)?shù)乩字Z應(yīng)力張量和熱通量的統(tǒng)計平均值后，將其插值到噪聲計算網(wǎng)格上，再根據(jù)這一統(tǒng)計平均結(jié)果對湍流進行人工重構(gòu)。該方法的優(yōu)勢在于噪聲求解器可以在各項同性更好的網(wǎng)格單元上計算，特別是在近壁面區(qū)域。

噪聲計算需要啟動非定常計算，取時間步長為0.000 05，流動求解時長為1 s，每個時間步長內(nèi)的最大迭代次數(shù)為4次，新的邊界被設(shè)置為吸收層邊界，遠場及衰減層數(shù)據(jù)由之前的RANS計算提供。由于非線性噪聲求解方法對于近壁面網(wǎng)格的要求放寬以及聲場計算域的縮小，使噪聲計算網(wǎng)格數(shù)量較RANS方法減少，從而降低了計算成本。圖1為 RANS計算方法與非線性噪聲求解方法的網(wǎng)格對比。

圖1 雷諾平均N-S方程(RANS)網(wǎng)格與噪聲網(wǎng)格比較Fig.1 Comparison of Reynolds-averaged N-S equations (RANS) and acoustics meshes

湍動能反應(yīng)了流場中不同區(qū)域的湍流脈動強弱。從圖2中可以看出，湍動能在彈艙空腔后緣附近的水平很高，表示此處為湍流脈動最為強烈的地方。

圖2 空腔對稱面湍動能云圖Fig.2 Turbulence kinetic energy contours in symmetry plane of cavity

RANS計算結(jié)束后，接下來需要進行湍流脈動量的重構(gòu)。重構(gòu)時，為更加精確地重構(gòu)亞格子源項，將傅里葉級數(shù)項設(shè)為200。選取空腔底部中心線上10%X/L～100%X/L(X/L為預(yù)測點縱向相對位置，其中，X為空腔底部監(jiān)測點與前緣的水平距離，L為空腔長度)共10個測點進行壓力脈動數(shù)據(jù)采集。由于脈動壓力在噪聲計算初始階段不具有周期性，本文只選取了0.5～1.0 s之間的脈動壓力值進行研究，這段時間內(nèi)的壓力脈動呈現(xiàn)出周期性的特點。將計算得到的隨時間變化的脈動壓力值進行傅里葉變換即可得到相應(yīng)的頻譜特性，本文采用聲壓級(SPL)表示空腔自激振蕩的振幅，并將其定義為

(12)

式中：pref=2.0×10-5Pa。

Rossiter[18]于1964年提出了預(yù)測空腔振蕩頻率的半經(jīng)驗公式。該公式是評價數(shù)值仿真結(jié)果的重要標(biāo)準(zhǔn)，被稱之為Rossiter公式。

(13)

式中：U∞為來流速度；Ma∞為來流馬赫數(shù)；n為模態(tài)數(shù)；γ為比熱比；κ為流渦和自由流的流速比；α為相位延遲量。

圖3為數(shù)值模擬特征頻譜與Rossiter公式計算值的對比。其中，虛線為Rossiter公式預(yù)測的前四階頻率值。從圖中可以看出，采用非線性噪聲求解方法計算得出的前兩階模態(tài)能夠被精確捕捉到，這與Rossiter公式計算出的典型頻譜相符。

圖3 數(shù)值模擬聲壓頻譜與Rossiter公式計算值對比Fig.3 Comparison of sound pressure spectrum calculated by numerical simulation and Rossiter formula

圖4 空腔底部聲壓級分布計算值與試驗值對比Fig.4 Comparison of sound pressure level (SPL) distribution between calculated values and experimental values

圖4為空腔底部不同站位處總聲壓級的計算值與試驗值的對比。顯然，空腔監(jiān)測點的頻譜曲線中存在幾個非常明顯的主頻，與圖3中所示內(nèi)容一致。此外，從圖4中還可以看出，所有的計算值和噪聲測點的聲壓級分布集中在155～180 dB之間，且上游聲壓級略低，而下游聲壓級略高。除了X/L=1處外，其他監(jiān)測點的計算值與試驗值誤差在5 dB以內(nèi)，且試驗值略高于計算值，這主要是由于風(fēng)洞噪聲試驗存在一定的背景噪聲。而空腔底面的最大聲壓已經(jīng)高于170 dB，因此必須對腔體噪聲進行控制。

3 氣簾噴流主動控制

本文采用氣簾噴流主動控制措施來抑制空腔噪聲。圖5和圖6分別為氣簾噴流發(fā)生器示意圖與計算網(wǎng)格示意圖。其中，Dsj為噴流縫隙，Lsj為噴流縫隙到彈艙前緣的距離。腔體幾何尺寸、計算條件和脈動壓力監(jiān)測點位置與空腔驗證計算保持一致。氣簾噴流壓力Pj與來流壓力p∞之比為

(14)

圖5 氣簾噴流發(fā)生器示意圖Fig.5 Sketch of jet screen

圖6 腔體表面網(wǎng)格示意圖Fig.6 Sketch of grid on cavity surface

3.1 超聲速流動氣簾控制

超聲速流動氣簾的控制工況為馬赫數(shù)Ma=1.5，施加主動控制前后的空腔對稱面馬赫數(shù)云圖分布如圖7所示。在加入主動控制前，來流脫離空腔前緣后未觸及空腔底部便流出腔口；剪切層與空腔后壁碰撞，使空腔后部壓力急劇上升，誘發(fā)激波產(chǎn)生；一部分氣流隨激波繼續(xù)向下流動，另一部分則反射回來在腔內(nèi)形成渦流，使腔內(nèi)中前部的低壓部分與后部高壓部分連通。氣體的流動特性及腔內(nèi)存在的激波、膨脹波和循環(huán)流動等復(fù)雜流動環(huán)境會導(dǎo)致空腔后部發(fā)生顫振現(xiàn)象，從而產(chǎn)生噪聲，而噪聲傳播到空腔前緣部分后又會導(dǎo)致剪切層的分離，使以上循環(huán)過程再次發(fā)生。當(dāng)滿足一定的相位條件時，聲波的反饋循環(huán)形成，從而發(fā)生空腔內(nèi)的自持振蕩，使腔內(nèi)壓力條件發(fā)生劇烈變化，導(dǎo)致更為強大的噪聲環(huán)境，如圖7(a)所示。加入主動控制之后，腔內(nèi)流動變得平緩，氣簾噴流位于空腔前緣之前，射出的噴流破壞了原本穩(wěn)定的剪切層，從而改變了附面層內(nèi)部瑞流相干結(jié)構(gòu)。從圖7(b)中可以看出，來流剪切層跨過整個空腔，減弱了與空腔后壁撞擊形成的噪聲，從而有效地降低了噪點，達到了對空腔流動實施控制的效果。

圖7 空腔對稱面馬赫數(shù)分布云圖Fig.7 Comparison of Mach contour in symmetry plane of cavity

圖8為空腔與帶氣簾噴流的渦量瞬時等值面，背景顏色代表流向速度。加入氣簾控制后，腔體特別是后緣處的渦結(jié)構(gòu)趨于簡單，撞擊后緣的氣流減小，不穩(wěn)定剪切層被抬高，從而降低了艙體的氣動噪聲。

圖8 渦量瞬時等值面Fig.8 Instantaneous streamwise vorticity isosurfaces

圖9為馬赫數(shù)Ma=1.5下氣簾控制加入前后，各監(jiān)測點總聲壓級的對比示意圖。從圖中可以看出，在加入氣簾控制后，腔體底部各點總聲壓級均降低了10 dB左右，噪聲抑制效果明顯，表明計算結(jié)果與分析一致。

圖10為原始空腔與帶氣簾控制在X/L=0.2、0.8這2個監(jiān)測點的聲壓頻譜曲線對比圖。從圖中可以看出，二者的變化趨勢基本相同，能量都集中在低頻區(qū)域；氣簾噴流的加入對艙體聲壓頻譜特性也產(chǎn)生了影響，使5 000 Hz頻率以下的噪聲都有所降低，特別是對艙體流激振蕩峰值頻率所對應(yīng)的噪聲有明顯的抑制效果。

圖9 原始與帶氣簾噴流空腔各監(jiān)測點聲壓級對比Fig.9 Comparison of cavity SPL distribution between original cavity and cavity with jet screen

圖10 原始與帶氣簾噴流空腔頻譜曲線對比Fig.10 Spectrum comparison between original cavity and cavity with jet screen

3.2 高超聲速流動氣簾控制

3.2.1 高超聲速空腔噪聲計算方法驗證

Batten提出的非線性噪聲求解方法的適用范圍為可壓縮理想氣體。馬赫數(shù)Ma=5時的真實氣體效應(yīng)很弱，實際上可以忽略不計，其流動基本符合理想氣體的運動規(guī)律，因此，理論上可應(yīng)用非線性噪聲求解方法進行計算。

由于國內(nèi)對于高超聲速空腔噪聲的研究起步較晚，相關(guān)試驗數(shù)據(jù)較少，為了驗證非線性噪聲求解方法對高超聲速空腔噪聲預(yù)測的精度，本節(jié)首先對Unalmis等[19]研究的空腔噪聲試驗進行了計算驗證。其中，空腔的長深比為L/H=6、寬深比為W/H=3，Ma=5，選取的計算點位置與試驗脈動壓力監(jiān)測點RW7一致，具體位置如圖11所示。

圖11 壓力傳感器位置Fig.11 Location of pressure sensor

圖12為RW7處0.2～0.3 s的脈動壓力P計算值。從圖中可以看出，在空腔后緣上部的RW7處，監(jiān)測壓力有強烈的脈動現(xiàn)象發(fā)生。圖13為RW7處的聲壓頻譜計算值與試驗結(jié)果的對比。從圖中可以看出，在低頻段時，非線性噪聲求解方法的計算值與試驗值吻合較好，能較準(zhǔn)確地預(yù)測出聲壓的頻譜特性；在高頻段時，計算值小于試驗值，且隨著頻率的增加，二者之間的差距有逐步擴大的趨勢。這主要是由以下2個原因造成：①由于計算時間步長為5×10-5，根據(jù)奈奎斯特最大采樣頻率理論，計算結(jié)果中可以分辨的最大頻率為fmax=1/2T=104Hz，而fmax/2=5 000 Hz以內(nèi)的聲壓級結(jié)果更具有可信性，所以計算值與試驗值相比，在5 000Hz處比在低頻時相差更大；而且隨著頻率的進一步增大，誤差也會增大；②由于文獻未提供風(fēng)洞試驗背景噪聲，在進行計算時沒有將這部分考慮在內(nèi)，所以計算值會略低于試驗值。

圖12 RW7處的脈動壓力計算值Fig.12 Calculated values of RW7 oscillating pressure

圖13 RW7處的頻譜計算值與試驗結(jié)果對比Fig.13 Comparison of spectrum between calculated values and experimental values

綜上所述，非線性噪聲求解方法能較好地預(yù)測Ma=5左右的空腔噪聲。因此，可將非線性超聲求解方法引入到脈動壓力主動控制效果評估中。

3.2.2 高超聲速空腔脈動壓力氣簾控制

高超聲速流場的突出特點是存在強激波，加入氣簾控制后，流場中除激波之間的相互干擾外，還存在旋渦以及激波與旋渦的干擾，進行數(shù)值模擬有一定難度。本節(jié)對高超聲速空腔流場中的噴流進行了研究。圖14～圖17分別為噴口附近的流場流線圖、流場等密度線、流場等壓強線以及流場等馬赫數(shù)線的分布。從圖14～圖17中可以看出，自由來流受到噴流阻礙后產(chǎn)生分離形成分離區(qū)，區(qū)內(nèi)包含分離渦。噴口內(nèi)高壓氣體在噴口兩側(cè)產(chǎn)生很強的膨脹波，波系的相互干擾導(dǎo)致噴流過度膨脹產(chǎn)生馬赫盤。同時，由于噴流的作用，噴口前緣處產(chǎn)生了二次分離。

圖14 噴口附近流場流線圖Fig.14 Streamline near jet screen

圖15 噴口附近流場等密度線Fig.15 Isodensity lines near jet screen

圖16 噴口附近流場等壓強線Fig.16 Isobaric lines near jet screen

圖17 流場等馬赫數(shù)線Fig.17 Mach contour of cavity flow field

隨后，對高超聲速條件下的空腔脈動壓力載荷及氣簾控制效果進行了評估。圖18為2種高超聲速條件下的空腔與帶氣簾噴流計算仿真的湍動能云圖。從圖中可以看出，在加入噴流控制后，艙體后緣附近的湍動能明顯降低，說明后緣點附近脈動壓力減小。

圖18 空腔對稱面湍流動能云圖Fig.18 Turbulence kinetic energy contours in symmetry plane of cavity

圖19為加入氣簾控制前后的空腔后緣點上壓力隨時間變化的曲線圖(V為來流速度)。顯然，加入氣簾控制之后，壓力降低9%左右，主動控制效果明顯。

圖19 原始與帶氣簾噴流空腔后緣點的壓力脈動變化Fig.19 Comparison of dynamic load on rear wall between original cavity and cavity with jet screen

4 結(jié) 論

1) 基于非線性噪聲求解方法，對三維空腔在超聲速和高超聲速條件下進行了噪聲預(yù)測。通過與試驗結(jié)果對比，非線性噪聲求解方法能夠準(zhǔn)確預(yù)測出空腔噪聲聲壓級，同時極大地減小了計算量。在現(xiàn)有計算條件下，與LES或DES等方法相比，具有更廣闊的應(yīng)用前景。

2) 氣簾噴流主動控制對高速飛行器的空腔噪聲抑制具有明顯作用。超聲速條件下，采用氣簾噴流主動控制后，空腔聲壓頻譜曲線的趨勢與不加控制時基本一致，但是聲壓峰值明顯降低，各預(yù)測點總聲壓級降低10 dB左右。高超聲速條件下，采用氣簾噴流主動控制后，脈動壓力降低9%左右。

3) 在實際工程應(yīng)用中，主動控制措施與本文模擬的氣簾噴流還有區(qū)別。后續(xù)將在此基礎(chǔ)上進一步細(xì)化主動控制研究，對非線性噪聲求解方法進行更為深入的探索，使其更好地應(yīng)用于高超聲速飛行器空腔噪聲研究。

[1] Wu J F， Luo X F， Fan Z L. Flow control method to improve cavity flow and store separation characteristic[J]. Acta Aeronautica et Astronautica Sinica， 2009， 30(10)： 1840-1845 (in Chinese). 吳繼飛， 羅新福， 范召林.內(nèi)埋式彈艙流場特性及武器分離特性改進措施[J]. 航空學(xué)報， 2009， 30(10)： 1840-1845.

[2] Feng Q， Cui X C. Study on integrated flow control for weapons bay of flying wing configuration aircraft[J]. Acta Aeronautica et Astronautica Sinica， 2012， 33(5)： 781-787 (in Chinese). 馮強， 崔曉春. 飛翼布局飛機武器艙綜合流動控制技術(shù)研究[J]. 航空學(xué)報， 2012， 33(5)： 781-787.

[3] Perng S W. Passive control of pressure oscillations in hypersonic cavity flow[D]. Austin： University of Texas at Austin， 1996.

[4] Barter J W. Prediction and passive control of fluct-uating pressure loads produced by shock-induced turbulent boundary layer separation[D]. Austin： University of Texas at Austin， 1995.

[5] Shaw L. Active control for cavity acoustics， A1AA-1998-2347[R].Reston： AIAA， 1998.

[6] Miehael J A.Control of cavity resonance through very high frequency forcing， AIAA-2000-1905[R]. Reston： AIAA， 2000.

[7] Smith B R， Welterlen T J， Maines B H， et al.Weapons bay acoustic suppression from rod spoilers，AIAA-2002-0662[R].Reston： AIAA， 2002.

[8] Alvi F S， Elavarasan R， Garg G， et al. Control of supersonic impinging jet using microjets， AIAA-2000-2236 [R]. Reston： AIAA， 2000.

[9] Bower W W， Kibens V， Cary A W， et al. High-frequency excitation active flow control for high-speed weapon release(HIFEX)， AIAA-2004-2513[R]. Reston： AIAA， 2004.

[10] Tam C K W. Computational aeroacoustics：an overview of computational challenges and applications[J]. International Journal of Computational Fluid Dynamics， 2004， 18(6)： 547-567.

[11] Long S L， Nie H， Xue C J， et al. Simulation and experiment on aeroacoustic noise characteristics of aircraft landing gear[J]. Acta Aeronautica et Astronautica Sinica，2012， 33(6)： 1002-1013 (in Chinese). 龍雙麗， 聶宏， 薛彩軍， 等. 飛機起落架氣動噪聲特性仿真與試驗[J]. 航空學(xué)報， 2012， 33(6)： 1002-1013.

[12] Hou Z X， Yi S H， Wang C Y. Numerical analysis of supersonic open cavity[J]. Journal of Propulsion Technology， 2001， 22(5)： 400-403 (in Chinese). 侯中喜， 易仕和， 王承堯. 超聲速開式空腔流動的數(shù)值模擬[J]. 推進技術(shù)， 2001， 22(5)： 400-403.

[13] Li X D， Liu J D， Gao J H. Numerical simulation of flow-induced oscillation and sound generation in a cavity[J]. Chinese Journal of Theoretical and Applied Mechanics， 2006， 38(5)： 599-604 (in Chinese). 李曉東， 劉靖東， 高軍輝. 空腔流激振蕩發(fā)聲的數(shù)值模擬研究[J]. 力學(xué)學(xué)報， 2006， 38(5)： 599-604.

[14] Batten P， Ribaldone E， Casella M， et al. Towards a generalized non-linear acoustics solver， AIAA-2004-3001[R]. Reston： AIAA， 2004

[15] Batten P，Goldberg U，Chakravarthy S. Reconstructed sub-grid methods for acoustics predictions at all Reynolds Numbers， AIAA-2002-2511[R]. Reston： AIAA， 2002.

[16] Silva C R I， Almeida O， Batten P. Investigation of an axi-symmetric subsonic turbulent jet using computational aeroacoustics tools，AIAA-2007-3656[R]. Reston： AIAA， 2007.

[17] Aflalo B S， Simoes L G C， Silva R G， et al. Comparative analysis of turbulence models for slat noise source calculations employing unstructured meshes， AIAA-2010-3838 [R]. Reston： AIAA， 2010.

[18] Rossiter J E.Wind tunnel experiments on the flow over rectangular cavities at subsonic and transonic speeds， RAE Technical Report No. 64037[R]. Hampshire： Royal Aircraft Establishment， 1964.

[19] Unalmis O H， Clemens N T， Dolling D S. Experimental study of shear-layer/acoustics coupling in Mach 5 cavity flow[J]. AIAA Journal， 2001， 39(2)： 242-252.

Tel： 025-84896031

E-mail： tongw@nuaa.edu.cn

張杰 男，博士，工程師。主要研究方向： 計算氣動聲學(xué)。

E-mail： jerry.zhang@visionstrategy.com.cn

*Corresponding author. Tel.： 025-84896031 E-mail: tongw@nuaa.edu.cn

Active control and nonlinear numerical simulation for oscillating pressure of high-speed aircraft cavity

WANG Yiding1,2, Guo Liang1,2, TONG Mingbo1,*, Zhang Jie3

1.CollegeofAerospaceEngineering,NanjingUniversityofAeronauticsandAstronautics,Nanjing210016,China2.ChengduAircraftDesignandResearchInstitute,Chengdu610091,China3.BeijingYinjingTechnologyCo.,Ltd.，Beijing100107，China

Prediction of oscillating pressure is a key technology for the weapon bay of high-speed aircraft cavity. Nonlinear numerical simulation is proposed as a new method to analyze noise recently. In order to evaluate the prediction performance of cavity noise, nonlinear numerical simulation solver is combined with Reynolds-averaged Navier-Stokes(RANS) equation. Firstly, the flow field around cavity is solved by RANS, and the average solution of initial turbulent statistics is obtained which contains the basic characteristics of average flow field and statistics description of turbulence fluctuation. After that, noise source is refactored and the spreading of pressure fluctuation is simulated precisely by the nonlinear acoustics solution. According to the comparison of the cavity noise calculation and experimental results underMa= 1.5 andMa= 5, it indicates that nonlinear numerical solution is able to well predict cavity flow noise at high speed. Based on that, the contribution to noise suppression made by active control such as adding jet screen at the leading edge of cavity underMa= 1.5 andMa= 5 is investigated. It is found that jet screen is suitable for suppression of oscillating pressure under supersonic condition as well as in hypersonic state.

oscillating; nonlinear; cavity; turbulence; active control

2014-07-02； Revised： 2014-09-10； Accepted： 2014-10-08； Published online： 2014-10-23 15:45

Aeronautical Science Foundation of China (2012ZC52035)

2014-07-02; 退修日期： 2014-09-10; 錄用日期： 2014-10-08; 網(wǎng)絡(luò)出版時間： 2014-10-23 15：45

www.cnki.net/kcms/detail/10.7527/S1000-6893.2014.0277.html

航空科學(xué)基金(2012ZC52035)

Wang Y D, Guo L, Tong M B, et al. Active control and nonlinear numerical simulation for oscillating pressure of high-speed aircraft cavity[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(1): 213-222. 王一丁， 郭亮， 童明波， 等． 高速飛行器空腔脈動壓力主動控制與非線性數(shù)值模擬[J]．航空學(xué)報, 2015, 36(1): 213-222

http://hkxb.buaa.edu.cn hkxb@buaa.edu.cn

10.7527/S1000-6893.2014.0277

V211.3

A

1000-6893(2015)01-0213-10

王一丁 男，博士研究生。主要研究方向： 計算氣動聲學(xué)，噪聲控制。

E-mail： wyding127@163.com

郭亮 男，博士研究生，工程師。主要研究方向： 飛機氣動綜合設(shè)計。

E-mail： brightgl@hotmail.com

童明波 男，博士，教授，博士生導(dǎo)師。主要研究方向： 飛行器總體設(shè)計。

*通訊作者.Tel.： 025-84896031 E-mail: tongw@nuaa.edu.cn

URL： www.cnki.net/kcms/detail/10.7527/S1000-6893.2014.0277.html