热颗粒和热辐射共同作用阴燃点燃松针燃料床：数值研究.pdf

收稿日期:2 0 2 2-1 0-1 2;修改日期:2 0 2 2-1 2-2 0基金项目:国家电网有限公司科技项目(5 2 0 0-2 0 2 1 2 6 0 9 9 A-0-0-0 0)作者简介:方伟(1 9 9 5-),中国科学技术大学火灾科学国家重点实验室博士研究生,研究方向为森林火灾,飞火引燃机理及风险。通讯作者:陈海翔,E-m a i l:h x c h e n u s t c.e d u.c n 第3 2卷 第1期2023年3月火 灾 科 学F I R ES A F E T YS C I E N C EV o l.3 2,N o.1M a r.2023文章编号:1 0 0 4-5 3 0 9(2 0 2 3)-0 0 0 1-1 1D O I:1 0.3 9 6 9/j.i s s n.1 0 0 4-5 3 0 9.2 0 2 3.0 1.0 1热颗粒和热辐射共同作用阴燃点燃松针燃料床:数值研究方 伟,张宏民,陈海翔*,彭志红(中国科学技术大学火灾科学国家重点实验室,合肥,2 3 0 0 2 6)摘要:建立了热颗粒和热辐射共同作用下阴燃点燃松针燃料床的二维数值模型,计算得到了热颗粒和热辐射单独点燃以及共同作用点燃的临界条件,分析了燃料床的点燃过程。与前人实验数据对比表明数值模型能够较好预测临界点燃条件。金属热颗粒单独作用时点火所需的临界温度与颗粒直径呈双曲线关系。金属热颗粒与热辐射共同作用时,临界辐射热通量呈现出随颗粒温度和直径的增加而显著减小的趋势。热辐射持续供热能有效维持表层炭氧化反应,两者共同作用下点燃危险性增加。研究同时发现,燃料含水率对共同作用下的临界辐射热通量有较大影响。关键词:金属热颗粒;热辐射;共同作用;阴燃点燃;数值研究中图分类号:X 9 3 2 文献标识码:A0 引言飞火是野火及森林-城镇交界域火灾火场加速扩展的重要机制之一,其难以预测,并可在较大的地理空间诱发大量新的火点。2 0 2 0年3月,四川省西昌大火中,火灾蔓延过程中伴有多处飞火,造成多处火线速燃态势。扑救过程中飞火切断了扑救队伍的撤退路线,最终导致1 9名扑救人员遇难1。2 0 2 0年5月,云南省安宁市发生森林火灾,受火场高温大风、植被易燃等不利因素的影响,火场多次出现飞火引燃现象,形成北、西南等多条火线,导致火灾扑救困难2。最具经济或社会灾难性的野火集中在森林-城镇交界域3,森林-城镇交界域火灾主要以三种方式进行蔓延:火焰直接接触、热辐射和飞火引燃4。研究指出在环境风驱动的大火中,飞火是导致火灾蔓延的主要机制5。有数据表明,至少5 0%的森林-城镇交界域火灾是由森林火灾产生的飞火颗粒引燃引发的6,一 些 火 灾 案 例 中,飞 火 引 燃 的 概 率 更高7。当前全球气候变化导致火灾季节延长8,人类向森林等区域的活动和迁徙导致交界域面积不断扩大,结构功能更趋复杂9,野火发生的可能性不断增长,飞火引燃的潜在危险性也不断增大。飞火引燃具体包括木质飞火颗粒引燃和输电线路故障产生的金属热颗粒引燃。前人对飞火颗粒及金 属 热 颗 粒 点 燃 交 界 域 燃 料 进 行 了 相 关 研 究。H a d d e n等1 0,Z a k等1 1,U r b a n等1 2,W a n g等1 3相继开展了不同实验条件下不同特性球形金属颗粒点燃野外燃料的实验。上述研究均表明点燃临界温度与颗粒直径呈相似的双曲线关系,即颗粒直径增加,点燃所需的临界温度降低。不同于金属热颗粒,飞火颗粒的种类、尺寸等特性以及颗粒本身的燃烧特性受大量因素的影响,因此飞火颗粒作用下的点燃过程更趋复杂。利用风洞,W a t e r m a n和T a k a-t a1 4、M a n z e l l o等1 5-1 7、P l u c i n s k i和A n d e r s o n1 8、G a n t e a u m e等1 9、E i l l s2 0,2 1进行了一系列试验,研究了环境风作用下飞火颗粒点燃可燃材料的模式。在无风条件下,V i e g a s等2 2实验研究了多种类型的飞火颗粒对地中海植被燃料的点燃模式。然而飞火颗粒和金属颗粒点燃可燃物是一个涉及固相和气相物理化学过程的复杂问题2 3,点燃过程受颗粒和燃料床特性以及外部环境条件等众多因素的影响1 8,当前仍然有许多问题亟待解决2 4。另一方面,关于交界域火灾中的热辐射点燃研究,主要围绕火灾产生的热辐射特征2 5,2 6以及森林-城 镇 交 界 域 火 灾 中 火 焰 辐 射 点 燃 可 燃 物 的 特性2 7,2 8来开展。B o o n m e e和Q u i n t i e r e2 7使用锥形热源研究了辐射热流作用时红木自燃的点燃特性,指出当辐射热通量低于4 0kW/m2时,样品在出现明火之前可能会先发生阴燃。M c A l l i s t e r等2 9选取两种地中海地区的代表性野外燃料(M a r i t i m ep i n en e e d l e sa n dK e r m e so a kl e a v e s),研究了燃料床在热辐射作用下的引燃特性,指出燃料床点燃时间的倒数(1/ti g)与入射辐射热流呈线性关系,表现出热薄型固体的热行为。同时他们研究了体积分数对点燃过程的影响,指出松针燃料床体积分数越大,点燃延滞时间越小。现有的研究主要关注于飞火颗粒和热辐射单独作用下点燃建筑构件或野外燃料。然而,研究表明存在飞火颗粒和热辐射共同作用下点燃建筑或野外燃料的实际火灾场景。R a m s a y等3 0指出,在1 9 9 4年的悉尼野火(S y d n e yw i l d l a n df i r e s)中,7 5%的房屋由飞火颗粒点燃,2 5%的房屋在飞火颗粒和热辐射共同作用下点燃。此外,B l a n c h i等7对2 0 0 3年堪培拉山火(C a n b e r r ab u s h f i r e s)的火灾调查报告进行分析并指出,超过9 0%的房屋在飞火颗粒单独袭击以及飞火颗粒和热辐射的共同作用下被点燃。回顾现有关于飞火点燃和热辐射点燃的相关研究,可以发现当前较为缺乏飞火颗粒和热辐射共同作用下点燃固体燃料的相关研究。最近,S u z u k i和M a n z e l l o3 1采用N I S Td r a g o n装置,实验研究了连续飞火颗粒和热辐射共同作用下的点燃过程。然而,该研究为实验统计性质,缺乏对耦合点燃物理过程的详细描述。鉴于飞火点燃的复杂性,以及飞火颗粒特性(尺寸、温度、材料等)描述的困难性,我们采用金属热颗粒代替飞火颗粒,实验研究了金属热颗粒和热辐射共同作用下阴燃点燃松针燃料床的过程和临界条件,分析了耦合点燃机理3 2。本文进一步建立了金属热颗粒和热辐射共同作用下阴燃点燃松针燃料床的二维数值模型,研究辐射热通量强度、金属热颗粒特性(尺寸和温度)和含水率对松针燃料床点燃规律的影响。作为基础研究,此模型的开发有助于揭示金属热颗粒和热辐射共同作用下的点燃机制。1 阴燃点燃数值模型1.1 模型计算域及假设基于C OM S O L M u l t i p h y s i c s多物理仿真软件建立图1所示二维数值模型,模型计算域由松针燃料床和金属热颗粒两部分组成。为简化研究,对模型进行如下假设:(1)松针燃料床是各向同性、均匀的连续多孔介质,热物性为常数3 3。燃料为灰体,且入射辐射热通量直接施加在燃料床的顶部3 4,作为边界条件的一项。(2)热颗粒完全包裹于燃料床内部,且靠近燃料上表面。金属颗粒对燃料床的传热作用仅通过接触热传导来完成,且忽略外部辐射对颗粒温度的影响。金属颗粒性质均一,内部无温度梯度(B i 0.1)3 5。(3)燃料热解产生的气体迅速离开燃料床表面,燃料床内部的气体密度始终为定值,满足理想气体公式。(4)燃料床底表面为绝缘表面,即不存在物质与热量的交换。(5)不考虑燃料体积变化,使用密度缩减模型表示燃料床的消耗情况3 6。图1 二维数值模型计算域示意图F i g.1 S c h e m a t i cd i a g r a mo f t h ec a l c u l a t i o nd o m a i no f t h e2-Dn u m e r i c a lm o d e l2火灾科学 F I R ES A F E T YS C I E N C E 第3 2卷第1期1.2 化学动力学本文建立的二维阴燃点燃模型采用水分蒸发(d r),松针热解(n p),燃料氧化(f o)和炭氧化(c o)四步化学反应机制3 7:F u e lH2OF u e l+vH2O,d rH2O (d r)F u e lv1C h a r+vg1G a s(n p)F u e l+vO2,f oO2v2C h a r+vg2G a s(f o)C h a r+vO2,c oO2v3A s h+vg3G a s(c o)(1)式中v表示各组分的化学计量数。模型假设燃料床由水分、原始燃料和炭三部分组成,基于假设(5),用各组 分 的 容 积 密 度 表 示 燃 料 床 的 构 成,则1=m1/VF B,2=m2/VF B,3=m3/VF B,式中VF B表示燃料床的总体积,为定值;mi和i分别表示各组分的质量和容积密度,下标1,2,3分别表示水分、干燥燃料和炭。每一步反应的反应速率用阿伦尼乌斯公式的形式表示为:1=(1-1)n1A1e(-E1/R T)/T3 4(2)2=(1-2)n2A2e(-E2/R T)(3)3=(1-3)n3A3e(-E3/R T)YO2/0.2 3 nO2,3(4)4=(1-4)n4A4e(-E4/R T)YO2/0.2 3 nO2,4(5)式(2)式(5)中,A和E分别表示各步反应的指前因子和活化能;下标1,2,3,4分别表示水分蒸发、热解、燃料氧化以及炭氧化;YO2表示空气中氧气的质量分数。i为各组分的转化率,表示为式(6),其中m0,i和0,i分别表示各组分的初始质量和实密度。模型的化学动力学参数如表1所示。i=(m0,i-mi)/(m0,i-m1,i)=(0,i-i)/(0,i-1,i)(6)表1 本模型中使用的化学动力学参数T a b l e1 T h ec h e m i c a lk i n e t i cp a r a m e t e r su s e d i nt h i sm o d e l反应i指前因子l gAi/l g(s-1)活化能Ei/k Jm o l-1反应级数ni反应级数nO2,iHi(M J/k g)d r5.7 8a4 9.5 2102.2 5n p6.3 43 81 0 53 80.8 73 803 80.4 1 83 4f o8.7 23 81 2 73 80.5 63 80.7 23 8-4.33 8c o6.5 53 81 2 43 80.6 33 80.6 83 8-1 1.93 8 a(l g(K1/2s-1)1.3 模型控制方程模型的固相总质量方程及各组分的质量守恒方程如式(7)式(1 0)所示:=1+2+3(7)1 t=-11(8)2 t=-22-32(9)3 t=v122+v232-43(1 0)燃料床内部气体的质量守恒方程为:O2YO2 t=O2DO2YO2 -vg,ii(1 1)其中,为燃料床的孔隙率,DO2为氧气的扩散系数。金属热颗粒的能量方程为:pcp Tp t=kpTp (1 2)燃料床能量方程为:c e f f Tf t=ke f fTf +Hii(1 3)式(1 2)和式(1 3)中,下标p和f分别表示金属颗粒和燃料,k为热传导系数,Hi(i=1,2,3,4)为各步反应的反应热,取值如表1所示。1.4 初始及边界条件为求解上述非稳态偏微分方程,模型相关初始条件及边界条件如下。式(1 4)和式(1 5)分别表示燃料床和金属热颗粒的初始温度;式(1 6)表示氧气的初始体积分数;式(1 7)式(1 9)表示组分的初始容积密度,其中MC 和b0分别表示燃料的含水率和初始容积密度。Tf(r,z,0)=2 9 3.1 5K(1 4)3V o l.3 2N o.1方 伟等:热颗粒和热辐射共同作用阴燃点燃松针燃料床:数值研究Tp(r,z,0)=Tp0(1 5)YO2(r,z,0)=0.2 3(1 6)1r,z,0 =MC b0(1 7)2r,z,0 =(1-MC)b0(1 8)3r,z,0 =0(1 9)基于假设(1),燃料床上表面接收辐射热通量,同时与外界环境通过对流和辐射进行热交换,上表面的温度边界条件如式(2 0)所示。式(2 1)式(2 4)分别表示燃料床侧面温度边界条件、燃料床底部温度边界条件、上表面氧浓度边界条件和燃料床底部氧浓度边界条件。ke f f Tf zz=0=q-q c o n v-q r a d=q-hc(Tf-T)-(T4fz=0-T4)(2 0)-ke f f Tf r=q r a d+q c o n v=-hc(Tf-T)-(T4f-T4)(2 1)-ke f f Tf zz=-5c m=0(2 2)YO2=0.2 3(2 3)YO2 zz=-5c m=0(2 4)上述 式 中q 为 燃 料 床 吸 收 的 辐 射 热 通 量;q c o n v=-hc(T-T)为燃料床与周围环境的对流换热,hc为对流换热系数,气体自然对流系数为02 0W/(m2K4);q r a d=(T4-T4)为燃料床对环境的 辐 射 热 损 失,为 史 蒂 芬-玻 尔 兹 曼 常 数,5.6 7 1 0-8W/(m2K4);和为燃料的吸收率和发射率,根据基尔霍夫定律,有=0.7 83 9;T为环境温度,模型取值为2 9 3.1 5K。模型中其他物理参数的取值如表2所示。模型采用自适应三角形网格,计算时间为6 0 0s,时间步长为0.1s。为对模型进行网格敏感性分析,对模型计算域网格进行细化,同时将时间步长缩短为0.0 1s。细化网格并缩短时间步长的模型计算时间比原来增加了1 0倍以上,而计算结果误差在1%以内。因此综合考虑模型计算准确性和计算时间经济性,本模型选取的时间步长和网格质量是合适的。模型考虑表3所示的三组工况,分别是金属热颗粒单独作用、热辐射单独作用以及金属热颗粒和热辐射共同作用下点燃松针燃料床。通过检查阴燃前沿能否从颗粒传播至少1 0mm并在该点显示为稳定或加速的 阴燃传播 来判定 阴 燃 点 燃 是 否 发生3 8,阴 燃 的 特 征 温 度 为5 0 0 1 3。金 属 颗 粒(S S 3 0 4)的熔点为13 9 914 5 5,故本文中金属颗粒温度不超过14 0 0。表2 模型其他相关参数取值T a b l e2 V a l u e so fo t h e rr e l e v a n tp a r a m e t e r s i nt h i sm o d e l物理意义参数符号数值单位对流换热系数hc1 5W/(m2K)燃料密度f6 3 03 4k g/m3燃料比热cf14 7 03 4J/(k gK)燃料热传导系数ke f f0.1 14 0W/(mK)理想气体常数R8.3 1 4J/(m o lK)金属颗粒密度p78 6 03 8k g/m3金属颗粒比热cp5 0 03 8J/(k gK)金属颗粒热传导系数kp2 1.53 8W/(mK)燃料孔隙率0.9 3-氧气质量扩散率Dg21 0-54 1m2/s燃料初始容积密度b04 0k g/m3表3 计算工况T a b l e3 C a l c u l a t i o nc o n d i t i o n s计算工况颗粒直径/mm辐射热通量/kWm-2含水率MC热颗粒单独作用61 4-热辐射单独作用-53 0(2.5)热颗粒和热辐射共同作用61 002 06%,1 5%2 结果与讨论2.1 热颗粒单独作用不同燃料含水率下,不同粒径金属热颗粒单独作用时阴燃点燃的临界温度如图2所示。图2中空心点为数值计算结果;实心误差棒为W a n g等1 3的实验数据,误差棒表示点燃概率为0.0 5和0.9 5所对应的点燃温度,空心方形点线和三角点线分别表示燃料含水率为6%和1 5%的实验工况下的5 0%点燃概 率 边 界。需 要 强 调 的 是,图2中 选 取 的W a n g等1 3的实验在2m/s风速下进行。与W a n g等1 3的研究相比,含水率为6%时,数4火灾科学 F I R ES A F E T YS C I E N C E 第3 2卷第1期值计算所得的临界温度略低于实验结果,但数值结果与实验结果表现出一致的变化趋势,且数值结果均处于误差棒之间,说明模型预测的临界值较为准确。燃料含水率为1 5%的工况下,当颗粒直径9mm,实验结果与数值预测值具有较好的一致性。同时数值结果表明,不同含水率下阴燃点燃所需的临界温度与颗粒直径呈双曲线关系,即颗粒直径增加,临界温度降低,这 种 关 系 对 小 颗 粒 更 加 明 显。这 一 结 论 与W a n g等1 3、U r b a n等3 8,4 2以及Z a k等1 1的研究结果相一致。金属颗粒点燃可燃物是一个涉及固相和气相物理化学过程的复杂问题2 3。金属颗粒能否点燃可燃物受颗粒、燃料床、外部环境条件以及颗粒在燃料床上的着陆特征等众多因素的影响1 8。数值结果对燃料含水率变化表现出更高的敏感性,与实验结果存在一定的误差,但反映出相同的变化规律,因此数值模型仍具有较高的合理性和准确性。图2 热颗粒单独作用时的临界点燃温度。空心方形点线和三角点线分别表示燃料含水率为6%和1 5%的实验工况下,5 0%点燃概率边界(数据来源:W a n g等1 3)F i g.2 T h ec r i t i c a l i g n i t i o nt e m p e r a t u r ew h e nah o tp a r t i c l ea c t sa l o n e.T h eh o l l o ws q u a r ep o i n t l i n ea n dt h eh o l l o wt r i a n g l ep o i n t l i n er e s p e c t i v e l yr e p r e s e n t t h e5 0%i g n i t i o np r o b a b i l i t yb o u n d a r yu n d e r t h e e x p e r i m e n t a lc o n d i t i o n sw i t hf u e lm o i s t u r ec o n t e n t so f 6%a n d1 5%(d a t a s o u r c e:W a n ge t a l.1 3)图3呈现了直径为8mm,初始温度分别为7 0 0和8 0 0的热颗粒单独作用于含水率为6%的松针燃料床时,燃料床不同位置的温度以及氧气质量分图3 直径为8mm,初始温度为(a)7 0 0(b)8 0 0的热颗粒单独作用时,含水率为6%的燃料床不同位置的温度和氧气质量分数随时间变化图F i g.3 T h e t e m p e r a t u r ea n do x y g e nm a s s f r a c t i o na td i f f e r e n tp o s i t i o n so f t h e f u e l b e du n d e r t h ee f f e c t o fah o tp a r t i c l ew i t ht h ed i a m e t e ro f8mma n dt h e i n i t i a lt e m p e r a t u r e so f(a)7 0 0,(b)8 0 0数随时间变化的关系图,图3中坐标表示为(r,z)。当热颗粒初始温度为7 0 0时(图3(a),燃料床在热颗粒的加热作用下先快速升温,热颗粒加热导致燃料发生热解和氧化,因此氧气质量分数呈现出与温度相反的变化趋势。由于颗粒能量较低,放出的能量不足以形成稳定的阴燃前沿,燃料床温度逐渐降低,氧气质量分数逐渐恢复至正常大气浓度。初始温度为8 0 0的热颗粒作用时(图3(b),燃料床不同位置的温度持续升高至阴燃特征温度,纵向各点的氧气质量分数不断下降,说明燃料氧化及炭氧5V o l.3 2N o.1方 伟等:热颗粒和热辐射共同作用阴燃点燃松针燃料床:数值研究化反应持续进行,阴燃前锋持续蔓延。热颗粒单独作用时,由于燃料床表面的对流冷却作用,颗粒在纵向表现出更强的加热能力。2.2 热辐射单独作用热辐射单独作用下,数值计算得到含水率为6%和1 5%的松针燃料床的点燃临界辐射热通量分别为1 6.7kW/m2和1 8.5kW/m2。M i n d y k o w s k i等4 3的实验研究指出,体积分数为0.0 20.0 8间的海岸松针(M a r i t i m ep i n en e e d l e s)燃料床的临界辐射热通量在1 2.5 0kW/m21 5.0 0kW/m2之间,体积分数为0.0 2 3的橡树叶(K e r m e so a k l e a v e s)的临界辐射热通量为1 7.5 0kW/m2。本数值模型使用的松针燃料床的体积分数为0.0 6,数值计算所得的临界辐射热 通量与前 人的实 验 研 究 结 果 较 为接近。图4所示为热辐射单独作用时,不同含水率燃料床表面最高温度随时间的变化。从图4中可以看出,施加热辐射的瞬间,燃料床的表面温度开始快速上升,随后温度变化趋于平缓。外加辐射热通量低于临界值时,温度基本维持不变,说明燃料表面达到了热平衡。当外加辐射热通量高于临界值时,温度逐渐升高达到点燃条件,燃料床被点燃,温度曲线出现向上的突变。辐射热通量越大,出现温度突变的时间越早,说明点燃延滞时间越短。数值结 果 同 时 表 明,高 辐 射 热 通 量 下(2 5kW/m2)燃料床含水率对点燃过程的影响较小,温度变化及点燃过程表现出相似的变化趋势;较低辐射热通量下(2 0kW/m22 2.5kW/m2),含水率为1 5%的燃料床的点燃延滞时间相较含水率为6%的工况有所延长。辐射热通量在点燃临界热通量以下时,燃料床表面达到热平衡时的温度与燃料含水率无关。2.3 热颗粒和热辐射共同作用热颗粒和热辐射共同作用时,不同温度热颗粒点燃不同含水率燃料床所需的临界辐射热通量与热颗粒直径的关系如图5所示。图5中空心划线和实心实线分别表示含水率为6%和1 5%的燃料床点燃所需的临界辐射热通量。6%含水率的工况中,共同作用时热颗粒能够显著降低点燃所需的临界辐射热通量。不同温度的热颗粒,临界辐射热通量随颗粒直径增加而减小,表现出一致的变化趋势。结合图2分析,热辐射的作用也能使热颗粒单独作用下未点燃的工况发生点燃,大大增加了火灾发生的可能图4 热辐射单独作用时,燃料床表面温度随时间的变化(a)M C=6%;(b)M C=1 5%F i g.4 T h e t e m p e r a t u r ee v o l u t i o no f t h e f u e l b e ds u r f a c eu n d e r t h e e f f e c t o f t h e r m a l r a d i a t i o nw i t hd i f f e r e n tr a d i a t i v eh e a t f l u x e s.T h e f u e lm o i s t u r ec o n t e n t sa r e(a)6%a n d(b)1 5%,r e s p e c t i v e l y性,严重威胁森林-城镇交界域的生命财产安全,因此伴随飞火和热辐射共同作用的火灾十分危险。同时燃料含水率对共同作用下的点燃临界辐射热通量有较大影响,含水率越大,点燃所需临界辐射热通量越高。2.4 共同作用下的点燃机理分析图6表示直径为8mm,初始温度为7 0 0 的热颗粒单独作用时,燃料床(MC=6%)不同位置的水分蒸发速率、热解速率、燃料氧化速率和炭氧化速率的变化。颗粒作用的初期阶段,燃料床的温度迅6火灾科学 F I R ES A F E T YS C I E N C E 第3 2卷第1期图5 热颗粒和热辐射共同作用时的临界辐射热通量与颗粒直径及温度的关系,其中空心划线和实心实线分别表示燃料含水率为6%和1 5%的工况F i g.5 T h ec r i t i c a l r a d i a n th e a t f l u xu n d e r t h ec o u p l e da c t i o no fah o tp a r t i c l ea n dt h e r m a l r a d i a t i o n.T h eh o l l o wl i n ea n dt h e s o l i d l i n er e p r e s e n t t h e f u e lm o i s t u r ec o n t e n t so f 6%a n d1 5%r e s p e c t i v e l y图6 初始温度为7 0 0,直径为8mm的热颗粒作用时,含水率为6%的燃料床不同位置的(a)水分蒸发速率、(b)热解速率、(c)燃料氧化速率和(d)炭氧化速率随时间的变化F i g.6 T i m e-d e p e n d e n t c h a n g e s o f(a)w a t e r e v a p o r a t i o nr a t e,(b)p y r o l y s i s r a t e,(c)f u e l o x i d a t i o nr a t e a n d(d)c h a ro x i d a t i o nr a t ea t d i f f e r e n t l o c a t i o n s o f t h e f u e l b e dw i t h t h ew a t e r c o n t e n t o f 6%u n d e r t h ea c t i o no f ah o t p a r t i c l e(Tp=7 0 0a n dd=8m m)7V o l.3 2N o.1方 伟等:热颗粒和热辐射共同作用阴燃点燃松针燃料床:数值研究速升高,颗粒周围发生较为剧烈的水分蒸发反应(图6(a),仅在点(0,-1 1)观察到缓慢的且持续时间较长的热解反应和燃料氧化反应(图6(b)和图6(c),仅在点(0,-1 1)处观察到相当缓慢的炭氧化反应(图6(d)。由于颗粒温度较低,热颗粒周围的燃料发生氧化反应放出的热量不足以克服冷却作用,最终燃料床未能被点燃(温度低于5 0 01 3)。图7所示为热颗粒和热辐射共同作用时(颗粒直径为8mm,初始温度为7 0 0;热辐射强度为2kW/m2),燃料床被成功点燃过程中不同位置的水分蒸发速率、热解速率、燃料氧化速率和炭氧化速率随时间变化的关系图。共同作用时,外界热辐射的持续加热抵消燃料床表面热损的同时,也使燃料床整体温度不断升高,提供了维持化学反应所需的能量,使燃料床的多步反应能够持续进行,最终成功点燃燃料床。结合图6以及图7分析得出,共同作用时燃料的热解、氧化以及炭氧化主要发生在燃料床表面。说明在热颗粒和热辐射共同作用下,先由金属热颗粒加热燃料床形成一定厚度的炭层,辐射提供大量热量使燃料床温度升高,并提供了维持反应所需的热量,最终导致阴燃前锋自维持蔓延,使燃料床点燃。图7 初始温度为7 0 0,直径为8mm的热颗粒与2k W/m2的辐射热通量共同作用时,含水率为6%的燃料床不同位置的(a)水分蒸发速率、(b)热解速率、(c)燃料氧化速率和(d)炭氧化速率随时间的变化F i g.7 T i m e-d e p e n d e n t c h a n g e so f(a)w a t e re v a p o r a t i o nr a t e,(b)p y r o l y s i s r a t e,(c)f u e l o x i d a t i o nr a t ea n d(d)c h a ro x i d a t i o nr a t ea td i f f e r e n t l o c a t i o n so f t h e f u e l b e dw i t ht h ew a t e rc o n t e n t o f6%u n d e r t h ec o u p l e da c t i o no fah o tp a r t i c l e(Tp=7 0 0,d=8mma n dt h e r m a l r a d i a t i o n2k W/m2)3 结论飞火和热辐射是交界域火灾蔓延的主要途径,现有的研究主要关注于两者独立作用下的点燃,当前缺乏两者共同作用下点燃交界域固体可燃材料的相关研究。由于飞火颗粒表征的复杂性,本文利用金属热颗粒代替飞火颗粒,采用水分蒸发(d r),松针热解(n p),燃料氧化(f o)和炭氧化(c o)四步化学8火灾科学 F I R ES A F E T YS C I E N C E 第3 2卷第1期反应机制,建立了热颗粒和热辐射共同作用以及热颗粒和热辐射单独作用点燃松针燃料床的二维数值模型,模拟了燃料床阴燃点燃过程。数值模拟计算得到了热颗粒和热辐射单独作用下点燃以及共同作用点燃的临界条件,分析了燃料床的点燃过程,通过与实验数据进行对比,说明数值模型能够较好地预测临界点燃条件。金属热颗粒与热辐射共同作用时,不同含水率燃料床的临界辐射热通量呈现出随颗粒温度和直径的增加而显著减小的趋势。热辐射持续供热能有效促进表层炭氧化反应维持,两者共同作用下点燃危险性增加。同时研究指出燃料含水率是影响点燃过程的重要因素。参考文献1汪龙华.森林火灾遇难1 9人系遭遇飞火切断退路N.中国青年报,2 0 2 0-0 4-0 3.2张魏桔.云南安宁森林火灾蔓延至楚雄境内3 0 0 0余人正全力扑救N.新华社,2 0 2 0-0 5-1 2.3B o wm a n D M,W i l l i a m s o n G J,A b a t z o g l o u J T,K o l d e nCA,C o c h r a n eM A,S m i t hA M.H u m a ne x p o-s u r e a n ds e n s i t i v i t y t og l o b a l l ye x t r e m ew i l d f i r e e v e n t sJ.N a t u r eE c o l o g y&E v o l u t i o n,2 0 1 7,1(3):0 0 5 8.4C a t o nSE,H a k e sRSP,G o r h a m DJ,Z h o u A X,G o l l n e r M J.R e v i e w o fp a t h w a y sf o rb u i l d i n gf i r es p r e a d i nt h ew i l d l a n du r b a n i n t e r f a c eP a r t I:E x p o s u r ec o n d i t i o n sJ.F i r eT e c h n o l o g y,2 0 1 6,5 3(2):4 2 9-4 7 3.5K o oE M,P a g n iPJ,W e i s eD R,W o y c h e e s eJP.F i r e b r a n d sa n ds p o t t i n g i g n i t i o n i nl a r g e-s c a l ef i r e sJ.I n t e r n a t i o n a l J o u r n a lo f W i l d l a n dF i r e,2 0 1 0,1 9(7):8 1 88 4 3.6M e l lW E,M a n z e l l oSL,M a r a n g h i d e sA,B u t r yD,R e h mRG.T h ew i l d l a n d-u r b a n i n t e r f a c e f i r ep r o b l e m-c u r r e n t a p p r o a c h e s a n dr e s e a r c hn e e d sJ.I n t e r n a t i o n a lJ o u r n a l o fW i l d l a n dF i r e,2 0 1 0,1 9(2):2 3 8-2 5 1.7B l a n c h iR,L e o n a r dJE,L e i c e s t e rRH.L e s s o n s l e a r n tf r o mp o s t-f i r e s u r v e y s a t t h eu r b a n i n t e r f a c e i nA u s t r a l-i aJ.F o r e s tE c o l o g ya n d M a n a g e m e n t,2 0 0 6,2 3 4:S 1 3 9.8J o l l yW M,C o c h r a n eM A,F r e e b o r nP H,H o l d e nZA,B r o w nTJ,W i l l i a m s o nGJ,B o wm a nDMJS.C l i-m a t e-i n d u c e dv a r i a t i o n si ng l o b a lw i l d f i r ed a n g e rf r o m1 9 7 9t o2 0 1 3J.N a t u r eC o mm u n i c a t i o n s,2 0 1 5,6(1):7 5 3 7.9R a d e l o f fVC,H e l m e r sDP,K r a m e rHA,M o c k r i nMH,A l e x a n d r eP M,B a r-M a s s a d aA,B u t s i cV,H a w-b a k e rTJ,M a r t i n u z z iS,S y p h a r dA D,S t e w a r tSI.R a p i dg r o w t ho f t h eU Sw i l d l a n d-u r b a ni n t e r f a c er a i s e sw i l d f i r er i s kJ.P r o c e e d i n g so f t h eN a t i o n a lA c a d e m yo fS c i e n c e so f t h eU n i t e dS t a t e so fAm e r i c a,2 0 1 8,1 1 5(1 3):3 3 1 4-3 3 1 9.1 0H a d d e nR M,S c o t tS,L a u t e n b e r g e rC,F e r n a n d e z-P e l l oAC.I g n i t i o no f c o m b u s t i b l e f u e l b e d sb yh o t p a r-t i c l e s:A ne x p e r i m e n t a la n dt h e o r e t i c a ls t u d yJ.F i r eT e c h n o l o g y,2 0 1 1,4 7(2):3 4 1-3 5 5.1 1Z a kCD,U r b a nJL,F e r n a n d e z-P e l l oC.C h a r a c t