薛寒寒，熊言義，姬 翔，李玉峰

(中國(guó)船舶集團(tuán)公司第七一三研究所，河南 鄭州 450015)

0 引 言

船舶作為水上可移動(dòng)的載物平臺(tái)，艙室內(nèi)設(shè)備越來越多，設(shè)備之間更加緊湊，艙內(nèi)的電纜發(fā)熱、設(shè)備正常運(yùn)行期間會(huì)向周圍空氣釋放大量的熱量，高溫環(huán)境下容易引發(fā)火災(zāi)，統(tǒng)計(jì)數(shù)據(jù)顯示火災(zāi)已成為船舶安全的主要威脅之一[1-2]。因此國(guó)內(nèi)外學(xué)者對(duì)船舶艙室火災(zāi)進(jìn)行了不少仿真和試驗(yàn)研究，Gouutk[3]實(shí)驗(yàn)研究了全尺寸艙室回燃現(xiàn)象，探討了回燃的主控參數(shù)；Rocket[4]對(duì)美國(guó)海軍艦船火災(zāi)增長(zhǎng)模型和煙氣蔓延模型進(jìn)行了研究和評(píng)價(jià)。鄒高萬等[5]對(duì)船舶機(jī)艙火災(zāi)熱流場(chǎng)特性進(jìn)行了數(shù)值仿真研究，用耗氧原理測(cè)量并得到不同通風(fēng)因子下轟燃臨界釋熱速率；董華等[6]對(duì)密閉艙室的火災(zāi)發(fā)展過程進(jìn)行了數(shù)值研究，胡隆華[7]對(duì)隧道火災(zāi)煙氣蔓延的熱物理特性研究。狹長(zhǎng)艙室火災(zāi)會(huì)有一些獨(dú)有特征，若能較好把握狹長(zhǎng)空間的火災(zāi)特性，有利于火災(zāi)初期的準(zhǔn)確自動(dòng)探測(cè)報(bào)警。本文使用流體仿真軟件通過計(jì)算機(jī)建模以仿真形式再現(xiàn)艦船狹長(zhǎng)艙室火災(zāi)場(chǎng)景，得到狹長(zhǎng)空間火災(zāi)發(fā)生后煙氣蔓延特性規(guī)律，艙內(nèi)溫度場(chǎng)分布以及隨時(shí)間的變化，不同位置的火災(zāi)特征參數(shù)變化的差異性，以及典型工況下特征參數(shù)的閾值，探討設(shè)備發(fā)熱對(duì)火災(zāi)熱流場(chǎng)的影響。通過分析火災(zāi)流場(chǎng)特性變化規(guī)律，為狹長(zhǎng)空間火災(zāi)探測(cè)設(shè)計(jì)提供參考，對(duì)于其他場(chǎng)所狹長(zhǎng)空間火災(zāi)規(guī)律研究也有重要意義。

1 數(shù)學(xué)模型

1.1 控制方程組

狹長(zhǎng)艙室內(nèi)的火災(zāi)熱流場(chǎng)主要為熱浮力驅(qū)動(dòng)的低馬赫數(shù)流動(dòng)過程，對(duì)于此類受重力作用的多組分理想氣體的通用控制方程組，根據(jù)火災(zāi)流場(chǎng)特性和計(jì)算求解速度的需要，對(duì)方程組進(jìn)行一些處理，具體過程在文獻(xiàn)[8-9]中已有詳細(xì)推導(dǎo)與論述。通過對(duì)控制方程組中速度散度項(xiàng) ▽·V及動(dòng)量方程的變形與簡(jiǎn)化，最終可得如下所需要求解的方程組：

1.2 濾波控制

為了得到大渦模擬技術(shù)處理的控制方程組，需要對(duì)與時(shí)間相關(guān)的N-S 方程組在傅里葉空間或物理空間進(jìn)行濾波。對(duì)任意一個(gè)流動(dòng)變量 ?(x)，可將其劃分為大尺度量(x) 和小尺度量?′(x)：

目前常用的濾波函數(shù)有高斯濾波、盒式濾波（又稱頂蓋濾波）和傅里葉階段濾波。其中比較常用的盒式濾波可以表示為：

