使用自行發(fā)展的MFlow求解器對第六屆AIAA阻力預(yù)測會議的翼身組合體外形進(jìn)行了數(shù)值模擬研究。在文章中著重對比了線性渦黏模型和考慮二次本構(gòu)關(guān)系的湍流模型的數(shù)值模擬結(jié)果,并與試驗結(jié)果進(jìn)行了比較。研究發(fā)現(xiàn),考慮二次本構(gòu)關(guān)系的湍流模型能夠比較準(zhǔn)確地預(yù)測翼身組合體的氣動特性,而線性渦黏模型從迎角3.5°開始,高估了翼身結(jié)合處的分離流動,導(dǎo)致預(yù)測的氣動力與試驗結(jié)果差異較大。對翼身結(jié)合處局部流動和壓力分布的詳細(xì)分析,展示了兩種模型預(yù)測的氣動特性差別很大的原因。 

【文章來源】：空氣動力學(xué)學(xué)報. 2020,38(05)北大核心CSCD

【文章頁數(shù)】：7 頁

【部分圖文】：

NASA通用研究模型

首先為了驗證計算程序和網(wǎng)格，使用逐漸加密的四套網(wǎng)格完成定升力CL=0.5狀態(tài)下的網(wǎng)格收斂性研究。計算的外形是迎角為2.75°下的靜氣動彈性變形模型。第二個算例是使用中等規(guī)模的網(wǎng)格完成靜氣動彈性影響的預(yù)測。迎角從2.5°到4°，以0.25°為間隔。使用的外形是在相應(yīng)迎角下的靜氣動彈性變形模型。計算狀態(tài)為馬赫數(shù)0.85，基于平均氣動弦長的雷諾數(shù)為5×106。2 數(shù)值方法

計算中使用SA一方程模型[18]，雷諾應(yīng)力分別由線性渦黏假設(shè)（后續(xù)結(jié)果中記為SA）和二次本構(gòu)關(guān)系（記為SA＿QCR）給出。本文的計算從自由來流開始，假定流動為全湍流。圖3為本文計算的殘差收斂和氣動力變化曲線，所有計算密度的殘差都下降了兩個量級以上，氣動力的振蕩在0.2%以內(nèi)。3 結(jié)果和討論


本文編號：3307518