王運濤, 孫巖, 孟德虹,*, 王光學
1.中國空氣動力研究與發(fā)展中心 計算空氣動力研究所, 綿陽 621000 2.中國空氣動力研究與發(fā)展中心 空氣動力學國家重點實驗室, 綿陽 621000
CRM翼/身/平尾組合體模型高階精度數(shù)值模擬
王運濤1, 孫巖2, 孟德虹1,*, 王光學1
1.中國空氣動力研究與發(fā)展中心 計算空氣動力研究所, 綿陽 621000 2.中國空氣動力研究與發(fā)展中心 空氣動力學國家重點實驗室, 綿陽 621000
基于五階空間離散精度的WCNS格式,開展了CRM翼/身/平尾組合體模型的高階精度數(shù)值模擬,以評估WCNS格式對復雜外形的模擬能力以及典型運輸機巡航構(gòu)型阻力預(yù)測方面的精度。首先依照DPW組委會提出的網(wǎng)格生成指導原則,利用ICEM軟件生成了粗、中、細3套網(wǎng)格,對應(yīng)的y+分別為1、2/3和4/9。通過對CRM模型的計算和分析,研究了網(wǎng)格規(guī)模對氣動特性、壓力分布和翼根后緣局部分離區(qū)的影響。通過與DPW4統(tǒng)計結(jié)果和部分實驗結(jié)果的對比,高階精度數(shù)值模擬結(jié)果表明,本文的氣動特性計算結(jié)果與統(tǒng)計平均結(jié)果吻合較好;網(wǎng)格密度對機翼內(nèi)側(cè)的展向站位壓力分布影響較小,對機翼外側(cè)展向壓力分布影響較大;網(wǎng)格密度對翼根后緣局部分離區(qū)略有影響。
RANS方程; WCNS格式; CRM模型; 流場模擬; 網(wǎng)格密度; 氣動特性
AIAA阻力預(yù)測會議DPW(Drag Prediction Workshop)由AIAA應(yīng)用空氣動力學技術(shù)委員會發(fā)起,旨在評估基于雷諾平均Navier-Stokes (RANS)方程的CFD方法在預(yù)測典型運輸機構(gòu)型氣動特性尤其是阻力預(yù)測方面的現(xiàn)狀,建立一個評估CFD可信度的國際平臺。DPW從發(fā)起到現(xiàn)在,已經(jīng)成功舉辦了5屆[1-5],持續(xù)了10多年時間,獲得了全世界近30個科研機構(gòu)和工業(yè)部門的廣泛參與,DPW系列會議已經(jīng)成為CFD驗證與確認發(fā)展歷程中最重要的一次國際合作。DPW系列會議的意義在于積累了豐富的CFD數(shù)據(jù)和實驗數(shù)據(jù),并引發(fā)了關(guān)于CFD技術(shù)的思考和討論,這對于評估CFD現(xiàn)狀,認識CFD當前的問題,指導CFD的發(fā)展具有重要意義。
第4屆DPW會議于2009年6月在美國德克薩斯州的圣地亞哥市召開,這次會議采用了一個全新的外形CRM (Common Research Model),CRM模型[6]代表了當代一類典型的高性能運輸機構(gòu)型。來自世界各地的17家研究機構(gòu)采用16種CFD軟件提供了30組計算結(jié)果[4],本文作者在前期工作中采用TRIP軟件也開展了CRM構(gòu)型的數(shù)值模擬[7]。以上這些基于RANS方程的計算結(jié)果均采用了二階精度的計算方法,采用三階以上差分格式的數(shù)值模擬結(jié)果尚未見公開報道。
高階精度格式一直是CFD領(lǐng)域的研究熱點,但其在復雜外形上的應(yīng)用才剛剛起步[8]。五階空間離散精度的WCNS(Weighted Compact Nonlinear Scheme)格式是鄧小剛和張涵信提出并發(fā)展起來的加權(quán)緊致非線性高階格式[9],具有耗散低、分辨率高的特點。近年來,通過在幾何守恒律方面持續(xù)不斷的研究工作[10],已經(jīng)成功應(yīng)用于典型運輸機構(gòu)型高階精度的數(shù)值模擬中[11-13],并取得了良好效果。
本文采用五階空間離散精度的WCNS格式對CRM翼/身/平尾(CRM-WBH)組合體模型開展了高階精度數(shù)值模擬,主要研究網(wǎng)格規(guī)模對高階精度數(shù)值模擬結(jié)果的影響,進一步確認WCNS格式模擬典型運輸機復雜構(gòu)型的能力。
CRM模型由NASA的亞聲速固定翼(Subsonic Fixed Wing, SFW)空氣動力技術(shù)研究小組和DPW的組織委員會合作設(shè)計開發(fā),主要目的是為CFD的驗證和確認工作提供基準外形[6],CRM模型的設(shè)計馬赫數(shù)Ma=0.85,升力系數(shù)CL=0.50。CRM-WBH計算構(gòu)型見圖1,平尾偏角為0°,計算外形的基本參數(shù)見表1。其中,Sref、Cref、xref、yref及zref分別表示計算外形的參考面積、平均氣動弦長以及3個坐標方向上的力矩參考點,b、ΛC/4、AR分別表示計算構(gòu)型的展長、1/4弦長處的后掠角以及展弦比。
圖1 CRM翼/身/平尾組合體外形Fig.1 CRM wing/body/horizontal tail configuration
表1 CRM-WBH模型幾何參數(shù)Table 1 Geometrical parameters of CRM-WBH model
注:1 in=2.54 cm。
本文研究采用有限差分方法離散任意坐標系下的RANS方程組,控制方程的對流項離散采用五階精度的WCNS格式,黏性項的離散采用六階精度中心格式,邊界及近邊界條件采用單邊四階精度離散,以上方法的詳細介紹見文獻[9-13];湍流模型采用Menter剪切應(yīng)力輸運(SST)兩方程模型[14],離散方程組的求解采用塊LU分解方法(Block Lower-Upper Symmetric Gauss-Seidel Scheme)[15-16]。
為了降低網(wǎng)格差異對計算結(jié)果分布的影響,DPW4組織委員會給出了網(wǎng)格生成的指導原則[4],對網(wǎng)格規(guī)模、物面第一層網(wǎng)格高度、邊界層網(wǎng)格增長率、遠場邊界距離等網(wǎng)格參數(shù)進行了約定,并提供了結(jié)構(gòu)、非結(jié)構(gòu)基準網(wǎng)格。
由于高階精度格式對計算網(wǎng)格質(zhì)量要求更高,本文研究中并沒有直接采用DPW4 組委會提供的計算網(wǎng)格,而是根據(jù)網(wǎng)格生成指導原則,采用ICEM軟件重新生成了不同規(guī)模的粗(Coarse)、中(Medium)、細(Fine)3套網(wǎng)格。為了提高邊界層流動模擬的精度,在機身、機翼、平尾上分別單獨包裹了一層O型貼體網(wǎng)格,粗網(wǎng)格的表面網(wǎng)格如圖2所示。3套網(wǎng)格的詳細信息見表2。其中,Nnum為總的網(wǎng)格節(jié)點數(shù);dy為第一層網(wǎng)格厚度;nBL、λBL分別為邊界層網(wǎng)格數(shù)量和網(wǎng)格增長率;nwing_chord、nwing_span和nhorizontal_chord、nhorizontal_span分別為機翼弦向、展向網(wǎng)格節(jié)點數(shù)和平尾弦向、展向網(wǎng)格節(jié)點數(shù);Lfar為計算網(wǎng)格遠場離模型的距離;Nblock為多塊網(wǎng)格包括的計算塊數(shù)量。
圖2 CRM-WBH模型表面網(wǎng)格(粗網(wǎng)格)Fig.2 Surface grid of CRM-WBH model (coarse grid)
表2 CRM-WBH模型網(wǎng)格參數(shù)Table 2 Grid parameters of CRM-WBH model
采用粗、中、細3套網(wǎng)格和高階精度計算方法,開展了固定升力系數(shù)CL下網(wǎng)格密度對計算結(jié)果的影響研究,主要包括氣動力系數(shù)、典型站位壓力系數(shù)和翼身結(jié)合部后緣局部分離區(qū)的大小等3個方面,計算采用全湍流方式。來流條件為:Ma=0.85,Re=5.0×106,CL=0.500±0.001。
1) 氣動力系數(shù)
表3給出了固定升力系數(shù)下(CL=0.500±0.001),采用粗、中、細3套網(wǎng)格得到的CRM-WBH組合體的迎角α、阻力系數(shù)CD、壓差阻力系數(shù)CDp、摩擦阻力系數(shù)CDf和俯仰力矩系數(shù)Cm。
表3CRM-WBH模型的氣動特性(CL=0.500±0.001)
Table3AerodynamiccharacteristicsofCRM-WBHmodel(CL=0.500±0.001)
Gridα/(°)CDCDpCDfCmCoarse2.2580.026810.014660.01215-0.0421Medium2.3110.026810.014520.01229-0.0390Fine2.3440.026910.014500.01241-0.0369
從表3可以看出,采用3套網(wǎng)格得到的固定升力系數(shù)下的迎角α、阻力系數(shù)分量、俯仰力矩系數(shù)均隨網(wǎng)格規(guī)模而單調(diào)變化,而總的阻力系數(shù)基本不變。Morrison[17]對DPW4提供的27組結(jié)果中的21組“核心解”(Core Solution)進行了統(tǒng)計分析,粗網(wǎng)格與密網(wǎng)格總阻力系數(shù)差量為-0.000 7,本文得到的總阻力系數(shù)差量為0.000 1。在統(tǒng)計平均的意義下,這說明采用高階精度方法得到的阻力系數(shù)對網(wǎng)格的依賴性要弱于DPW4中采用的各類二階精度方法。
表4給出了RE(Richardson Extrapolation)方法[4]得到的網(wǎng)格無關(guān)性結(jié)果,DPW組委會統(tǒng)計27組計算結(jié)果得到的相應(yīng)物理量平均值及其標準差,以及NASA蘭利國家跨聲速設(shè)備(NTF)的風洞試驗結(jié)果[4,18]。從表4可以看出,本文計算的網(wǎng)格無關(guān)性結(jié)果均位于文獻[4]的統(tǒng)計結(jié)果分布范圍內(nèi)。Rivers等[19-20]采用非結(jié)構(gòu)計算軟件USM3D 6.0研究了風洞試驗?zāi)P偷闹螜C構(gòu)和模型的靜氣動彈性對數(shù)值模擬結(jié)果的影響,研究表明:計算模型中考慮風洞模型支撐裝置和靜氣動彈性變化可以顯著提高計算結(jié)果與試驗結(jié)果的吻合程度。以2°迎角為例,計算模型中考慮支撐裝置和靜氣動彈性影響使得阻力系數(shù)下降0.003 7,俯仰力矩系數(shù)增加0.035 9,升力系數(shù)下降0.026 0。除了模型支撐裝置和靜氣動彈性外,導致氣動力計算結(jié)果與試驗結(jié)果之間差異的其他原因還包括轉(zhuǎn)捩位置和洞壁影響等。
表4 CRM-WBH模型網(wǎng)格無關(guān)性結(jié)果Table 4 Grid-converged results of CRM-WBH model
2) 表面壓力系數(shù)
圖3給出了粗、中、細3套不同網(wǎng)格規(guī)模下,CRM-WBH模型機翼3個典型展向位置剖面的壓力系數(shù)分布曲線以及相應(yīng)的試驗測壓結(jié)果。從圖3可以看出,從靠近翼根(展向站位系數(shù)η=0.115)的位置直到機翼中部(η=0.502)的位置,網(wǎng)格規(guī)模的變化對壓力系數(shù)的分布基本沒有影響,計算結(jié)果與試驗結(jié)果之間吻合良好。從機翼中部(η=0.502)的位置直到翼梢(η=0.950),網(wǎng)格規(guī)模的變化對壓力系數(shù)的影響主要體現(xiàn)在上翼面激波的位置,計算結(jié)果與試驗結(jié)果的差別主要是由于試驗結(jié)果包含了靜氣動變形的影響,而計算模型采用了剛性模型。
圖3 CRM-WBH模型典型站位壓力系數(shù)分布 Fig.3 Pressure coefficients distribution at typical spanwise locations of CRM-WBH model
3) 局部分離區(qū)
圖4給出了不同網(wǎng)格規(guī)模下翼身結(jié)合位置物面流線的變化。表5給出了不同網(wǎng)格規(guī)模下分離泡的尺寸變化及利用Richardson外插方法得到的網(wǎng)格無關(guān)結(jié)果,其中ΔBL、ΔFS分別為分離泡的展向?qū)挾扰c縱向?qū)挾?。從圖4和表5可以看出,隨著網(wǎng)格規(guī)模的增加,預(yù)測的分離泡尺寸也逐漸增加,但分離泡尺寸增加的速度逐漸減慢。
圖4 CRM-WBH模型表面極限流線(局部) Fig.4 Streamlines on upper surface of CRM-WBH model (local)
表5 不同網(wǎng)格下分離泡尺寸Table 5 Separation bubble dimension with different grids
1) 在固定升力系數(shù)下,采用粗、中、細3套網(wǎng)格和高階精度計算方法得到了具有網(wǎng)格收斂性的氣動力結(jié)果。
2) 氣動特性高精度數(shù)值模擬結(jié)果的外差值與DPW4的統(tǒng)計結(jié)果平均值吻合較好。
3) 網(wǎng)格規(guī)模對機翼外側(cè)上表面的激波位置和翼身結(jié)合部后緣局部分離區(qū)大小略有影響。
計算結(jié)果與試驗結(jié)果之間的差別需要從支撐干擾、靜氣動彈性影響等方面進一步開展研究。
致 謝
感謝中國空氣動力研究與發(fā)展中心張玉倫、洪俊武、張書俊、李偉、楊小川等同志在高階精度格式程序?qū)崿F(xiàn)方面所作的研究工作,感謝中國航空研究院白文博士在數(shù)據(jù)分析方面提供的幫助。
[1] LEVY D W, VASSBERG J C, WAHLS R A, et al. Summary of data from the first AIAA CFD drag prediction workshop[J]. Journal of Aircraft, 2003, 40(5): 875-882.
[2] LAFLIN K R, VASSBERG J C, WAHLS R A, et al. Summary of data from the second AIAA CFD drag prediction workshop[J]. Journal of Aircraft, 2005, 42(5): 1165-1178.
[3] VASSBERG J C, TINOCO E N, MANI M, et al. Abridged summary of the third AIAA CFD drag prediction workshop[J]. Journal of Aircraft, 2008, 45(3): 781-798.
[4] VASSBERG J C, TINOCO E N, MANI M, et al. Summary of the fourth AIAA computational fluid dynamics drag prediction workshop[J]. Journal of Aircraft, 2014, 51(4): 1070-1089.
[5] LEVY D W, LAFLIN K R, TINOCO E N, et al. Summary of data from the fifth computational fluid dynamics drag prediction workshop[J]. Journal of Aircraft, 2014, 51(4): 1194-1213.
[6] VASSBERG J C, DEHAAN M A, RIVERS S M, et al. Development of a common research model for applied CFD validation studies: AIAA-2008-6919[R]. Reston: AIAA, 2008.
[7] 王運濤, 張書俊, 孟德虹. DPW4翼/身/平尾組合體的數(shù)值模擬[J]. 空氣動力學學報, 2013, 31(6): 739-744.
WANG Y T, ZHANG S J, MENG D H. Numerical simulation and study for DPW4 wing/body/tail[J]. Acta Aerodynamica Sinica, 2013, 31(6): 739-744 (in Chinese).
[8] SLOTNICK J, KHODADOUST A, ALONSO J, et al. CFD vision 2030 study: A path to revolutionary computational aerosciences: NASA/CR-2014-218178[R]. Washington, D.C.: NASA, 2014.
[9] DENG X G, ZHANG H X. Developing high-order werighted compact nonlinear schemes[J]. Journal of Computational Physics, 2000, 165: 24-44.
[10] DENG X G, MIN R B, MAO M L, et al. Further studies on geometric conservation law and application to high-order finite difference scheme with stationary grid[J]. Journal of Computational Physics, 2013, 239: 90-111.
[11] 王運濤, 孫巖, 王光學, 等. DLR-F6翼身組合體的高階精度數(shù)值模擬[J]. 航空學報, 2015, 36(9): 2923-2929.
WANG Y T, SUN Y, WANG G X, et al. High-order numerical simulation of DLR-F6 wing-body configuration[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(9): 2923-2929 (in Chinese).
[12] WANG Y T, ZHANG Y L, LI S, et al. Calibration of aγ-Reθtransition model and its validation with high-order numerical method[J]. Chinese Journal of Aeronautics, 2015, 28(3): 704-711.
[13] 王運濤, 孟德虹, 孫巖, 等. DLR-F6/FX2B翼身組合體構(gòu)型高階精度數(shù)值模擬[J]. 航空學報, 2016, 37(2): 484-490.
WANG Y T, MENG D H, SUN Y, et al. High-order numerical simulation of DLR-F6/FX2B wing-body configuration[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(2): 484-490 (in Chinese).
[14] MENTER F R. Two-equation eddy-viscosity turbulence models for engineering application[J]. AIAA Journal, 1994, 32(8): 1598-1605.
[15] CHEN R F, WANG Z J. Fast, block lower-upper symmetric Gauss-Seidel scheme for arbitrary grids[J]. AIAA Journal, 2000, 38(12): 2238-2245.
[16] 王光學, 張玉倫, 王運濤, 等. BLU-SGS方法在WCNS高階精度格式上的數(shù)值分析[J]. 空氣動力學學報, 2015, 33(6): 733-739.
WANG G X, ZHANG Y L, WANG Y T, et al. Numerical analysis of BLU-SGS method in WCNS high-order scheme[J]. Acta Aerodynamica Sinica, 2015, 33(6): 733-739 (in Chinese).
[17] MORRISON J H. Statistical analysis of CFD solutions from the fourth AIAA drag prediction workshop: AIAA-2010-4673[R]. Reston: AIAA, 2010.
[18] RIVERS M B, DITTBERNER A. Experimental investigation of the NASA common research model(invited): AIAA-2010-4218[R]. Reston: AIAA, 2010.
[19] RIVERS M B, HUNTER C A. Support system effects on the NASA common research model: AIAA-2012-0707[R]. Reston: AIAA, 2012.
[20] RIVERS M B, HUNTER C A, CAMPBELL R L. Further investigation of he support system effects and wing twist on the NASA common research model: AIAA-2012-3209[R]. Reston: AIAA, 2012.
High-orderprecisionnumericalsimulationofCRMwing/body/horizontaltailmodel
WANGYuntao1,SUNYan2,MENGDehong1,*,WANGGuangxue1
1.ComputationalAerodynamicsInstitute,ChinaAerodynamicsResearchandDevelopmentCenter,Mianyang621000,China2.StateKeyLaboratoryofAerodynamics,ChinaAerodynamicsResearchandDevelopmentCenter,Mianyang621000,China
High-ordernumericalsimulationoncommonresearchmodel(CRM)wing/body/horizontaltailmodelispresentedwiththefifth-orderWCNSschemetoassesstheabilityofhigh-orderWCNSschemeoncomplexconfigurationsimulationandtheprecisionofpredicatingcruisedragoftransonicconfiguration.Threegrids(coarse,mediumandfine)arecreatedwithsoftwareICEMaccordingtothegriddingguidelinesprovidedbyDragPredictionWorkshop(DPW)organizingcommittee,andthey+ofthethreegridsare1,2/3and4/9.Computationandanalysisonthreegridsarecarriedouttoinvestigatethegrideffectonaerodynamiccharacteristics,pressuredistributionandthelocalseparationbubbleatthewingroottrailingedge.ComparedtothestatisticresultsfromDPW4andsomeexperimentaldatafromNTF,thehigh-ordernumericalresultsshowthatthecomputationalresultsofaerodynamiccharacteristicsagreewellwithstatisticdata;thegriddensityhaslittleinfluenceonpressuredistributionontheinboardstations,whereasithassomeinfluenceontheoutboardstations;thegriddensityhassomeinfluenceonthesizeofthelocalseparationbubbleatthewingroottrailingedge.
RANSequations;WCNSscheme;CRMmodel;flowfieldsimulation;griddensity;aerodynamiccharacteristics
2016-01-20;Revised2016-02-15;Accepted2016-03-14;Publishedonline2016-03-221541
URL:www.cnki.net/kcms/detail/11.1929.V.20160322.1541.002.html
NationalKeyResearchandDevelopmentProgram(2016YFB0200700)
2016-01-20;退修日期2016-02-15;錄用日期2016-03-14; < class="emphasis_bold">網(wǎng)絡(luò)出版時間
時間:2016-03-221541
www.cnki.net/kcms/detail/11.1929.V.20160322.1541.002.html
國家重點研發(fā)計劃 (2016YFB0200700)
*
.Tel.:0816-2463062E-mailmdh157@163.com
王運濤, 孫巖, 孟德虹, 等.CRM翼/身/平尾組合體模型高精度數(shù)值模擬J. 航空學報,2016,37(12):3692-3697.WANGYT,SUNY,MENGDH,etal.High-orderprecisionnumericalsimulationofCRMwing/body/horizontaltailmodelJ.ActaAeronauticaetAstronauticaSinica,2016,37(12):3692-3697.
http://hkxb.buaa.edu.cnhkxb@buaa.edu.cn
10.7527/S1000-6893.2016.0080
V211.7
A
1000-6893(2016)12-3692-06
王運濤男, 博士, 研究員, 博士生導師。主要研究方向: 計算空氣動力學。Tel.: 0816-2463037E-mail: ytwang@skla.cardc.cn
孫巖男, 博士研究生。主要研究方向: 計算流體力學。Tel.: 0816-7067915E-mail: supersunyan@163.com
孟德虹男, 碩士, 助理研究員。主要研究方向: 計算空氣動力學。Tel.: 0816-2463062E-mail: mdh157@163.com
王光學男, 碩士, 副研究員。主要研究方向: 計算空氣動力學和大攻角分離流。Tel.: 0816-2463062E-mail: wgx111@sina.com
*Correspondingauthor.Tel.:0816-2463062E-mailmdh157@163.com