FDS在原木火災研究中的應用探討

薛 偉，耿志偉 (東北林業大學工程技術學院，黑龍江哈爾濱150040)

貯木場具有木材生產、貯存和銷售三大特性，在木材的產銷之間起調節和緩沖作用。因木材生產季節性強，貯木場內長時間存放大量原木、原條和木材加工產品，一旦發生火災，經濟損失非常慘重。因此，為預防火災的發生并在火災發生后能夠及時控制火情，一系列針對原木及其楞堆燃燒規律的研究也相應展開。目前主要的研究方法為試驗模擬與數值模擬，而FDS(Fire Dynamics Simulator，簡稱FDS)作為重要的場模擬軟件之一，在原木以及楞堆火災的數值模擬中發揮了重要的作用。FDS數值模擬在國內的使用與研究不斷發展，其主要的應用領域已涉及到機場、隧道、劇院、核電廠及大型廠房等的火災［1］。但是關于FDS在原木火災研究的應用并沒有系統的探討及總結。筆者從FDS軟件的應用出發，歸納FDS在原木火災研究中的應用現狀，并對FDS在原木火災數值模擬應用過程中的問題進行探討，以促進原木火災發展規律研究的開展以及FDS在原木火災研究中的廣泛運用。

1 FDS的基本介紹

FDS模擬是美國國家標準研究所(NIST)建筑與火災研究實驗室(BFRL)開發的三維計算流體動力學軟件。結合Smokeview這一可視化程序，對火災湍流流動過程進行模擬并直觀地展示模擬結果。FDS的主要使用流程如圖1所示。

場模擬中處理湍流的方法可分為直接數值模擬(Direct Numerieal Simulation，簡稱DNS)、雷諾平均模擬(Rynld Averaged Navier-stokes Simulation，簡稱 RANS)以及大渦模擬(Large Eddy Simulation，簡稱LES)。DNS可得到所有時間尺度和空間尺度上的湍流特征量，但其計算量巨大，因此一般用于解決簡單的小尺度的流動問題。RANS求解的最小時間尺度和空間尺度是和流動特征值可比的，對網格劃分要求不高，計算量相對較小，可用于復雜物理場景的模擬。LES直接求解表征流體流動特性的含能大渦，而對隨機性較強的耗散小渦則采用亞網格模型進行模擬，既克服了RANS不能反映流體瞬時特性的缺點，又減少了數值模擬的計算量，目前正處于發展的活躍階段。

原木及其楞堆的火災多發生在貯木場這種開放環境中，計算區域較大，物理場景也較為復雜，且原木火災是集煙氣流動、傳熱、輻射、燃燒、炭化及其相互作用的復雜燃燒過程，模擬計算量大，而LES不僅可以滿足場景模擬要求而且計算量相對較小。因此，目前在原木火災的數值模擬中主要采用的是LES方法。

2 FDS在原木火災研究中的應用探討

FDS自2000年公開發布以來就受到普遍的關注。據統計，該模型大約有50%應用于煙氣控制系統，50%應用于居民和工業建筑火災后的重建和修復［2］。隨著貯木場消防安全管理的發展，原木及其楞堆燃燒的理論有待進一步深入研究的情況下，原木火災的實驗研究與數值模擬也逐漸展開。

2.1 FDS在原木火災研究中的應用現狀 最初卞偉等［3－4］在森工貯木場楞堆火災模擬模型分析與評價中分析討論了FDS場模型的可行性，之后FDS在原木火災研究中的應用逐漸多起來。總體來看，FDS在原木火災研究中的應用主要在以下方面。

2.1.1 原木楞堆火災溫度場的數值模擬。張光俊等［5－6］通過室外原木燃燒模擬試驗，驗證數值模擬程序的通用性與可靠性。試驗對象為直徑在0.10～0.12 m，長度為1 m的云杉原木，其平均含水率為20%，楞堆由下至上依次為10－9－8－7－6根原木排列，環境溫度13℃，平均風速3 m/s。

運用LES方法對大空間的原木楞堆燃燒進行模擬，模擬結果分析包括楞堆的端面溫度與楞堆的燃燒周圍溫度場2個方面。模擬燃燒過程采用固定熱釋放速率。所建立的適用于原木火災研究的數學物理模型，包括基本方程組、輻射模型、燃燒模型。對比分析的結果表明，在輸入參數合理的情況下FDS可以對貯木場楞堆燃燒進行合理的預測，而且模擬曲線比實測溫度變化平穩且規律性更強，因此將FDS應用于貯木場防火領域是可行的，且對于研究原木火災的發展規律具有現實意義。

