消防安全工程 第6部分：烟气层的计算要求 GB/T31593.6-2015

master

前言

中华人民共和国国家标准
消防安全工程 第6部分：烟气层的计算要求
Fire Safety engineering——Part 6；Requirements governing algebraic equaUons for smoke layers
GB/T 31593.6-2015
发布日期：2015年06月02日
实施日期：2015年08月01日
发布部门：中华人民共和国国家质量监督检验检疫总局
中国国家标准化管理委员会

    GB/T 31593《消防安全工程》分为以下九个部分：
    ——第1部分：计算方法的评估、验证和确认；
    ——第2部分：所需数据类型与信息；
    ——第3部分：火灾风险评估指南；
    ——第4部分：设定火灾场景和设定火灾的选择；
    ——第5部分：火羽流的计算要求；
    ——第6部分：烟气层的计算要求；
    ——第7部分：顶棚射流的计算要求；
    ——第8部分：开口气流的计算要求；
    ——第9部分：人员疏散评估指南。
    本部分为GB/T 31593的第6部分。
    本部分按照GB/T 1.1-2009给出的规则起草。
    本部分使用重新起草法修改采用ISO 16735:2006《消防安全工程 计算公式控制要求 烟气层》。
    本部分与ISO 16735:2006相比在结构上有较多调整，附录A中列出了本部分与ISO 16735:2006的章条编号对照一览表。
    本部分与ISO 16735:2006的技术性差异及其原因如下：
    ——关于规范性引用文件，本部分做了具有技术性差异的调整，以适应我国的技术条件，调整的情况集中反映在第2章“规范性引用文件”中，具体调整如下：
    ● 用GB/T 5907（所有部分）代替了ISO 13943（见第3章）；
    ● 用修改采用国际标准的GB/T 31593.1代替了ISO/TR 13387-3（见5.1、7.3）；
    ● 删除了ISO 16735:2006当中仅被资料性附录引用的ISO 16734:2006、ISO 16737。
    ——将国际标准中1.1～1.4改为列项，以符合我国标准的编写习惯（见第1章）。
    为了方便使用和符合我国相关标准编写要求，本部分还做了下列编辑性修改：
    ——删除了国际标准的前言，重新起草了前言；
    ——修改了国际标准的引言，将其作为本部分的引言；
    ——将国际标准的“本国际标准”一词改为“本部分”；
    ——将国际标准中某些标点符号修改为符合汉语习惯的标点符号；
    ——调整了部分公式的编排格式；
    ——增加了附录A（资料性附录），给出了本部分与国际标准的章条编号对照情况。
    本部分由中华人民共和国公安部提出。
    本部分由全国消防标准化技术委员会建筑消防安全工程分技术委员会（SAC/TC 113/SC 13）归口。
    本部分起草单位：公安部天津消防研究所、公安部四川消防研究所、中国科学技术大学、中国建筑科学研究院。
    本部分主要起草人：姚松经、韩伟平、毕少颖、智会强、阚强、张玉贤、陆守香、胡忠日、张向阳、邓松华、郑巍、郭歌。

引言

    本部分所讨论的计算公式适用于设定火灾场景的量化分析，专业人员无须进行繁杂的数值计算，就能够依据这些公式快速判断初步消防安全设计是否需要调整以及如何调整，以满足预期的性能指标要求。这些公式在以下方面已经得到应用：
   ——测定火羽流传热，包括对流和辐射；
   ——预测顶棚射流特性以调节探测器的响应时间；
   ——计算通风口的排烟能力；
   ——分析分隔空间内的烟气传递和轰燃等火灾危害。
   对于烟气层，计算公式经常用来估算烟气充满一个分隔空间所需的时间以及烟气层的温度和浓度。
   针对不同类型的烟气层，本部分的附录B和附录C分别给出了符合本部分要求的计算公式示例。附录B给出了烟气层的一般信息和守恒关系，附录C给出了计算烟气层特性的具体公式及其计算示例。
   GB/T 31593.1规定了运用数学公式进行火灾动力学计算的通用要求，本部分是这些通用要求的具体实践。如果采用数学模型计算火灾发展过程及其产生的后果，本部分所讨论的计算公式可用于对其结论进行校验。
   本部分适合从事消防安全工程的专业人员使用，包括消防设计人员、消防救援人员、消防标准规范制定者、保险商、消防安全管理人员，标准使用者需要充分理解本部分给出的计算方法中相关参数的含义和应用条件。

