基于Fluent的底部排氣彈減阻特性數(shù)值模擬

朱越亭,馬鐵華,張紅艷,沈大偉

(中北大學(xué) 儀器科學(xué)與動態(tài)測試教育部重點實驗室， 太原 030051)

底部排氣減阻增程是從曳光彈射擊試驗中得到啟發(fā)的，其原理如下：彈丸飛行時，在彈丸后部形成低壓區(qū)，產(chǎn)生底部阻力，在彈丸底部附加一個排氣裝置，向彈底低壓區(qū)排入質(zhì)量和能量，通過提高底壓來減阻增程。底排減阻的空氣動力學(xué)研究已經(jīng)有將近70年的歷史，底部排氣彈的研制也有將近50年[1]。

目前，數(shù)值模擬方法是國內(nèi)外研究底排減阻技術(shù)的主要手段之一。國內(nèi)方面，史金光等[2]運用滑移網(wǎng)格技術(shù)進(jìn)行底部排氣彈三維流場的數(shù)值模擬；卓長飛等[3-4]對超聲速底部排氣彈底排真實氣體流場以及底部排氣彈底排真實氣體底部二次燃燒現(xiàn)象進(jìn)行了數(shù)值模擬；駱曉臣等[5]采用計算空氣動力學(xué)和質(zhì)點外彈道學(xué)耦合的方法求解底部排氣彈的飛行彈道；謝利平等[6]建立了底排內(nèi)彈道和外彈道模型，分析了底排裝置結(jié)構(gòu)參數(shù)和藥劑燃速系數(shù)等對底部排氣彈減阻增程效果的影響；羅盟等[7]采用標(biāo)準(zhǔn)k-ω湍流模型和有限速率化學(xué)反應(yīng)燃燒模型對二維軸對稱雷諾平均N-S方程進(jìn)行了數(shù)值計算。國外方面，Subbareddy P等[8]和Shin J R等[9]先后采用分離渦數(shù)值模擬的方法對超聲速飛行的底部排氣彈進(jìn)行了數(shù)值模擬；Charles J Nietubica[10]采用層流流動模型和層流燃燒有限速率基元反應(yīng)模型，模擬了M864彈的二維軸對稱流場；LEE Y K等[11]通過數(shù)值模擬帶孔后體模型在馬赫數(shù)為2.47時底部排氣的流動特征，發(fā)現(xiàn)對應(yīng)于最大的底排壓力有一個使底阻最小的底排條件；Petri Kaurinkoski[12]和Jeong-Yeol Chor[13]分別采用k-ε、k-ω湍流流動模型和層流燃燒有限速率基元反應(yīng)模型，模擬了155 mm彈的三維、二維流場；Choi J Y等[14]對含外部燃燒和中心推進(jìn)噴射的底部排氣彈的減阻特性進(jìn)行了數(shù)值模擬。

區(qū)別于以往數(shù)值計算所采用的軟件和方法，本研究以155mm彈為物理模型，運用前處理軟件Gambit建立彈丸模型并進(jìn)行網(wǎng)格劃分，運用Fluent軟件進(jìn)行仿真計算，對有無底排的彈丸流場及其氣動力參數(shù)進(jìn)行分析，研究不同來流馬赫數(shù)和攻角下彈丸的阻力系數(shù)和升力系數(shù)的變化規(guī)律，為底部排氣彈的研究提供理論參考。

1 模擬方法

以155 mm彈為物理模型，彈丸由尖拱形母線的彈頭部、圓柱彈身以及船尾組成，運用前處理軟件Gambit建立彈丸模型和劃分網(wǎng)格，彈丸表面網(wǎng)格的分布如圖1所示。

圖1 彈丸表面網(wǎng)格分布

彈丸的計算域網(wǎng)格分布如圖2所示，在運用Gambit劃分網(wǎng)格時，Equisize Skew通過單元的大小計算歪斜度，絕大多數(shù)網(wǎng)格控制在0.4以內(nèi)，以滿足計算要求。設(shè)置壁面、對稱面和壓力遠(yuǎn)場的邊界條件：其中壁面要求選擇的面為彈體的壁面；對稱面選擇為計算域的截面和彈體的截面；壓力遠(yuǎn)場選擇大計算域除了界面的另外三個面。

圖2 計算域網(wǎng)格分布

網(wǎng)格模型確定后，將其導(dǎo)入Fluent軟件中，數(shù)值模擬的具體相關(guān)參數(shù)均在該軟件中修改。材料屬性選擇理想氣體ideal-gas，薩蘭德(sutherland)定律計算黏性。計算域外邊界采用壓力遠(yuǎn)場的邊界條件，大氣壓為101 325 Pa，馬赫數(shù)以及攻角的余弦和正弦值應(yīng)根據(jù)具體模型來確定。湍流參數(shù)中的湍流黏度比(Turbulent Viscosity Ratio，TVR)保持為10。在大馬赫數(shù)下柯朗數(shù)(Courant Number)為0.1-0.5，在小馬赫數(shù)下為默認(rèn)值；Flux Type在小馬赫數(shù)情況下為Roe-FDS，在大馬赫數(shù)的情況下為AUSM。選擇二階迎風(fēng)格式。

