泵噴推進(jìn)器水動力噪聲數(shù)值模擬方法研究現(xiàn)狀?

白 日 龔京風(fēng) 趙 勝 尤明星 王 晴 劉 春

（1.武漢科技大學(xué)汽車與交通工程學(xué)院 武漢 430065）（2.中國艦船研究設(shè)計中心 武漢 430064）

1 引言

在泵噴推進(jìn)器設(shè)計制造初期，由于艦艇航速有限以及機(jī)械噪聲和螺旋槳噪聲的掩蓋作用，水動力噪聲并未得到設(shè)計制造者的足夠重視。隨著機(jī)械噪聲和螺旋槳噪聲得到有效控制，高航速帶來的水動力噪聲在水下艦艇輻射噪聲中的貢獻(xiàn)逐步凸顯，探索高效準(zhǔn)確的水動力噪聲數(shù)值模擬方法日益迫切。而水動力噪聲研究的關(guān)鍵就是明確噪聲源特性及選擇恰當(dāng)?shù)难芯糠椒ā?/p>

水動力噪聲是指由于流體流過槳葉等部件表面與固體邊界產(chǎn)生粘性作用，以及流體內(nèi)部湍流、渦脫落等而引起的輻射噪聲［1］。泵噴作為噪聲源的輻射噪聲中，根據(jù)噪聲物理特性的不同，可分為以下四類；1）空泡噪聲；2）低頻離散譜噪聲；3）低頻連續(xù)譜噪聲；4）高頻連續(xù)譜噪聲。在一般情況下，潛艇中高速航行時，低頻噪聲是泵噴推進(jìn)器水動力噪聲的主要貢獻(xiàn)者，也是學(xué)者研究的重點(diǎn)。

在確定了水動力噪聲機(jī)理和特性之后，最為關(guān)鍵的就是研究方法的探究。無論是過去還是現(xiàn)在，計算水動力的噪聲的基本思路都是先計算泵噴推進(jìn)器水動力性能，而后進(jìn)行聲場求解預(yù)報。本文也將從這兩個方面進(jìn)行回顧和評述，從而為研究人員提供參考。雖然眾多學(xué)者嘗試了各種數(shù)值模擬方法，但由于泵噴推進(jìn)器結(jié)構(gòu)復(fù)雜，三維流場模擬又十分棘手，許多文章對其作用機(jī)理解釋不清，也沒有很好地將各種數(shù)值方法進(jìn)行對比分析，多年來還未在國內(nèi)外見過相關(guān)的綜述性文獻(xiàn)資料，因此，本文對泵噴推進(jìn)器水動力噪聲的研究工作具有現(xiàn)實(shí)意義。

2 水動力性能數(shù)值模擬方法

目前，泵噴推進(jìn)器水動力性能數(shù)值模擬方法主要有面元法和CFD方法，此外，也有學(xué)者試圖結(jié)合勢流理論計算簡便、計算周期短和CFD方法精細(xì)真實(shí)的優(yōu)勢，采用勢流和粘流耦合的方法對泵噴推進(jìn)器進(jìn)行水動力性能預(yù)報。

2.1 面元法

面元法是一種比較成熟的旋轉(zhuǎn)機(jī)械設(shè)計和數(shù)值分析方法。根據(jù)不同的分類標(biāo)準(zhǔn)，面元法可以分為不同的類型，各種類型都有其應(yīng)用特點(diǎn)和適用性。目前主流的面元法分類有：基于速度的面元法和基于速度勢的面元法、低階面元法和高階面元法。

基于速度的面元法在其整個計算過程中只利用了速度求解，而速度勢本身在整個計算過程中沒有在實(shí)質(zhì)上引入［2］?；谒俣葎莸拿嬖ǖ挠嬎阍硎抢脕砹鳁l件事先確定物體表面上的源匯分布情況，然后得到求解模型表面偶極子強(qiáng)度的求解方程［3～4］。學(xué)者對二者的適用性進(jìn)行了對比研究，對于較薄的旋轉(zhuǎn)機(jī)械部件，基于速度勢的面元法計算精度更高［5］。由于泵噴推進(jìn)器的轉(zhuǎn)子也具有較薄葉片的特點(diǎn)，因此采用基于速度勢的面元法對泵噴推進(jìn)器流場特性分析更為合適。就目前研究現(xiàn)狀而言，基于速度勢的面元法正逐步成為主流的面元法。

