趙明陽, 王 革，韓萬之，張 琦

(哈爾濱工程大學(xué) 航天與建筑工程學(xué)院， 哈爾濱 150001)

【航天工程】

基于接觸判斷方法的尾噴管堵蓋分離運動分析

趙明陽, 王 革，韓萬之，張 琦

(哈爾濱工程大學(xué) 航天與建筑工程學(xué)院， 哈爾濱 150001)

為研究堵蓋在初始分離及后續(xù)過程中與燃?xì)獾南嗷プ饔?，?yīng)用CFD和動網(wǎng)格技術(shù)，結(jié)合接觸判斷及多孔介質(zhì)模型，建立了堵蓋接觸分離及后續(xù)運動的物理模型。研究結(jié)果表明，堵蓋與高壓燃?xì)獾南嗷プ饔?，使堵蓋獲得較高的沖擊動能；同時相較于非接觸式預(yù)留縫隙的處理方法，接觸判斷耦合多孔介質(zhì)模型的方法，更準(zhǔn)確的模擬出堵蓋的分離延遲效應(yīng)。

尾噴管堵蓋；接觸分離運動；接觸判斷；多孔介質(zhì)

超壓破蓋技術(shù)應(yīng)用廣泛，且日趨完善。點火過程中經(jīng)由尾噴管堵蓋的阻塞作用，實現(xiàn)高壓燃?xì)獾目焖俳?，從而獲得對預(yù)設(shè)承壓標(biāo)準(zhǔn)的堵蓋的開啟。過程中的壓力匹配問題，高沖擊動能的堵蓋運動問題等，都是采用尾噴管超壓破蓋技術(shù)亟待解決的問題。堵蓋打開過程是整個運動過程的開始段，對堵蓋的初始運動影響明顯，因此有必要對打開過程做一個合適的處理。

尾噴管超壓破蓋技術(shù)的研究，國外起步較早，并有一些成熟的應(yīng)用，但大多集中于貯運發(fā)射箱的破蓋過程。國內(nèi)對此的研究相對滯后，分析過程以實驗與工程經(jīng)驗為主，尤其對于破蓋過程中的壓力匹配，缺乏定量的分析，嚴(yán)重影響了設(shè)置的精度。同時，對于堵蓋打開過程大都采用預(yù)留縫隙加約束力的粗狂處理方式，也影響了計算的精度。

本文通過數(shù)值計算，分別模擬了非接觸式預(yù)留縫隙方法及接觸判斷處理方法作用下尾噴管堵蓋受燃?xì)庾饔玫拇蜷_及飛出過程，通過檢測堵蓋前的壓力變化研究打開過程中沖擊波的作用，同時由堵蓋的運動規(guī)律，獲得其不同時刻的沖擊動能，為堵蓋的撞擊強(qiáng)度分析提供依據(jù)。

1 相關(guān)理論

1.1 流動控制方程

燃?xì)饬鲃涌刂品匠滩捎枚S非定常方程組，湍流方程采用標(biāo)準(zhǔn)k-ε方程，考慮燃?xì)馀c空氣的混合，添加組分輸運方程:

(1)

式中：Ri為化學(xué)反應(yīng)帶來的組分i凈增量，Si為源項。

在動網(wǎng)格計算區(qū)域內(nèi)，在任意控制體V中任意標(biāo)量的積分形式的控制方程可以表示為：

(2)

式中：φ為通量變量；ug為網(wǎng)格移動速度；Γ為耗散系數(shù)。

1.2 運動控制方程

物體的運動主要是由其初速度和加速度決定的，而加速度是由物體的受力決定的。對于所要研究的運動體有：

(3)

(4)

(5)

1.3 接觸判斷與多孔介質(zhì)模型

為了解決此類問題常規(guī)做法(預(yù)留縫隙)中的弊端，在此使用接觸判斷與多孔介質(zhì)模型相結(jié)合的方法實現(xiàn)計算域中運動體與固壁縫隙處流動狀態(tài)的轉(zhuǎn)換。