1.3 燃燒模型

計(jì)算中將燃燒過程簡(jiǎn)化處理，認(rèn)為燃料和氧氣混合燃燒速率無限快，燃燒過程僅和二者混合分?jǐn)?shù)有關(guān)，可以簡(jiǎn)化計(jì)算量，適用于較大空間內(nèi)火災(zāi)的模擬，其對(duì)于單步瞬時(shí)反應(yīng)假設(shè)燃料與氧氣混合后按如下化學(xué)反應(yīng)進(jìn)行：

1.4 熱輻射模型

熱輻射在火災(zāi)流場(chǎng)計(jì)算中對(duì)計(jì)算結(jié)果有較大影響，考慮到熱輻射的復(fù)雜性，要完全模擬難度較大，一般對(duì)模型合理簡(jiǎn)化處理。輻射控制方程為：

2 仿真模型

2.1 幾何模型

艦船艙室狹長(zhǎng)空間設(shè)備較多，結(jié)構(gòu)比較復(fù)雜，圖1為仿真計(jì)算時(shí)簡(jiǎn)化的狹長(zhǎng)艙室空間。簡(jiǎn)化后縱切面為半梯形結(jié)構(gòu)，右側(cè)壁為發(fā)熱設(shè)備，左側(cè)為V 型空間，上甲板整體尺寸（長(zhǎng)×寬）為：36 m×1.6 m，下底板整體尺寸（長(zhǎng)×寬）為：36 m×0.4 m，高為1.5 m，艙內(nèi)其他設(shè)備進(jìn)行同樣簡(jiǎn)化，簡(jiǎn)化后計(jì)算網(wǎng)格如圖2 所示。為準(zhǔn)確捕捉火場(chǎng)信息，火源及其附近區(qū)域網(wǎng)格尺寸為0.02 m，其他區(qū)域采用邊長(zhǎng)不大于0.03 m 的網(wǎng)格進(jìn)行劃分，整個(gè)計(jì)算空間的總體網(wǎng)格數(shù)約為350 萬，艙室表面固體均假設(shè)為絕熱邊界，沒有熱量向外散失。

圖1 艙室狹長(zhǎng)空間三維模型Fig. 1 Three-dimensional model of ship's long-narrow space

計(jì)算過程中，基于Deardorff 亞格子模型的大渦模擬（LES）技術(shù)求解湍流模型，用有限體積法（FVM）計(jì)算熱輻射，以顯性二階預(yù)測(cè)/校正機(jī)制迭代變量的時(shí)間微分離散，用CFL 條件保證迭代的穩(wěn)定。

2.2 火源位置和大小

根據(jù)狹長(zhǎng)艙室內(nèi)設(shè)備實(shí)際情況，結(jié)合相關(guān)文獻(xiàn)中記錄的電纜火災(zāi)發(fā)生的原因[8-9]，同時(shí)考慮不同火源位置對(duì)探測(cè)報(bào)警裝置的影響，分別模擬火源位于狹長(zhǎng)艙室邊緣上部、邊緣下部2 個(gè)典型位置發(fā)生火災(zāi)時(shí)的熱流場(chǎng)，火源設(shè)置在距離端面4 m 處，具體位置如圖1所示。

實(shí)際的火源情況非常復(fù)雜，在盡可能真實(shí)地反映燃燒過程基礎(chǔ)上，將火源抽象并簡(jiǎn)化，結(jié)合國(guó)際標(biāo)準(zhǔn)ISO/TS16733，最普遍使用的 αt2火災(zāi)模型。不同火災(zāi)增長(zhǎng)系數(shù) α定義了4 種標(biāo)準(zhǔn)t2火災(zāi)，t2火災(zāi)特性發(fā)展曲線圖中給出了慢速火、中速火、快速火和超快速火的火源功率達(dá)到1MW 的時(shí)間，它們分別是75 s，150 s，300 s，600 s。參考SFPE 消防手冊(cè)Electric cable trays部分Lee 數(shù)據(jù)，電纜（PE/PVC cable）的HRR 值為1071 kW/m2。塑料泡沫火災(zāi)屬于快速火，增長(zhǎng)系數(shù)為0.0469 kW/s2，著火面積為2 m×0.09 m=0.18 m2，經(jīng)計(jì)算需要55 s 達(dá)到最大1071×2×0.09=192.8 kW，所以程序設(shè)置火源增長(zhǎng)時(shí)間TAU_Q=?64，HRR=1071 kW/m2。在電纜熱釋放增長(zhǎng)速率如圖2 所示。

