二、要明确的几个问题

1 、外部辐射系数（ External Emissivity ） VS 内部辐射系数（Internal Emissivity）。FLUENT 中采用这样两个相似的名字有它一定的道理，它们都是用来计算辐射的时候要在总辐射能量的前面用到的一个系数。但同时这样的命名也给理解造成了一定麻烦，很容易混淆。要是从物理概念上理解这两个参数就不会弄混了。外部辐射考虑当在流场外有一个辐射源向流场辐射热量的时候而用到的参数，也就是说只有你选择辐射壁或者混合壁的时候这个参数才出现，要根据流场外的辐射源来确定这个参数。内部辐射系数，是在你考虑辐射模型的条件下才出现，例如在你选择P1、DO 等计算热辐射的时候，所以这是一个根据壁面固体材料特性选择确定的参数。

特别要注意的是，不同的固体材料差别很大，具体应用的时候查材料特性手册得到。我要强调的是大家是做流体计算的，好多人都忽略了固体材料的事情，用默认的万能的铝，什么都不改，有时这是很成问题的！ 2、壁面热产生率一定要和壁面厚度配合使用。

3 、分清壁面条件中给定的都是什么温度，用到的温度有三个： Temperature 、Free Stream Temperature 和External Radiation

Temperature，特别要注意，只有第一个是给定了壁面的温度，后面两个分别给出自由流的参考温度和外部辐射源的参考温度。

___________________________________________ 湍流的数值模拟

目前采用的数值计算方法可以大致分为以下三大类：

2.1直接模拟

直接模拟就是用三维的非稳态的纳维—斯托克斯方程对湍流进行直接数值计算的方法。要对高度复杂的湍流运动进行直接的数值计算，必须采用很小的时间与空间步长，才能分辨出湍流的详细的空间结构以及变化剧烈的特性。因此，湍流的直接数值模拟对计算机内存空间和计算速度的要求非常高，目前还无法用于工程数值计算。只有少数使用超级计算机的研究者才能从事这一类研究和计算。

2.2大涡模拟

按照湍流的涡旋学说，湍流的脉动与混合主要由大尺度的涡旋造成。大尺度的涡从主流中获得能量，他们是高度的非各向同性，而且随流动的情形而异。大尺度的涡通过互相作用把能量传递给小尺度的涡。小尺度涡的主要作用是耗散能量，它们几乎是各向同性的，而且不同的流动中的小尺度的涡有许多共性。关于

17

涡旋的上述认识就导致了大尺度涡模拟的数值解法。这种方法旨在用非稳态的纳维-斯托克斯方程来直接模拟大尺度涡，但不直接计算小尺度涡，小涡对大涡的影响通过近似的模型来考虑，这种大涡模拟对计算机内存以及计算速度的要求虽然仍比较高，但远低于直接模拟的方法对计算机资源的要求，在工作站上甚至个人电脑上都可以进行一定的研究工作，因而近年来的研究与应用日趋广泛。

2.3应用Reynolds时均方程的模拟方法

在这类方法中将非稳态的控制方程对时间做平均，在所得的关于时均物理量的控制方程中包含了脉动量乘积的时均值等物理量，于是方程的个数将小于未知量的个数，方程组不封闭。要使方程组封闭，必须做出建设，即建立模型。 在Reynolds时均方程法中，又有Reynolds应力方程法及湍流粘性系数法两大类。其中湍流粘性系数法是目前工程流动与数值计算中应用最广泛的方法。

________________________________________________

多相流模型和离散相模型的区别

两相流：通常把含有大量固体或液体颗粒的气体或液体流动称为两相流；其中含有多种尺寸组颗粒群为一个“相”，气体或液体为另一“相”，由此就有气—液，气—固，液—固等两相流之分。

两相流的研究：对两相流的研究有两种不同的观点：一是把流体作为连续介质，而把颗粒群作为离散体系；而另一是除了把流体作为连续介质外，还把颗粒群当作拟连续介质或拟流体。

引入两种坐标系：即拉格朗日坐标和欧拉坐标，以变形前的初始坐标为自变量称为拉格朗日Langrangian 坐标或物质坐标；以变形后瞬时坐标为自变量称为欧拉Eulerian 坐标或空间坐标。

1.离散相模型

??FLUENT在求解连续相的输运方程的同时，在拉格朗日坐标下模拟流场中离散相的第二相；

??离散相模型解决的问题：煤粉燃烧、颗粒分离、喷雾干燥、液体燃料的燃烧等；

??应用范围：FLUENT中的离散相模型假定第二相体积分数一般说来要小于10-12%（但颗粒质量承载率可以大于10-12%，即可模拟离散相质量流率等/大于连续相的流动）；不适用于模拟在连续相中无限期悬浮的颗粒流问题，包括：搅拌釜、流化床等；

