胡 坤， 張建華, 黃海峰

(海軍潛艇學(xué)院， 山東 青島 266199)

當(dāng)潛艇在水下航行過程中出現(xiàn)大縱傾及艙室破損等緊急情況時，如果只靠車、舵已無法控制險情，為保證潛艇安全，必須及時、準(zhǔn)確地向主壓載水艙釋放高壓氣，排出水艙內(nèi)的壓載水，以迅速提供足夠的正浮力，挽回潛艇的深度和縱傾。[1]

在潛艇應(yīng)急吹除過程中，高壓氣閥柱是高壓氣從氣瓶流向水艙全過程中途經(jīng)的第一個部件，對高壓氣起著匯聚和再分配的重要作用。因此，在對吹除過程進(jìn)行理論研究時應(yīng)先建立相應(yīng)的數(shù)學(xué)模型。然而，由于閥柱內(nèi)的流場形態(tài)過于復(fù)雜，為便于研究，目前，國內(nèi)的研究學(xué)者[2-4]在對高壓氣應(yīng)急吹除主壓載水艙過程進(jìn)行數(shù)學(xué)建模時，均將其忽略，這樣的簡化處理方法是否合理尚有待考證。另外，對于閥柱內(nèi)部的流場形態(tài)如何、其存在對吹除速率產(chǎn)生什么影響以及在不同吹除組合方式下其內(nèi)部流場形態(tài)會有何不同等問題的研究目前也處于空白。因此，采用計算流體動力學(xué)(Compwtational Fluid Dynamics, CFD)數(shù)值模擬的方法，對簡化的閥柱模型進(jìn)行仿真分析，旨在摸清閥柱內(nèi)部的流場形態(tài)及其變化規(guī)律，為建立完整的高壓氣吹除主壓載水艙理論數(shù)學(xué)模型提供依據(jù)。

1 高壓氣吹除主壓載水艙基本原理

潛艇高壓氣吹除主壓載水艙基本原理示意見圖1。高壓氣從氣瓶釋放之后，首先匯集至高壓氣站，然后再由高壓氣站的高壓氣閥柱對高氣壓進(jìn)行再分配，最終進(jìn)入需要吹除的主壓載水艙中并將水艙中的海水通過通海閥吹除至舷外。高壓氣應(yīng)急吹除主壓載水艙時，可選擇采用多組氣瓶同時向1個水艙供氣，也可采用1組氣瓶同時向多個水艙供氣，或是多組氣瓶同時向多個水艙供氣，這一重要功能實現(xiàn)的關(guān)鍵部件就是高壓氣閥柱。

圖1 高壓氣吹除主壓載水艙示意

2 控制方程與湍流模型

由于高壓氣瓶內(nèi)的氣體壓力遠(yuǎn)大于主壓載水艙處的背壓，氣體將以射流的形式自進(jìn)氣管路進(jìn)入閥柱并迅速擴(kuò)散，而可實現(xiàn)k-ε湍流模型能較精確地預(yù)測平面和圓形射流擴(kuò)散作用，且對于旋轉(zhuǎn)流動、流動分離和二次流有很好的表現(xiàn)[5]。因此,數(shù)值仿真時采用Realizablek-ε湍流模型來封閉時均化的Navier-Stokes方程。

2.1 連續(xù)性方程

(1)

2.2 動量方程

(2)

2.3 能量方程

(3)

式(3)中：cp為比熱容；k為流體的導(dǎo)熱系數(shù)；ST為黏性耗散項；T為溫度。

2.4 湍流模型

Realizablek-ε模型的湍動能k與耗散率ε的輸運方程[6-7]為

(4)

(5)

式(4)～式(5)中：σk=1.0；σε=1.2；C2=1.9；Cμ=0.09；μt、Gk、Gb、YM等參數(shù)的意義及其表達(dá)式與標(biāo)準(zhǔn)模型相同，具體見文獻(xiàn)[6]。

3 數(shù)值仿真方法

3.1 試驗方案與模型建立

為研究閥柱內(nèi)部壓力與氣瓶組壓力之間的關(guān)系和選擇不同的水艙進(jìn)行吹除時閥柱內(nèi)部流場形態(tài)間的差異，同時考慮到網(wǎng)格劃分與運算時間等方面的數(shù)值計算成本，建立包含有2個入口和2個出口的簡化系統(tǒng)物理模型,見圖2。其中：氣體入口管路直徑、出口管路直徑和閥柱的內(nèi)徑分別為21 mm、34 mm和52 mm；入口管路長度和出口管路長度均有250 mm和500 mm兩種規(guī)格。

圖2 高壓氣閥柱簡化幾何模型

為便于比較，制訂4組試驗方案，具體見表1。

表1 閥柱內(nèi)流場數(shù)值模擬試驗方案