1 范围

    GB/T 31593的本部分规定了火灾烟气层特征值计算公式的应用要求，提供了与烟气层计算公式应用相关的下列通用要求：
    a） 物理现象的描述；
    b） 计算书；
    c） 计算公式的局限性； 
    d） 计算公式的输入参数；
    e） 计算公式的适用范围。
    本部分适用于建设工程消防性能化设计和评估中烟气层的计算。

2 规范性引用文件

    下列文件对于本文件的应用是必不可少的。凡是注日期的引用文件，仅注日期的版本适用于本文件。凡是不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。
    GB/T 5907（所有部分） 消防词汇
    GB/T 31593.1 消防安全工程 第1部分：计算方法的评估、验证和确认（GB/T 31593.1-2015，ISO 16730:2008，MOD）
    GB/T 31593.5 消防安全工程 第5部分：火羽流的计算要求（GB/T 31593.5-2015，ISO 16734:2006，MOD）
    GB/T 31593.8 消防安全工程 第8部分：开口气流的计算要求（GB/T 31593.8-2015，ISO 16737:2006，MOD）

3 术语和定义

    GB/T 5907界定的以及下列术语和定义适用于本文件。
3.1 边界面 boundary
    确定室内空间范围的表面。
3.2 火羽流 fire plume
    由燃烧所产生的浮力形成的向上湍流流动，通常包括下部的燃烧区域。
3.3 火焰 flame
    火羽流（3.2）的发光区域。
3.4 界面位置 interface position
    烟气层高度
    烟气层（3.7）相对于参考高度的高度。
    注：参考高度通常为室内空间的最低边界面（3.1）高度。
3.5 准稳态 quasi-steady state
    假定火源热释放速率变化引起的所有影响在流场中能立刻反映出来的状态。
3.6 烟气 smoke
    烟
    物质高温分解或燃烧时产生的固体和液体微粒、气体，连同夹带和混入的部分空气形成的气流。
3.7 烟气层 smoke layer
    热上层 hot upper layer
    热气层 hot-gas layer由火灾引发，在室内空间的最高边界面（3.1）下面形成并聚集，相对均匀的一定量的烟气（3.6）。
3.8  烟气层界面 smoke layer interface 
    烟气层（3.7）与其下面无烟层的水平界面。
3.9 开口 vent 
    在室内空间边界面（3.1）上，空气和烟气（3.6）通过自然或机械方式而流通的孔洞。
3.10 开口气流 vent flow
    通过室内空间边界面（3.1）开口（3.9）的烟气（3.6）或空气流。

4 物理现象的描述

4.1 在室内空间中，火源生成浮动烟气层是一种复杂的热物理现象。它可能是高瞬变状态，也可能是接近稳态。在有焰燃烧区域和未发生燃烧的区域，都可能出现烟气层。除了受浮力影响，烟气层还会受到通风排烟的影响。 
4.2 应采用图表描述火源类型、室内空间流动边界面条件和其他场景因素。
4.3 应明确识别需要计算的烟气层特征参数及其适用范围，适当时应包括由相关计算量推导出来的特征参数（例如通过质量和能量守恒推导出的烟气浓度与升高的气体温度之间的关系）以及其他与烟气层对物体和人员的热作用相关的特征参数。
4.4 应明确识别具体计算公式适用的物理现象（如单纯的蓄烟，机械排烟等）。
4.5 不同的计算公式描述不同的烟气层特性（见4.3）或适用于不同的火灾场景（见4.4），当有多种方法可用于计算同一个给定量的值时，应明示其结果与选用的计算方法无关。

5 计算书

5.1 计算书的一般要求见GB/T 31593.1。
5.2 计算步骤应由一系列的计算公式表述。
5.3 每个计算公式应由独立的条款表述，其内容应包含公式输出的详细描述，以及该公式解释性说明和限定条件。
5.4 应明确定义计算公式中的各变量，给出适用的SI单位；计算公式应优先选用量关系式。
5.5 应酌情通过引用公认手册、科技文献或通过推导等方法给出公式的科学依据。
5.6 应给出公式的应用实例，演示如何使用符合第4章要求的输入参数和具体计算过程。

6 计算公式的局限性

6.1 应给出直接应用计算公式计算输出参数的定量限制条件，并符合第4章所描述场景的要求。
6.2 应给出在更通用的计算方法中引用计算公式的注意事项，包括检查与计算方法中用到的其他关系式的一致性以及所采用的计算方案。

7 计算公式的输入参数

