聚光阴燃点火的特性：实验和模拟

林少润，王思言，孙培艺，袁 晗，黄鑫炎

聚光阴燃点火的特性：实验和模拟

林少润1, 2，王思言1，孙培艺1，袁 晗1，黄鑫炎1, 2

(1. 香港理工大学消防安全工程研究中心，香港 999077；2. 香港理工大学深圳研究院，深圳 518057)

阳光经汇聚后形成的光斑可点燃常见的可燃物，是生活中潜在的火灾风险，威胁人们的生命和财产安全．本文首先通过小尺度实验研究了多层薄纸在汇聚阳光下的阴燃点火现象．实验将直射的太阳光通过一个直径为150mm、焦距为108mm的玻璃球汇聚后作用在燃料上，并在焦距范围内通过改变薄纸的位置以控制光斑的直径大小(1.5～20.0mm)．实验发现：随着光斑直径减小，燃料内部热传导散热的影响增强，点火所需的临界辐射热通量增大．然后，研究进一步基于开源代码Gpyro建立了阴燃点火和火蔓延的二维模型，成功再现实验现象并验证实验结果．本文量化了太阳聚光后的火灾风险，以期对未来的火灾防护工作提供一定的参考意义和科学 指导．

阴燃；临界热通量；光斑；加热面积；火灾

俗话说，“万物生长靠太阳”，太阳光是一切生命的源泉．然而，阳光经汇聚后会形成拥有巨大能量的光斑，可点燃生活中常见的可燃物，进而造成巨大的经济损失和人员伤亡．古希腊的阿基米德曾利用一面巨大的凹面镜汇聚阳光并用其点燃了罗马的战船，在著名的叙拉古战役中立下大功．这段历史虽被后代学者证实为夸大之词，但聚光镜能通过光学原理聚光并引燃周围的可燃物却是不争事实[1]．然而，人们对于“光与镜”相遇时的火灾风险常常不以为意．反光镜、玻璃球、圆形鱼缸等这类日常生活中不起眼的物品，本身虽不易燃，却很可能因为不当摆放成为现实版的“阿基米德之镜”，成为火灾的点火源，对人类生命和财产安全造成巨大威胁．从2010年到2015年，英国报道了125起由常见物品聚光点火而引发的火灾事故．如2009年，一个由玻璃球制成的门把手汇聚阳光并点燃了屋内的可燃物，最终造成了触目惊心的火灾悲剧．2014年，在英国伦敦，一颗被放在窗边的玻璃球汇聚阳光后点燃了屋内的电视机，最终烧毁整栋房屋．不只是建筑火灾，太阳聚光也常常是导致野火的罪魁祸首．在2019年，一场森林火灾袭击了我国的广西省．事后调查发现，火灾的起因正是装水的塑料瓶汇聚的阳光点燃了森林中的可燃物．近日在上海，一对父子进行放大镜聚光实验时意外引燃了一块巨大的生态园草坪，给园区造成了巨大损失．

可燃物的点火对火灾的发生和发展起着至关重要的作用[2]．与其他点火方式不同，太阳的聚光点火主要发生在较小的加热区域，可视为点状点火源．点状点火源是现实火灾中常见的起火因素，如飞火颗 粒[3]、金属热颗粒[4]、高温熔融液滴[5]、激光[6]和闪 电[7]等．在过去几十年，大量的文献研究了飞火颗粒在森林或城林交界域的点火机理[3]，而针对其他点状点火源引发火灾的机理则较少研究．飞火颗粒点火通常包含传热传质过程[8-9]，而太阳聚光点火只包含能量的传递，该过程与激光点火相似．因此，激光点火的实验现象可以为太阳聚光点火提供一定的参考．Kashiwagi[10-11]发现，在加热直径为2～3cm的激光束下，红橡木和PMMA塑料发生明火自燃(flaming autoignition)的临界热通量分别为90kW/m2和160kW/m2．即使将激光束的加热直径增大到3.5cm，其所需的临界热通量(80～90kW/m2)仍远远高于锥形量热仪下标准测试(加热直径为10cm)所测的值(25～50kW/m2)[12]．