低階面元法是將面元上的奇點(diǎn)散布形式假設(shè)為均勻散布，從而減小了計算難度，使實(shí)際問題更加簡化，但與實(shí)際求解模型有較大誤差。而高階面元法是將面元上的奇點(diǎn)散布假設(shè)為線性分布或曲線分布，與實(shí)際模型相差較小，精度更高［6］。針對低階面元法計算量小，計算效率高的優(yōu)勢，諸多學(xué)者采用低階面元法對各種推進(jìn)器進(jìn)行數(shù)值仿真研究。Baltazar［7］等采用低階面元法對Ka4-70 導(dǎo)管槳的敞水性能進(jìn)行了模擬，分析了尾流模型、導(dǎo)管邊界層以及尾渦螺距等因素對計算結(jié)果的影響。Takinaci［8］考慮了邊界層的影響，采用低階面元法求解了三維矩形的薄翼型結(jié)構(gòu)的流場速度分布，通過與實(shí)驗(yàn)數(shù)據(jù)相比較，發(fā)現(xiàn)預(yù)報的速度分布和實(shí)驗(yàn)值基本一致。隨著研究的進(jìn)一步深入，對泵噴推進(jìn)器水動力性能求解精度有了更加嚴(yán)格的要求，學(xué)者開始探索高階面元法的應(yīng)用，雖然高階面元法精度高，但是它的計算量也大。目前，學(xué)者采用的解決方法是對其模型進(jìn)行簡化，其中，對導(dǎo)管和轉(zhuǎn)子之間的相互影響進(jìn)行了時域簡化，對導(dǎo)管和葉輪、葉輪和定子葉片之間的誘導(dǎo)速度計算采用周向混合平均化處理，降低模型的復(fù)雜性，從而減少計算量［9］。最終通過實(shí)驗(yàn)驗(yàn)證發(fā)現(xiàn)，簡化后的模型對計算結(jié)果影響不大。隨著泵噴推進(jìn)器水動力計算精度要求的進(jìn)一步提升，高階面元法很有可能成為下一代主流的面元法。

2.2 CFD法

計算流體力學(xué)（CFD）是一種基于粘性理論的流體仿真方法，目前主要的CFD求解方法有直接模擬（DNS）法、大渦模擬（LES）法和雷諾時均（RANS）法。

DNS是直接求解封閉的N-S方程，不對流動情況做任何模型假設(shè)，計算耗時長，占用資源大，受計算機(jī)計算能力限制，在泵噴推進(jìn)器水動力應(yīng)用很少。RANS方法求解時均化的N-S方程，其基本原理是把湍流運(yùn)動看成是一個時均運(yùn)動加上一個脈動量，這樣時均化的N-S方程引入了新的湍流脈動值項(xiàng)，需要建立新的應(yīng)力表達(dá)式才能封閉方程組，因而引入了各種湍流模型，在泵噴推進(jìn)器水動力性能計算方面的應(yīng)用較為廣泛。LES的基本原理是假設(shè)動量、能量、質(zhì)量等主要由大渦輸運(yùn)，而流動特性只由大渦體現(xiàn)，隨后只對大渦進(jìn)行直接求解，而小渦則被過濾掉利用模擬的方法考慮其對輸運(yùn)的影響，但大渦模擬對小渦的網(wǎng)格節(jié)點(diǎn)劃分會非常密集，極大的增加了計算量，對于泵噴推進(jìn)器三維復(fù)雜流場的模擬仍有很大困難，目前應(yīng)用較少，但LES作為計算量介于DNS和RANS之間的數(shù)值模型，兼具二者優(yōu)勢，相信隨著LES理論和軟件技術(shù)的不斷完善，LES將具有很好的應(yīng)用前景?？偟膩碇v，RANS方法是目前泵噴推進(jìn)器水動力性能數(shù)值模擬方法中最為成熟的CFD方法。

在RANS方法中SST k-ω湍流模型比較適合于泵噴推進(jìn)器水動力計算并且已經(jīng)得到了廣泛應(yīng)用。該模型的最大的優(yōu)點(diǎn)是考慮了湍流剪切應(yīng)力對流場特性的影響，從而避免了對渦流粘度的過渡預(yù)測，結(jié)合了k-ε模型湍流區(qū)模擬精確和k-ω模型在粘性底層模擬精確的優(yōu)勢［10］。劉占一［11］、孫瑜［12］、姜漢［13］、饒志強(qiáng)［14］、潘光［15］等利用 SST k-ω模型對泵噴推進(jìn)器水動力性能進(jìn)行了深入研究，其中鹿麟［16］、劉登成［17］等采用SST k-ω模型研究了推進(jìn)器的間隙流動，導(dǎo)管槳或者泵噴推進(jìn)器中轉(zhuǎn)子和導(dǎo)管間隙非常小，網(wǎng)格節(jié)點(diǎn)的布置非常困難，是泵噴推進(jìn)器水動力數(shù)值模擬的難點(diǎn)之一，但間隙流場對泵噴推進(jìn)器流場特性影響卻非常大，不容忽視，此外，間隙流場也對聲場產(chǎn)生了重要影響［16］，可見采用SST k-ω模型可以捕捉到比較精細(xì)的流場信息，對泵噴推進(jìn)器水動力性能計算的適用性很好。