7.1 应明确识别计算公式的输入参数，如热释放速率或几何尺寸。
7.2 应识别或明确提供输入参数的数据来源。
7.3 应按GB/T 31593.1的规定列出输入参数的有效范围。

8 计算公式的适用范围

8.1 应通过一组或多组测量数据确定公式的适用范围。这些数据应通过文件化的程序或标准（参见GB/T 6379）进行评估并达到一定的质量水平（如重复性、再现性）。
8.2 应按照计算方法的评估、验证和确认原则，通过与8.1测量数据的比对确定计算公式的适用范围。
8.3 应明确识别用于第4章特定场景的计算公式限定条件造成的潜在误差。

附录A （资料性附录）本部分与ISO 16735：2006的章条编号对照情况

    表A.1给出了本部分章条编号与ISO 16735:2006的章条编号的对照一览表。
 

附录B （资料性附录）烟气层计算概述

B.1 物理现象的描述
    B.1.1 计算方法概述
    本附录旨在描述室内空间中火灾的界面位置、平均温度和边界面以下烟气层的特定化学物平均浓度的计算方法。这些计算方法基于质量、物质和能量守恒定律，适用于与热力学参考空间相同的烟气层。
    燃烧时，烟会在室内空间的上部积聚。当烟气层的温度和物质浓度比较均匀时，根据烟气层的质量、物质和能量守恒定律，可计算温度平均值、烟浓度和界面位置。火羽流和开口气流的详细计算要求见GB/T 31593.5—2015和GB/丁31593.8。
B.1.2 所需计算的烟气层特征值
    所需计算的烟气层特征值为平均烟气层温度、物质浓度和界面位置。
B.2 计算书
    B.2.1 概述
    在室内空间的火源上方产生的烟气层如图B.1所示。B.3.2～B.2.4规定了质量、能量和特定的化学物守恒。

B.2.2 质量守恒
    烟气层中质量守恒应考虑选择如图B.1虚线所示的适当参考空间。通过各界面流入该参考空间（流出该空间的为负值）的质量流速应等于该空间烟气层的质量累积速率。必要时应考虑火羽流、开口气流和其他流动。
B.2.3 能量守恒
    烟气层的能量守恒采用与质量守恒相同的方法。通过烟气层界面流入（流出的为负值）的能量流速应等于该空间烟气层的能量累积速率。除了火羽流和开口气流外，还应考虑室内空间边界面发生的辐射损失和热量吸收。
    注：当无法确定火焰的辐射热损失时，质量流速可近似用附录B提供的热释放速率代替。
B.2.4 特定化学物的守恒 
    特定化学物的质量守恒采用与总的质量守恒相同的方法。如果烟气层中发生气相化学反应，应适当考虑反应速率。
B.2.5 通过界面的火羽流的质量流速
    在界面（烟气层的底面）上，火羽流的质量流速应作为火的热释放速率和火源底面与烟气层界面间的垂直距离的函数给出。GB/T 31593.5—2015给出了具有这种火羽流特性的公式示例。
B.2.6 通过开口的烟的质量流速
    通过开口的质量流速应作为烟气层的温度、烟气层与相邻空间的压差、开口宽度和开口高度的函数给出。GB/T 31593.8给出了具有这种开口特性的公式示例。
B.2.7 状态方程
    状态方程表达了烟的温度和密度的关系。烟近似于理想气体，其性质等同于空气。 

附录C （资料性附录）满足附录B要求的烟气层计算公式

C.1 本附录所用的符号说明
C.2 物理现象的描述
C.3 计算书
C.4 计算公式的依据
C.5 计算公式的局限性
C.6 公式的输出参数
C.7 计算公式的输入参数
C.8 计算公式的适用范围

C.1 本附录所用的符号说明

    本附录所用的符号说明见表C.1。

 

  

C.2 物理现象的描述

C.2.1 概述
    本附录的这些计算方法适用于计算平均温度、烟浓度和不同火灾场景下的界面位置。如果能够确认和验证，也可用适用的其他方法计算这些量。
C.2.2 公式计算法适用的火灾场景
    这些公式适用于静态环境下火源上方的烟气层。如果与火无关的湍流比较明显，则这些公式不适用。例如，如果采暖通风与空调系统或外部通风产生的空气流比较明显，则应注意其影响。如果洒水喷头等动作的灭火系统对烟气层具有明显的影响，则这些公式不适用。
    火源必须足够小，以保证主要火焰高度低于界面位置和火羽流的特征宽度比室内空间窄（具有火羽流特征的公式所需的附加条件）。
    烟气层条件的计算方法有两个限定阶段。一个是在火初发期（t2型火），控烟设备尚未启动时，单一的室内空间蓄烟阶段；另一个是烟的产出率等于烟气层流出率的准稳态排放阶段。本附录不涉及中间阶段（即排烟系统启动后，继续充烟的阶段）。