2 數(shù)值模擬結(jié)果與分析

2.1 無底排氣動特性數(shù)值模擬

無底排狀態(tài)下，來流馬赫數(shù)分別為1.5、2.0、2.5、3.0、3.5，攻角分別為0，2，4，6(單位度)共計20種狀態(tài)。馬赫數(shù)為2.0，攻角為0的彈丸密度和壓力等值線圖如圖3所示。

圖3 Ma=2.0，α=0°的等值線

彈丸超音速飛行時，氣流在彈頭部產(chǎn)生斜激波，氣流經(jīng)過激波后沿彈身向下流動，當(dāng)氣流到達(dá)彈丸船尾角時，由于彈尾物面的折轉(zhuǎn)，氣流開始膨脹，產(chǎn)生膨脹波。當(dāng)氣流達(dá)到彈丸底部時，物面再一次迅速折轉(zhuǎn)，氣流再一次膨脹，產(chǎn)生劇烈的膨脹波。氣流繼續(xù)向下游流動，并逐漸恢復(fù)來流方向，在彈丸的尾跡區(qū)由于氣流的壓縮產(chǎn)生再壓縮激波。

由密度和壓力的大小分布得知，氣流通過激波后壓力與密度迅速增大，圖示彈頭部呈現(xiàn)紅色。氣流繼續(xù)朝彈丸的底部流動,由于彈底部氣體的低速回流，并且無外流的補充，產(chǎn)生一個低壓區(qū)，低壓區(qū)的存在導(dǎo)致彈丸底阻的形成。

圖4為馬赫數(shù)1.5、3.5，攻角0°時的等壓線圖。隨著來流馬赫數(shù)的增加，氣流壓縮變得更加強烈，彈體表面與激波面夾角逐漸減小，導(dǎo)致激波會逐漸貼近彈體表面。

當(dāng)來流有攻角時，以來流馬赫數(shù)2.0，攻角6°時為例，流場變化如圖5所示。彈丸迎風(fēng)面的激波面與物面夾角隨攻角的增大而減小，彈丸背風(fēng)面的激波面與物面夾角隨攻角的增大而增大。此外，彈丸底部流場在攻角變化時，迎風(fēng)面的尾跡區(qū)受到的擠壓隨攻角增大而越明顯，背風(fēng)面的尾跡區(qū)向下擠壓迎風(fēng)面尾跡區(qū)的程度隨攻角增大而增強。

彈丸的阻力系數(shù)(Cd)和升力系數(shù)(Cl)隨來流攻角及馬赫數(shù)變化曲線如圖6所示。

圖4 不同馬赫數(shù)的等壓線

圖5 Ma=2.0，α=6°的等壓線

圖6 阻力系數(shù)和升力系數(shù)變化曲線

隨來流攻角的增大，彈丸的阻力系數(shù)呈拋物線性增長趨勢。在同一攻角下，彈丸的阻力系數(shù)隨來流馬赫數(shù)的增大不斷減小，呈拋物線性減小趨勢。

彈丸升力系數(shù)隨著攻角的增大呈線性增長的趨勢。在同一來流攻角下：當(dāng)攻角為0°時，彈丸升力系數(shù)基本不隨來流馬赫數(shù)的變化而變化，維持為0。當(dāng)彈丸有攻角且攻角比較小的情況下如α=2°時，升力系數(shù)隨馬赫數(shù)的增加呈減小趨勢，但變化不是很明顯，當(dāng)攻角相對較大的情況下如α=6°時，其升力系數(shù)隨來流馬赫數(shù)的減小更為明顯。

2.2 有底排氣動特性數(shù)值模擬

在無底排的基礎(chǔ)上修改邊界條件，簡化物理模型，通過在彈底排出速度v=25 m/s的氣體，分析對底部流場、氣動力系數(shù)和減阻的影響。

以來流馬赫數(shù)2.0，攻角0°時繞彈丸流場的壓力和密度等值線圖為例來說明有底部排氣時的彈丸流場特性，如圖7。

圖7 Ma=2.0，α=0°的等值線