2.1.2 原木楞堆火災煙氣溫度變化的數值模擬。辛穎等［7］采用FDS軟件對原木楞堆在長方體單室結構的煙氣變化規律進行了數值模擬。單室(10 m×10 m×7 m)內有一個可燃原木楞堆，頂部設置一個矩形通風口，右側設置一個2 m×4 m的門，墻壁厚度為0.025 m;可燃單根原木尺寸為0.25 m×0.25 m ×1.00 m，位于單室內中央。

通過提取楞堆中心線所在平面的溫度切片，原木在燃燒過程時會在楞堆上方形成連續火焰區、間斷火焰區和無火焰熱氣流區。楞堆火源上方形成高溫煙氣羽流，并在流動過程中與沿程區域不斷進行熱交換，使得楞堆周圍環境溫度不斷升高。通過提取不同高度Z=2、4、6 m煙氣溫度，原木楞堆火災經過一段時間的燃燒后會進入穩定燃燒的狀態(該試驗條件下是300 s)，此時會有比較穩定的煙氣分層現象，而且在原木楞堆燃燒的垂直往上的空間環境中，火災煙氣的溫度與水平方向上距離的關系大致呈現出雙Lorentzian曲線的規律，且距離原木楞堆越近，曲線的峰值會越大，但是并沒有能夠得出溫度隨垂直高度變化的規律公式。

2.1.3 細水霧作用下滅火過程的數值模擬。李雙玉等［8］對細水霧作用下原木楞堆燃燒溫度場進行了試驗研究與數值模擬，以落葉松為試驗原料，選用長度為1 m、直徑0.25～0.30 m的原木，其含水率在18%～24%。原木堆碼方式為逐層減一的間壓式，最底層8根，最上層3根，共6層33根原木。

模擬過程中采用熱釋放速率隨時間變化而改變的火源。考慮到原木楞堆從初起到蔓延需1 100 s左右，試驗在1 200 s時施加細水霧，并計算溫度下降速率。在比照試驗過程對細水霧滅火過程進行模擬時，重點分析細水霧施加30 s內原木火災的變化，并動態展示了原木火災的滅火過程，計算不同狀況下的溫度下降速率。同時，運用FDS分析不同霧滴直徑、不同噴霧速度和噴頭距離原木楞堆頂面不同距離等情況下熱煙氣溫度的變化規律，深化了對細水霧滅火機理的認識，為細水霧化消防技術用于貯木場火災提供了科學指導和理論依據。

2.2 FDS在原木火災研究應用中的探討 火災研究最直接最有效的方式就是火災場景的再現，但大尺度火災試驗也是破壞性的和消耗性的，尤其是原木楞堆火災的模擬試驗，不僅要消耗大量木材與測試元件，而且火災場景的再現也不完美。但借助于FDS數值模擬結合Smokeview可視化程序，不僅可以將所求解的原木火災溫度場、煙氣場以形象、直觀地方式表現出來，而且降低了試驗成本。但統觀研究發現在原木火災的FDS數值模擬過程中存在以下問題。

2.2.1 FDS模擬原木火災范圍受限。實際中原木是圓柱形，但FDS只給出長方體的構造方法，且過于細致的構造原木又極大地增加了編寫程序的難度與運行程序時的計算量。因此在原木火災的數值模擬過程中都是將原木進行簡化，以長方體來代替，這樣的模型模擬出來的結果必然存在一定的誤差。隨其輔助軟件的發展，如Firesim建模軟件，可以構造較為復雜的場景，因此在將來運用FDS模擬原木火災可以適當采用輔助軟件，以提高其模擬精度與可靠性。

2.2.2 FDS對室外原木火災環境的模擬不夠精確。FDS通過時間平均方程求解，對流場的紊流解析存在誤差，無法精確模擬強通風或火旋風等現象，而燃燒時間與燃燒速度、風速大小及方向等都是影響結果的重要因素。在原木楞堆火災溫度場的數值模擬中就發現受傳熱方式、通風條件等不確定條件的影響，模擬曲線的最高溫度比實際測量要高，最高溫度差異最大可達100℃，所以在計算區域、網格劃分、邊界條件的設置中要檢驗其合理性，使得FDS壁面條件設定盡量滿足對室外之一特定的空間環境的模擬要求。

