1.2 數(shù)值方法
采用NUMECA 商用軟件包FINE/TURBO 進(jìn)行數(shù)值計算。該軟件采用時間相關(guān)法求解雷諾平均的N-S 方程,中心節(jié)點(diǎn)的有限體積法離散, 顯式Runge-Kutta 法求解,全多重網(wǎng)格初場處理和多重網(wǎng)格迭代加速,以及低速流動的預(yù)處理技術(shù)等。
結(jié)構(gòu)化網(wǎng)格由Auto Grid 模塊自動生成。圖2 給出計算域網(wǎng)格及翼型周圍的局部放大圖,計算域尺度為60 倍翼型弦長,翼型周圍網(wǎng)格點(diǎn)數(shù)為275,網(wǎng)格總數(shù)約15萬,Y+小于3。

計算域外邊界設(shè)為遠(yuǎn)場條件,給定來流速度,靜壓和靜溫等,壁面為絕熱無滑移條件。計算全局殘差收斂到5 個量級以上。
湍流模型采用AGS 模型。
2. 計算結(jié)果及分析
2.1 計算結(jié)果確認(rèn)

對DF-003-300-DU、DF-002-350-DU、DF-001-400-DU3個翼型在不同攻角下進(jìn)行數(shù)值模擬,計算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對比, 如圖3 所示, 其中“exp”表示翼型的實(shí)驗(yàn)數(shù)據(jù),翼型名稱表示模擬結(jié)果,以下各圖相同。
從圖3 可以看出,3 個翼型數(shù)值計算的升力系數(shù)與阻力系數(shù)與實(shí)驗(yàn)值趨勢一致,大部分區(qū)域吻合良好,確認(rèn)了線性段數(shù)值計算的可靠性;隨著相對厚度的增大,數(shù)值模擬的結(jié)果與實(shí)驗(yàn)值的吻合度越來越低,這是由于選取了相對厚度較大的3 個翼型進(jìn)行數(shù)值模擬,而大厚度翼型在CFD 計算又是一個難點(diǎn),計算很難收斂,本文選用AGS 模型進(jìn)行計算。在小攻角的工況下,AGS 模型能夠得到與實(shí)驗(yàn)值較為一致的計算結(jié)果,這是由于在附著流中,轉(zhuǎn)捩的經(jīng)驗(yàn)計算式能夠較為準(zhǔn)確地預(yù)測轉(zhuǎn)捩;隨著攻角的增加,流動變?yōu)榉蛛x流,由于現(xiàn)有的轉(zhuǎn)捩預(yù)測模型對分離流預(yù)測的并不理想,因此,數(shù)值計算值與實(shí)驗(yàn)值誤差較大。同時,由于實(shí)驗(yàn)翼型表面粗糙度以及實(shí)驗(yàn)空氣流場粘度、密度的影響等,使模擬結(jié)果與實(shí)驗(yàn)值之間會產(chǎn)生差別。