圖2 熱釋放速率曲線Fig. 2 Heat release rate curve

2.3 測(cè)點(diǎn)布置

為分析各種工況下溫度變化規(guī)律，在計(jì)算過程中對(duì)特征參數(shù)進(jìn)行記錄，在火源中心截面平面上布置溫度測(cè)點(diǎn)，具體測(cè)點(diǎn)位置如圖3 所示。采樣周期為0.2 s，頻率為5 Hz。

圖3 截面和監(jiān)測(cè)點(diǎn)位置示意圖Fig. 3 Diagram of sections and monitor points

2.4 計(jì)算方法

本文采用CFD 數(shù)值模擬方法，根據(jù)狹長(zhǎng)艙室火災(zāi)流場(chǎng)特性，計(jì)算過程描述如下：1）設(shè)置艙室初始環(huán)境條件包括時(shí)間步和各變量值和；2）計(jì)算下一時(shí)刻的參數(shù)預(yù)測(cè)值；3）根據(jù)速度預(yù)測(cè)值求與第2 步得到的速度值判定收斂與否；4）結(jié)合CFL 判定是否穩(wěn)定；5）求解參數(shù)校正值；6）與在第5 步中求得的速度散度值比較，判定收斂與否；7）返回第2 步繼續(xù)計(jì)算；8）反復(fù)迭代完成計(jì)算。

2.5 計(jì)算工況

假設(shè)狹長(zhǎng)艙室內(nèi)電纜因?yàn)楦邷匦纬苫馂?zāi)，在電纜保持發(fā)熱的前提下，分別考慮甲板69 ℃高溫和發(fā)熱設(shè)備以80 ℃溫度發(fā)熱獨(dú)立和耦合作用時(shí)，以電纜（固體）蔓延火為火源，模擬火源位于狹長(zhǎng)艙室艙邊緣上部典型位置發(fā)生火災(zāi)時(shí)的熱流場(chǎng)。計(jì)算時(shí)間為1000 s，其中800 s 為艙內(nèi)環(huán)境達(dá)到穩(wěn)態(tài)的時(shí)間，200 s 為火災(zāi)發(fā)生初期所形成的熱流場(chǎng)。通過模擬結(jié)果判斷在火源工況下狹長(zhǎng)艙室內(nèi)各測(cè)點(diǎn)處的溫度隨火災(zāi)發(fā)展的變化情況。模擬設(shè)備發(fā)熱和不發(fā)熱2 種工況，如表1 所示。

表1 計(jì)算工況Tab. 1 Conditions of calculation

3 計(jì)算結(jié)果和分析

圖4 為火源位于艙室邊緣上層時(shí)不同時(shí)刻煙氣蔓延分布。在t=30 s 時(shí)位于艙室邊緣處火源產(chǎn)生的煙氣在火源處向上運(yùn)動(dòng)，撞擊到艙室上壁面后開始向周圍橫向蔓延，但由于艙室一端是封閉狀態(tài)，因此左端壁面對(duì)火災(zāi)煙氣的擴(kuò)散蔓延起到了一定的影響，表現(xiàn)在火災(zāi)發(fā)生90 s 時(shí)煙氣碰到左端壁面時(shí)煙氣向艙室內(nèi)下部蔓延，在火災(zāi)發(fā)生180 s 時(shí)艙室左側(cè)下部邊緣處空間已基本都被火災(zāi)煙氣覆蓋，而艙室另一端基本未受煙氣侵?jǐn)_（煙氣長(zhǎng)度約25 m，有11 m 左右沒有被煙氣蔓延），這是由于狹長(zhǎng)空間結(jié)構(gòu)導(dǎo)致煙氣在水平方向阻力大于向下的浮力阻力，煙氣不足以繼續(xù)向前運(yùn)動(dòng)而改為向下運(yùn)動(dòng)。而普通空間的火災(zāi)一般煙氣在整個(gè)空間的頂部完全彌散后才向空間下部蔓延。