通過對比方案1與方案2的數(shù)值模擬結(jié)果可研究氣瓶組的壓力大小對閥柱內(nèi)流場分布的影響；通過對比方案2與方案3的數(shù)值模擬結(jié)果可研究氣瓶組的數(shù)量對閥柱內(nèi)流場分布的影響；通過對比方案3與方案4的數(shù)值模擬結(jié)果可研究氣瓶組壓力相差懸殊但均大于出口背壓時的氣體流動規(guī)律。

3.2 網(wǎng)格劃分

為獲得較高的計算效率和計算精度，采用六面體結(jié)構(gòu)網(wǎng)格單元對計算域進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分；為能更好地捕捉流場信息，對管路與閥柱相貫處流動較劇烈的區(qū)域進(jìn)行網(wǎng)格加密；為有效模擬近壁面處的流動，邊界層區(qū)域內(nèi)合理布置網(wǎng)格并選擇合適的網(wǎng)格尺度；為減少計算工作量，加快計算速度，僅對計算域的1/2進(jìn)行網(wǎng)格劃分，見圖3。

圖3 閥柱內(nèi)流場結(jié)構(gòu)化網(wǎng)格

3.3 邊界條件與數(shù)值方法

1) 入口1、入口2、出口1和出口2：根據(jù)不同的試驗方案，具體邊界類型設(shè)置如表1所示，其中:入口1與入口2壓力隨時間的變化采用用戶定義函數(shù)(User Defined Function，UDF)技術(shù)來實現(xiàn)。

2) 對稱面：設(shè)置為對稱邊界條件，該面上所有流動變量的通量為零。

3) 其他面：除入口、出口和對稱面之外的其它面均設(shè)為無滑移固壁邊界條件。

采用有限體積法對控制方程進(jìn)行離散，采用SIMPLE方法計算壓力-速度關(guān)聯(lián)方程，可實現(xiàn)k-ε兩方程湍流模型封閉控制方程組，壓力項采用標(biāo)準(zhǔn)格式進(jìn)行離散，其他項均采用二階迎風(fēng)格式進(jìn)行離散。

4 計算結(jié)果與分析

高壓氣進(jìn)氣管路進(jìn)入閥柱，被閥柱內(nèi)壁阻擋后再發(fā)生旋轉(zhuǎn)，并與其他進(jìn)氣管路的氣體發(fā)生混合，這導(dǎo)致閥柱內(nèi)的氣體流場形態(tài)要比在管道內(nèi)復(fù)雜得多，因為在一定的主流速度和幾何邊界條件下，閥柱內(nèi)的高壓氣會產(chǎn)生二次流，這也是除幾何尺寸、離心力、流動方向等因素之外，高壓氣閥柱產(chǎn)生壓降的一個重要原因。在高壓氣應(yīng)急吹除主壓載水艙過程中，影響主壓載水艙排水速率的2個關(guān)鍵因素為噴入水艙的氣體質(zhì)量流量和通海閥處的排水背壓，而根據(jù)氣體動力學(xué)相關(guān)理論[8]，氣體質(zhì)量流量的大小由閥柱內(nèi)氣體壓力與排水背壓共同決定。因此，該仿真結(jié)果重點關(guān)注閥柱內(nèi)的壓力變化規(guī)律和各出流支管的氣體質(zhì)量流量變化規(guī)律，見圖4～圖7。

a) 入口總壓與閥柱內(nèi)總壓

b) 入口與出口質(zhì)量流量

a) 入口總壓與閥柱內(nèi)總壓

b) 入口與出口質(zhì)量流量

a) 入口總壓與閥柱內(nèi)總壓

b) 入口與出口質(zhì)量流量

a) 入口總壓與閥柱內(nèi)總壓

b) 入口與出口質(zhì)量流量

由圖4～圖7可知：閥柱內(nèi)總壓變化規(guī)律與入口總壓變化規(guī)律一致，均呈線性變化，但前者的斜率要小于后者。所有入口的氣體質(zhì)量流量之和等于所有出口的氣體質(zhì)量流量之和，滿足質(zhì)量守恒定律，各試驗方案所進(jìn)行的數(shù)值計算具有良好的收斂性(說明：在數(shù)值計算時，F(xiàn)LUENT將流入系統(tǒng)的質(zhì)量流量記為正值，將從系統(tǒng)流出的質(zhì)量流量記為負(fù)值，為便于分析，此處將方案1～方案3中的出口質(zhì)量流量也記為正值)。

由圖4與圖5可知：閥柱內(nèi)氣體總壓和管路中總的氣體質(zhì)量流量均與入口總壓密切相關(guān)。當(dāng)入口總壓升高時，閥柱內(nèi)總壓亦隨之升高，同時各管路中氣體質(zhì)量流量也相應(yīng)變大。但對比圖4b和5b中入口1的氣體質(zhì)量流量可發(fā)現(xiàn)：入口2的壓力升高后，入口1的氣體質(zhì)量流量有所減小，這是閥柱內(nèi)的壓力升高導(dǎo)致入口管路1兩端的壓差變小所致；另外，無論入口壓力如何變化，出口管路2始終比出口管路1的氣體質(zhì)量流量大，這是因為出口管路2的長度較管路1短，這表明當(dāng)潛艇在進(jìn)行高壓氣應(yīng)急吹除時，距離閥柱越近的水艙進(jìn)氣量越多，其排水速率相應(yīng)也會越大。因此，在實艇進(jìn)行全部主壓載水艙應(yīng)急吹除時，由于艏部和舯前部的主壓載水艙距離高壓氣閥柱較近，管路較短，從而排水速率會略大于艉部主壓載水艙。