此外，Grishin等[13]发现了多孔森林可燃物在钨丝灯产生的光束下发生明火点火所需的临界热通量随着加热直径的增大而减小．在2014年，Warren[1]专门对常见聚光物品产生的光斑进行了火灾风险的定性研究，他将一个小型球状玻璃鱼缸装满水后放在太阳光下，随后对其建立光学模型算出大致的光斑辐照热流，并证明了太阳光聚光后点燃印刷纸的可能性．最近，美国桑迪亚国家实验室[14-15]在一个聚光太阳能塔中使用了数百个反射镜，产生了大于2000kW/m2(直径约为0.5m)的辐射热流，并在几秒钟内点燃PMMA板和松树．然而，目前学术圈对太阳聚光点火尚无系统的研究，其点火的机理仍未得到很好的理解与总结．

本文利用直径为150mm、焦距为108mm的玻璃球汇聚的太阳光，探索了太阳聚光对多层再生餐巾纸的阴燃点火机理．在玻璃球的焦距内，实验所选取的光斑直径为1.5～20mm，其辐射强度可高达780kW/m2．实验对不同光斑大小和辐射强度下的阴燃点火实验和临界热通量进行了定量分析，从而确定太阳聚光点火的临界条件．最后，本文建立了一个二维计算模型对实验结果进行验证．

1 实验设计

1.1 实验装置与燃料

实验采用未经漂白的薄纸巾作为燃料(图1(a))，因其常见于室内，并常作为垃圾被遗弃于野外．实验开始前，薄纸先在70℃的恒温烤箱内充分干燥48h，然后存储在电子恒湿恒温箱内，以避免其再次吸收空气中的水分．薄纸的干密度约为98kg/m3，厚度约为1mm．

本实验利用直径为150mm、折射率为1.53的透明玻璃球汇聚阳光，其焦距和后焦距分别为108mm和33mm．太阳光经玻璃球汇聚后会在燃料表面形成圆形光斑并持续加热燃料(图1)．为了使实验装置与太阳光平行(保证光束垂直作用于燃料表面)，装置的底座上安装了一个由两个铝块组成的定位装置(靠近太阳的铝板正中间开有一个1mm的小孔)．

1.2 光斑辐射强度的量化

图1 实验装置图及光学模拟的聚光系数与测量值的对比

为了验证模拟的结果，利用常规的热流计测量的较大光斑的辐射强度并计算出其聚光系数．如图1(b)所示，计算结果与辐射热流的测量值基本一致，验证了光学模拟计算结果的准确性．由于存在球面像差，穿过玻璃球的阳光束会形成焦散区，即视图外部逐渐变得模糊．因此，当光束超过焦距后，光斑的边缘将难以辨别．为了保证实验的精度，实验中只采用了后焦距内边界清晰的光斑．为了研究光斑大小对阴燃点火临界条件的影响，本实验选取了直径分别为20mm、9.0mm、5.5mm和1.5mm的4种光斑，其距离玻璃球表面的直线距离分别为3mm、13mm、19mm和33mm．

1.3 实验流程

实验前，通过轨道调节样品架至所需光斑大小所对应的位置．然后，用一块隔热板将光斑隔离，再将样品插入样品架内．当太阳辐射计的读数相对稳定后，移开隔热板使光斑作用在试样上一段时间．加热后，将试样移出并放置在稳定无风的环境内观察5min，以判断阴燃点火是否成功．在相同光斑下的不同辐射强度测试了超过100组实验以确定其阴燃点火的临界条件．然后，通过调节样品架的位置以研究光斑直径的大小对点火临界条件的影响．本研究在阳光充足的室外环境下进行了超过600组实验．太阳的初始辐射强度在0.2～1.6kW/m2之间．实验过程中，大气温度为(29±2)℃，相对湿度为(82±10)%，大气压强约为101kPa．