假設(shè)多孔介質(zhì)是各向同性的，則多孔介質(zhì)模型的應(yīng)用可以表達(dá)為在動量方程中添加了如下源項：

(6)

式中，D和C分別為黏性阻力系數(shù)和慣性阻力系數(shù)。

通過添加人工的粘性阻力和慣性阻力，可以使流體速度接近為零，作為近似的固體區(qū)域。

在實現(xiàn)計算的過程中，需要分別實現(xiàn)以下三個方面：

1) 采用接觸判斷的方法，以運動體與固壁的距離作為縫隙處流動狀態(tài)轉(zhuǎn)換的判據(jù)，即：

d≤dcon

(7)

其中，d為堵蓋壁面與噴管壁面的最小距離；dcon為預(yù)設(shè)的接觸距離。

2) 滿足判據(jù)條件，將縫隙區(qū)域處理為固體區(qū)，即令D、C為一無窮大數(shù)；

3) 判據(jù)條件失效，將縫隙區(qū)域處理為流體區(qū)，不再施加源項。

2 仿真計算

2.1 計算模型

計算模擬燃?xì)馄颇：蟠蜷_堵蓋及推動其飛出噴管的整個過程，關(guān)注堵蓋飛出后的沖擊特性，計算模型包括發(fā)動機(jī)尾噴管，堵蓋以及足夠堵蓋運動的環(huán)境區(qū)域(2 m×20 m)，如圖1所示。根據(jù)堵蓋的設(shè)計外形，將其簡化為帶凹坑的圓臺，并按質(zhì)量均勻分布來處理。由于堵蓋運動的時間短，速度大，將其運動簡化為單方向的平動。考慮模型的對稱性，采用二維軸對稱模型來減少計算量。

圖1 仿真計算物理模型

模型網(wǎng)格劃分如圖2所示，采用結(jié)構(gòu)化網(wǎng)格，在噴管區(qū)域進(jìn)行加密處理，縫隙處預(yù)留距離為4 mm，網(wǎng)格總數(shù)13萬。

圖2 仿真計算網(wǎng)格模型

2.2 計算條件

對實驗所測得的發(fā)動機(jī)壓力變化曲線進(jìn)行數(shù)值離散擬合出計算模型入口的總壓變化曲線如圖3所示，其它燃?xì)鈪?shù)見表1。環(huán)境邊界采用壓力遠(yuǎn)場邊界，溫度293 K，壓強(qiáng)為0.101 325 MPa。

根據(jù)實驗測得的破膜壓力為3.058 2 MPa，經(jīng)積分面積得到的堵蓋受力為68 000 N。因此仿真計算中初始時刻將防潮膜前定義為燃?xì)鈪^(qū)，膜后定義為空氣區(qū)，并將堵蓋受力達(dá)到68 000 N的時刻設(shè)置為堵蓋運動的初始時刻。

圖3 入口總壓變化曲線

項目單位數(shù)值燃?xì)獗葻岜取?．16188燃?xì)舛▔罕葻酛/(kg·K)1954．34燃?xì)鈿怏w常數(shù)J/(kg·K)272．32燃?xì)庀鄬Ψ肿淤|(zhì)量kg/kmol30．53燃?xì)饪倻豄3486

3 計算結(jié)果分析

3.1 破蓋過程結(jié)果分析

圖4所示為非接觸式預(yù)留縫隙方法與接觸判斷方法作用下，尾噴管堵蓋啟動前的壓力云圖。從圖中可以明顯的看出，非接觸式預(yù)留縫隙方法對于接觸縫隙的處理，使得燃?xì)庠诙律w啟動前已發(fā)生泄漏；而接觸判斷方法則避免了這一情況。

圖4 堵蓋啟動前壓力云圖

圖5所示為兩種方法下尾噴管堵蓋運動前的受力變化。從圖中可以清晰地看到兩種方法的受力變化趨勢一致，但接觸判斷方法的引入，使得破蓋時間加長，并且同時刻的受力明顯下降。

