曾清珣 毕海权 苟红松 王 菁 刘方吉
(1.西南交通大学机械工程学院 成都 610031;2.中铁二院工程集团有限责任公司 成都 610031)
中国公路网不断向崇山峻岭、离岸深水区延伸,公路隧道总量和建设规模持续增大。近40年来,中国公路隧道由单洞2 车道逐步向双洞4 车道、6 车道、8 车道隧道发展[1]。随着所建隧道的长度及宽度日趋增大,隧道内发生火灾时的危害也随之增大。隧道内通风排烟系统多采用自然排烟[2]、纵向通风排烟[3]及顶部排烟[4]方式。但侧向排烟方式具有节省隧道净高、降低造价等优点,因此这一排烟方式在高度方向受限制、两孔一管廊隧道中被广泛应用[5]。
热烟气的浮力效应为垂直向上[6],顶部设置排烟口有利于烟气的排出,而当排烟口位于隧道侧壁面时,排烟口开口方向与烟气浮力方向不同,在侧向排烟方式下隧道内的烟气特性将不同于以往研究。目前关于侧向排烟的研究主要是基于单洞3 车道隧道,研究排烟量[7]、排烟口设置[8]、纵向风速[9]等影响因素的最佳设置以及其对人员疏散安全的影响,而对采用侧向排烟系统的隧道在断面宽度改变后的烟气特性方面研究不足。
本文通过数值模拟研究隧道断面为不同宽度时,各行车道处拱顶处、人员高度处烟气温度的分布规律以及侧向排烟口的排烟效率,研究结果以期为火灾探测、消防救援、隧道结构设计等提供参考。
火灾的发生、发展为非稳态、非定常变化的复杂的三维过程,但是流体的流动都遵循质量、动量、能量及组分守恒等基本定律,在FDS 中的基本控制方程如下[10]:
质量守恒方程:
动量守恒方程:
组分守恒方程:
能量守恒方程:
气体状态方程:
式中,ρ为密度,kg/m3;t为时间,s;为速度矢量,m/s;p为环境压力,Pa;g为重力加速度,m/s2;f为除去重力的外力,N;τij为粘性应力张量;im′ 为单位体积内第i种组分的质量生成率,kg/(m3·s);Di为i组分质量扩散系数,m2/s;Yi为i组分体积分数;qr辐射热通量向量,W/m2;hs为显焓,J/kg;T为热力学温度,K;k为导热系数,W/(m·K)。R为气体常数,J/(mol·K);为气体混合物的分子量,kg/mol。
本文根据某隧道实际尺寸建立计算模型,为了保证在模拟中隧道内气流能达到充分的发展,建立的隧道长度为600m,模拟时间设置为500s。隧道出入口设置为Open 边界,排烟道两端开口设置为Exhaust 边界,环境温度设置为20℃。隧道壁面材料为混凝土,该材料导热系数为1.8W/(m·K),比热为1.04kJ/(kg·K)。以隧道正中位置为中心,在两侧对称位置各布置两个排烟口,排烟口尺寸为3m×2m,排烟口间距为60m,排烟口底部距行车道的垂直高度为4.6m。隧道模型示意图(以单洞4车道隧道为例)如图1所示。
图1 隧道模型示意图(以单洞4 车道隧道为例)Fig.1 Schematic diagram of the tunnel model(a single-hole four-lane tunnel as an example)
在进行数值计算时,合理的网格尺寸设置可以在保证计算精度的同时节约计算时间。前人研究指出,通过1/4 到1/16 倍的火源特征直径D*来确定着火点附近的网格尺寸可以与试验结果很好的吻合[11],并且相关研究表明,火源区、浮力羽流区、烟气层区对于网格密度的敏感程度依次递减,加密火源区的网格可以在控制计算时间的同时有效提高计算精度[12]。火源特征直径D*由公式(6)计算得出:
式中,D*为火源特征直径,m;Q为总热释放速率,kW;ρ∞为环境空气密度,此处取1.204kg/m3;Cp为环境空气比热,此处取1.005,kJ/kg·K;T∞为环境空气温度,K。
计算出火源功率30MW 所对应的火源特征直径为3.7m。考虑到火源附近温度梯度较大,本文中火源附近纵向50m 区域内采用0.2m 的网格进行加密,沿隧道纵向其他区域的网格尺寸为0.4m。
在同一时间,考虑隧道内只发生一次火灾,并且发生火灾的位置均位于隧道正中位置。火源功率为30MW,尺寸为8m×2.4m,使用软件自带的甲烷燃烧反应进行隧道火灾模拟。由于本文主要研究火灾发展到最大并达到稳定时隧道内的烟气特性,所以采用恒定火源。
在隧道中,拱顶及人员高度处的温度分布、排烟口处的排烟效率反映了排烟系统的排烟效果、系统性能,关乎发生火灾时人员的安全疏散。因此,本文对这些关键参数进行研究,以分析在不同断面宽度的隧道中采用侧向排烟时,隧道内的烟气特性。温度测点在拱顶下方以及人员高度处沿隧道纵向布置,在隧道横向方向上,温度测点布置在第一行车道的内边沿以及各个行车道的正中位置。对于排烟效率,本文选取CO2作为参考依据,计算CO2的生成量和测定排烟口处的CO2排出量,以此得到排烟口的排烟效率。各个排烟孔的排烟效率为单位时间内该排烟孔的烟气排出量占火源产生烟气总量的百分比,排烟系统总排烟效率为单位时间内所有排烟孔的烟气排出量占火源产生烟气总量的百分比,即所有排烟孔排烟效率之和。
式中:η为排烟效率,%;me为排烟孔的烟气排出量,kg/s;mei为第i个排烟孔的烟气排出量,kg/s;mp为烟气产生量,kg/s。
通过计算得到30MW 火灾在单位时间甲烷燃烧热值为52043.68kJ/kg,对应的CO2生成量为1.56kg/s。
为研究不同隧道宽度采用侧向排烟时隧道内的烟气特性,对单洞2 车道、3 车道、4 车道、5车道、6 车道隧道进行了计算模型建立。模型隧道宽度分别10.5m、14.4m、17m、23.1m 和28m,高度均为7m,如图2所示。
图2 不同隧道宽度数值计算模型横截面Fig.2 Cross-sections of numerical calculation models for different tunnel widths
如表1所示为不同隧道宽度数值模拟工况设置。
表1 不同隧道宽度数值计算工况Table 1 Numerical calculation conditions of different tunnel widths
如图3所示为不同隧道宽度下各行车道上方拱顶烟气温度对比。从图中可以看出,同一侧向排烟量下,隧道越宽,同一行车道上方的拱顶烟气温度越低。特别是在两个排烟口之间,温差明显,在y=60m 处,单洞2 车道隧道与6 车道隧道拱顶烟气温度相差50℃左右。当烟气经过第二个排烟口后,温差减小,在y=120m 处,温差降至20℃左右。行车道1 边沿(x=1m)及行车道1(x=2.8m)离排烟口较近,两个排烟口附近发生了温度的陡降(见图a、b),这是由于排烟口不断抽吸走高温烟气。单洞2 车道、3 车道隧道在第二个排烟口(y=90m)处温度陡降较其余工况更为明显,并且在单洞2 车道隧道中,烟气在经过第二个排烟口后,温度回升,呈现出v 型温度曲线。
图3 各行车道上方拱顶烟气温度纵向分布Fig.3 Longitudinal distribution of flue gas temperature on the vault above each carriageway
如表2所示为拱顶烟气温度在火源所在位置处横向方向上的分布情况,表中L与R分别代表在火源左侧及右侧各行车道上方的温度测点,其中火源左侧(L)为排烟道与火源之间的横向位置。下标代表了与火源横向距离,下标越大,离火源横向距离越远,以4 车道隧道为例,L1、L2及R1、R2位置示意图如图4所示。由表2 可知,随着横向位置远离火源,拱顶烟气温度逐渐降低。由于左侧测点离排烟道更近,受侧向排烟作用更显著,温度下降程度较右侧更大,火源处的拱顶烟气温度横向分布呈“左低右高”的非对称分布形态。
表2 火源处横向方向上拱顶烟气温度Table 2 The vault smoke temperature in the transverse direction at the fire source
图4 L1、L2 及R1、R2 位置示意Fig.4 Location of L1,L2,R1,R2
如图5所示为人员高度处温度纵向分布,从图中可以看出,在两个排烟口之间,人员高度处温度随着隧道宽度的增加而降低,在经过第二个排烟口后,规律相反。由于排烟口附近的抽吸力扰动了烟气分层,加强了卷吸,以及在隧道纵向方向上烟气蔓延的最远端,上层烟气与下层冷空气运动方向相反,剪切力导致卷吸剧烈,使人员高度处温度在这两个位置有较明显的增幅,对于人员疏散而言,危险性较高。
图5 人员高度处温度纵向分布Fig.5 Longitudinal distribution of temperature at personnel height
如图6所示为不同宽度隧道采用侧向排烟方式时的排烟效率。从图6 中可以看出,在该排烟方式下,排烟效率随着隧道宽度增加而降低。从单洞2 车道隧道至6 车道隧道,排烟效率从89.2%减小为61.7%,降低了31%。从2 车道隧道至3 车道隧道及从3 车道隧道至4 车道隧道,排烟效率的衰减率均为8%,在此之后,衰减率随着隧道宽度的增加而增加。
图6 排烟效率Fig.6 Smoke exhaust efficiency
本文采用数值计算方法对采用侧向排烟系统的不同宽度隧道内的烟气特性进行了数值模拟研究,分析各行车道上方拱顶处、人员高度处的温度分布及排烟口处排烟效率,主要结论如下:
(1)增加隧道宽度有利于降低拱顶烟气温度、人员高度处温度。但是增加隧道宽度也导致侧向排烟效率下降,特别当隧道行车道增加至4 车道及以上时,排烟效率的衰减率也逐渐增大。
(2)侧向排烟系统使火源处横向方向上拱顶温度呈非对称分布。朝向排烟道方向的拱顶温度较远离排烟道方向更低。距离排烟风机较近的排烟口抽吸力较强,在经过该排烟口后,不同隧道断面宽度下,同一行车道上方拱顶烟气温度之间的温差减弱。对于人员疏散安全来说,在烟气与冷空气卷吸严重的排烟口附近及隧道纵向方向上烟气蔓延的最远端,温度存在较明显的增幅,危险性较大。