朱亞迪，陳峰, 3，王子甲, 3

?

基于格子Boltzmann方法的地鐵站內(nèi)火災(zāi)仿真

朱亞迪1, 2，陳峰1, 2, 3，王子甲1, 2, 3

(1. 北京交通大學(xué)土木建筑工程學(xué)院，北京，100044；2. 北京交通大學(xué)北京市軌道交通線路安全與防災(zāi)工程技術(shù)研究中心，北京，100044；3. 北京交通大學(xué)軌道交通安全協(xié)同創(chuàng)新中心，北京，100044)

引入格子Boltzmann方法，構(gòu)建以半拉格朗日方法求解的煙流擴散理論，真實重現(xiàn)煙流溫度演化過程。采用C++語言，實現(xiàn)該理論方法。在此基礎(chǔ)上，建立二維數(shù)值模型并仿真；將其與實驗數(shù)據(jù)進行對比。研究結(jié)果表明：該方法能較好地重現(xiàn)地鐵站內(nèi)火災(zāi)煙流的擴散過程，關(guān)鍵節(jié)點的變化趨勢與實驗結(jié)果一致，距離火源越近仿真結(jié)果與實驗結(jié)果吻合程度越高。

交通工程；火災(zāi)仿真；格子Boltzmann；半拉格朗日方法；地鐵車站

城市軌道交通是一種大運量公共交通方式，軌道交通車站是其客流集散的場所，客流總量大，客流密度高。據(jù)統(tǒng)計北京地鐵2013年最大車站高峰小時進站客流達1.4萬人次[1]。如此大的客流，加上封閉的車站環(huán)境，發(fā)生重大事故后快速有效的疏散措施成為保證站內(nèi)乘客生命安全的關(guān)鍵。據(jù)統(tǒng)計，重大地鐵事故中有57%是火災(zāi)事故[2]。這也使得地鐵站內(nèi)火災(zāi)疏散成為研究熱點。由于火災(zāi)疏散實驗研究成本較高，目前國內(nèi)外在地鐵站內(nèi)火災(zāi)疏散方面的研究主要通過計算機仿真實現(xiàn)。然而，在地鐵站內(nèi)火災(zāi)仿真領(lǐng)域，專用的仿真軟件較少。早期，地鐵站內(nèi)火災(zāi)仿真主要利用Fluent等通用流體軟件進行[3?4]。2000年美國國家標(biāo)準(zhǔn)研究所和建筑火災(zāi)研究實驗室發(fā)布了以數(shù)值方法求解納維斯托克方程的FDS火災(zāi)模擬器，隨后被眾多學(xué)者用于地鐵站內(nèi)火災(zāi)疏散的仿真研究[2]。為了建立更適合地鐵車站環(huán)境的仿真平臺，王英輝[5]基于PHONEICS利用POWER?BUILDER和C++進行二次開發(fā)，建立了SIMU_DDCZRHJ平臺，可以對地鐵車站內(nèi)熱環(huán)境控制以及火災(zāi)條件熱環(huán)境的發(fā)展進行仿真，但是僅對站內(nèi)熱環(huán)境控制效果進行了仿真驗證；關(guān)于熱環(huán)境模擬的核心計算方法仍然采用PHONEICS的算法。在火災(zāi)演化與人員疏散協(xié)同仿真領(lǐng)域，目前研究主要是利用FDS，BRI2002和Smartfire等軟件進行火災(zāi)演變仿真，再將導(dǎo)出數(shù)據(jù)與元胞自動機模型、社會力模型等行人疏散仿真模型結(jié)合實現(xiàn)[6?8]；并且研究主要集中在人員疏散模型的修正完善。這使得火災(zāi)仿真和人員疏散綜合評價較為困難；尤其對建模思路有差異的2種模型融合，還會產(chǎn)生數(shù)據(jù)傳遞錯誤[6]。采用該方法進行城市軌道交通疏散安全性評價時，一般僅根據(jù)規(guī)范標(biāo)準(zhǔn)進行，與考慮車站環(huán)境和環(huán)控系統(tǒng)功能得出的容許疏散時間差別較大[9]。因此，為了實現(xiàn)火災(zāi)演變與人員疏散仿真的融合，需要尋求新的火災(zāi)仿真理論，以期能與人員仿真方法進行結(jié)合。格子Boltzmann方法[10?11]作為一種從格子氣自動機演化而來的介觀流體動理學(xué)仿真方法，與行人仿真領(lǐng)域所采用的格子氣自動機[8]以及元胞自動機[10?14]等仿真理論有一定相似性；因此，利用該方法進行火災(zāi)煙流仿真實現(xiàn)可以更好地融合人員疏散仿真。