2.3 勢流/粘流耦合法

粘流/勢流耦合法是一種結(jié)合面元法和CFD方法二者優(yōu)勢的綜合方法。目前，粘勢/流耦合方法在國內(nèi)外已經(jīng)取得了一定的研究成果。黃璐［18］利用粘流/勢流耦合法，將CFD方法計算得到的尾渦模型提取出來作為面元法的尾渦模型進(jìn)行了流場計算，并于試驗(yàn)值進(jìn)行了對比。試驗(yàn)結(jié)果表明，粘流/勢流耦合法相對于傳統(tǒng)的面元法精度有一定提高，但對轉(zhuǎn)矩系數(shù)計算結(jié)果的改進(jìn)并不明顯。Black［19］采用粘性流與高階面元法耦合方法，進(jìn)行了吊艙推進(jìn)器的結(jié)構(gòu)設(shè)計和流場特性預(yù)報的研究分析。Warren［20］等在混合網(wǎng)格技術(shù)的基礎(chǔ)上，探索發(fā)展了粘流/勢流結(jié)合的方法，利用二者優(yōu)勢，對泵噴推進(jìn)器流場特性進(jìn)行了分析。粘流/勢流耦合方法雖然能夠在一定程度上提高計算精度，但目前研究較少，仍存在諸如精度提高不明顯等問題，而且目前粘流/勢流耦合方法在水動力性能計算方面的應(yīng)用較少，該方法目前處于探索階段，技術(shù)原理仍不夠成熟。

3 噪聲預(yù)報方法

在目前已公開的研究中，學(xué)者們普遍將泵噴推進(jìn)器噪聲預(yù)報分為旋轉(zhuǎn)葉輪和靜止部件（定子、導(dǎo)管）兩部分，在考慮聲反射情況下，將兩部分的聲場疊加即得到泵噴推進(jìn)器聲場。靜止部件噪聲計算是基于聲類比方程求解的。對于轉(zhuǎn)子聲場的求解，目前較為成熟的模型是扇聲源模型（點(diǎn)源模型）或基于聲類比方程求解。其中，聲類比方程主要應(yīng)用于轉(zhuǎn)子的自由聲場計算，不能考慮泵噴推進(jìn)器導(dǎo)管對轉(zhuǎn)子聲場的輻射和折射效應(yīng)［21～22］，而扇聲源模型可以考慮固體邊界的聲散射效應(yīng)［23～24］。

對于扇聲源模型，轉(zhuǎn)子輻射壓力在軸向和切向的貢獻(xiàn)量［25］為

在徑向的貢獻(xiàn)量為

其中：m是諧波數(shù)；B是轉(zhuǎn)子葉片的數(shù)量；Ω是轉(zhuǎn)速；R是場點(diǎn)到軸系原點(diǎn)距離；c0是聲速；Fs是轉(zhuǎn)子葉片上脈動壓力的傅里葉級數(shù)；M是旋轉(zhuǎn)馬赫數(shù)；如圖1所示。

圖1 風(fēng)扇示意圖［25］

上述表達(dá)式中，壓力場由轉(zhuǎn)子葉片上的載荷計算得到，葉片載荷可通過流場分析獲得。將泵噴推進(jìn)器非定常水動力性能計算得到的脈動壓力映射到聲學(xué)網(wǎng)格。流場求解得到的非定常脈動壓力可以看做是一個分布在聲學(xué)網(wǎng)格上的旋轉(zhuǎn)點(diǎn)力源，通過對整個旋轉(zhuǎn)運(yùn)動聲場的力源聲場積分便可得到整個旋轉(zhuǎn)運(yùn)動聲場［24］。扇聲源模型將聲學(xué)網(wǎng)格上的點(diǎn)力源等效為緊致聲源，大大減少了計算量。