C.2.3 采用公式计算的烟气层特性
    采用公式可计算气体温度、物质浓度和界面位置。
C.2.4 公式适用的烟气层条件
    给定的公式适用于没有排烟以及在机械或自然排烟达到准稳态的情况下，室内空间中蓄烟的短暂阶段。
C.2.5 公式的协调一致性
    公式具有自身一致性。
C.2.6 采用这些公式的标准和其他文件
    无。

C.3 计算书

C.3.1 公式的范围
    本附录提供了四组公式。一组用于火的初发期单一的室内空间中蓄烟过程；其余三组用于稳态下通过机械排风或自然开口进行的烟气控制。
C.3.2 室内空间的蓄烟过程
C.3.2.1 公式在过程中的应用
    如图C.1所示，烟在室内空间的上部累积，直到烟气层界面下移到垂直开口的上边缘为止。由于热膨胀，过量的空气被挤出室内空间。
    注：这个前提在烟气层的底部位于开口上界面的以上位置时是成立的。当烟气层降低到开口上界面以下位置时，随着新鲜空气的进入，烟气会流出室内空间。

    热释放速率见式（C.1）：

     当n＝0时，表示稳定燃烧的火；n＝2时，表示正在增大的火与时间成平方关系。当辐射释放的份数为X，则对流热释放速率见式（C.2）：

     参考文献[4]中给出了火源上方高度为z处的火羽流质量流速公式，见式（C.3）：

   注：在GB/T 31593.5—2015的附录B中，该公式可近似用作火羽流公式。该公式只适用于平均火焰高度以上的情况；如果界面低于平均火焰高度，则计算结果可能不正确。
C.3.2.2 界面位置
    计算界面的位置[见式（C.4）和式（C.5）]，以便在具有均匀密度的上层累积火羽流的质量流量。

     注：为了计算界面位置，必须设定烟密度。在实际应用中，对大体积室内空间初期的蓄烟过程取ρ_s＝1.0会得到保守的结果（见参考文献[5]）。在稍后的蓄烟过程中，具有明显的热膨胀。在这种情况下，由参考文献[6]和[4]给出的下列公式适用于t²型即＝atⁿ）。

    其中：

    烟气层温度T_s由式（C.9）计算得出。