2 实验结果

2.1 阴燃点火与燃烧现象

图2(a)展示了纸张在加热直径为1.5mm、辐射热流为600kW/m2的光斑下的阴燃点火和火蔓延现象．当光斑作用在纸张上时，可观察到大量的白烟从加热区域溢出．随着加热的进行，加热区域表面的纸张开始炭化并破碎，使下层的纸张直接受到光斑的加热．在加热8s后，纸张从样品架取出，并放置在无风的环境下观察5min．如图2(a)所示，炭化区域开始向外扩散，并最终将整份纸张烧烬．图2(b)展示了纸张在加热直径为5.5mm、辐射热流为350kW/m2的光斑下的实验现象．在加热的初始阶段，可以观察到跟图2(a)相同的热解和炭化趋势．然而，在持续加热8s之后，炭化区域没有向外扩散，因此可认为点火失败．

图2 纸张在不同光斑下阴燃点火与火蔓延现象

2.2 临界点火热通量

如公式(5)所示，点火所需的临界热通量随着光斑直径的增大而减小，并逐渐趋近于一维模型下所得的最小点火辐射热通量(12kW/m2)，这与实验的现象完全相符．

图3 不同光斑直径下的阴燃点火临界热通量

3 数值模拟

3.1 二维瞬态守恒方程和阴燃动力学参数

为了更好地解释实验现象，本文基于开源代码Gpyro[22-24]建立了二维模型，以预测光斑大小对点火临界条件的影响，其数值模型如图4所示．利用试样几何结构的对称性，模型选取了真实尺寸的一半进行计算，以节约计算成本和时间．

该模型求解了二维固相和气相瞬态守恒方程，并假设了气相和固相组分之间的热平衡．此处仅列出模型中的主要守恒方程，其他的细节可以从手册中查阅[22–24]．

固相质量守恒方程：

固相组分守恒方程：

固相能量守恒方程：

气相质量守恒方程：

气相组分守恒方程：

气相动量守恒方程[25]：

图4 太阳聚光点火的二维计算模型

由于纸张主要由木浆制成，本文采用与木材阴燃相似的4步异相化学反应模型：包括半纤维素、纤维素和木质素的热解及炭的氧化，见式(12)～式(15)：

表1 薄纸阴燃4步反应动力学参数

Tab.1 4-step smoldering kinetics of paper

表2 凝聚相物种的热物理参数

Tab.2　　Thermophysical properties of condensed-phase species

3.2 模拟结果

为了更好地展示瞬态特性，图5于对称轴左、右两边分别展示了图5(a)点火成功与图5(b)点火失败两个算例的二维温度和和热释放速率分布图．如图5(a)所示，在25kW/m2的辐射热流下，当加热光斑的直径为20mm时，在加热50s后，加热区域附近可观察到一个释放出强烈热量的高温阴燃区，并随着时间的推移逐渐向外端蔓延．在500s时，该燃烧区域蔓延至试样的末端，并最终将样品烧尽．相反地，保持辐射热流的强度不变，当加热光斑的直径减小到5mm，在加热了200s后，加热区域温度升至约400℃并释放出少量热量，如图5(b)所示．然而，随着时间的推移，加热区域的温度和热释放逐渐降低，并最终熄灭．即使将加热的时间延长至1000s，点火仍不能发生．这些结果均与实验观察到的现象一致．

（a）点火成功(D＝20mm) （b）点火失败(D＝5mm)

图6总结了模型预测的不同光斑下阴燃点火临界辐射热通量与实验值的对比．总的来说，预测结果与实验数据吻合较好，进一步验证了模型的准确性．同时，该模型还较好地预测了加热光斑直径对点火临界热通量的影响；即点火临界热通量随光斑直径的增大而减小．例如，当光斑直径从1mm增大到20mm时，点火的临界辐射热通量从290kW/m2急降至13kW/m2，所有的这些预测均与实验测量一致．