由圖5a與圖6a比較可知：出口數(shù)量對閥柱內(nèi)的壓力具有重要影響。在方案2中，出口管路數(shù)量為2條，而在方案3中出口管路數(shù)量僅為1條，在其他條件均相同的情況下，后者閥柱內(nèi)總壓約為前者的1.5倍。這是因為當(dāng)出口管路為2條時，在穩(wěn)定狀態(tài)下，入口與出口處的質(zhì)量流量應(yīng)按照圖5b所示的規(guī)律進(jìn)行變化，且流入閥柱的氣體質(zhì)量流量等于流出閥柱的氣體質(zhì)量流量。當(dāng)關(guān)閉1條出口管路時(即方案3)，則相當(dāng)于本該從該管路排出的氣體卻無法排出，從而在閥柱中被“堆積”起來，導(dǎo)致閥柱內(nèi)壓力升高。閥柱內(nèi)壓力可看做是入口管路的背壓，同時也可看做是出口管路的初始總壓。根據(jù)氣體動力學(xué)理論，其值變大會導(dǎo)致入口氣體質(zhì)量流量減小(當(dāng)流動為亞聲速流時)，而出口氣體質(zhì)量流量增大如圖6b所示；根據(jù)質(zhì)量守恒定律，進(jìn)入與流出閥柱的氣體質(zhì)量流量相等，系統(tǒng)最終會達(dá)到一個新的平衡狀態(tài)，即按照圖6所示的規(guī)律進(jìn)行變化。由此可推斷出以下結(jié)論：若參與供氣的氣瓶組和壓力固定，則被吹除的水艙越多，氣瓶組壓降就越快；同理，若被吹除的水艙固定，則參與供氣的氣瓶組越多，閥柱內(nèi)壓力就越高，且越接近于氣瓶內(nèi)壓力，進(jìn)入水艙的氣體質(zhì)量流量也越大，排水速率也越快。因此，當(dāng)潛艇出現(xiàn)損失浮力等險情時，應(yīng)盡可能讓更多的氣瓶參與供氣，不建議同時吹除所有水艙，應(yīng)首選艏部和舯部水艙，這將有利于潛艇形成有利于抗沉的艉傾姿態(tài)。

由圖6與圖7比較可知：當(dāng)閥柱有多個入口且各入口壓力值不同時，則閥柱內(nèi)的壓力可能低于所有入口的壓力。但如果各入口壓力值相差懸殊，則閥柱內(nèi)壓力也可能出現(xiàn)高于其中1個或多個入口壓力的情況，壓力值低于閥柱內(nèi)壓力的入口則變成“出口”，氣體將通過該截面從系統(tǒng)中流出(如圖7b所示)。然而，在實艇高壓氣吹除主壓載水艙過程中不會出現(xiàn)此類現(xiàn)象，因為控制氣瓶組參與供氣的閥為單向閥，氣體只能從氣瓶中向閥柱單向流通。由此可知：當(dāng)高壓氣應(yīng)急吹除主壓載水艙時，所有與閥柱連通的氣瓶組中，只有壓力高于閥柱內(nèi)壓力的氣瓶組才實際參與供氣排水。

5 結(jié)束語

通過CFD數(shù)值仿真結(jié)果可知：潛艇高壓氣系統(tǒng)在正常工作時，閥柱對于氣瓶內(nèi)的壓力變化和主壓載水艙的排水速率均具有重要影響，在進(jìn)行有關(guān)理論研究時，不可簡單將其忽略。

1) 由于舯艉部的主壓載水艙距離高壓氣閥柱較遠(yuǎn)，因此，當(dāng)進(jìn)行全部主壓載水艙吹除時，艏部的主壓載水艙排水速率較快，將在一定程度上有助于潛艇產(chǎn)生艉傾。

2) 在相同吹除條件下，被吹除的水艙越多，氣瓶組壓降將越快。

3) 在相同吹除條件下，參與供氣的氣瓶組越多，進(jìn)入水艙的氣體質(zhì)量流量越大，排水速率也越快。

4) 在實際高壓氣應(yīng)急使用時，應(yīng)盡可能讓更多的氣瓶參與供氣，且相比較于同時吹除全部主壓載水艙，僅吹除艏部和舯部主壓載水艙的效果可能會更好，因其能使得潛艇形成有利于抗沉的艉傾姿態(tài)。