圖6所示為兩種方法下，圖1監(jiān)測點在相應(yīng)破蓋時間內(nèi)的壓力變化曲線。從圖中可以看出，兩種方法的壓力變化趨勢基本一致，但接觸判斷方法的引入，使得同時刻的壓力有所升高。

圖5 堵蓋運動前受力變化

圖6 堵蓋運動前監(jiān)測點壓力變化

分析圖4、圖5及圖6可知：接觸判斷方法的引入，克服了非接觸式預(yù)留縫隙方法中，燃?xì)庑孤┫扔诙律w運動的情況；通過構(gòu)造人工阻力項的方法，抵消了堵蓋側(cè)壁所受燃?xì)庾饔昧Φ挠绊?，這就造成了同時刻堵蓋受力略小，需要更長時間來實現(xiàn)超壓破蓋的過程。同時，由于人工阻力項的加入，阻礙了非接觸式預(yù)留縫隙方法中堵蓋與固壁縫隙處的流動，這也造成了同時刻堵蓋前壓力的增高。

3.2 堵蓋飛出過程結(jié)果分析

圖7，圖8以及圖9所示分別為兩種處理方法下堵蓋啟動位置，2.5 m位置以及5 m位置處的馬赫數(shù)及壓力云圖。從圖中可以明顯的看出，非接觸式預(yù)留縫隙方法對于縫隙的處理，使得啟動時刻的堵蓋后方已經(jīng)建立起一定尺度的壓力；同時，兩種處理方法在啟動時刻帶來的差異作用到后續(xù)結(jié)果中，表現(xiàn)為接觸判斷處理方法在噴管出口位置形成的激波串結(jié)構(gòu)減弱，馬赫盤基本消失。

圖7 堵蓋啟動位置的馬赫數(shù)及壓力云圖

圖8 2.5 m位置的馬赫數(shù)及壓力云圖

圖 9 5 m位置處的馬赫數(shù)及壓力云圖

圖10，圖11分別為堵蓋運動后，前5 m的受力與速度變化曲線，從圖中可以明顯看出，兩種處理方法對于堵蓋在噴管內(nèi)的運動影響較小；由于非接觸式預(yù)留縫隙方法作用下激波串結(jié)構(gòu)的產(chǎn)生，使得堵蓋在1.2 m以及2.2m位置處分別產(chǎn)生了兩次加速運動；接觸判斷處理方式作用下的堵蓋運動，受力更為平穩(wěn)，近似為一勻加速運動。

圖10 堵蓋的受力變化

圖11 堵蓋運動的速度變化

分析堵蓋啟動后的運動數(shù)據(jù)可知：接觸判斷方法的引入，使得堵塊的啟動狀態(tài)發(fā)生變化，數(shù)值上表現(xiàn)為流場結(jié)構(gòu)以及堵塊運動規(guī)律的差異。分析其成因，由于縫隙位置處理方式的不同，使得兩種方法得到的啟動時堵蓋后方物理場不同；同時由于堵蓋運動的延遲時間不同，使得堵蓋運動耦合物理場的時間節(jié)點發(fā)生偏移。

4 結(jié)論

1) 燃?xì)馄颇ば纬勺矒舳律w的激波，經(jīng)堵蓋的反射作用，快速的建立起燃?xì)鈪^(qū)啟動壓強(qiáng)；接觸判斷結(jié)合多孔介質(zhì)模型的處理方法，避免了縫隙位置的燃?xì)庑孤?，更?zhǔn)確的模擬出堵蓋的啟動壓強(qiáng)以及運動的延遲時間。

2) 堵蓋啟動后，伴隨堵蓋的運動，燃?xì)饨?jīng)堵蓋與噴管的縫隙流出并向外擴(kuò)散；相較于非接觸式預(yù)留縫隙的方法，建立在接觸判斷方法基礎(chǔ)上的堵蓋與燃?xì)獾鸟詈献饔酶鼮榭尚拧?/p>

