半横向通风下公路隧道火灾工况的数值模拟.docxVIP

半横向通风下公路隧道火灾工况的数值模拟 1 隧道全面调查,为疫情防控工作服务 隧道的特性非常重要。外围是土壤或岩石,空间发生火灾时热烟气不能及时排除,热量聚集,温度上升快,可能较早出现轰燃,火灾一旦发生就威胁人们的安全和造成巨大经济损失。较长的公路隧道由于交通事故或汽车发动机故障、汽车油路、油箱跑冒、滴、漏造成火灾的机率较高,并且隧道长度长,纵向坡度大,单纯采用水消防,消防水流至道路面,即可迅速向低处扩散流淌。由于油轻于水且不溶于水,使迅速流向低处的水即成为着火油料的载体, 并继续在水面上燃烧,加剧火势,因此对较长的隧道采用通风防火有着重要意义。 2 火灾模型 隧道火灾通风的研究通常采用试验法和计算机流体力学(CFD)模拟。试验耗费大量的人力物力, CFD模拟具有缩短周期,节省费用等方面的优点,能够直观反映隧道内各处烟气的温度、速度、压力等参数。火灾发生时隧道内气流为三维稳态流动过程,可用下列连续、动量、能量、组分等基本守恒方程描述。 连续方程: ?(ρui)?xi=0(1)?(ρui)?xi=0(1) 动量方程: ?(ρujui)?xj=-?Ρ?xi+??xj(μ?ui?xj-ρˉu′iu′j)+ρgi+Sv(2)?(ρujui)?xj=??P?xi+??xj(μ?ui?xj?ρu′iu′jˉˉˉˉˉˉˉ)+ρgi+Sv(2) 能量方程: ?Τ?t+?(Τuj)?xj=??xj(kρCp?Τ?xj)+˙qcCp(3)?T?t+?(Tuj)?xj=??xj(kρCp?T?xj)+q˙cCp(3) 组分质量守恒方程: ?(ρujφ)?xj=??xi(Γ?φ?xi-ρˉu′iφ′)+˙Μ(4)?(ρujφ)?xj=??xi(Γ?φ?xi?ρu′iφ′ˉˉˉˉˉˉˉ)+M˙(4) 由于浮力作用,紊流模型k-ε方程的k、ε的变化遵循下列方程: ?(ρujk)?xj=??xi(μt?kσk?xi)+Sk(5)?(ρujk)?xj=??xi(μt?kσk?xi)+Sk(5) ?(ρujε)?xj=??xi(μt?εσε?xi)+Sε(6)?(ρujε)?xj=??xi(μt?εσε?xi)+Sε(6) μt=ρCμk2ε(7) Sk=Gk+GB-ρε(8) Sε=C1εk(Gk+GB)(1+C3Rf)-C2ρε2k(9) Gk=μt(?uj?xi+?ui?xi),GB=-βgiμtσt?Τ?xi(10) 流体密度和温度的关系遵循理想气体方程: P=ρRT(11) 式中:ρ为密度;u为速度;S为源项;Cp为定压比热;qc为单位体积热源;Γ为广义扩散系数;φ为组分的质量分数;σ为Prandtl数;β为体积膨胀系数;gi为i方向重力加速度分量(i=x,y,z);T为温度;C为常数。 本文选取体积热源模型(VHS)作为火灾模型,它是最简单的燃烧模型。该模型不考虑化学反应过程,而是设定一个与稳态火源相当的热释放率˙qc,烟气质量释放率˙Μ,将火源模拟成具有固定体积的热源,按受限空间中可能出现的最大释热率计算,它代表了该空间中可能发生的最严重的火灾情况,认为火势不会蔓延。 热释放率: 隧道内火灾规模按中型20MW设计,根据PIARC(1987),20MW的火灾规模取火源尺寸为6m×2m。 烟气质量释放率: ˉΜ=0.2Ly32(12) 式中:L为火源的周长,y为烟层的净高。 3 隧道内通风方案 计算模型取某实际越江隧道的一段,全长600 m,坡度为1%,单洞单向,横截面积如图1所示。一旦发生火灾,隧道采用射流风机纵向送风,竖井内的轴流风机排风,局部开启吊顶风口,从排烟道排除烟气的半横向通风方式。隧道、竖井、排烟道各横断面尺寸严格执行施工图要求。 坐标系选取隧道长度方向从左至右为X轴正向,竖井方向为Y轴的正向,隧道宽Z轴沿X轴对称分部,坐标原点在(0,0,0)处。火源中心位于(40,0,0)处,火源上游截面洞口视为通风入口,风速矢量垂直于入口断面。根据设计参考:吊顶风口大小取1.2 m×1.5 m,间距为30 m。在此情况下研究吊顶风口间距和开启个数的合理性,建立三种通风方案:A.吊顶风口间距30 m且开启火区下游最近的2个风口;B.吊顶风口间距30 m且开启火区下游最近的3个风口;C吊顶风口间距15 m且开启火区下游最近的3个风口。根据方案建立3种几何模型,图2为根据B方案建立的几何模型。 4 结论分析 4.1 风机补给的空气 图3给出了3m/s的通风速度下,火灾发生时,射流风机纵向送风抑止热烟气回流,并且开启吊顶风口排除火区下游烟气的半横向通风方案的X,Y方向合成速度场。由于通风的风量充足,使得风速在X方向为正,无回流现象产生,但是射流风机补给的空气也使火灾沿隧道纵向蔓延的速度加快。同时可注意到,火区气流的速度场在Y轴方向先为正,后为负,如此不断