地铁火灾烟气蔓延FDS数值模拟

1、地铁火灾烟气蔓延FDS数值模拟研究韩凯旋 任神河 何娣（1. 长安大学 电子与控制工程学院，陕西 西安 710064； 2.长安大学 信息工程学院，陕西 西安 710064；3. 长安大学 电子与控制工程学院，陕西 西安 710064） 摘 要: 为了有效解决地铁隧道火灾时烟雾分布对人员疏散的影响问题，以西安地铁2号线为研究对象，针对火灾列车停留在隧道中的火灾工况，重点研究不同规模火灾条件下隧道温度、烟雾蔓延范围、可见度等参数的分布情况及变化规律。根据该隧道特定的内部几何构造，建立FDS(Fire Dynamics Simulator，火灾动态模拟仿真模型。利用该软件对隧道开展数值模拟研究，获

2、得了隧道火灾发展及烟气蔓延的一般性规律。关键词：地铁隧道；人员疏散；FDS数值模拟；烟气蔓延FDS numerical simulation on fire spread in subway tunnelHan Kaixuan Ren Shenhe He Di(School of electronic and control engineering, Changan University, xian 710064，China Modern engineer training center, Changan University, xian 710064，China。School of elec

3、tronic and control engineering, Changan University, xian 710064 Abstract: The purpose of this Paper is to effectively solve the problems concerning the smoke distribution, which would influence the passengers safe evacuation, during the fire in subway tunnne1. The research object is the Subway Line

4、Two in xian. According to hypothesis condition in the emergency when the burned subway train stops in the middle of the tunne1, It focuses on the study of the distribution and variation of temperature, smoke spreading scope, visibility and other parameters under different scales fire. Contrary to sp

5、ecific internal geometry of the tunnel, we built the FDS (Fire Dynamics Simulator simulation model of the tunnel. We use the software to execute the numerical simulation study. Lastly, we obtain the general rules of the fire development and the smoke spread rules in the tunnels. Key Words: subway tu

6、nnels; human evacuation; FDS numerical simulation; smoke spread0 引 言鉴于地铁隧道火灾的危害性，国内外学者试图通过研究找出火灾发生的规律，制定一套隧道火灾的预防措施和救援方法。本文利用计算流体动力学软件FDS ( Fire Dynamics Simulator 对西安地铁2号线进行火灾仿真模拟，以Navier-Stokes方程为基础，引入浮力修正的湍流模型、湍流燃烧模型和辐射换热模型，建立了适用于描述地铁隧道内烟气温度分布和气体流动的计算流体动力学模型，实现了对地铁隧道内火灾发生时温度场的数值模拟分析，获取火灾参数。1 公路隧道

7、热释放速率依据瑞典国家测试研究所Ingason.H的火灾热释放理论,现行采用的火灾热释放率数学模型主要有以下几种:1）线性增长模型增长阶段采用线性增长、稳定燃烧阶段保持恒定、下降阶段为线性下降。2）平方增长模型增长阶段采用平方增长、稳定燃烧阶段保持恒定、下降阶段采用指数模型，数学模型函数如表1.1: 表1.1 平方增长模型火灾热释放率平方增长模型热释放率数学模型时间HRR=at20 max HRR=at2maxtmax d HRR=Qmaxea(t- td td 其中:tmax是火灾达到最大热释放率的时间；td维持最大热释放率的时间；Qmax火灾最大热释放率；HRR火灾的热释放率。3）指数增长

8、模型Ingason.H采用一个指数函数来描述火源热释放率的变化，燃料控制的火源热释放率模型依据Numajiri和Furukawa的建议，给出以下数学模型: （1.1）其中:最大热释放率；r，k根据实际条件定出的变量；n选取的变量，无物理意义。以上几种数学模型中，线性模型直观明了，反应了火灾变化的整个过程，但线性热释放率对应的火荷载(热释放率函数对时间求积分即热释放率曲线和坐标轴围成的面积称之为火荷载与实际的火荷载偏差较大。指数模型给出了不同控制条件下热释放率的数学模型，但是这些模型比较复杂，且模型里未定的参数较多，不易选定。Ingason.H给出的平方增长模型，即:增长段采用平方函数，稳定段为

9、常值、衰减段为指数函数，形式简单，容易确定。2 地铁隧道火灾数值模拟理论基础2.1 基本方程FDS以低马赫数的LES方程式来描述受火灾浮力驱动的气体流动现象，其方程式如下: （2.1.1）FDS根据boussinesq approximation将温度、密度与压力区分为空间平均项与振动项，其形式如下所示: （2.1.2）其中 气体密度(kg/m3描述公路隧道火灭发展过程的数学模型建立在N-S方程基础上，在一般坐标系表示为如下形式:质量守恒方程: （2.1.3）组分守恒方程: （2.1.4）动量守恒方程: （2.1.5）能量守恒方程: （2.1.6）气体状态方程: （2.1.7）其中:方程(2.

10、1.5中流体受到的外力f可以包括水喷淋作用时，液滴对流体的阻力作用及除重力外的其他外力；q表示流体因燃烧反应放出的热量。即表示压力项的物质导数。综合上述，FDS由(2.1.6式的能量方程式，(2.1.5式的动量方程式以及(2.1.2的空间平均温度、密度与压力方程式联立求解计算区域的速度、温度、密度与压力。在方程式的数值方法方面，FDS对空间坐标的微分项采用二阶中央差分法，时间的微分项则以显性二阶Runge-Kutta法离散化。上述方程组描述了一般形式下的火灾动力学演化过程，如果不是直接模拟求解，它是不封闭的。若要对特定的火灾场景进行模拟计算，必须对上述方程中表示湍流、燃烧、辐射传热等基本物理过

11、程进行正确的模化，同时还必须给出正确的初始条件和边界条件。2.2 燃烧模型火灾过程几乎都是湍流燃烧过程，火灾中的燃烧可能是扩散燃烧也可能是预混燃烧。在火灾动力学模拟中，采用的湍流燃烧模型有混合组分燃烧模型和有限化学反应速度模型。混合组分燃烧模型假定系统为:燃烧单步不可逆反应的简单化学反应系统。即燃烧反应可以简单表示成: （2.2.1）如果只关心火灾过程的热效应，该模型是简单实用的；若需要研究了解火灾过程中污染物和有毒有害气体的产生，则需要引入包含这些物质产生机理和速率的有限化学反应模型。对于一般碳氢化合物燃烧反应可表示为: （2.2.2） 相应的化学反应速率为: （2.2.3）3 西安地铁2号

12、线火灾动态烟气蔓延数值模拟3.1火灾场景及模型参数设定本文将燃烧火源处理成一个燃烧面积固定的火源。另外，由于研究的是隧道内烟气的动态蔓延过程，故不考虑火灾燃烧物的构成比例以及化学反应引起烟气成分的变化。隧道采用入口纵向通风，风速为2.5m/s。以下参考Ingason.H给出的几种火源热释放率模型，并结合各国所给出不同火灾规模对应的热释率火灾规模值，本文的仿真实验将对火灾场景分别设为30MW、100MW的燃烧进行动态模拟，所取的值应是是合理的。线性模型直观明了，但线性热释放率对应的火荷载与实际的火荷载偏差较大。指数模型给出了不同控制条件下热释放率的数学模型，但是这些模型比较复杂，且模型里未定的参

13、数较多，不易选定。Ingason.H给出的平方增长模型，即:增长段采用平方函数，稳定段为常值、衰减段为指数函数，形式简单，容易确定。3.2 30MW和100Mw热释放率火灾烟气蔓延过程表3.1 不同火源在相同风速下持续时间火源风速增长模型持续时间30MW2.5m/s平方增长模型6.5min100MW2.5m/s平方增长模型9min30MW火源平方增长模型:100MW火源平方增长模型:30MW 100MW图3.1热释放率曲线 30MW 100MW图3.2截面温度场变化情况30MW 100MW图3.3右洞中心线截面烟气可见度场分布结论：在纵向通风风速为2.5m/s的条件下，火灾规模越大，隧道内各点

14、的温度越高，温度场的扩散范围越大；火灾时，隧道内温度有一个急剧增加的过程；不同规模的火灾，隧道内所产生的最高温度从200至1000以上不等，30MW的火灾在火源处火焰的最高温度可达200左右，100MW的火灾在火源处壁顶的最高温度可达1000以上。 隧道内纵向温度分布特点表现为:火源温度最高，随着远离火区温度逐渐降低；随着时间增长，火源附近上游区域的温度高于下游区域的温度；竖直方向呈上层高底层低。参考文献1.覃文清、李风.隧道火灾与防范J.消防科学与技术.2004(1:54.2.兰利敏(译，交通隧道的防火:大规模试验的初步研究成果M.隧道译，1994(23.马积薪(译.公路隧道的防火设备M.隧

15、道译丛.1989(8.4.柴永模.隧道内发生火灾时的温度分布规律初探N.消防技术与产品信息.2002.3.5.王明年、杨其新等.公路隧道火灾温度场的分布规律研究J.地下空间.2003，23(3:317-322.6.王明年、杨其新、袁雪戡、杨忠，公路隧道火灾温度场的分布规律研究J.地下空间，2003(3: 317-3237.史聪灵，钟茂华，涂旭炜等深埋地铁车站火灾实验与数值分析M北京:科学出版社20098.虞利强城市公路隧道建设中火灾数值模拟与灭火救援体系研究D.硕士学位论文.上海:同济大学20089.杨涛公路隧道火灾热释放率及通风方式研究D.硕士学位论文.西安:长安大学200910.蒋玲火灾模拟软件在性能化防火设计中的应用研究D.硕士学位论文.大连:大连理工大学2009文章标题地铁火灾烟气蔓延FDS数值模拟研究作者单位陕西省西安市长安大学作者姓名韩凯旋研究方向隧道通风E-mail hankaixuan7816 联系电话邮 编710064手 信地址陕西省西安市长安大学雁塔校区投稿栏目通信与信息技术文章标题地