2.2.3 火源熱釋放速率描述曲線的準確度不夠。FDS軟件中火源設置采用了體積熱源法，以降低模擬的難度。如果將火源設計為一個有限大小的固定熱源，則是將火源熱釋放速率設計為穩態火源或按t2規律變化。但是不恰當的火源性質以及熱釋放速率曲線也會使得模擬結果誤差很大，或者不能夠精確描述火場分布。火源的熱釋放速率變化規律影響火災的蔓延過程，同時輻射換熱量超過一定的限度時也會加速火災的蔓延，甚至引起轟燃現象。在模擬過程中，應當對火源設定的準確性予以檢驗，以滿足原木火災計算區域尺度大、各控制方程耦合求解的要求，同時也能注意火源的燃燒及蔓延過程對于煙氣流動的影響。

2.2.4 網格劃分單一。在FDS模擬過程中根據模擬對象的特征尺度來劃分網格。為保證計算的精度，一般是選取比特征尺度小1～2個數量級的網格維度。因此在原木火災的數值模擬中，網格大小均小了1個數量級。但在實際的模擬中，不同位置的計算區域所起的作用不同，受網格大小影響也不同。比如火源周圍及其他關鍵部位的網格尺寸應該足夠小以確保計算的正確性;而遠離火源受到影響較小的區域可以適當地放大網格尺寸，以減少計算時間。因此，應當根據實際情況，采用多種網格劃分的方式來減少模擬的計算量。

2.2.5 FDS軟件中木材性質的種類少。木材的種類、含水率等木材性質都對模擬結果有影響，但FDS軟件只給出了一種木材以供參考，所以有必要改寫其屬性，使它的化學性質更加接近于試驗所選原木的燃燒特性，使模擬試驗更加接近于實際燃燒情況。因此，在固定火源燃燒時間及燃燒溫度的情況下，合理設置木材燃燒性質也是實現正確模擬試驗的關鍵。

2.2.6 FDS在原木火災的數值模擬中應用范圍有待拓展。FDS在原木火災研究中的應用主要在原木楞堆火災溫度場、煙氣溫度變化、細水霧滅火過程幾方面，主要觀察其溫度場、煙氣溫度變化、煙氣層分布以及一氧化碳等物質的分布情況，對風的考慮較少并且目前還沒有對原木火災速度場變化的研究，如質量流率、卷吸速度等。將來的研究中可以引入不同的氣流組織方式，對在不同風速下的原木火災進行模擬，求解其速度場的分布。

3 結語

FDS軟件在原木火災研究中的應用正在逐漸發展與完善，不僅借助Smokeview可視化軟件實現了原木火災場景的再現，將所求解的溫度場以直觀的方式表現出來，同時降低了原木火災試驗研究的消耗性，并得到更多的數據以供參考。FDS成本低、速度快，能夠適應各種情況和計算條件，可以從可能的設計方案快速地選擇出最優方案，但FDS也存在建模不靈活、火源設置不確切等問題，因此在分析建模、設置火災場景時，一定要以科學負責的態度對待，重視關鍵性因素與計算細節。在今后的研究中不斷修正原木火災的FDS模型，充分發揮FDS的優勢，為貯木場原木火災理論的豐富、先進消防技術的運用提供了一定的科學指導和理論依據。

［1］叢北華，廖光煊，韋亞星.計算機模擬在火災科學與工程研究中的應用［J］.防災減災工程學報，2003，23(2):63 －68.

［2］趙澤文.常用火災數值模擬軟件FDS的應用探討［J］.山西建筑，2010，36(13):364－365.

［3］卞偉，薛偉.森工貯木場楞堆火災模擬模型分析與評價［J］.林業機械與木工設備，2006(4):36－38.

［4］FLOYD J E，MCGRATTAN K B，HOSTIKKA S.CFD fire simulation using mixture fraction combustion and finite volume radiative heat transfer［J］.Journal of Fire Protection Engineering，2003(13):11 －36.

［5］薛偉，張光俊.FDS 火災模擬與應用［J］.吉林林業科技，2006，35(6):18－20.

［6］張光俊.貯木場楞堆燃燒的數值模擬［D］.哈爾濱:東北林業大學，2007.

［7］辛穎，薛偉.原木楞堆火災煙溫變化數值模擬研究［J］.消防科學與技術，2011，30(3):183 －185.

［8］李雙玉.細水霧作用下原木楞堆燃燒過程實驗研究［D］.哈爾濱:東北林業大學，2013.