圖4 火源位于艙室邊緣上層時(shí)不同時(shí)刻煙氣蔓延Fig. 4 Smoke spread over time when fire is on the edge of upper layer

圖5 為z=0.34 m 時(shí)平均剖面溫度在35 ℃～55 ℃范圍內(nèi)的變化。t=30 s 時(shí)，溫度差異并不明顯，整體溫度在47 ℃左右，t=60 s 時(shí)，火源下方處的溫度明顯升高，并逐步向兩側(cè)延伸，由于煙氣在艙室邊緣處聚集，在火災(zāi)發(fā)生180 s 時(shí)，高溫區(qū)域擴(kuò)大，而且高溫?zé)煔夥植疾痪鶆驅(qū)е聹囟确植汲霈F(xiàn)較大差異。z=0.78 m處平均剖面溫度在35 ℃～65 ℃，高于艙室下部截面z=0.34 m，在t=30 s 時(shí)剖面溫度在53 ℃，比z=0.34 m 溫度高6 K，在火災(zāi)發(fā)生180 s 時(shí)，在CD 艙長(zhǎng)度方向上溫度分布不再是連續(xù)升高，而是局部區(qū)域出現(xiàn)溫度降低又升高現(xiàn)象，這是由于在該位置處煙氣發(fā)生沉降而發(fā)生局部溫度較高。z=1.11m 處平均剖面溫度在35 ℃～111 ℃范圍內(nèi)變化，t=30 s 時(shí)，火源上方靠側(cè)壁局部溫度遠(yuǎn)高于其他區(qū)域，溫度分布很不均勻，t=180 s 時(shí)，以火源位置為對(duì)稱軸到艙室左右兩側(cè)溫度值均在111 ℃，沿艙室延伸方向溫度逐漸降低。

圖5 上邊電纜著火時(shí)不同時(shí)刻不同高度截面的溫度場(chǎng)分布Fig. 5 Temperature field distribution at different heights and time

圖6為火災(zāi)發(fā)生后10 個(gè)溫度探測(cè)器T1～T10溫度隨時(shí)間變化結(jié)果?？炕鹪锤浇黅10測(cè)點(diǎn)最先開始升溫，其余測(cè)點(diǎn)變化趨勢(shì)相同。在火源正上方T7～T10，隨著高度的增加溫度變小，而在遠(yuǎn)離火源的T1～T6隨著高度的增大，溫度變大。其中遠(yuǎn)離火源且處于下層測(cè)點(diǎn)的T4～T6的溫度變化幅度小，火源上方的測(cè)點(diǎn)T7～T10溫度整體高于遠(yuǎn)離火源的測(cè)點(diǎn)T1～T6溫度。在40 s 后上層測(cè)點(diǎn)T1和T7溫度曲線基本一致，隨著煙氣的蔓延在160 s 后上層測(cè)點(diǎn)T1溫度又低于火源上方測(cè)點(diǎn)T7的溫度。整個(gè)火災(zāi)過程在t=55 s 時(shí)達(dá)到火源上方距離最近的T10溫度最高達(dá)到747 ℃，而此時(shí)靠近艙室下層的T6只有62 ℃。

圖6 各測(cè)點(diǎn)測(cè)點(diǎn)溫度隨時(shí)間的變化Fig. 6 The variation in monitor points' temperature over time

圖7為溫度梯度。在火災(zāi)充分發(fā)展之前，也即前50 sT10～T7溫度梯度逐漸變小，T6～T1溫度梯度逐漸增大。50 s 之后溫度梯度變化最大的是T2測(cè)點(diǎn)，最大梯度達(dá)到了40 ℃/s，這是由于在T2位置煙氣渦旋運(yùn)動(dòng)更劇烈導(dǎo)致。

圖7 各測(cè)點(diǎn)溫度梯度隨時(shí)間的變化Fig. 7 Monitor points′ temperature gradient over time

