火灾过程动力学特性分析2011-08-04.doc

第6章 多层多室建筑火灾烟气蔓延过程计算模拟 6.1 计算模拟方法介绍 本章采用就算流体力学、计算燃烧学与传热传质学的方法，对五层模型楼火灾烟气扩散的影响进行数值模拟研究，给出火灾条件下不同楼层内的烟气温度、有毒气体浓度、可见度、速度场等特征参数的分布情况。 计算采用模拟软件 FDS version 3 进行数值模拟，FDS （Fire Dynamics Sinulator 火灾动力学模拟）是一种火灾驱动流体流动的计算流体动力学软件，是有美国国家标准技术局开发（李元洲，2001：史聪灵，2005）。其默认湍流模型采用Smagorinsky形式 大涡模型，燃烧模型采用的是混合分数模型，辐射传热通过对非散射灰体近似的气体利用有限体积方法求解其辐射传输方程，大约为100个离散角。这一软件可以通过对几何参数和网格数目以及燃烧过程的一些设置，如燃烧反应模型、火源功率与燃烧时间、边界热交换属性、通通风口大小与流量等：然后模拟计算可得到的一些输出数据，如某点、某条线或某个面上的温度、密度、压力及混合组分在活在过程中随时间的变化等，然后再利用一些数据处理工具就可以获得研究所需的数据以及曲线等。FDS原理是火灾场模拟计算，场模拟是利用计算机求解火灾过程中状态参数的空间及其随时间变化的模拟方式，场是指状态参数如加速度、温度、各组分的浓度等的空间分布。场模拟的理论依据是自然界普遍成立的质量守恒、动量守恒、能量守恒以及化学反应的定律等。火灾过程中状态参数的变化也遵循着这些规律，因而可以用场模拟方法求解火灾过程。FDS通过大涡模型对连续方程、动量、能量方程以及压力收敛方程进行求解，可得到温度、压力、气体成分、可见度等参数的空间分布（刘方，2002）。 6.2 计算流动、燃烧、传热模型 6.2.1 基本控制方程求解过程 计算求解的控制方程（史聪灵等，2003）如下。 质量守恒方程： （6.1） 式中，为密度;t为时;为速度矢量。 组分守恒方程： （6.2） 式中，下标l为第l种组分;为第l种组分的浓度; 为单位空间内第l种组分的质量产生的速率；为第l种组分的扩散系数。 动量方程： （6.3） 式中，为外部施加的矢量；为黏性力张量；P为压力。 能量方程： （6.4） 式中，为焓；为热辐射通量；T为温度。 状态方程： (6.5) 式中，R为气体常数，M为混合气体分子质量，下标为第种组分。 通过耦合连续方程、能量方程、状态方程，可以得到速度散度方程： （6.6） 式中，为定压比热容；为第种组分的定压比热容；为参考压力。 进一步对散度方程进行简化： （6.7） 式中，为单位空间内能量产生的速率。 对整个区域进行积分，得到参考压力方程： (6.8) 式中，为整个区域，V为体积。 定义，经过推导，动量方程最终求解形式为 （6.9） 6.2.2 湍流流动模型 湍流模型采用Smagorinsky形式的大涡模型LES。动量方程中的粘性应力张量为 (6.10) 式中，I为辐射强度；应变张量为 （6.11） 式中，、、为速度。 根据Smagorinsky的分析，湍流黏性系数可以表示为 （6.12） 式中，为经验系数，为网格特征尺度。并且，应变张量与弥散方程有关系： （6.13） 在LES模拟中，热扩散和物质扩散均与湍流黏性系数的关系为 ， (6.14) 式中，为普朗特殊数，为施密特数。 6.2.3 燃烧模型 燃烧模型 采用的是混合分数模型，定义混合分数： ，