如前所述，泵噴低頻噪聲是泵噴推進(jìn)器整體輻射噪聲中的主要貢獻(xiàn)者。目前對泵噴推進(jìn)器噪聲預(yù)報的研究也主要集中于低頻噪聲領(lǐng)域。付建［26］結(jié)合分塊離散方法將轉(zhuǎn)子葉片上的偶極子等效為緊致聲源，預(yù)報泵噴推進(jìn)器低頻噪聲，大大降低了計算量，數(shù)值計算結(jié)果與實(shí)驗(yàn)值誤差較小。結(jié)果表明，靜止部件聲場和旋轉(zhuǎn)部件聲場合成后的泵噴推進(jìn)器總聲場特性與螺旋槳聲場存在較大差異，不能直接套用求解螺旋槳聲場的方式求解泵噴推進(jìn)器的總聲場。劉敏［27］結(jié)合扇聲源模型和有限元方法，通過多組數(shù)值仿真結(jié)果的比較得出導(dǎo)管在高頻段對轉(zhuǎn)子噪聲降噪效果明顯的結(jié)論。盧丁?。?8］結(jié)合扇聲源模型和邊界元法求解某泵噴模型聲場，以葉頻處的聲場特性為對象作進(jìn)一步分析，結(jié)果表明導(dǎo)管對轉(zhuǎn)子聲場的徑向分布有較大影響，如圖2所示。

圖2 BPF處的聲壓分布云圖［29］

考慮泵噴推進(jìn)器流噪聲預(yù)報涉及到流場和聲場的耦合問題，由于流體流動能量只有很少一部分轉(zhuǎn)化為聲能在介質(zhì)中傳播，流動特征尺度和聲擾動特征尺度差別很大，因此，數(shù)值模擬中的微小誤差可能導(dǎo)致噪聲預(yù)報較大偏差。因此流致噪聲的數(shù)值模擬除了一般流體模擬的要求外，在兩個方面需要特殊考慮：一是流致噪聲數(shù)值模擬要有準(zhǔn)確的邊界條件，一般計算域很難將非線性區(qū)域全部考慮到，這需要邊界條件能吸收邊界以外的流體擾動，避免產(chǎn)生額外的聲反射。二是非定常水動力性能計算得到的脈動壓力映射到聲學(xué)網(wǎng)格數(shù)據(jù)映射精度的問題。對于第一個方面，楊瓊方［30］在其國家自然科學(xué)基金項(xiàng)目中分析得出對于水下旋轉(zhuǎn)聲源，只要噪聲分析頻率不超過30kHz，其線性效應(yīng)就可以完全忽略不計。Virtual.Lab在邊界層應(yīng)用PML或AML技術(shù)可以避免產(chǎn)生額外的聲發(fā)射。對于第二個方面，Virtual.Lab對于網(wǎng)格數(shù)據(jù)的映射提供了三種算法以滿足不同情況，從而降低數(shù)值誤差，且目前有學(xué)者通過Virtual.Lab13版本嘗試對泵噴推進(jìn)器流噪聲一次性整體預(yù)報，從而避免了數(shù)據(jù)映射帶來的誤差［31］。

4 結(jié)語

目前，主流的泵噴推進(jìn)器水動力性能研究方法是面元法和CFD方法。面元法有自己的優(yōu)勢：但面元法無法捕捉流場的精細(xì)流動，且在面元法計算過程中引入了假定的尾渦模型，降低了計算精度，因此相較于CFD方法計算精度仍有待提高。為保證計算結(jié)果的可靠性，多數(shù)學(xué)者采用基于粘流理論的CFD方法作為泵噴推進(jìn)器水動力性能的研究方法。其中采用SST k-ω湍流模型已成為一個較為成熟的研究方法。

粘流/勢流耦合法理論上是較為理想的方法，但將其用于泵噴推進(jìn)器水動力研究并不普遍，計算精度及穩(wěn)定性易受到黏流求解器設(shè)置和勢流程序的影響，該方法有待于進(jìn)一步發(fā)展。

扇聲源模型能夠很好地預(yù)報泵噴推進(jìn)器轉(zhuǎn)子低頻噪聲。已公開的研究結(jié)果表明由于導(dǎo)管的存在，泵噴推進(jìn)器轉(zhuǎn)子聲場與螺旋槳自由聲場在徑向方向差異較大。在聲學(xué)軟件Virtual.Lab中能考慮到聲音在導(dǎo)管中的反射，預(yù)報結(jié)果更符合泵噴推進(jìn)器的實(shí)際情況。