加入底部排氣后，彈丸前部流場結(jié)構(gòu)變化不大，彈丸底部折轉(zhuǎn)處的膨脹波明顯減弱。這主要是由于底排氣體的補充，剪切層被抬高，氣流流經(jīng)彈底時膨脹波減弱，折轉(zhuǎn)角減小。此外，底排高能氣體填充底部低壓回流區(qū)，將彈丸底部尾跡區(qū)明顯后推，導(dǎo)致再壓縮激波位置后移。密度等值線中，由于補充進(jìn)入的氣體，彈丸底部低壓區(qū)的剪切層被抬高，氣流流經(jīng)尾部時的折轉(zhuǎn)角明顯減小。彈丸底部補充的氣體改變了彈丸底部的流場結(jié)構(gòu)。

來流馬赫數(shù)變化時，以馬赫數(shù)1.5和2.5為例，彈丸流場變化情況如圖8所示。氣流受到的壓縮作用隨著馬赫數(shù)的增加而更加強烈，產(chǎn)生的激波角慢慢減小，尾跡區(qū)的再壓縮激波和彈尾的膨脹波的變化趨勢也相同。彈丸周圍流場的壓力和密度相較低馬赫數(shù)也有較大提高。彈丸底部受到底排氣體的影響，尾跡區(qū)都向后推移，根據(jù)等值線圖可以看出，在高馬赫數(shù)的情況下，彈丸底排氣體于尾跡區(qū)的后推效果要好于低馬赫數(shù)的情況。

當(dāng)來流攻角變化時，以來流馬赫數(shù)2.0，攻角6°時的流場圖為例，如圖9所示。彈丸迎風(fēng)面的激波角減小，背風(fēng)面的激波角增大，在彈丸有攻角時彈丸尾跡區(qū)流場結(jié)構(gòu)不再對稱，迎風(fēng)面尾跡區(qū)域被背風(fēng)面尾跡區(qū)域擠壓明顯，背風(fēng)面底部流場由于受到外流的擠壓讓整個彈丸底部排氣區(qū)域向迎風(fēng)面的底部區(qū)域靠近。背風(fēng)面彈丸尾跡區(qū)后推作用相對于迎風(fēng)面被削弱。

圖8 不同馬赫數(shù)的等壓線

圖9 Ma=2.0，α=6°的等壓線

彈丸的阻力系數(shù)(Cd)和升力系數(shù)(Cl)隨來流攻角及馬赫數(shù)的變化曲線如圖10所示。

圖10 阻力系數(shù)和升力系數(shù)變化曲線

彈丸的阻力系數(shù)隨攻角的增大呈拋物線性增長，隨著來流馬赫數(shù)的增大呈拋物線性減小。與無底部排氣時彈丸阻力系數(shù)的變化趨勢相同。

彈丸升力系數(shù)隨攻角增大呈線性增長，0°攻角時升力系數(shù)不隨馬赫數(shù)增大而改變，均維持在0左右。攻角增大后，升力系數(shù)隨馬赫數(shù)的增加呈現(xiàn)出線性減小趨勢，但變化趨勢不是很明顯。

2.3 有無底排情況計算結(jié)果對比

圖11為任取同一來流馬赫數(shù)下，加入底排前后彈丸氣動力諸系數(shù)隨來流攻角的變化。

彈丸升力系數(shù)都隨攻角的增大呈線性增長的趨勢，有無底排前后大小幾乎沒有發(fā)生變化，這是因為彈丸的底部阻力是由于彈丸底部所受到壓力遠(yuǎn)低于來流導(dǎo)致的，并沒有直接作用于彈體上,對于彈丸的法向力和軸向力的影響很小。彈丸的全彈阻力系數(shù)有較大改變，加入底排后，彈丸的阻力系數(shù)大幅下降，說明底排減阻效果明顯。

作者還研究了有無底排狀態(tài)下的阻力系數(shù)。將加入底排前后的阻力系數(shù)整理比較之后，得到阻力系數(shù)的最高減阻率達(dá)到34.8%，最低達(dá)到14.8%。上面兩個數(shù)據(jù)可以說明，底排氣體減小了彈丸的底阻，使得彈丸的整體阻力下降，射程增加。由此可見加入底排氣體確實可以降低彈丸阻力，從而實現(xiàn)增程的目的，驗證了理論的正確性。

圖11 底排前后的阻力系數(shù)和升力系數(shù)

3 結(jié)論

1) 超音速來流下，阻力系數(shù)隨來流攻角的增大呈拋物線性增大，隨來流馬赫數(shù)的增大呈拋物線性減小。

2) 超音速來流下，升力系數(shù)隨來流攻角的增大線性增大，隨來流馬赫數(shù)的增大而減小，攻角越大減小幅度越明顯。

3) 彈丸在加入底排前后，升力系數(shù)無明顯變化，阻力系數(shù)大幅下降，達(dá)到減阻效果，且阻力系數(shù)的最高減阻率達(dá)34.8%，結(jié)果符合空氣動力學(xué)規(guī)律，驗證了底部排氣彈減阻增程的理論依據(jù)，為底部排氣彈的研究提供參考。