C.3.2.3 烟气层温度
    计算烟气层的温度[见式（C.8）），使得火释放的热量对体积为A（H－z）的烟气层加热。此时，可忽略室内空间界面的热量吸收。

    注1：符号λ表示室内空间界面吸收的热量的份数。除非火羽流、烟气层和室内空间之间的辐射热交换相互叠加，否则，推荐设定λ＝0，即所有的热量均用来加热烟气层。
    注2：在实际应用中，对大体积室内空间初期的蓄烟过程取ρ_s＝1.0会得到可接受的结果。
    注3：在稍后的蓄烟过程中，当烟气层的热膨胀比较明显时，t²型火的烟气层温度由式（C.9）计算：

    其中Λ和X由式（C.6）和式（C.7）计算得到。
C.3.2.4 特定的化学物的浓度
    计算特定的化学物的浓度[见式（C.10）]，使得生成物的质量在烟气层中均匀分布。

C.3.2.5 计算举例
    ——火源＝O.05t2（α＝0.05kW/s²，n＝2，D＝1mm）位于图C.1所示的室内空间内；
    ——室内空间的底面面积A为100m²，
    ——室内空间高度H为8m，门的开口高度H_u为2m；
    ——假设辐射热所占份数X为0.333；
    ——忽略室内空间界面的热吸收（λ＝0）；
    ——CO₂的产率η为7.61×10^-5kg/kJ；
    ——计算界面高度、温度和60s时CO₂的浓度。
    采用式（C.4），界面高度为：

    将结果带入式（C.8）和式（C.10），则界面温度和CO₂的浓度为：

    火焰高度必须低于界面的高度，火羽流的式（C.3）才适用。此例中，火焰的平均高度低于界面高度，这样：

    该计算的依据为GB/T 31593.5-2015的附录B。
    同样，烟气层高度、温度和CO₂浓度计算示例见图C.2。为使公式有效，烟气层的底面必须高于火焰的平均高度和门的开口上边缘高度。该例在126s时平均火焰高度和烟气层高度近似相等，此时，

    在142s，烟气层高度与门的开口上边缘高度近似相等，此时：

    因此，该公式适用于126s之前的阶段。
    图C.2中的计算有效范围为：A＝100m²，H＝8m，＝0.05t²，X＝0.333，λ＝0.0。粗线是根据式（C.4）、式（C.8）和式（C.10）的计算结果。细线是考虑烟气层的热膨胀而根据式（C.5）、式（C.9）和式（C.11）得出的计算结果。

C.3.3 通过机械排气系统进行烟气控制的稳态
C.3.3.1 公式适用的过程
    如图C.3所示，在烟气控制状态下，烟通过机械排气系统排放。烟气层的性质通过产烟和排烟达到平衡的准稳态进行计算。假设室内空间边界面在较低的部分具有足够的开口能使空气容易流动并且公式中的热释放速率为常数，式（C.3）给出了火羽流的质量流速。把体积流速作为试验设计参数，排烟的质量流速用式（C.18）计算：

 

    计算界面高度，使排烟的质量流速等于火羽流的质量流速。

C.3.3.2 界面位置
    计算界面位置，使火羽流的质量流速等于排烟的质量流速。

    注：为了计算界面位置，需要知道烟气层的密度（也就是温度），可以通过适当推测获得，或者结合下面提供的公式进行计算。
C.3.3.3 烟气层密度
    烟气层密度由下式计算：

    注：对大多数工程计算，烟气层近似为理想气体。烟气层的温度由下条的公式计算。
C.3.3.4 烟气层温度
    计算烟气层温度，使烟气层的热流量等于排放和室内空间表面吸收的热量损失之和。

C.3.3.5 有效传热系数
    有效传热系数的计算取决于室内空间边界面的建筑材料。传热既可近似为厚壁热行为（近似为半无限体）也可近似为薄壁热行为（近似薄材料稳态下的温度分布）。

    注；特征时间t_c通常取1000s。
C.3.3.6 特定化学物的浓度
    计算特定化学物的浓度，使其产出速率等于排放速率。

C.3.3.7 计算示例
    火源位于图C.3所示的室内空间中心，相关计算参数为：
    ——室内空间底面面积A为100m²（10m×10m）；
    ——室内空间高度H为8m，火源的热释放速率为300 kW；
    ——火源的热辐射份数x为0.333；
    ——火源直径D为1.0m；
    ——机械排气的体积流速为4m³/s；
    ——室内空间边界面由100mm厚的混凝土板制成；
    ——混凝土的导热系数取k＝0.0015 kW/（m·K），热容量c_v＝2026kJ/（m³·K）；
    ——空气的参考温度T₀为20℃（293 K），相对应的参考密度ρ₀为1.205 kg/m³。
    界面位置与温度的公式相关联，采用迭代法求解。计算出界面位置与温度后，可以直接计算物质的浓度。
    1） 假设界面高度是室内空间高度的50％，即：

    2） 由式（C.3）计算界面高度的火羽流质量流速：
 

    3）由式（C.23）计算有效热传导系数：
    假设室内空间边界面具有厚壁热行为，即：

    因此，有效热传导系数为：

    4） 由式（C.22）计算烟气层温度：

 

    5） 由式（C.21）计算烟密度：

    6） 由式（C.18）计算机械排气系统的质量流速：

    7） 由式（C.20）校正界面高度，以便火羽流的质量流速等于排放的质量流速：

    8） 重复步骤2）～7），直到火羽流的质量流速与排放的质量流速一致。此例中，经过3次迭代，得出以下满意结果：

    9） 为了满足火羽流式（C.3），平均火焰高度必须小于界面高度。此例中，依据GB/T 31592.5-2015的附录B进行计算，结果符合要求。

    10） 将已知值带入式（C.24）计算物质浓度。在通风良好的条件下，二氧化碳的产率η＝7.61×10^-5>kg/kJ。

C.3.4 水平开口稳态下的烟气控制
C.3.4.1 公式适用的过程
    烟气由图C.4所示的自然开口排放。假设新鲜空气可流入室内空间底部，通过准稳态下热量和质量平衡来计算烟气层特性。质量流速的平衡由式（C.39）给出：

    在这个公式中，假设热释放速率始终为常量。火羽流的质量流速在式（C.3）中给出，通过开口的质量流速按照GB/T 31593.8计算。

C.3.4.2 界面位置
    采用式（C.20）计算界面位置。
    注：排放的质量流速由式（C.40）计算。
C.3.4.3 烟气层密度
    采用式（C.21）计算烟气层密度。
C.3.4.4 烟气层温度
    采用式（C.22）计算烟气层温度。
C.3.4.5 有效传热系数
    采用式（C.23）计算有效传热系数。
C.3.4.6 通过水平开口排烟的质量流速
    排烟的质量流速由开口流速的常规公式计算。

C.3.4.7 地板面的压差
    地板面的压差由适用于较低开口流速的常规公式计算。

C.3.4.8 化学物的浓度
    采用式（C.24）计算化学物的浓度Y。
C.3.4.9 计算举例
    火源位于图C.4所示的室内空间内，相关计算参数为：
    ——室内空间地板面积A为100m²（10m×10m）；室内空间高度H为8m；
    ——水平开口的面积A_vent为2m²；
    ——较低的进气开口面积A_open为4m²；
    ——火源的热释放速率为300kW；
    ——火源的热辐射份数x为0.333；
    ——CO₂的产率为η为7.51×10^-5kg/kJ；
    ——燃料直径D为1.0m；
    ——参考温度T₀为20℃（293 K）。
    室内空间界面由等同于示例C.3.3.7的混凝土制成。
    界面位置与温度的公式相关联，采用迭代法求解。计算出界面位置与温度后，可以直接计算物质的浓度。
    1） 假设界面高度是室内空间高度的50％，即：

    2） 由式（C.3）计算界面高度的火羽流质量流速：

    3） 采用和C.3.3.7的3）相同的步骤由式（C.23）计算有效热传导系数：

    4） 由式（C.22）计算烟气层温度：

    5） 由式（C.21）计算烟密度：

    6） 由式（C.41）计算参考高度的压差：

    7） 由式（C.40）计算通过水平开口的质量流速：

    8） 由式（C.20）校正界面高度，以使火羽流的质量流速等于排放的质量流速：

    注：为了数值的稳定性，在迭代计算中，式（C.20）中的排放质量流速m。由代替。收敛后，依然保持质量平衡。这一变更对量终结果没有影响。
    9） 重复步骤2）～8），直到火羽流的质量流速与排放的质量流速一致，即。此例中，经过4次迭代，得出以下满意结果：

    10） 为了满足火羽流式（C.3），火焰高度必须小于界面高度。此例中，火焰的平均高度为1.28m，与C.3.3.7的9）相同。
    11） 将上述7）中的结果带入式（C.24）计算物质浓度：

C.3.5 垂直开口稳态下的烟气控制
C.3.5.1 公式适用的过程 
    在烟气控制状态，烟气通过图C.5所示的垂直开口排放。假设新鲜空气从开口的下部流入，而烟气从开口的上部流出，烟气层的特性通过准稳态下烟/热的产生和排放/流出速率的平衡进行计算。在该公式中，假设热释放速率始终为常量。式（C.3）给出了火羽流的质量流速，通过开口的质量流速按照GB/T 31593.8计算。

C.3.5.2 界面位置
    采用式（C.20）计算界面位置。
    注：该公式中，界面位置的计算是隐含的，它要满足稳态的质量平衡

，详见本条下面的示例。
C.3.5.3 烟气层密度
    采用式（C.21）计算烟气层密度。
C.3.5.4 烟气层温度
    采用式（C.22）计算烟气层温度。
C.3.5.5 有效传热系数
    采用式（C.23）计算有效传热系数。
C.3.5，6 开口下部水平方向上的压差
    计算开口下部水平方向上的压差，使物质排出速率等于火羽流的质量流速。

C.3.5.7 通过开口下部的空气质量流速
    将式（C.55）计算出的压差带入，计算通过开口下部的空气质量流速。

C.3.5.8 特定化学物浓度
    采用式（C.24）计算特定化学物浓度。
C.3.5.9 计算步骤和示例
    火源位于图C.5所示的室内空间内，该空间具有垂直开口，相关计算参数为：
    ——室内空间地板面积A为100m²。；
    ——垂直开口的高度为5m，在地板上方1m处（H₁＝1m，H_u＝6m）；
    ——开口宽度B为4m；
    ——火源的热释放速率为300 kW；
    ——火源的热辐射份数X为0.333；
    ——CO₂的产率η为7.51×10^-5kg/kJ；
    ——火源直径D为1.0m。
    室内空间界面由等同于示例C.3.3.7的混凝土制成。
    界面位置与温度的公式相关联，采用迭代法求解。计算出界面位置与温度后，可以直接计算物质的浓度。
    1） 假设界面高度的下限z₁是整个开口高度的1/3：

    2） 由式（C.3）计算界面高度的火羽流质量流速：

    3） 采用和C.3.3.7的3）相同的步骤由式（C.23）计算有效热传导系数：

    4） 由式（C.22）计算烟气层温度：

    5） 由式（C.21）计算烟密度：

    6） 由式（C.55）计算开口下部水平方向上的压差：

    7） 由式（C.56）计算通过开口下部的空气质量流速：
 

    8） 计算质量流速的误差：

    9） 火羽流的质量流速应等于空气的质量流速。为了得到正确的结果，将界面高度取上限z₂重复步骤2）～8）的计算。初算时，z₂取开口高度的2/3：

    采用上述界面高度，计算出以下结果：

    这个假设的界面高度过高，致使压差为负值。这种情况下，没有空气进入室内空间：

    注：对质量达到平衡的稳态，这种情况是不存在的.在迭代的过程中，当压差为负值时，取零值，以满足质量平衡。
    质量流速的误差为：

    10） 在界面高度两次取值之间插入以下估算新值：

    注：计算步骤如图C.6的虚线所示。
    11） 重复步骤2）～8），计算z₃值下质量流速的误差：

    12） 由于为负值，将式（C.71）中的z₁值换成z₃，重复步骤2）～11），直到足够小。
    在迭代过程中，如果为正值，则用z₃代替z₂进行计算。此例中，共进行了7次迭代才得出如下所示的收数值：

    注：如图C.6所示，最终结果为（z）和（z）的交差点开圆。
    13） 为了满足火羽流式（C.3），火焰高度必须小于界面高度。此例中，火焰的平均高度为1.28m，与C.3.3.7的9）相同。
    14） 将步骤12）中计算结果带入式（C.24），计算物质浓度：

 

C.4 计算公式的依据

    公式基于附录B给出的通用的热和质量守恒定律。蓄烟过程的研究可追溯到Turner等人（参见参考文献[5]）对具有明显密度分界形成的室内空间所做的流体动力学基础研究。针对火灾初期蓄烟过程，Zukoski（参见参考文献[6]、[7]）提出了类似的具体公式。在Zukoski的理论研究之后，Mulholland 等人（参见参考文献[8]）进行了实验研究，验证了小的密度变化下的假设。Tanaka等人（参见参考文献[2]）的实验工作发现如果将小的密度变化带入计算，则测量数据具有很好的复现性。[Delichatsios（参见参考文献[3]、[4]）]

C.5 计算公式的局限性

C.5.1 火羽流
    本部分的质量流速公式与GB/T 31593.5-2015中有关火源、火焰尺寸、边界、气流扰动等的限定条件相同。
C.5.2 烟气层的均匀性
    公式的前提是烟气层具有均匀的特性。相对于平均值，如果烟气层特性具有重大变化，例如竖长、杆状的室内空间和长廊，则不推荐使用这些公式。

C.6 公式的输出参数

    公式的输出参数为界面位置、烟气层温度和物质浓度。另外，公式还给出了有关火羽流的质量流速和开口流量信息。

C.7 计算公式的输入参数

C.7.1 火的热释放速率
    参数是特定环境条件下火的实际热释放速率。该参数可通过热量计测量，即测量收集到的气体产物而确定氧气、二氧化碳和一氧化碳的生成速率；也可采用其他指定的方法测量。该参数通常可由设定火灾场景中获得。

C.7.2 热释放的辐射份数
    参数x取决于燃料的种类，典型的取值范围为0.3～0.4。更多的信息可参见GB/T 31593.5-2015的B.6.2。
C.7.3 室内空间边界面初期吸收热量份数
    参数λ是计算蓄烟过程必需的。对于稳态下烟气控制（机械通风、屋顶或侧面开口）的计算，可采用不同的方法计算室内空间的热量吸收。参数的值主要取决于室内空间边界面的建筑。对具有大的热惯性（如混凝土结构）或轻型、非隔热的建筑物（如玻璃房）形成的室内空间，热量吸收可能是重要的。为了计算精确值，需要详细的室内空间的热辐射传导资料。在实际应用中，推荐取λ＝0。
C.7.4 机械排烟的效果
    在C.3.3中，假设机械排烟系统只排出烟。然而，如果烟气层不十分厚，较低的空气层将被排烟系统排出（见参考文献[12]），这种情况下，应根据夹带空气的比率而减小排放速率的值（参见参考文献[13]）。
C.7.5 物质的产率
    参数η取决于燃料的种类和燃料/空气的比率。由于公式中假设火相对于室内空间比较小，η可取通风良好的值（参见参考文献[14]）。

C.8 计算公式的适用范围

    日本建筑研究所的Tanaka和Yamana见参考文献[2]在大尺寸的中庭对公式已经进行了一系列的比对实验。实验房间的地板面积为720㎡，室内净高为26.3m，热释放速率约为130 kW。另外，Karlsson等人（参见参考文献[15]）与Hägglund等人（参见参考文献[16]）对蓄烟过程的公式进行了比对，房间的尺寸为5.62m×5.62m，净高6.15m（有效高度为5.95m），采用煤油油盘火源，热释放速率为186 kW。两个比对均证明公式具有适当的准确度。对多房间的烟扩散，可用一组相似的公式求解（参见参考文献[17]、[18]、[19]）。

参考文献

    [1] GB/T 6379（所有部分），测量方法与结果的准确度
    [2] ISO/TR 13387-3:1999 Fire safety engineering-Part 3: Assessment and verification of mathematical fire models
    [3] ISO 13943:2008 Fire safety-Vocabulary
    [4] ZUKOSKI,E.E.,Mass Flux in Fire Plumes,Fire Safety Science-Proc.Fourth Int.Symposium （ed.T.Kashiwagi）,International Association for Fire Safety Science,pp.137-147,1994.
    [5] TANAKA,T.and YAMANA,T.,Smoke Control in Large Scale Spaces （Part 1/Part 2）,Fire Science and Technology,5,No.1,Science University of Tokyo,1985.
    [6] DELICHATSIOS,M.A.,Closed form approximate solutions for smoke filling in enclosures including the volume expansion term,Fire Safety Journal,38,pp.97-101,2003.
    [7] DELICHATSIOS,M.A.,Tenability Conditions and Filling Times for Fires in Large Spaces,Fire Safety Journal,39,pp.643-662,2004.
    [8] BAINES,W.D.and TURNER,J.S.,Turbulent Buoyant Convection from a Source in a Confined Region,Journal of Fluid Mechanics,37（1）,pp.51-80,1969.
    [9] ZUKOSKI,E.E.,Development of a Stratified Ceiling Layer in the Early Stages of a Closed room Fire,Fire and materials,2（2）,pp.54-62,1978.
    [10] ZUKOSKI,E.E.,Two-Layer Modeling of Smoke Movement in Building Fires,Fire and materials,4（1）,pp.17-27,1980.
    [11] MULHOLLAND,G.,HANDA,T.,SUGAWA,O.and YAMAMOTO,H.,Smoke Filling in an Enclosure,20th Joint ASME/AIChE National Heat Transfer Conference,81-HT-B,1981.
    [12] SPRATT,D.and HESELDEN,A.J.M.,Efficient Extraction of Smoke from a Thin Layer under a Ceiling,Fire Research Note,No.1001,Building Research Establishment,1974.
    [13] NII,D.,NITTA,K.,HARADA,K.and YAMAGUCHI,J.,Air Entrainment into Mechanical Smoke Vent on Ceiling,Proceedings of the 7th International Symposium on Fire Safety Science,（Worcester,June 2002）,pp.729-740,2003.
    [14] TEWARSON,A.,Generation of Heat and Chemical Compounds in Fires,SFPE Handbook of Fire Protection Engineering,2nd ed.,National Fire Protection Association,1995.
    [15] KARLSSON,B.and QUINTIERE,J.,G.,Conservation Equations and Smoke Filling,Chapter 8 in Enclosure Fire Dynamics,CRC Press LLC,2000.
    [16] HÄGGLUND B.,JANSSON R.and NIREUS K.,Smoke Filling Experiments in 6x6x6 meter Enclosure,FOA Report C20585-D6,National Defence Research Establishment,Sweden,1985.
    [17] TANAKA,T.,A Model of Multiroom Fire Spread,NBSIR 83-2718,US Department of Commerce,National Bureau of Standards,1983.
    [18] TANAKA,T.,A Model on Fire Spread in Small Scale Buildings,BRI Research Paper,No.79 and 84 ,Building Research Institute,1978,1980.
    [19] PEACOCK,R.,JONES,W.,W.,FORNEY,G.,P.,RENEKE,P.,A.,BUKOWSKI,R.and W.,KLOTE J.,H.,An Update Guide for HAZARD I version 1.2,NISTIR 5410 ,Building and Fire Research Laboratory,National Institute of Standards and Technology,1994.