隔離艙室污染氣體擴散的大渦模擬研究

劉琦，謝承利，李偉光，劉文川，龍新平

1 中國艦船研究設計中心，湖北 武漢 430064

2 煤礦災害動力學與控制國家重點實驗室，重慶 400044

3 重慶大學 資源與安全學院，重慶 400044

4 武漢大學 水射流理論與新技術湖北省重點實驗室，湖北 武漢 430072

0 引 言

船舶艙室是一個相對密閉的室內環境，其空氣品質受污染源、顆粒物、溫度、濕度等諸多因素影響。日本“鉆石公主”號（Diamond Princess）郵輪新冠病毒感染事件使船舶系統設計人員更加清晰地意識到，必須設計可靠的船舶空氣調節和通風系統、建立合理的氣流組織，以有效保障艙室的空氣品質，從而降低病毒通過空氣傳播導致交叉感染的風險。

為確保船舶傳染病隔離病房內含有病毒的污染空氣不侵入周圍的潔凈空間，通常需要設計負壓系統以維持病房內部較低的負壓環境。美國疾病控制中心認為，控制氣流進出房間的最小壓差僅為0.25 Pa，但考慮門開閉和熱浮升力效應等影響，推薦隔離房間負壓大于2.5 Pa[1]。我國《傳染病醫院建筑設計規范》[2]規定，負壓隔離病房與其相鄰及相通的緩沖間、走廊壓差應保持不小于5 Pa的負壓差。此類準則和規范為相應的負壓系統設計提供了重要參考，但是無法準確評估病房內壓力波動等實際因素對其內部氣流組織和污染物分布的影響，同時缺少針對船舶隔離艙室負壓系統設計的相關研究。近年來，計算流體動力學（computational fluid dynamics，CFD）技術已在室內污染物擴散分布研究領域得到廣泛應用[3-7]。研究人員基于雷諾平均Navier-Stokes方程（RANS）方法并結合動網格技術，對人員走動和門開閉過程的流場特征進行動態模擬，結果表明，室內物體移動誘導的氣流運動會對室內流場和污染物的擴散分布產生重要影響[8-10]。由于物體運動與流場特性高度耦合，固體表面誘導的流體運動具有極強的瞬態特性，大渦模擬（large eddy simulation，LES）方法較RANS方法可以更好地捕捉此類流場特性和污染物運動規律[11-12]。

壓差是隔離艙室負壓系統的重要設計輸入，通常可參考借鑒相關標準規范，但目前相關評估指標都是基于靜態參數，缺少對隔離艙室門開啟的過程中污染物擴散規律的動態評估。試驗研究方法具有費用高、周期長、靈活性差等局限性，CFD方法是解決此類問題的有效手段。因此，本文將以某型應急醫療救援船隔離艙室的負壓系統設計為研究對象，對其典型隔離艙室開展三維實尺度污染物擴散大渦模擬研究，并對比分析空間離散精度以及時間步長對模擬結果帶來的影響；重點對不同壓差工況下，模擬分析隔離艙室門動態開啟過程中的流場特性以及污染物擴散規律；以期為傳染病隔離艙室安全性評估及壓差設計提供參考依據。

1 模擬計算方案

1.1 房間三維模型

本文計算采用的三維模型如圖1所示，其中隔離艙室的長（L1）、寬（W1）、高（H1）分別為4，3.9和2.6 m，緩沖間長（L2）、寬（W2）、高（H2）分別設置為2.9，3.9和2.6 m，中間艙室門的寬度（WD）為0.9 m、高度（HD）為2.0 m，高度方向即為z軸方向。初始時刻，隔離艙室壓力值設定為P1，緩沖室壓力設定為P2，則隔離艙室壓差為 ΔP=P1-P2。由于本文壓差都為負值，下文統一取其絕對值進行表述。艙室門初始時刻處于關閉位置，即艙室門將隔離艙室和緩沖間隔開。計算開始后，艙室門以角速度ω = π/2 rad/s轉向緩沖間一側，在t=1 s時刻停止運動并保持在90°全開位置。

圖 1 房間三維模型示意圖Fig.1 Schematic diagram of three-dimensional room model

1.2 數學模型

將隔離艙室在壓差作用下門開啟過程中污染物的輸運和擴散過程視為不可壓縮流動，且不考慮溫度效應以及重力的影響。因此，采用盒式濾波函數對瞬態控制方程進行濾波后，可以得到如下LES控制方程

本文將污染氣體視為連續介質，采用CO2作為示蹤氣體。初始時刻t= 0 s時，僅在隔離艙室存在一定濃度的均勻分布的CO2，而緩沖間內只有潔凈空氣。在艙室門開啟過程中，空氣和CO2在房門誘導和房間壓差的共同作用下開始進行對流擴散。因此，采用組分輸運模型描述CO2的運動如下：

2 結果及分析

2.1 網格尺度及時間步長敏感性分析