3) 堵蓋在飛出過程中，受噴管內(nèi)高壓燃?xì)獾淖饔?，速度急劇增加，?dāng)移動到噴管出口位置外，受力趨于穩(wěn)定，速度變化趨穩(wěn)，在計算區(qū)間內(nèi)，未見明顯減速。

[1] 徐強(qiáng),李開明,張福祥.封閉式導(dǎo)彈發(fā)射箱內(nèi)燃?xì)饬魈匦詫嶒炑芯縖J].彈道學(xué)報,1995,7(2):52 -56.

[2] 徐強(qiáng),李軍.燃?xì)馍淞髌鹗紱_擊波形成機(jī)理的實驗研究[J].推進(jìn)技術(shù),2000,21(3):16-18.

[3] 劉琦,傅德彬,姜毅.貯運發(fā)射箱內(nèi)燃?xì)馍淞鞯姆嵌ǔ_擊波流場數(shù)值模擬[J].彈箭與制導(dǎo)學(xué)報,2005,25(S4):382 -384.

[4] 傅德彬,姜毅.某導(dǎo)彈易碎蓋的開啟過程[J].固體火箭技術(shù),2007,30(4):275 -277.

[5] 苗佩云, 袁曾鳳.同心發(fā)射筒燃?xì)忾_蓋技術(shù)[J].北京理工大學(xué)學(xué)報,2004,24(4):283 -285.

[6] ARNAB C, ABDELLAH H. Numerical simulation of transient supersonic nozzle flows [C]//50th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition. Nashville, TN:AIAA,2012.

[7] 于邵禎,姜毅,周笑飛,等.耦合尾噴管堵蓋運動的發(fā)射箱內(nèi)流場研究[J].兵工學(xué)報,2014,35(11):1806-1812.

[8] 張玉東,紀(jì)楚群.多體分離非定常氣動特性數(shù)值模擬[J]. 空氣動力學(xué)學(xué)報,2006,24(1):1-4.

[9] 劉君,白曉征,郭正,等.有相對運動多體動力學(xué)系統(tǒng)流動計算方法若干問題討論[J].空氣動力學(xué)學(xué)報, 2008, 26(1): 14-19.

[10] 王飛,朱立新,顧文彬. 基于ALE算法的空氣沖擊波繞流數(shù)值模擬研究[J].工程爆破,2002,8(2):13-16.

AnalysisofSeparationMovementofNozzleClosureBasedonContactDetectionMethod

ZHAO Mingyang， WANG Ge， HAN Wanzhi， ZHANG Qi

(College of Architecture and Aerospace Engineering, Harbin Engineering University, Harbin 150001, China)

In order to study the interaction between the closure and the gas in the initial separation from the nozzle and subsequent process, a simulation model is established by using a numerical method combined with the Contact Detection and Porous Media model, in which a dynamic mesh technique is used to update the meshes of flow field. The research results show that the closure obtains a high impact kinetic energy from the high-pressure gas,. And compared with the practice of reserved gap, Contact Detection and Porous Media model can simulate the delayed effect of the Separation.

nozzle closure; separation movement; contact detection; porous media

2017-09-28；

2017-10-19

趙明陽(1992—)，男，碩士研究生，主要從事火箭發(fā)動機(jī)燃燒、流動及多物理場計算研究。

王革(1966—)，男，博士，教授，主要從事火箭發(fā)動機(jī)燃燒、流動及內(nèi)彈道計算研究。

10.11809/scbgxb2017.12.046

本文引用格式：趙明陽, 王革，韓萬之，等.基于接觸判斷方法的尾噴管堵蓋分離運動分析[J].兵器裝備工程學(xué)報，2017(12):207-210,252.

formatZHAO Mingyang， WANG Ge， HAN Wanzhi, et al.Analysis of Separation Movement of Nozzle Closure Based on Contact Detection Method[J].Journal of Ordnance Equipment Engineering,2017(12):207-210,252.

V435；TJ768

A

2096-2304(2017)12-0207-04

(責(zé)任編輯楊繼森)