P ω 增強型k-ω湍流模型在三角翼旋渦流動的應(yīng)用
發(fā)布時間:2021-09-03 03:38
在三角翼旋渦繞流數(shù)值模擬中,標準Wilcox k-ω湍流模型生成項未考慮旋度的影響而導(dǎo)致預(yù)測的旋渦強度較弱。通過引入探測因子區(qū)分剪切層和渦核,在旋渦流動的高旋度區(qū)域增加ω方程生成項的方法,基于結(jié)構(gòu)化網(wǎng)格上的RANS求解器,加入了Pω增強型k-ω湍流模型,對繞尖前緣三角翼亞聲速和跨聲速旋渦流場進行了數(shù)值模擬。計算結(jié)果與NASA的NTF風洞和DLR的DNW-TWG風洞試驗數(shù)據(jù)進行了對比分析,結(jié)果表明:不論在亞聲速還是跨聲速自由來流條件下,Pω增強型k-ω湍流模型計算的壓力分布、渦破裂位置均與試驗數(shù)據(jù)吻合良好,準確地預(yù)測出了三角翼上翼面的主渦、二次渦結(jié)構(gòu),特別是跨聲速條件下激波干擾導(dǎo)致的渦破裂的臨界迎角及渦破裂位置,表明Pω增強型k-ω湍流模型在繞三角翼旋渦流動數(shù)值模擬中具有良好的適用性。
【文章來源】:空氣動力學(xué)學(xué)報. 2016,34(04)北大核心CSCD
【文章頁數(shù)】:8 頁
【部分圖文】:
圖165°后掠三角翼模型[12]Fig.1Geometryofthe65°sweptleadingedgedeltawing
表1計算狀態(tài)對比Table1CFDcaseschosenforsummarycomparison計算狀態(tài)Ma∞Re∞(×106)α(°)12亞聲速0.40.46618.523.03456跨聲速0.850.850.850.80666318.523.024.626.03.1網(wǎng)格影響分析本文網(wǎng)格采用了H-O型網(wǎng)格拓撲。為了分析網(wǎng)格影響,相對粗網(wǎng)格,對三角翼前、后緣,支架頂點和物面法向區(qū)域網(wǎng)格進行加密處理。圖2為密網(wǎng)格的拓撲結(jié)構(gòu)、物面網(wǎng)格和空間網(wǎng)格切片,部分網(wǎng)格分布參數(shù)見表2。圖2三角翼拓撲結(jié)構(gòu)及計算網(wǎng)格示意圖(finegrid)Fig.2Meshandgridtopologyaroundthedeltawing表2計算網(wǎng)格量對比Table2Comparisonofcomputationalgridsize網(wǎng)格類型網(wǎng)格量(×106)壁面距離y+三角翼上網(wǎng)格分布展向流向法向coarse≈7.01.0×10-5cr2.1129165129fine≈13.90.5×10-5cr1.0219229161463第4期張冬云等:Pω增強型k-ω湍流模型在三角翼旋渦流動的應(yīng)用
進行網(wǎng)格影響分析選用的計算狀態(tài)為:Ma∞=0.85,α=23.0°,Re∞=6×106。在文獻[12]中,NASANTF風洞試驗數(shù)據(jù)顯示三角翼上翼面并未發(fā)生主渦破裂現(xiàn)象,該迎角與發(fā)生渦破裂的臨界迎角接近,而文獻[12]的CFD方法預(yù)測出了渦破裂現(xiàn)象。圖3為選用Pω增強型k-ω湍流模型,粗、密網(wǎng)格的CFD計算結(jié)果對比,表明粗、密網(wǎng)格均預(yù)測出了三角翼上翼面的主渦峰值和二次渦峰值,前三個站位處,密網(wǎng)格預(yù)測的主渦和二次渦峰值均高于粗網(wǎng)格,更接近于試驗數(shù)據(jù)。x/cr=0.8站位處,粗網(wǎng)格預(yù)測的主渦峰值及站位與試驗結(jié)果吻合良好,但沿展向外側(cè)壓力分布均低于試驗數(shù)據(jù),二次渦峰值不明顯;密網(wǎng)格預(yù)測的主渦和二次渦峰值均高于試驗數(shù)據(jù)。x/cr=0.95站位處,粗網(wǎng)格預(yù)測的主渦峰值與試驗數(shù)據(jù)吻合較好,但主渦形態(tài)較為平滑,密網(wǎng)格預(yù)測的主渦峰值偏高,但主渦形態(tài)與試驗數(shù)據(jù)接近(即沿展向方向,在主渦外側(cè)的壓力系數(shù)曲線更為陡峭)。上述對比表明,兩套網(wǎng)格均較為準確地預(yù)測出了主渦和二次渦結(jié)構(gòu),網(wǎng)格密度達到研究要求。為了更為準確地預(yù)測繞三角翼旋渦流動中的分離渦、剪切層、激波與渦的干擾現(xiàn)象等,本文后續(xù)的研究均采用密網(wǎng)格。(a)x/cr=0.2(b)x/cr=0.4(c)x/cr=0.6(d)x/cr=0.8(e)x/cr=0.95(f)spanwisestations圖3不同網(wǎng)格密度在不同展向站位物面壓力系數(shù)曲線對比Fig.3Comparisonofsurfacepressuredistribu
【參考文獻】:
期刊論文
[1]繞跨聲速三角翼的激波/渦干擾流場數(shù)值模擬[J]. 李喜樂,楊永,張強,夏貞鋒. 航空學(xué)報. 2013(04)
[2]高階精度格式WCNS在三角翼大攻角模擬中的應(yīng)用研究[J]. 王光學(xué),鄧小剛,劉化勇,王運濤. 空氣動力學(xué)學(xué)報. 2012(01)
本文編號:3380325
【文章來源】:空氣動力學(xué)學(xué)報. 2016,34(04)北大核心CSCD
【文章頁數(shù)】:8 頁
【部分圖文】:
圖165°后掠三角翼模型[12]Fig.1Geometryofthe65°sweptleadingedgedeltawing
表1計算狀態(tài)對比Table1CFDcaseschosenforsummarycomparison計算狀態(tài)Ma∞Re∞(×106)α(°)12亞聲速0.40.46618.523.03456跨聲速0.850.850.850.80666318.523.024.626.03.1網(wǎng)格影響分析本文網(wǎng)格采用了H-O型網(wǎng)格拓撲。為了分析網(wǎng)格影響,相對粗網(wǎng)格,對三角翼前、后緣,支架頂點和物面法向區(qū)域網(wǎng)格進行加密處理。圖2為密網(wǎng)格的拓撲結(jié)構(gòu)、物面網(wǎng)格和空間網(wǎng)格切片,部分網(wǎng)格分布參數(shù)見表2。圖2三角翼拓撲結(jié)構(gòu)及計算網(wǎng)格示意圖(finegrid)Fig.2Meshandgridtopologyaroundthedeltawing表2計算網(wǎng)格量對比Table2Comparisonofcomputationalgridsize網(wǎng)格類型網(wǎng)格量(×106)壁面距離y+三角翼上網(wǎng)格分布展向流向法向coarse≈7.01.0×10-5cr2.1129165129fine≈13.90.5×10-5cr1.0219229161463第4期張冬云等:Pω增強型k-ω湍流模型在三角翼旋渦流動的應(yīng)用
進行網(wǎng)格影響分析選用的計算狀態(tài)為:Ma∞=0.85,α=23.0°,Re∞=6×106。在文獻[12]中,NASANTF風洞試驗數(shù)據(jù)顯示三角翼上翼面并未發(fā)生主渦破裂現(xiàn)象,該迎角與發(fā)生渦破裂的臨界迎角接近,而文獻[12]的CFD方法預(yù)測出了渦破裂現(xiàn)象。圖3為選用Pω增強型k-ω湍流模型,粗、密網(wǎng)格的CFD計算結(jié)果對比,表明粗、密網(wǎng)格均預(yù)測出了三角翼上翼面的主渦峰值和二次渦峰值,前三個站位處,密網(wǎng)格預(yù)測的主渦和二次渦峰值均高于粗網(wǎng)格,更接近于試驗數(shù)據(jù)。x/cr=0.8站位處,粗網(wǎng)格預(yù)測的主渦峰值及站位與試驗結(jié)果吻合良好,但沿展向外側(cè)壓力分布均低于試驗數(shù)據(jù),二次渦峰值不明顯;密網(wǎng)格預(yù)測的主渦和二次渦峰值均高于試驗數(shù)據(jù)。x/cr=0.95站位處,粗網(wǎng)格預(yù)測的主渦峰值與試驗數(shù)據(jù)吻合較好,但主渦形態(tài)較為平滑,密網(wǎng)格預(yù)測的主渦峰值偏高,但主渦形態(tài)與試驗數(shù)據(jù)接近(即沿展向方向,在主渦外側(cè)的壓力系數(shù)曲線更為陡峭)。上述對比表明,兩套網(wǎng)格均較為準確地預(yù)測出了主渦和二次渦結(jié)構(gòu),網(wǎng)格密度達到研究要求。為了更為準確地預(yù)測繞三角翼旋渦流動中的分離渦、剪切層、激波與渦的干擾現(xiàn)象等,本文后續(xù)的研究均采用密網(wǎng)格。(a)x/cr=0.2(b)x/cr=0.4(c)x/cr=0.6(d)x/cr=0.8(e)x/cr=0.95(f)spanwisestations圖3不同網(wǎng)格密度在不同展向站位物面壓力系數(shù)曲線對比Fig.3Comparisonofsurfacepressuredistribu
【參考文獻】:
期刊論文
[1]繞跨聲速三角翼的激波/渦干擾流場數(shù)值模擬[J]. 李喜樂,楊永,張強,夏貞鋒. 航空學(xué)報. 2013(04)
[2]高階精度格式WCNS在三角翼大攻角模擬中的應(yīng)用研究[J]. 王光學(xué),鄧小剛,劉化勇,王運濤. 空氣動力學(xué)學(xué)報. 2012(01)
本文編號:3380325
本文鏈接:http://sikaile.net/kejilunwen/lxlw/3380325.html
最近更新
教材專著