1 格子Boltzmann方法

格子Boltzmann方法(LBM)是一種相對較新的計算流體力學(xué)和計算傳熱學(xué)數(shù)值模擬方法，從精細(xì)程度上說屬于介觀方法。該方法將流體離散成流體微團，同時物理區(qū)域也離散成一系列網(wǎng)格，時間同樣也被離散化。離散的流體微團在網(wǎng)格節(jié)點上按照離散的控制方程進行碰撞遷移，最后根據(jù)網(wǎng)格節(jié)點上流體微團的特性值計算出流體的宏觀物理量。

LBM具有易于實現(xiàn)、程序健壯性好以及并行性能好等優(yōu)點，因此該方法在誕生的20多年時間里得到了迅速的發(fā)展，并在火災(zāi)煙流的擴散仿真中逐步應(yīng)用。其中，馬明等[15?18]將LBM用于地下礦巷道中火災(zāi)演變以及瓦斯運移的仿真中；羅金輝[19]將LBM用于煤層氣的擴散和滲透。ALIM等[20?21]為了優(yōu)化煙流擴散的效果建立了D3bQ15模型，并加入因密度引起的浮力項和渦力項，取得了較好的實現(xiàn)效果。德國的Christoph通過融合多松弛LBM和Smagorinsky亞格子模型進行了室內(nèi)空氣湍流的模擬[22]。

1.1 基本理論

氣體系統(tǒng)由大量自由運動的氣體粒子組成，系統(tǒng)的宏觀物理量由大量氣體粒子的統(tǒng)計結(jié)果決定；每個氣體粒子的運動是隨機的，但是平衡狀態(tài)下總體服從麥克斯韋?玻爾茲曼分布，利用該分布計算出每個粒子的狀態(tài)概率，進而計算出來氣體系統(tǒng)的宏觀物理量?；谝陨纤枷胪茖?dǎo)出Boltzmann方程[10]：

式中：為某一時刻粒子位置和動量的分布函數(shù)；為粒子速度；為系統(tǒng)時間；為該時刻粒子的位置向量；為粒子受外力引起的加速度，等式右端項代表碰 撞項。

Boltzmann方程中的碰撞項較為復(fù)雜，求解比較困難；因此，眾多學(xué)者不斷尋求簡單形式來代替碰撞項，其中BGK近似[23]是最為常用的一種近似形式。用BGK近似替換碰撞項后，式(1)就化為Boltzmann?BGK方程：

為實現(xiàn)數(shù)值求解，對Boltzmann?BGK方程進行離散，將粒子速度和分布函數(shù)離散到有限維的空間中；最終得到格子Boltzmann?BGK方程[10]：

由于邊界節(jié)點上部分分布函數(shù)是未知的，因此要采用不同于式(3)的邊界處理格式，常用的有周期性邊界處理格式，反彈格式，非平衡外推格式等[11]。由于非平衡外推格式對于雙分布函數(shù)模型具有較大優(yōu)勢，因此本文選擇非平衡外推格式處理邊界問題。