??颗粒-颗粒之间的相互作用、颗粒体积分数对连续相的影响未考虑；

18

??湍流中颗粒处理的两种模型：Stochastic Tracking，应用随机方法来考虑瞬时湍流速度对颗粒轨道的影响；Cloud Tracking，运用统计方法来跟踪颗粒围绕某一平均轨道的湍流扩散。通过计算颗粒的系统平均运动方程得到颗粒的某个“平均轨道”

FLUENT 提供五种雾化模型：

??（1）平口喷嘴雾化(plain-orifice atomizer) ??（2）压力－旋流雾化(pressure-swirl atomizer) ??（3）靶式雾化(flat-fan atomizer)

??（4）气体辅助雾化(air-blast/air-assisted atomizer) ??（5）气泡雾化(effervescent/flashing atomizer)

用户可以在Set Injection Properties 面板中选择喷嘴类型及其相应参数下面就介绍各个喷嘴模型：

概述：所有的模型都是用喷嘴的物理及尺寸参数（例如喷口直径、质量流率）来计算初始颗粒尺寸、速度、位置。

对于实际的喷嘴模拟来说，无论是颗粒的喷射角度还是其喷出时间都是随机分布的。但对FLUENT 的非雾化喷射入口来说，液滴都是在初始时刻以一个固定的轨道喷射出去（到流场中去）。喷雾模型中使用随机选择模型得到液滴的随机分布。

随机选择轨道表明初始液滴的喷射方向是随机的。所有的喷嘴模型中都要设第初始喷射角（范围），颗粒通过随机的方法在这个范围内得到一个初始喷射方向。这种方法提高了由喷射占主导地位流动的计算精度。在喷嘴附近，液滴在计算网格内的分布趋向于更加均匀，这样，通过气相作用于液滴上的曳力就加强了气相－液滴之间的耦合作用。

1. 平口喷嘴雾化(plain-orifice atomizer)模型

平口喷嘴是最常见也是最简单的一种雾化器。但对于其内部与外部的流动机制却很复杂。液体在喷嘴内部得到加速，然后喷出，形成液滴。这个看似简单的过程实际却及其复杂。平口喷嘴可分为三个不同的工作区：单相区、空穴区、以

19

及回流区。不同工作区的转变是个突然的过程，并且产生截然不同的喷雾状态。喷嘴内部区域决定了流体在喷嘴处的速度、初始颗粒尺寸、以及液滴分散角。 2. 压力－旋流雾化喷嘴模型

另一种重要的喷嘴类型就是压力－旋流雾化喷嘴。气体透平工业的人把它称作单相喷嘴（simplex atomizer）。这种喷嘴，然后流体通过一个称作旋流片的喷头被加速后，进入中心旋流室。在旋流室内，旋转的液体被挤压到固壁，在流体中心形成空气柱，然后，液体以不稳定的薄膜状态从喷口喷出，破碎成丝状物及液滴。在气体透平、燃油炉、直接喷射点火式汽车内燃机的液体燃料燃烧中，压力－旋流雾化喷嘴使用很广泛。液体从内部流到完全雾化的过程可分为三个步骤：液膜形成、液膜破碎及雾滴形成。 3. 空气辅助雾化模型：

为了加速液膜的破碎，喷嘴经常会添加上辅助空气。液体通过喷座的作用形成液膜,空气则直接冲击液膜以加速液膜的破碎。这种雾化被称为空气辅助雾化或气泡（air-blast）迸裂雾化（依赖于空气量及其速度）。通过辅助空气的作用，可以得到更小的雾滴。这种提高雾化质量的作用机制并不清楚。一般的看法是，辅助空气加剧了液膜的不稳定性。同时，空气有助于液滴的分散，防止液滴间的碰撞。空气雾化喷嘴同压力－旋流雾化喷嘴一样被广泛应用，尤其是用在要求雾化粒径很小的场合。

FLUENT 中的空气雾化模型由压力－旋流雾化模型所衍生。有个不同点是，在空气辅助雾化模型里，用户需要直接设定液膜厚度。因为空气雾化喷嘴的液膜形成机制很多，所以必须要有此设定。这样，在空气辅助雾化喷嘴模型里就没有压力－旋流雾化模型里所具有的液膜形成方程。

用户还必须设定液膜与空气间的最大速度差。尽管这个量可以计算出来，但设定一个值之后用户就可以不必计算喷嘴的内部流动了。这个特点对大区域（喷嘴相对很小）的流动模拟很方便。

其他方面的设定与压力－雾化喷嘴模型的一样。用户必须设定质量流率和喷射角度。液膜离开喷口之后，它的初始轨道沿着设定的喷射角。如果初始液膜的轨道指向中心线，那么，喷射角度为负值。用户还需要设定喷口处液膜的内外半径。

20