图6　　不同光斑直径下预测的阴燃点火临界热通量与测量值的对比

4 结 论

(1)本文通过小尺度实验研究了多层薄纸在汇聚的阳光下阴燃点火现象．实验中，自然的阳光通过一个直径为150mm、焦距为108mm的玻璃球汇聚后作用在燃料上，并通过改变玻璃球和燃料的距离以研究光斑的大小对阴燃点火临界热通量的影响．为了量化聚光后的辐射热通量，本实验利用光学软件Tracepro模拟了不同光斑下的辐射热通量，并通过实测值验证模型的可靠性．

(2)实验发现，随着太阳辐射增大或光斑直径减小，聚光后的辐射热流均增大．随着光斑直径的增大，阴燃点火的临界热通量逐渐降低并趋近于在锥形量热仪下测到的最小值(11kW/m2)．研究通过简化的传热分析揭示了燃料内部传导散热的影响将随着光斑的变小而逐渐增强．最后，本文基于开源代码Gpyro建立二维模型再现实验现象并验证实验结果．本项工作量化了阳光经汇聚后潜在的火灾风险，并揭示了燃料内部传导散热对小光斑加热的影响．未来的工作将利用实验和模拟进一步探讨不同燃料在汇聚的太阳光下的火灾风险．

［1］ Warren D W. Optical modeling of fire hazards arising from sunlight focused by water[J].2014，50：1327-1334.

［2］ Babrauskas V.[M]. Issaquah：Society of Fire Protection Engineers，2003.

［3］ Manzello S L，Suzuki S，Gollner M J，et al. Role of firebrand combustion in large outdoor fire spread[J].2020，76：100801.

［4］ Wang S，Huang X，Chen H，et al. Ignition of low-density expandable polystyrene foam by a hot particle[J].2015，162(11)：4112-4118.

［5］ Sun P，Lin S，Huang X. Ignition of thin fuel by thermoplastic drips：An experimental study for the dripping ignition theory[J].2020，115：103006.

［6］ Chen J，Taniguchi M，Narato K，et al. Laser ignition of pulverized coals[J].1994，97(1)：107-117.

［7］ Zhang H，Qiao Y，Chen H，et al. Experimental study on flaming ignition of pine needles by simulated lightning discharge[J].2021，120：103029.

［8］ Fernandez-Pello A C. Wildland fire spot ignition by sparks and firebrands[J].2017，91：2-10.

［9］ Wang S，Huang X，Chen H，et al. Interaction between flaming and smouldering in hot-particle ignition of forest fuels and effects of moisture and wind[J].2017，26(1)：71-81.

［10］ Kashiwagi T. Effects of sample orientation on radiative ignition[J].，1982，44(1-3)：223-245.

［11］ Kashiwagi T. Experimental observation of radiative ignition mechanisms[J].1979，34：231-244.

［12］ Yang L Z，Zhou Y P，Wang Y F，et al. Autoignition of solid combustibles subjected to a uniform incident heat flux：The effect of distance from the radiation source[J].2011，158(5)：1015-1017.

［13］ Grishin A M，Golovanov A N，Medvedev V V. On the ignition of a layer of combustible forest materials by light radiation[J].1999，35：618-621.

［14］ Engerer J D，Brown A L，Christian J M. Ignition and damage thresholds of materials at extreme incident radiative heat flux[C]//12Atlanta，GA，USA，2018，AIAA 2018-3764.

［15］ Brown A L，Engerer J D，Ricks A J，et al. Datasets for material ignition from high radiant flux[J].，2021，121：103131.

［16］ Altair Engineering Inc. TracePro. Software for design and analysis of illumination and optical systems[EB/OL]. https://lambdares.com/tracepro/，2021.