根據(jù)式(3)確定合適的平衡態(tài)分布函數(shù)，計算得到不同時刻的分布函數(shù)，利用式(4)計算得到系統(tǒng)的宏觀密度和速度。

對于等溫不可壓流體的模擬，可以忽略式(3)中的外力項進行求解；對于本文研究的自然對流系統(tǒng)來說，通常采用Boussinesq假設(shè)，建立溫度場的演化方程，并尋找合適的外力項使溫度場與速度場耦合，建立雙分布函數(shù)模型。

本文中采用溫度場演化方程[10]：

采用的外力項為[24]

1.2 煙流擴散理論

為了真實地模擬煙霧擴散，需要引入煙流擴散的演化機理。FEDKIW等[25]假定煙流密度的對流擴散僅與煙流的運動速度有關(guān)：

另一方面，由于煙流的密度大于空氣密度，會反過來影響流體的速度；因此，需要引入外力項[25]，見式(8)，將其耦合到速度場演化方程中。

式中：為反映密度對速度影響的正常數(shù)；為力的作用方向，通常為(0, 0, 1)。

為了能夠真實地反映火災(zāi)煙流的湍流特點，同時保證算法的穩(wěn)定性；火災(zāi)煙流擴散演化方程的算法實現(xiàn)采用氣象學(xué)中預(yù)測云團運動的無條件穩(wěn)定的半拉格朗日方法[25?26]。該方法的主要思想是該時刻的云團是由上一時刻某一位置的云團遷移而來，因此，

式中：為云團的物理量函數(shù)；為Δ時間粒子運動的距離。

因此，半拉格朗日方法的關(guān)鍵是確定的數(shù)值，

對于的計算根據(jù)節(jié)點位置的格子速度采用二階龍格庫塔方法進行確定，再根據(jù)格子節(jié)點的物理量運用雙線性插值方法得出該時刻云團的對應(yīng)物理量。

2 算法實現(xiàn)

LBM實現(xiàn)的關(guān)鍵在于根據(jù)需要模擬的模型尋求合適的平衡態(tài)分布函數(shù)。以二維仿真為例，選取標(biāo)準(zhǔn)D2Q9模型，采用的平衡態(tài)密度分布函數(shù)[10]和平衡態(tài)溫度分布函數(shù)[24]分別如式(11)和(12)所示。

式中：s為格子聲速。

結(jié)合式(3)~(6)，按照圖1所示的流程就可以實現(xiàn)基本模型進行不可壓非等溫流體的對流仿真。

圖1 LBM實現(xiàn)流程圖