從式（3）可以發現，大渦模擬中亞格子渦黏系數與網格尺度相關，較細密的網格尺寸理論上可以求解更多尺度的流體運動，但是其計算量也會相應增加。因此，需要針對網格尺寸對模擬結果的影響進行研究，這里以房間壓差ΔP=2.5 Pa的工況為例進行說明。如圖2所示，3套網格Mesh 1，Mesh 2和Mesh 3的節點數分別約為148.66萬、341.9萬和570.9萬，并對門框中心線處5個不同時刻的速度分布值進行對比。如圖3所示，模擬結果在z軸方向（艙室高度方向）上呈現出相似的分布規律：在門洞上部區域速度值變化較大；在門洞下部區域速度分布較均勻。同時可以發現，基于Mesh 2的模擬結果與Mesh 1的預測結果有較為明顯的差異；而隨著網格的進一步加密，Mesh 3與Mesh 2的預測結果更加趨于一致。

本文重點研究門開啟過程中隔離艙室污染物進入緩沖間的運動規律，因此，可定義污染物總擴散量M作為模擬評估指標，其計算公式為

圖 2 模型表面網格節點分布Fig.2 Grids at the room model surface

圖 3 不同網格及不同時刻中線處速度分布對比Fig.3 Velocities along the centerline at different time for different meshes

圖 4 不同網格污染物總擴散量隨時間的變化Fig.4 Time histories of total contaminants dispersion for different meshes

因此，本文將基于Mesh 2針對時間步長對LES模擬的影響開展進一步研究。在LES模擬中，為了確保模擬結果的可靠性以及求解過程的穩定性，時間步長Δt應該滿足CFL準則，即

式中：u，v和w分別為計算域中流體速度在x，y，z方向上的分量；Δx, Δy和Δz分別為網格在x，y，z方向的尺度。可根據特征速度和特征長度初步確定LES模擬中的時間步長，因此，本文選取Δt=0.005 s作為參考基準，并設置0.01和0.002 5 s這2種時間步長進行對比分析。

圖5所示為采用3種時間步長設置的LES模擬預測的中線處不同時刻速度分布情況。在門洞上部區域，3種時間步長的預測結果一致性較好；但在下部區域，在t=0.9 s和t=1.2 s時刻，0.002 5 s時間步長的預測結果與其他2種時間步長的模擬結果明顯偏離。圖6所示為采用3種時間步長預測污染物總擴散量隨時間的變化規律，可見Δt=0.01 s計算結果出現明顯偏差，而Δt=0.005 s和Δt=0.002 5 s這2種時間步長的預測結果能夠較好吻合。

由此可見，網格尺度和時間步長都會對LES模擬結果產生一定影響，盡管加密網格以及減小時間步長理論上可以求解更小尺寸的流體流動特征，但也會使得計算量大幅增加，這不利于LES在復雜工程問題上的應用。因此，后續工況的計算基于本節的網格尺寸和時間步長分析結果。

2.2 壓差對污染氣體擴散的影響分析

選取4種壓差工況，即ΔP= 0.25, 2.5, 5.0和10.0 Pa，基于LES對污染物在0～5 s時間段內的擴散規律分別進行模擬計算。計算過程中，分別對圖7所示的4個監測點的速度值隨時間變化規律進行監測，其中點1距離房間地面高度為0.25 m，點2高度為0.5 m，點3高度為1.5 m，點4高度為1.75 m。

圖 5 不同時間步長及不同時刻中線處速度分布對比Fig.5 Velocities along the centerline at different time for different time steps

圖 6 不同時間步長污染物總擴散量隨時間的變化Fig.6 Time histories of total contaminants dispersion for different time steps

圖 7 監測點示意圖Fig.7 Illustration of the monitoring points

圖8～圖11所示分別為4個監測點處在不同壓差工況下的速度變化規律。由圖可見，4個監測點處的速度分布呈現出相同規律：在ΔP= 0.25 Pa工況下，速度先緩慢增大，在t=1 s時（即門停止運動時），速度迅速降低，隨后逐漸降低至最小值并緩慢回升，最終速度值趨于平緩；在ΔP= 2.5 Pa工況，速度值在0～1 s時間段內逐漸降低，隨后速度開始回升并趨于平緩；在ΔP= 5.0和10.0 Pa工況，速度值在0～0.5 s之內達到最小值，隨后便開始迅速增大。值得注意的是，ΔP= 2.5，5.0和10.0 Pa這3個工況下，在所監測的4個點處都出現了明顯的速度峰值，其中ΔP=10.0 Pa時速度峰值出現在t= 1 s附近，隨著壓差減小，速度峰值也開始減小，且對應的峰值時間開始延后。由此可見，在門開啟過程的誘導作用與房間內外不同壓差的配合下，流體運動出現了截然不同的規律，這對污染物的擴散也將產生影響。

