3.1气体泄漏扩散的模拟方法
目前在研究气体扩散领域应用较多的模拟方法主要有三种，即：物理模拟方法、 数学模拟方法和CFD 数值模拟方法。当然在实际的模拟仿真过程中，经常是两种或是三种方法同时使用，以此来验证模拟的准确性。
3.1.1物理模拟方法
物理模拟是模拟的基础方法，[31]指在不同与实体的规模上将某一过程再现，并 分析其物理特性和线性尺度对实体的影响，进而对所研究实体或过程进行直接实 验。将实际地形物理按比例的缩小模型置于实验体(如风洞、水槽等)内，在满足基本相似条件(主要包括几何、运动、热力、动力和边界条件相似)的基础上，模拟真实过程的主要特征，如空气动力规律和扩散规律。
物理模型建立的理论基础是相似理论。进行进行物理模拟研究，必须解决如何 设计和制作模型以及将模型实验的结论在实体上应用等问题。相似原理是研究、支配力学相似系统的性质及如何用模型实验解决实际问题的一门科学，是进行模型实验研究的依据。
根据相似理论，物理模型若能与原型保持相似，则由物理模型经过实验得到的 规律，原型也同样适用。建立物理模型要遵循很多相似条件，如几何相似、运动相似、动力相似及热相似等。在建立模型时，由于所有相似条件不可能完全满足，所以针对研究的具体要求，要适当做出取舍，恰当选取相似参数是实现物理模拟的关键。物理模拟主要用于数值计算模式难于处理的复杂地形以及受到建筑物影响时的扩散研究。与现场实验相比，特别是复杂条件下的现场试验相比，物理模拟实验条件易控制、可重复，且可节省人力、物力，可进行较全面和规律性实验，是大气扩散研究的重要手段。
3.1.2数学模拟方法
数学模拟方法是解决简单扩散问题的常用方法，此方法是[31]通过用数学模型、 在一定条件下来研究一个物理或化学过程，或通过模型描述一个复杂的物理或化学过程的某些特点。此种方法所借助的数学模型的方式没有固定限制，可以是一系列代数式或微分、积分方程，也可以简化为一个关系式。
其中常见的数学模型：高斯模型、箱及相似模型、浅层模型、Sutton 模型以及唯象模型。
3.1.3CFD 数值模拟方法
CFD 模拟是一种数值模拟方法，用此方法解决流体运动问题于数学方法的研究思路不同，此种方法对扩散的研究不必依赖偏微分方程的求解，在解决问题时，根据具体研究的要求，不是去求解析解，而是运用有限元的思想对具体问题建模，并通过相应的软件技术对模型进行模拟仿真计算，使对具体的流动过程的分析和研模拟。
采用这种数值模拟方法进行模拟有一定的程序。第一步，根据泄漏介质的特点和泄漏条件建立基本守恒方程，包括质量方程、动量方程、能量方程以及组分方程等；第二步，判断和选择初始和边界条件，对扩散中的各种场函数进行模拟，这些场主要有流场、温度场、浓度场等。第三步，对各种描述结果进行分析，完成模拟目的。[40]此种方法是在借助计算机的基础上完成的，模拟过程中不需要对
气体的密度性质进行判断，对于扩散的初始状态也不用进行详细的阶段分析，因此用这种方法进行模拟就克服了箱及相似模型中的许多缺点，对于箱及相似模型模拟是遇到的问题都可以不考虑，如辨识和模拟重气的下沉、扩散过程中空气的卷吸、等许多问题。另外，该方法对模拟对象的要求不高，[41]既可以模拟平坦均匀地形情况的扩散问题，又有模拟更复杂情形流动的能力。这种模拟方法是在Navier-Stokes 方程(以下简称N-S方程)的基础上产生的，是一种完全三维的流体力学模拟方法，至少在原理上，此种方法可以模拟所有重要的物理过程。 物体的流动有层流和湍流之分。通常情况下，大气的流动都为湍流。[42]相对于 层流来说，湍流过程更为复杂。经过长时间的研究和总结工作，目前关于湍流的工程模式和计算机数值模拟方法有许多，但总结来看，已经普遍采用只有直接数值模拟方法、大涡模拟方法和湍流统观模拟方法这三种。
3.2 Fluent 软件介绍
尽管天然气在室内的扩散满足微观质量守恒定律、动量守恒定律和能量守恒定 律，但是由于受到各种因素内外界因素的影响，[28]天然气在运动过程中由于多种原因会发生很大的形变，这种形变的方式采用数学方法很难进行求解。即使利用数学方法得到了偏微分方程，并通过近似等方式得到了方程的解析解，但这种方式不能很好的反应天然气的扩散状态。所以对具体的扩散问题只有借助计算机技术来完成，运用数值求解方法来对扩散问题进行研究，计算流体力学(CFD)是进行此类研究的基础。
Fluent 是CFD 软件中应用很广泛的一种，它[29]在美国的市场占有率已经超过 60%，是目前国际上比较流行的CFD 软件包，可以模拟跟流体、热传递及化学反应等有关的很多问题。此软件使用方便、功能强大。其主要特点是：物理模型全面、数值方法先进、及的前、后处理功能强大，目前在很多领域都有着广泛的应用，并且取得了很好的效果。特别是航空航天研究、汽车设计、石油天然气泄漏扩散等方面。
Fluent 软件设计是基于CFD 软件群的思想，针对各种复杂流动的现象，从用户需求角度出发，陈对模拟对象的特点，采用不同的离散格式，利用不同的数值求解方法，最终达到模拟实验的目的，并要求实验中各种参数在自身领域内的稳定性。同时也要对计算速度和精度进行设置，这样可以高效的结局复杂的流动问题，而且可以解决各个领域的问题。基于上述思想，Fluent 开发了各种模拟软件，这些软件适用于各个领域、多种流动的模拟。这些软件能够模拟流体流动、化学反应、物质传热和其它多种复杂的物理现象。为了方便用户，Fluent 在建立网格时采用了统一的生成技术，而且应用相同的图形界面，使用起来很方便。究变得很方便。本章就应用CFD 软件Fluent 对天然气室内扩散过程进行仿真
3.2.1 Fluent 用于气体泄漏扩散的适用性
目前在泄漏扩散模拟仿真方面已经应用Fluent 做过的研究有化学试验烟羽轨迹跟踪，冷却塔烟羽扩散模式研究等。另外，袋虑室、旋风除尘器、有毒气体排放、烟气脱硫设备、SO2 洗涤、工业污染物排放等过程和设备的模拟等都可以应用Fluent进行模拟。通过以上的成功经验证明，Fluent 对于居室内天然气的泄露扩散问题的模拟具有适用性。
3.3天然气泄漏的仿真过程
本小节是在分析天然气室内泄漏扩散特点的基础上，确定基本控制方程，利用计 算流体力学(CFD)方法，对室内燃气泄漏扩散过程的状态进行模拟和仿真。
3.3.1天然气泄漏扩散模型
1）物理模型的建立
本文选择普通居民居住的房间作为研究对象。几何尺寸:长8m ,宽3m .在通风的条件下，天然气在管道口发生泄漏
在此条件下做以下假设：
a.把空气质点的平均运动看作是不可压缩流体的运动；
b.在数值计算中忽略空气粘性的影响；
c.本文使用湍流粘性系数来表征大气的湍流性质，并采用双方程湍流模型，即k- 湍流模型求解；
d.假设泄漏时温度不发生变化，并忽略可燃性气体与空气的温度差异。
2）几何模型的建立
Fluent 建立几何模型主要有两种方式：[49]一是在前处理器Gambit 中直接建立；另一种是利用专业的软件工具建立模型并将其导入Fluent 中。Gambit 只适用于创建简单的几何体，对于复杂形体而言，其绘图功能是远远不够的，这时Gambit 允许引入一些其他软件创建的文件，常用的有Autocad 创建的ASCI 形式的文件。本文的模型是在前处理器Gambit 中直接建立。为了更直观、更形象的反映天然气泄漏扩散的过程，本文采用的是二维建模的方式。
3）网格的划分
在使用CFD 软件的工作中，有大约80%的时间是花费在划分网格上，可以说 网格划分能力的高低是决定工作效率的主要因素之一。
[50]Fluent 软件采用非结构与适应性网格相结合的方式进行网格划分。与结构化网格和分块网格相比，非结构网格划分便于处理复杂外形的网格划分，而适应性网格则便于计算流场参数变化剧烈、梯度很大的流动，同时这种划分方式也便于网格的细化或粗化，使得网格划分更加灵活、简便。
Fluent 划分网格的途径有两种：一种是用Fluent 提供的专用网格划分软件
Gambit 进行网格划分，另一种则是由其他的制图软件，如CAD 软件等完成造型工作，再导入Gambit 中生成网格。还可以用其他网格生成软件生成与Fluent 兼容的网格用于计算。
本文是在Fluent 前处理器Gambit 中进行网格划分，按照由线到面、由面到体 的原则，采用自动划分的方式对网格进行划分，共划分网格12996个，网格划分情况如图3-3 所示。
3.3.2基本控制方程
Fluent 的求解是基于各种方程的，[51]本文涉及到的基本方程主要有连续性方程、动量方程、能量方程和组分方程四种。
1）连续性方程
对于湍流的求解要借助连续性方程来完成。根据物质扩散的各种性质不同，可 将连续性方程可分为以下[52]4 类，第一类为不可压缩流体模型(incompressible fluidmodel,IFM)；第二类为非弹性流体模型(anelastic fluid model,AFM)；第三类为变密度流体模型(variabledensity fluid model,VDFM)；第四类为变密度湍流流体模型(variabledensity turbulent fluidmodel,VDTFM)。
不可压缩流体
∂ui=0 (3.1) ∂xi
非弹性流体
∂ρui=0 (3.2) ∂xi
变密度流体
∂ρ∂(ρui)+=0 (3.3) ∂t∂xi
变密度湍流流体
∂ρ∂(ρui)∂(ρ'ui')++=0 (3.4) ∂t∂xi∂xi
对于本课题研究的气体射流问题，考虑到天然气泄漏为射流喷口，且室内泄漏 过程中忽略温度的变化影响，且其密度与空气的不同，在扩散过程中其密度也在不断的发生变化。根据这些特点，天然气室内泄漏过程的基本控制方程为如式3.3 所示，式中，ρ 为泄漏物质的密度，ui为三个方向的速度(u，v，w )。
2）动量方程
气体扩散同样遵循动量守恒方程，可得出扩散过程中x、y、z 三个方向上的动 量守恒方程，通式为：
∂uj∂(ρui)∂∂ui∂p∂∂+(ρuiuj)=-+(μt)+(μt)+(ρ-ρa)gi(3.5) ∂t∂xj∂xi∂xj∂xj∂xj∂ui
∂(ρu)∂∂∂+(ρu2)+(ρuv)+(ρuw)=∂t∂t∂y∂z
x 方向：-∂ρ∂∂u∂∂u∂∂u+(μt)+(μt)+(μt)+ (3.6) ∂x∂x∂x∂y∂y∂z∂z
∂∂u∂∂v∂∂w(μt)+(μt)+(μt)∂x∂x∂y∂x∂t∂x
∂(ρv)∂∂∂+(ρuv)+(ρv2)+(ρuv)=∂t∂t∂y∂z
y 方向：-∂p∂∂v∂∂v∂∂v+(μt)+(μt)+(μt)+ (3.7) ∂y∂x∂x∂y∂y∂z∂z
∂∂u∂∂v∂∂w(μt)+(μt)+(μt)∂x∂y∂y∂y∂z∂y
∂(ρw)∂∂∂+(ρuw)+(ρvw)+(ρw2)=∂t∂t∂y∂z
∂p∂∂w∂∂w∂∂wz 方向：-+(μt)+(μt)+(μt)+ (3.8) ∂z∂x∂x∂y∂y∂z∂z
∂∂u∂∂v∂∂w(μt)+(μt)+(μt)∂x∂z∂y∂z∂z∂z
式中，u 、v 、w 分别为流体扩散过程中在x 、y 、z 三个方向上的速度分量， 为流体的湍流粘度， g 为重力加速度， p 为扩散开始的绝对压力。
3）能量方程
扩散过程中的能量方程如下所示：
cp-cpaμt∂w∂T∂(ρT)∂1∂∂T (3.9) +(ρujT)=(kt)+v[()]∂t∂xjcp∂xj∂xjcpσc∂xj∂xj
式中，T 为流体的温度，kt为流体的湍流导热系数，可通过流体的湍流普朗特 (Prandtl)数将σt与μt联系起来：ktPrμ=()t。σT取为常数，一般为 0.9~1， kσTμ
Pr=μcp/k。其中：cp、cpv分别为混合物流体的定压比热和为泄漏物质的定压比
热，cpa为空气的定压比热， cp=ωcpv+(1-ω)cpa。
4）组分方程
由组分质量守恒定律可得出组分方程：
∂(ρω)∂∂∂ω+(ρujω)=(ρDt) ∂t∂xj∂xj∂xj
式中，ω 为各成分的组分质量分率，Dt为流体扩散的湍流扩散系数，此系数可通过
流体的湍流施密特（Schmidtl）数将 σc 与 ut联系起来：DtSμ=(c)t。σc取Dσcμ
为常数，一般可取为1，而Sc=
μ∂(ρω)∂∂∂ω于是有：+(ρujω)=(μt)。
ρD∂t∂xj∂xj∂xj
3.3.3Fluent 数值计算
将网格划分后，进进入Fluent 计算阶段，其基本计算步骤是：(1)定义流畅的几何参数；(2)启动相关的求解器；(3)将划分好的网格输入到Fluent 中；(4)对网格进行检查；(5)选择求解器格式；(6)选择求解器所用的基本方程；(7)定义物质属性；(8)定义边界条件；(9)调整解的控制参数；(10)初始化流场；(11)开始求解。以下就对其中的几个关键步骤进行介绍。
a.确定求解器
Fluent 提供三种计算方式：分离方式、耦合隐式和耦合显式。这三种计算方式 都可以给出精确的计算结果，只是针对某些特殊问题时，某种计算方式可能比其他两种方式更快一些。分离计算和耦合计算的区别在于求解连续、动量、能量和组元方程的方法有所不同。分离方式是分别求解上面的几个方程，最后得到全部方程的解；耦合方式则是用求解方程组的方式，同时进行计算并最后获得方程的解。两种方式的共同特点是，在求解附带的标量方程时，都是采用单独求解的方式，就是先求解控制方程，再求解湍流模型方程或是辐射方程。显式和隐式的区别在于对方程的线化方式有所不同。分离方式一般用于不可压缩流或弱可压缩流的计算。耦合方式则通常用于高速可压缩流的计算。由本文第一章中的假设条件可知，本文假设天然气及空气都为不可压缩流，所以本文选择分离式求解器进行求解。
b.确定湍流模型
泄漏物质扩散过程中浓度变化主要受泄漏物质自身性质和气流场两方面因素
的影响。而天然气在室内的泄漏是一个复杂的过程。通过物理模拟和数值模拟的方法都不能方便准确的得到天然气扩散过程中浓度的变化，所以本研究采用CFD 数值模拟方法。天然气泄漏后在居室内的运动并不稳定，一般是以湍流流动为表现形式。湍流本身很复杂，所以系统还要遵守附加的湍流输运方程，模拟要选择好相应的湍流模型以考虑湍动量的影响因素。
常用的湍流模型有涡粘性模式和Reynolds 应力模式，[53]涡粘模型中常用的有 零方程模型，一方程模型和两方程k −ε 、k −ω 模型；Reynolds 应力模式中常用的有SST、SSG 模型等。Sklavounos 等比较了k −ε 、k −ω 、SST、SSG4种模型在重气输运数值模拟的结果，发现k −ε 及 SST 两种模型的效果较好，本文计算采用标准k −ε模型(standard k −ε model)。
c.设置边界条件
设置边界条件是模拟效果好坏的关键因素之一。在Fluent 中所谓边界条件，是 指在模拟过程中，在模型求解域的边界上所求解的变量和参数。只有给定了合理准确的边界条件，才可能准确计算出流场中各种参数的解。因此，只有根据模型的要求以及模拟的目的准确的确定了边界条件，要计算的问题才会有定解，模拟才会达到想要的效果。
泄漏口的初始速度3m/s ,压力直径0.02m ，甲烷的质量分数 1
空气入口的速度0.3m/s ,压力直径0.2m ，氧气的质量分数 0.23 出口压力直径0.2m
d.设置初始化条件
在开始迭代前，必须先进行初始化流场。所谓初始化，就是给各流场参数赋初 始值，以便迭代计算可以有一个起点。初始化流场的方式有两种：一种是初始化流场边界；一种是直接初始化流场变量。
本文假设从泄漏口开始计算，迭代的时间步长为1s，天然气的初始浓度为0，泄漏的初始速度为3m/s。
3.3.4 Fluent的仿真结果
50S
100s
150s
200s
3.1气体泄漏扩散的模拟方法
目前在研究气体扩散领域应用较多的模拟方法主要有三种，即：物理模拟方法、 数学模拟方法和CFD 数值模拟方法。当然在实际的模拟仿真过程中，经常是两种或是三种方法同时使用，以此来验证模拟的准确性。
3.1.1物理模拟方法
物理模拟是模拟的基础方法，[31]指在不同与实体的规模上将某一过程再现，并 分析其物理特性和线性尺度对实体的影响，进而对所研究实体或过程进行直接实 验。将实际地形物理按比例的缩小模型置于实验体(如风洞、水槽等)内，在满足基本相似条件(主要包括几何、运动、热力、动力和边界条件相似)的基础上，模拟真实过程的主要特征，如空气动力规律和扩散规律。
物理模型建立的理论基础是相似理论。进行进行物理模拟研究，必须解决如何 设计和制作模型以及将模型实验的结论在实体上应用等问题。相似原理是研究、支配力学相似系统的性质及如何用模型实验解决实际问题的一门科学，是进行模型实验研究的依据。
根据相似理论，物理模型若能与原型保持相似，则由物理模型经过实验得到的 规律，原型也同样适用。建立物理模型要遵循很多相似条件，如几何相似、运动相似、动力相似及热相似等。在建立模型时，由于所有相似条件不可能完全满足，所以针对研究的具体要求，要适当做出取舍，恰当选取相似参数是实现物理模拟的关键。物理模拟主要用于数值计算模式难于处理的复杂地形以及受到建筑物影响时的扩散研究。与现场实验相比，特别是复杂条件下的现场试验相比，物理模拟实验条件易控制、可重复，且可节省人力、物力，可进行较全面和规律性实验，是大气扩散研究的重要手段。
3.1.2数学模拟方法
数学模拟方法是解决简单扩散问题的常用方法，此方法是[31]通过用数学模型、 在一定条件下来研究一个物理或化学过程，或通过模型描述一个复杂的物理或化学过程的某些特点。此种方法所借助的数学模型的方式没有固定限制，可以是一系列代数式或微分、积分方程，也可以简化为一个关系式。
其中常见的数学模型：高斯模型、箱及相似模型、浅层模型、Sutton 模型以及唯象模型。
3.1.3CFD 数值模拟方法
CFD 模拟是一种数值模拟方法，用此方法解决流体运动问题于数学方法的研究思路不同，此种方法对扩散的研究不必依赖偏微分方程的求解，在解决问题时，根据具体研究的要求，不是去求解析解，而是运用有限元的思想对具体问题建模，并通过相应的软件技术对模型进行模拟仿真计算，使对具体的流动过程的分析和研模拟。
采用这种数值模拟方法进行模拟有一定的程序。第一步，根据泄漏介质的特点和泄漏条件建立基本守恒方程，包括质量方程、动量方程、能量方程以及组分方程等；第二步，判断和选择初始和边界条件，对扩散中的各种场函数进行模拟，这些场主要有流场、温度场、浓度场等。第三步，对各种描述结果进行分析，完成模拟目的。[40]此种方法是在借助计算机的基础上完成的，模拟过程中不需要对
气体的密度性质进行判断，对于扩散的初始状态也不用进行详细的阶段分析，因此用这种方法进行模拟就克服了箱及相似模型中的许多缺点，对于箱及相似模型模拟是遇到的问题都可以不考虑，如辨识和模拟重气的下沉、扩散过程中空气的卷吸、等许多问题。另外，该方法对模拟对象的要求不高，[41]既可以模拟平坦均匀地形情况的扩散问题，又有模拟更复杂情形流动的能力。这种模拟方法是在Navier-Stokes 方程(以下简称N-S方程)的基础上产生的，是一种完全三维的流体力学模拟方法，至少在原理上，此种方法可以模拟所有重要的物理过程。 物体的流动有层流和湍流之分。通常情况下，大气的流动都为湍流。[42]相对于 层流来说，湍流过程更为复杂。经过长时间的研究和总结工作，目前关于湍流的工程模式和计算机数值模拟方法有许多，但总结来看，已经普遍采用只有直接数值模拟方法、大涡模拟方法和湍流统观模拟方法这三种。
3.2 Fluent 软件介绍
尽管天然气在室内的扩散满足微观质量守恒定律、动量守恒定律和能量守恒定 律，但是由于受到各种因素内外界因素的影响，[28]天然气在运动过程中由于多种原因会发生很大的形变，这种形变的方式采用数学方法很难进行求解。即使利用数学方法得到了偏微分方程，并通过近似等方式得到了方程的解析解，但这种方式不能很好的反应天然气的扩散状态。所以对具体的扩散问题只有借助计算机技术来完成，运用数值求解方法来对扩散问题进行研究，计算流体力学(CFD)是进行此类研究的基础。
Fluent 是CFD 软件中应用很广泛的一种，它[29]在美国的市场占有率已经超过 60%，是目前国际上比较流行的CFD 软件包，可以模拟跟流体、热传递及化学反应等有关的很多问题。此软件使用方便、功能强大。其主要特点是：物理模型全面、数值方法先进、及的前、后处理功能强大，目前在很多领域都有着广泛的应用，并且取得了很好的效果。特别是航空航天研究、汽车设计、石油天然气泄漏扩散等方面。
Fluent 软件设计是基于CFD 软件群的思想，针对各种复杂流动的现象，从用户需求角度出发，陈对模拟对象的特点，采用不同的离散格式，利用不同的数值求解方法，最终达到模拟实验的目的，并要求实验中各种参数在自身领域内的稳定性。同时也要对计算速度和精度进行设置，这样可以高效的结局复杂的流动问题，而且可以解决各个领域的问题。基于上述思想，Fluent 开发了各种模拟软件，这些软件适用于各个领域、多种流动的模拟。这些软件能够模拟流体流动、化学反应、物质传热和其它多种复杂的物理现象。为了方便用户，Fluent 在建立网格时采用了统一的生成技术，而且应用相同的图形界面，使用起来很方便。究变得很方便。本章就应用CFD 软件Fluent 对天然气室内扩散过程进行仿真
3.2.1 Fluent 用于气体泄漏扩散的适用性
目前在泄漏扩散模拟仿真方面已经应用Fluent 做过的研究有化学试验烟羽轨迹跟踪，冷却塔烟羽扩散模式研究等。另外，袋虑室、旋风除尘器、有毒气体排放、烟气脱硫设备、SO2 洗涤、工业污染物排放等过程和设备的模拟等都可以应用Fluent进行模拟。通过以上的成功经验证明，Fluent 对于居室内天然气的泄露扩散问题的模拟具有适用性。
3.3天然气泄漏的仿真过程
本小节是在分析天然气室内泄漏扩散特点的基础上，确定基本控制方程，利用计 算流体力学(CFD)方法，对室内燃气泄漏扩散过程的状态进行模拟和仿真。
3.3.1天然气泄漏扩散模型
1）物理模型的建立
本文选择普通居民居住的房间作为研究对象。几何尺寸:长8m ,宽3m .在通风的条件下，天然气在管道口发生泄漏
在此条件下做以下假设：
a.把空气质点的平均运动看作是不可压缩流体的运动；
b.在数值计算中忽略空气粘性的影响；
c.本文使用湍流粘性系数来表征大气的湍流性质，并采用双方程湍流模型，即k- 湍流模型求解；
d.假设泄漏时温度不发生变化，并忽略可燃性气体与空气的温度差异。
2）几何模型的建立
Fluent 建立几何模型主要有两种方式：[49]一是在前处理器Gambit 中直接建立；另一种是利用专业的软件工具建立模型并将其导入Fluent 中。Gambit 只适用于创建简单的几何体，对于复杂形体而言，其绘图功能是远远不够的，这时Gambit 允许引入一些其他软件创建的文件，常用的有Autocad 创建的ASCI 形式的文件。本文的模型是在前处理器Gambit 中直接建立。为了更直观、更形象的反映天然气泄漏扩散的过程，本文采用的是二维建模的方式。
3）网格的划分
在使用CFD 软件的工作中，有大约80%的时间是花费在划分网格上，可以说 网格划分能力的高低是决定工作效率的主要因素之一。
[50]Fluent 软件采用非结构与适应性网格相结合的方式进行网格划分。与结构化网格和分块网格相比，非结构网格划分便于处理复杂外形的网格划分，而适应性网格则便于计算流场参数变化剧烈、梯度很大的流动，同时这种划分方式也便于网格的细化或粗化，使得网格划分更加灵活、简便。
Fluent 划分网格的途径有两种：一种是用Fluent 提供的专用网格划分软件
Gambit 进行网格划分，另一种则是由其他的制图软件，如CAD 软件等完成造型工作，再导入Gambit 中生成网格。还可以用其他网格生成软件生成与Fluent 兼容的网格用于计算。
本文是在Fluent 前处理器Gambit 中进行网格划分，按照由线到面、由面到体 的原则，采用自动划分的方式对网格进行划分，共划分网格12996个，网格划分情况如图3-3 所示。
3.3.2基本控制方程
Fluent 的求解是基于各种方程的，[51]本文涉及到的基本方程主要有连续性方程、动量方程、能量方程和组分方程四种。
1）连续性方程
对于湍流的求解要借助连续性方程来完成。根据物质扩散的各种性质不同，可 将连续性方程可分为以下[52]4 类，第一类为不可压缩流体模型(incompressible fluidmodel,IFM)；第二类为非弹性流体模型(anelastic fluid model,AFM)；第三类为变密度流体模型(variabledensity fluid model,VDFM)；第四类为变密度湍流流体模型(variabledensity turbulent fluidmodel,VDTFM)。
不可压缩流体
∂ui=0 (3.1) ∂xi
非弹性流体
∂ρui=0 (3.2) ∂xi
变密度流体
∂ρ∂(ρui)+=0 (3.3) ∂t∂xi
变密度湍流流体
∂ρ∂(ρui)∂(ρ'ui')++=0 (3.4) ∂t∂xi∂xi
对于本课题研究的气体射流问题，考虑到天然气泄漏为射流喷口，且室内泄漏 过程中忽略温度的变化影响，且其密度与空气的不同，在扩散过程中其密度也在不断的发生变化。根据这些特点，天然气室内泄漏过程的基本控制方程为如式3.3 所示，式中，ρ 为泄漏物质的密度，ui为三个方向的速度(u，v，w )。
2）动量方程
气体扩散同样遵循动量守恒方程，可得出扩散过程中x、y、z 三个方向上的动 量守恒方程，通式为：
∂uj∂(ρui)∂∂ui∂p∂∂+(ρuiuj)=-+(μt)+(μt)+(ρ-ρa)gi(3.5) ∂t∂xj∂xi∂xj∂xj∂xj∂ui
∂(ρu)∂∂∂+(ρu2)+(ρuv)+(ρuw)=∂t∂t∂y∂z
x 方向：-∂ρ∂∂u∂∂u∂∂u+(μt)+(μt)+(μt)+ (3.6) ∂x∂x∂x∂y∂y∂z∂z
∂∂u∂∂v∂∂w(μt)+(μt)+(μt)∂x∂x∂y∂x∂t∂x
∂(ρv)∂∂∂+(ρuv)+(ρv2)+(ρuv)=∂t∂t∂y∂z
y 方向：-∂p∂∂v∂∂v∂∂v+(μt)+(μt)+(μt)+ (3.7) ∂y∂x∂x∂y∂y∂z∂z
∂∂u∂∂v∂∂w(μt)+(μt)+(μt)∂x∂y∂y∂y∂z∂y
∂(ρw)∂∂∂+(ρuw)+(ρvw)+(ρw2)=∂t∂t∂y∂z
∂p∂∂w∂∂w∂∂wz 方向：-+(μt)+(μt)+(μt)+ (3.8) ∂z∂x∂x∂y∂y∂z∂z
∂∂u∂∂v∂∂w(μt)+(μt)+(μt)∂x∂z∂y∂z∂z∂z
式中，u 、v 、w 分别为流体扩散过程中在x 、y 、z 三个方向上的速度分量， 为流体的湍流粘度， g 为重力加速度， p 为扩散开始的绝对压力。
3）能量方程
扩散过程中的能量方程如下所示：
cp-cpaμt∂w∂T∂(ρT)∂1∂∂T (3.9) +(ρujT)=(kt)+v[()]∂t∂xjcp∂xj∂xjcpσc∂xj∂xj
式中，T 为流体的温度，kt为流体的湍流导热系数，可通过流体的湍流普朗特 (Prandtl)数将σt与μt联系起来：ktPrμ=()t。σT取为常数，一般为 0.9~1， kσTμ
Pr=μcp/k。其中：cp、cpv分别为混合物流体的定压比热和为泄漏物质的定压比
热，cpa为空气的定压比热， cp=ωcpv+(1-ω)cpa。
4）组分方程
由组分质量守恒定律可得出组分方程：
∂(ρω)∂∂∂ω+(ρujω)=(ρDt) ∂t∂xj∂xj∂xj
式中，ω 为各成分的组分质量分率，Dt为流体扩散的湍流扩散系数，此系数可通过
流体的湍流施密特（Schmidtl）数将 σc 与 ut联系起来：DtSμ=(c)t。σc取Dσcμ
为常数，一般可取为1，而Sc=
μ∂(ρω)∂∂∂ω于是有：+(ρujω)=(μt)。
ρD∂t∂xj∂xj∂xj
3.3.3Fluent 数值计算
将网格划分后，进进入Fluent 计算阶段，其基本计算步骤是：(1)定义流畅的几何参数；(2)启动相关的求解器；(3)将划分好的网格输入到Fluent 中；(4)对网格进行检查；(5)选择求解器格式；(6)选择求解器所用的基本方程；(7)定义物质属性；(8)定义边界条件；(9)调整解的控制参数；(10)初始化流场；(11)开始求解。以下就对其中的几个关键步骤进行介绍。
a.确定求解器
Fluent 提供三种计算方式：分离方式、耦合隐式和耦合显式。这三种计算方式 都可以给出精确的计算结果，只是针对某些特殊问题时，某种计算方式可能比其他两种方式更快一些。分离计算和耦合计算的区别在于求解连续、动量、能量和组元方程的方法有所不同。分离方式是分别求解上面的几个方程，最后得到全部方程的解；耦合方式则是用求解方程组的方式，同时进行计算并最后获得方程的解。两种方式的共同特点是，在求解附带的标量方程时，都是采用单独求解的方式，就是先求解控制方程，再求解湍流模型方程或是辐射方程。显式和隐式的区别在于对方程的线化方式有所不同。分离方式一般用于不可压缩流或弱可压缩流的计算。耦合方式则通常用于高速可压缩流的计算。由本文第一章中的假设条件可知，本文假设天然气及空气都为不可压缩流，所以本文选择分离式求解器进行求解。
b.确定湍流模型
泄漏物质扩散过程中浓度变化主要受泄漏物质自身性质和气流场两方面因素
的影响。而天然气在室内的泄漏是一个复杂的过程。通过物理模拟和数值模拟的方法都不能方便准确的得到天然气扩散过程中浓度的变化，所以本研究采用CFD 数值模拟方法。天然气泄漏后在居室内的运动并不稳定，一般是以湍流流动为表现形式。湍流本身很复杂，所以系统还要遵守附加的湍流输运方程，模拟要选择好相应的湍流模型以考虑湍动量的影响因素。
常用的湍流模型有涡粘性模式和Reynolds 应力模式，[53]涡粘模型中常用的有 零方程模型，一方程模型和两方程k −ε 、k −ω 模型；Reynolds 应力模式中常用的有SST、SSG 模型等。Sklavounos 等比较了k −ε 、k −ω 、SST、SSG4种模型在重气输运数值模拟的结果，发现k −ε 及 SST 两种模型的效果较好，本文计算采用标准k −ε模型(standard k −ε model)。
c.设置边界条件
设置边界条件是模拟效果好坏的关键因素之一。在Fluent 中所谓边界条件，是 指在模拟过程中，在模型求解域的边界上所求解的变量和参数。只有给定了合理准确的边界条件，才可能准确计算出流场中各种参数的解。因此，只有根据模型的要求以及模拟的目的准确的确定了边界条件，要计算的问题才会有定解，模拟才会达到想要的效果。
泄漏口的初始速度3m/s ,压力直径0.02m ，甲烷的质量分数 1
空气入口的速度0.3m/s ,压力直径0.2m ，氧气的质量分数 0.23 出口压力直径0.2m
d.设置初始化条件
在开始迭代前，必须先进行初始化流场。所谓初始化，就是给各流场参数赋初 始值，以便迭代计算可以有一个起点。初始化流场的方式有两种：一种是初始化流场边界；一种是直接初始化流场变量。
本文假设从泄漏口开始计算，迭代的时间步长为1s，天然气的初始浓度为0，泄漏的初始速度为3m/s。
3.3.4 Fluent的仿真结果
50S
100s
150s
200s