胡 坤， 張建華, 黃海峰

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

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

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

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

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

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

2 控制方程與湍流模型

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

2.1 連續(xù)性方程

(1)

2.2 動(dòng)量方程

(2)

2.3 能量方程

(3)

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

2.4 湍流模型

Realizablek-ε模型的湍動(dòng)能k與耗散率ε的輸運(yùn)方程[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 試驗(yàn)方案與模型建立

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

圖2 高壓氣閥柱簡(jiǎn)化幾何模型

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

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

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

3.2 網(wǎng)格劃分

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

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

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

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

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

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

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

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

高壓氣進(jìn)氣管路進(jìn)入閥柱，被閥柱內(nèi)壁阻擋后再發(fā)生旋轉(zhuǎn)，并與其他進(jìn)氣管路的氣體發(fā)生混合，這導(dǎo)致閥柱內(nèi)的氣體流場(chǎng)形態(tài)要比在管道內(nèi)復(fù)雜得多，因?yàn)樵谝欢ǖ闹髁魉俣群蛶缀芜吔鐥l件下，閥柱內(nèi)的高壓氣會(huì)產(chǎn)生二次流，這也是除幾何尺寸、離心力、流動(dòng)方向等因素之外，高壓氣閥柱產(chǎn)生壓降的一個(gè)重要原因。在高壓氣應(yīng)急吹除主壓載水艙過程中，影響主壓載水艙排水速率的2個(gè)關(guān)鍵因素為噴入水艙的氣體質(zhì)量流量和通海閥處的排水背壓，而根據(jù)氣體動(dòng)力學(xué)相關(guān)理論[8]，氣體質(zhì)量流量的大小由閥柱內(nèi)氣體壓力與排水背壓共同決定。因此，該仿真結(jié)果重點(diǎn)關(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ì)量守恒定律，各試驗(yàn)方案所進(jìn)行的數(shù)值計(jì)算具有良好的收斂性(說明：在數(shù)值計(jì)算時(shí)，F(xiàn)LUENT將流入系統(tǒng)的質(zhì)量流量記為正值，將從系統(tǒng)流出的質(zhì)量流量記為負(fù)值，為便于分析，此處將方案1～方案3中的出口質(zhì)量流量也記為正值)。

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

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

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

5 結(jié)束語

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

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

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

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

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