［17］ Lin S，Huang X，Urban J，et al. Piloted ignition of cylindrical wildland fuels under irradiation[J].2019，5：54.

［18］ Babrauskas V.[M]. London：Springer，2016.

［19］ Incropera F P.r[M]. London：John Wiley，2007.

［20］ Holman J P.[M]. 10th Edition. New York：McGraw-Hill，2008.

［21］ Lin S，Huang X. Quenching of smoldering：Effect of wall cooling on extinction[J].2021，38(3)：5015-5022.

［22］ Lautenberger C. Gpyro-A Generalized Pyrolysis Model for Combustible Solids：Technical Reference[R]. Berkeley：University of California，Berkeley，2014.

［23］ Yuan H，Huang X，Zhao G，et al. Gpyro Workbook on Pyrolysis & Smoldering Problems[R]. London：Imperial College London，2019.

［24］ Lautenberger C. A Generalized Pyrolysis Model for Combustible Solids[R]. Berkeley：University of California，Berkeley，2007.

［25］ 黄鑫炎，林少润，刘乃安. 林火中的阴燃现象：研究前沿与展望[J]. 工程热物理学报，2021，42：512-528.

Huang Xinyan，Lin Shaorun，Liu Naian. A review of smoldering wildfire：Research advances and prospects[J].2021，42(2)：512-528(in Chinese).

［26］ Huang X，Rein G. Thermochemical conversion of biomass in smouldering combustion across scales：The roles of heterogeneous kinetics，oxygen and transport phenomena[J].2016，207：409-421.

［27］ Li K，Huang X，Fleischmann C，et al. Pyrolysis of medium-density fiberboard：Optimized search for kinetics scheme and parameters via a genetic algorithm driven by Kissinger’s method[J].2014，28(9)：6130-6139.

［28］ Huang X，Rein G. Smouldering combustion of peat in wildfires：Inverse modelling of the drying and the thermal and oxidative decomposition kinetics[J].2014，161(6)：1633-1644.

［29］ Huang X，Li K，Zhang H. Modelling bench-scale fire on engineered wood：Effects of transient flame and physicochemical properties[J].2017，36(2)：3167-3175.

Smoldering Ignition Characteristics by Concentrated Solar Irradiation：Insight from Experimental and Numerical Study

Lin Shaorun1, 2，Wang Siyan1，Sun Peiyi1，Yuan Han1，Huang Xinyan1, 2

(1. Research Centre for Fire Safety Engineering，The Hong Kong Polytechnic University，Hong Kong 999077，China；2. Shenzhen Research Institute，The Hong Kong Polytechnic University，Shenzhen 518057，China)

The concentrated solar irradiation spot can ignite common fuels，which is a potential fire risk that threatens the safety of people’s life and properties. In this work，small-scale experiments were performed to investigate the smoldering ignition of multi-layer tissue paper using concentrated solar irradiation spots. The sunlight was concentrated through a transparent glass sphere with a diameter of 150mm and a focus length of 108mm. The diameter of solar spot applied here(1.5—20.0mm)was obtained by adjusting the position of the sample within the focus length. It was found that，as the diameter of solar spot decreases，the internal conductive heat loss increases，and the critical radiant heat flux required for ignition increases. Afterwards，a 2-D numerical model was established based on the open-source code Gpyro to reproduce the experimental phenomena and validate the experimental results. This work quantifies the fire risk of concentrated solar spots and provides some guidance for future fire protection strategies.

smoldering；critical heat flux；solar spot；heating area；fire

10.11715/rskxjs.R202305024

TK11

A

1006-8740(2023)04-0390-07

2022-03-07．

国家自然科学基金资助项目(51876183).

林少润(1994—  )，男，博士，flynn.lin@connect.polyu.hk.

黄鑫炎，男，博士，副教授，xy.huang@polyu.edu.hk.

(责任编辑：隋韶颖)