依據(jù)圖1，首先建立了封閉方腔自然對流的物理模型，分別對雷諾數(shù)=104~106范圍的流體進行求解，以溫度變化收斂后的冷壁面的努塞爾數(shù)()為評判標(biāo)準(zhǔn)與基準(zhǔn)解[27]進行對比，驗證基本模型的正確性，見表 1?？梢姡合鄬φ`差較小，基本模型能夠很好地模擬非等溫流體的自然對流。

表1 數(shù)值解與基準(zhǔn)解的比較

Table 1 Comparison of numerical solution and benchmark solution

注：max為最大努塞爾數(shù)；y為最大努塞爾數(shù)出現(xiàn)的位置；ave為努塞爾數(shù)平均值。

根據(jù)上述方法，并結(jié)合式(8)計算出煙流產(chǎn)生的浮力；將其耦合到碰撞過程中的演化方程上。

3 實例分析

采用VC++ 2008實現(xiàn)了以上理論模型，對地鐵站內(nèi)火災(zāi)煙流的擴散演變進行了模擬。為了驗證仿真結(jié)果，結(jié)合文獻[2]中對地鐵站模型進行的燃燒實驗，采用本程序建立簡化的二維模型。根據(jù)計算機仿真結(jié)果，與實驗相關(guān)數(shù)據(jù)進行對比，分析本文所建立的模型在仿真地鐵火災(zāi)煙流擴散演化時的可靠性。

3.1 模型信息

本案例涉及實際全尺寸系統(tǒng)、縮尺實驗?zāi)Ｐ拖到y(tǒng)及其對應(yīng)的LBM仿真模型系統(tǒng)3個系統(tǒng)。其中，縮尺模型系統(tǒng)是根據(jù)模型實驗相似性原理，針對全尺寸系統(tǒng)以1:10為比例建立的小尺寸模型；LBM仿真模型系統(tǒng)是對縮尺模型系統(tǒng)進行網(wǎng)格化而建立。

3.1.1 縮尺模型

縮尺模型是三維模型，如圖2(a)所示；圖中黑色區(qū)域為實驗起火位置及蓄煙腔所處位置，模型的尺寸及相關(guān)參數(shù)見表2。

(a) 實驗?zāi)Ｐ秃唸D；(b) LBM模型

表2 實驗?zāi)Ｐ蛥?shù)

Table 2 Parameters of experiment model

根據(jù)相似性原理，縮尺模型系統(tǒng)溫度m、時間m，與實際全尺寸系統(tǒng)溫度r、時間r關(guān)系為

3.1.2 LBM仿真模型

取圖2(a)中陰影剖面，建立LBM二維仿真模型，如圖2(b)所示，將模型按文獻[28]轉(zhuǎn)換到格子系統(tǒng)，取每個格子尺寸為1，相關(guān)參數(shù)見表 3。LBM模型系統(tǒng)溫度L、時間L，與縮尺模型系統(tǒng)溫度m、時間m關(guān)系為

表3 LBM模型的格子參數(shù)

Table 3 Lattice parameters of LBM model

溫度變化按照式(15)輸入。

（15）

3.2 結(jié)果及分析

火災(zāi)演變一般經(jīng)歷發(fā)展、穩(wěn)定和衰減3個階段，本文僅分析發(fā)展階段。為使表達清晰，以下參量單位均采用全尺寸系統(tǒng)單位描述。

利用本文建立模型，對火災(zāi)發(fā)展階段475 s內(nèi)演變情況進行計算，得到站臺層空間不同時刻的等溫圖如圖 3所示。從圖3可以看出：隨著溫度的增加，空氣呈現(xiàn)非穩(wěn)定流態(tài)，出現(xiàn)紊流渦，與實際的變化情況較為相似。對比文獻[2]中的實驗場景：95 s時站臺煙氣溫度比較低，煙氣比較彌散；并且煙氣在190 s到達扶梯口1處，在360 s到達扶梯口2處。圖 3中呈現(xiàn)的對應(yīng)時間節(jié)點的煙流溫度與實驗描述基本一致，60 ℃的高溫?zé)熈鞣謩e于190 s和360 s到達扶梯口1和2。可見：本文模型能很好地模擬地鐵站內(nèi)火災(zāi)條件下煙流擴散的總體情況。

tr/s：(a) 32；(b) 95；(c) 190；(d) 285；(e) 360；(f) 475

依據(jù)消防安全疏散中關(guān)于煙層臨界危險高度的定義，進一步分析煙層臨界危險高度處溫度的變化情況。本案例系統(tǒng)的危險高度為2 m，因此提取火源正上方以及火源左側(cè)12 m上方，高度2 m處的溫度數(shù)據(jù)，繪制溫度隨時間的變化關(guān)系圖；并與文獻[2]中實驗數(shù)據(jù)進行對比，如圖 4所示。

圖4 仿真與實驗所得溫度隨時間?變化關(guān)系

從圖4可以看出：火源正上方的溫度較高，且仿真所得的溫度變化與實驗結(jié)果在相對誤差允許范圍內(nèi)吻合較好；在火源發(fā)展結(jié)束時刻，即475 s時數(shù)值解溫度為73 ℃，實驗溫度為73 ℃。然而，距離火源水平距離12 m處危險高度點的溫度，隨著時間的推移仿真結(jié)果與實驗結(jié)果呈現(xiàn)出較大差異，但是發(fā)展趨勢基本相同；在規(guī)范規(guī)定疏散極限時間6 min內(nèi)，仿真所得的平均溫度為54 ℃，實驗溫度為36 ℃，相差較大。

圖 4中仿真結(jié)果與實驗結(jié)果在火源左側(cè)12 m處出現(xiàn)明顯差異主要是由于實驗是在三維空間中進行，并且賦予了不同的邊界條件，比如開啟站臺排煙系統(tǒng)、存在樓扶梯接口等，使得煙霧自然擴散受到一定影響，距離火源較遠(yuǎn)的地方溫度下降較快；而本文模型是從二維層面實現(xiàn)的，并且沒有考慮復(fù)雜的邊界條件，煙霧的運動是自然擴散的結(jié)果。另外，本文的煙霧溫度傳播機制仍待進一步完善，使其能夠符合實際情況。

4 結(jié)論

1) 利用本文完善的格子Boltzmann方法可以在二維層面上很好地重現(xiàn)高溫火災(zāi)煙流的總體擴散規(guī)律。

2) 利用格子Boltzmann方法進行數(shù)值求解得到的關(guān)鍵節(jié)點溫度與模型實驗溫度的發(fā)展規(guī)律趨勢吻合較好，并且離火源越近吻合程度越高；但是由于仿真條件的簡化使得仿真結(jié)果與實驗結(jié)果出現(xiàn)一定差異。

3) 本文模型能夠有效地模擬火災(zāi)條件下地鐵站內(nèi)環(huán)境溫度的發(fā)展演變過程。

[1] 閻欣雨, 劉春瑞. 1小時1.4萬人天通苑北如何安檢?[EB/OL]. 2013?12?20. http://www.bjnews.com.cn/news/2013/12/20/298333. html.YAN Xinyu, LIU Chunrui. 1 hour 14 thousand persons how Tiantongyuan north Station conduct safety inspection?[EB/ OL]. 2013?12?20. http://www.bjnews.com.cn/news/2013/12/20/298333.html.

[2] 史聰靈, 鐘茂華, 涂旭煒, 等. 深埋地鐵車站火災(zāi)實驗與數(shù)值分析[M]. 北京: 科學(xué)出版社, 2009: 1?4, 131?159. SHI Zhongling, ZHONG Maohua, TU Xuwei, et al. Fire experiment and numerical analysis of deep-depth subway station[M]. Beijing: Science Press, 2009: 1?4, 131?159.

[3] 李兆文. 地鐵站火災(zāi)煙氣擴散及控制的研究[D]. 南京:南京工業(yè)大學(xué)城市建設(shè)與安全工程學(xué)院, 2005: 20?21. LI Zhaowen. The research on fire smoke spreading and controlling in subway station[D]. Nanjing: Nanjing University of Technology. College of Urban Construction and Safety Engineering, 2005: 20?21.

[4] 趙軼. 地鐵火災(zāi)煙氣數(shù)值模擬及通風(fēng)控制[D]. 大連:大連交通大學(xué)環(huán)境與化學(xué)工程學(xué)院, 2009: 7?8. ZHAO Yi. Numerical simulation of fire smoke and ventilation controlling in subway[D]. Dalian: Dalian Jiaotong University. College of Environmental and Chemical Engineering, 2009: 7?8.

[5] 王英輝. 基于通用CFD的地鐵車站熱環(huán)境模擬系統(tǒng)的開發(fā)與應(yīng)用[D]. 天津:天津大學(xué)環(huán)境科學(xué)與工程學(xué)院, 2004: 46?63. WANG Yinghui. The development and appliance of subway station’s heat environmental emulators based on all-purpose CFD software[D]. Tianjin: Tianjin University. School of Environmental Science and Engineering, 2004: 46?63.

[6] FANG Zhiming, SONG Weiguo, ZHANG Jun, et al. A multi-grid model for evacuation coupling with the effects of fire products[J]. Fire Technology, 2012, 48(1): 91?104.

[7] KAKEGAWA S, NOTAKE H, SEKIZAWA A, et al. Evacuation and smoke movement interactive simulation model[J]. Biochemical Journal, 1995, 312(1): 159?162.

[8] HAO Shuqing, HUANG Hongwei, YUAN Yong. Railway station pedestrian simulation on fire smoke based on Smartfire[C]//International Forum on Information Technology and Applications (IFITA 2009). Los Alamitos: IEEE Computer Society, 2009: 455?457.

[9] 何利英. 中庭式地鐵車站火災(zāi)安全疏散仿真研究[J]. 地下空間與工程學(xué)報, 2010, 6(4): 861?866. HE Liying. Evacuation simulation of atrium metro station in fire accidents[J]. Chinese Journal of Underground Space and Engineering, 2010, 6(4): 861?866.

[10] 何雅玲, 王勇, 李慶. 格子Boltzmann 方法的理論及應(yīng)用[M]. 北京: 科學(xué)出版社, 2009: 33?52. HE Yaling, WANG Yong, LI Qing. Lattice Boltzmann method: theory and applications[M]. Beijing: Science Press, 2009: 33?52.

[11] 郭照立, 鄭楚光. 格子Boltzmann方法的理論及應(yīng)用研究[M]. 北京: 科學(xué)出版社, 2009: 29?46. GUO Zhaoli, ZHENG Chuguang. Theory and applications of Lattice Boltzmann method[M]. Beijing: Science Press, 2009: 29?46.

[12] 胡明偉, 史其信. 行人交通仿真模型與相關(guān)軟件的對比分析[J]. 交通信息與安全, 2009(4): 122?127. HU Mingwei, SHI Qixin. Comparative study of pedestrian simulation model and related software[J]. Journal of Transport Information and Safety, 2009(4): 122?127.

[13] 孫劍, 李克平. 行人運動建模及仿真研究綜述[J]. 計算機仿真, 2008(12): 12?16. SUN Jian, LI Keping. A review of pedestrian modelling and simulation[J]. Computer Simulation, 2008(12): 12?16.

[14] 郭謹(jǐn)一, 劉爽, 陳紹寬, 等. 行人運動仿真研究綜述[J]. 系統(tǒng)仿真學(xué)報, 2008(9): 2237?2242. GUO Jinyi, LIU Shuang, CHEN Shaokuan, et al. Review of pedestrian movement simulation studies[J]. Journal of System Simulation, 2008(9):2237?2242.

[15] 馬明. 基于格子Boltzmann方法的火災(zāi)巷道煙流滾退仿真研究[J]. 科技信息, 2010(3): 595?596, 609. MA Ming. The backflow of smoke model research on fire tunnel based on Lattice Boltzmann method[J]. Science and Technology Information, 2010(3): 595?596, 609.

[16] 馬明. 地下礦巷道火災(zāi)風(fēng)流紊亂現(xiàn)象仿真方法研究[D]. 西安:西安建筑科技大學(xué)材料與礦資學(xué)院, 2010: 8?43. MA Ming. Research on airflow turbulence simulation in fire accident on underground mine tunnel[D]. Xi’an: Xi’an University of Architecture and Technology. College of Materials and Mineral Resources, 2010: 8?43.

[17] 周楠. 基于Boltzmann原理煙霧蔓延的仿真[D]. 西安: 西安建筑科技大學(xué)管理學(xué)院, 2010: 30?48. ZHOU Nan. Research of new method of computer simulation on smoke-flute based on Lattice Boltzmann[D]. Xi’an: Xi’an University of Architecture and Technology. College of Management, 2010: 30?48.

[18] 陸秋琴. 地下煤礦瓦斯運移數(shù)值模擬及積聚危險性評價研究[D]. 西安: 西安建筑科技大學(xué)管理學(xué)院, 2010: 41?118. LU Qiuqin. Research on numeric simulation of gas migration and risk evaluation of gas aggregation in underground coal mine[D]. Xi’an: Xi’an University of Architecture and Technology. College of Management, 2010: 41?118.

[19] 羅金輝. 煤層氣運移LBM模型與井間干擾模擬研究[D]. 北京: 中國礦業(yè)大學(xué)資源與地球科學(xué)學(xué)院, 2012: 53?103. LUO Jinhui. Modelling of CBM migration and well interference via LBM[D]. Beijing: China University of Mining and Technology. School of Resources and Geoscience, 2012: 53?103.

[20] ALIM U R, ENTEZARI A, MOLLER T. The Lattice?Boltzmann method on optimal sampling lattices[J]. IEEE Transactions on Visualization and Computer Graphics, 2009, 15(4): 630?641.

[21] ALIM U R. The Lattice?Boltzmann model for the visual simulation of smoke[D]. Rochester: University of Rochester. Department of Computer Science, 2001: 31?32.

[22] VAN?TREECK C, RANK E, KRAFCZYK M, et al. Extension of a hybrid thermal LBE scheme for large-eddy simulations of turbulent convective flows[J]. Computers & Fluids, 2006, 35(8/9): 863?871.

[23] BHATNAGAR P L, GROSS E P, KROOK M A. A model for collision processes in gases, I: small amplitude processes in charged and neutral one-component systems[J]. Physical Review, 1954, 94: 511?525.

[24] MOHAMAD A A. Lattice?Boltzmann method[M]. London: Springer, 2011: 96?97.

[25] FEDKIW R, STAM J, JENSEN H W. Visual simulation of smoke[C]//Proceedings of the 28th Annual Conference on Computer Graphics and Interactive Techniques.New York: ACM, 2007: 15?22.

[26] STANIFORTH A, COTE J. Semi?Lagrangian integration schemes for atmospheric models: a review[J]. Monthly Weather Review, 1991, 119(9): 2206?2223.

[27] HORTMANN M, PERIC? M, SCHEUERER G. Finite volume multigrid prediction of laminar natural convection: benchmark solutions[J]. International Journal for Numerical Methods in Fluids, 1990, 11(2): 189?207.

[28] LATT J. Choice of units in Lattice?Boltzmann simulations[EB/OL]. 2008?04?01. http://lbmethod.org/\_media/ howtos: lbunits.pdf.

(編輯 趙俊)

Fire simulation in subway station based on Lattice?Boltzmann method

ZHU Yadi1, 2, CHEN Feng1, 2, 3, WANG Zijia1, 2, 3

(1. School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China;2. Research Center of Beijing Rail Transit Line Security and Disaster-resistance Technology,Beijing Jiaotong University, Beijing 100044, China;3. Collaborative Innovation Center for Rail Transit Safety, Beijing Jiaotong University, Beijing 100044, China)

Lattice?Boltzmann method was introduced, and a smoke diffusion theory solved by semi?Lagrangian method was built to show the evolutionary process more truly. Then, this method was implemented by C++. On this basis, a simulation based on a model experiment was conducted. The results show that this method can simulate smoke diffusion process under fire condition in subway station very well, and the key nodes’ temperature evolution has the same tendency with experiment results, and the closer to the fire source, the more the simulation results are in agreement with the experiment results.

transportation engineering; fire simulation; Lattice?Boltzmann; semi?Lagrangian method; subway station

10.11817/j.issn.1672-7207.2016.06.046

U231.96

A

1672?7207(2016)06?2160?06

2015?06?04；

2015?09?22

國家自然科學(xué)基金青年科學(xué)基金資助項目(51408029)(Project(51408029) supported by the National Natural Science Youth Foundation of China)

王子甲，博士，講師，從事城市軌道交通的研究；E-mail：zjwang@bjtu.edu.cn