實(shí)際工程中，火災(zāi)發(fā)生有一定的隨機(jī)性，探測(cè)器不在火源正上方是大概率事件，探測(cè)要盡早盡準(zhǔn)，因此重點(diǎn)分析非火源正上方的T1～T6測(cè)點(diǎn)在火災(zāi)發(fā)生初期的溫度和溫度梯度隨時(shí)間變化。一般來說探測(cè)器達(dá)到某一設(shè)定的閾值溫度或者溫度梯度時(shí)開始報(bào)警，或者同時(shí)達(dá)到溫度值和溫度梯度時(shí)開始報(bào)警。本次仿真溫度閾值取80 ℃，溫度梯度取3 ℃/s 時(shí)。T1～T6測(cè)點(diǎn)安溫度和溫度梯度達(dá)到閾值的時(shí)間如圖8 所示。

圖8 溫度和溫度梯度達(dá)到閾值的時(shí)間Fig. 8 Time for temperature and temperature gradient to reach threshold

圖8 顯示T1～T5五個(gè)測(cè)點(diǎn)的溫度閾值達(dá)到80 ℃時(shí)間都大于溫度梯度達(dá)到3 ℃/s 所需時(shí)間。而T6測(cè)點(diǎn)在240 s 達(dá)到溫度閾值，整個(gè)火災(zāi)過程溫度梯度都沒有達(dá)到閾值。T1～T2無論是溫度還是溫度梯度達(dá)到閾值的時(shí)間均在15 s 以下，溫度梯度綜合考慮了溫度和時(shí)間2 個(gè)因素，可以實(shí)現(xiàn)及早報(bào)警。當(dāng)然無論是溫度還是溫度梯度只要閾值設(shè)置的足夠低，就能實(shí)現(xiàn)及早報(bào)警。當(dāng)閾值足夠低時(shí)，可能不是火災(zāi)卻達(dá)到了報(bào)警閾值，從而錯(cuò)誤報(bào)警。因此無論是溫度還是溫度梯度結(jié)合火災(zāi)特性都要設(shè)置一個(gè)合理的閾值。探測(cè)器布置在頂層附近，對(duì)艙室其他設(shè)備的干涉影響較小，此外也能實(shí)現(xiàn)較早的報(bào)警。

從圖9 可以看出，在火災(zāi)發(fā)生初期，有設(shè)備發(fā)熱時(shí)，T1，T6，T7和T10升溫速度更快。而在火災(zāi)充分發(fā)展階段，設(shè)備發(fā)熱對(duì)T1，T7和T10三個(gè)測(cè)點(diǎn)的溫度的升高作用不明顯，由于新的熱源加入，高溫?zé)煔獾倪\(yùn)動(dòng)更紊亂，3 個(gè)測(cè)點(diǎn)溫度的波動(dòng)范圍更大。而下層測(cè)點(diǎn)T6溫度變化顯示，有設(shè)備發(fā)熱時(shí)溫升速率更快，溫升更大。

圖9 有/無發(fā)熱設(shè)備對(duì)火災(zāi)各測(cè)點(diǎn)溫度的影響Fig. 9 Influence of heating equipment on measuring points′ temperature

4 結(jié) 語

通過以上計(jì)算和分析得出以下結(jié)論：

1）對(duì)于狹長(zhǎng)空間，煙氣在火源處向上運(yùn)動(dòng)，漂浮至上壁面后開始向周圍橫向蔓延，不同于普通空間煙氣在頂部完全彌散，在180 s 艙室一端被煙氣完全覆蓋，另一端未受煙氣侵?jǐn)_，這是由于狹長(zhǎng)空間結(jié)構(gòu)導(dǎo)致煙氣在水平方向阻力持續(xù)增大，直至大于向下的阻力導(dǎo)致。

2）當(dāng)溫度閾值80 ℃，溫度梯度閾值為3 ℃/s 的時(shí)間，上層測(cè)點(diǎn)T1和T2，無論溫度還是溫度梯度達(dá)到閾值的時(shí)間均在15s 之內(nèi)。

3）當(dāng)設(shè)備發(fā)熱時(shí)，上層測(cè)點(diǎn)溫升不明顯，溫度波動(dòng)范圍更大；下層測(cè)點(diǎn)溫升速率更快，溫升更大。