圖12～圖15所示為4個不同壓差工況下污染物擴散三維視圖，以及對應時刻距離地面高度z=1 m平面上的污染物擴散俯視圖。污染物濃度采用隔離艙室內初始濃度進行無量綱處理，深藍色表示初始濃度的污染空氣，白色表示潔凈空氣。如圖12所示，ΔP= 0.25 Pa工況下，門開啟過程誘導出一團污染物，隨著時間推演污染物逐漸向艙室外部移動并逐漸減小，在t= 4.0 s時仍有一部分污染物停留在艙門前部區域，若此時房門關閉，該部分污染物將有滯留在隔離艙室外部的風險。將圖13所示ΔP= 2.5 Pa工況的污染物擴散規律與圖12對比，可見ΔP=0.25 Pa 時，污染物在艙室外部的分布范圍明顯減小。在圖14和圖15所示的ΔP= 5.0 和 10.0 Pa工況下，門誘導出的污染物體積明顯減小，且污染物向艙室外部擴散的范圍也得到有效控制。本文模擬結果較好地捕捉到如圖16所示試驗所得的門開啟誘導空氣流動的典型擴散特征[16]。

圖 8 監測點1處不同壓差工況下速度隨時間變化規律Fig.8 Time histories of velocity under different differential pressures at monitoring point 1

圖 10 監測點3處不同壓差工況下速度隨時間變化規律Fig.10 Time histories of velocity under different differential pressures at monitoring point 3

圖 11 監測點4處不同壓差工況下速度隨時間變化規律Fig.11 Time histories of velocity under different differential pressures at monitoring point 4

圖 13 Δ P=2.5Pa時污染物擴散規律Fig.13 Evolution of contaminants dispersion under ΔP=2.5Pa

圖 14 Δ P=5.0Pa時污染物擴散規律Fig.14 Evolution of contaminants dispersion under ΔP=5.0Pa

圖 15 Δ P=10.0Pa時污染物擴散規律Fig.15 Evolution of contaminants dispersion under ΔP=10.0Pa

圖 16 門開啟誘導空氣流動的試驗結果[16]Fig.16 Smoke visualization of door opening induced flow[16]

圖 17 不同壓差工況下污染物總擴散量隨時間的變化Fig.17 Time histories of total contaminants dispersion under different differential pressures

圖17所示為基于式（5）計算的不同壓差工況下進入隔離艙室外部空間的污染物總擴散量隨時間的變化規律。ΔP= 0.25 Pa壓差下，進入緩沖間的污染物在t= 3.3 s達到最大值，隨后開始緩慢降低，在t= 5 s時仍滯留在外部空間的污染物總量仍處于一個較高水平。ΔP= 2.5 Pa壓差下，污染物總擴散量在t= 1.1 s時達到最大值，隨后也呈現下降趨勢，在t= 5 s時緩沖間內尚有少量污染氣體。ΔP= 5.0 和 10.0 Pa工況下，污染物總擴散量都在t= 1 s之內達到最大值，且其峰值明顯小于其他2種工況。在ΔP= 5.0 Pa工況下，t= 2 s時污染物總擴散量為0。在t= 3 s時刻，ΔP= 5.0 Pa工況下擴散出的污染物在壓差作用下也已經全部返回隔離病房。

3 結 論

本文基于大渦模擬（LES）方法結合示蹤氣體組分輸運方程對隔離艙室門開啟過程中污染物的對流擴散規律進行了數值模擬研究。通過對空間離散尺度和時間離散步長開展敏感性分析，并對不同壓差工況下隔離艙室門動態開啟過程中流場速度分布特性及污染物擴散規律進行對比分析，得到如下結論：

1） 網格尺度和時間步長都會對LES模擬結果產生一定影響，以污染物總擴散量為指標對網格尺寸和時間步長進行敏感性分析，可確保本文LES模擬結果的可重復性。

2） 在門開啟過程的誘導與不同壓差的協同作用下，流體呈現出不同的運動規律，速度峰值的大小和出現峰值的時間有明顯差異，壓差越大，速度峰值越大，且出現峰值所需時間越短。

3） 不同壓差工況下，污染物擴散速度和擴散范圍都呈現出較大差異。當壓差 ΔP=0.25Pa時，污染物總擴散量在t=5 s時刻仍維持較高水平；當壓差 ΔP=2.5Pa時，污染物總擴散量峰值大幅下降，但污染物在t=5 s時刻仍未完全消失；當壓差ΔP=5.0Pa時，污染物總擴散量在t=3 s時即變為0；當壓差 ΔP=10.0Pa時，污染物總擴散量在t=2 s時變為0。因此，提高隔離病房負壓值一方面可以減小污染物擴散峰值，同時，也可以縮短已擴散污染物在壓差作用下返回隔離病房的時間。

4 致 謝

本文的數值計算工作得到了武漢大學超級計算中心的支持和幫助。