离散相笔记(重要).docVIP

Fluent 笔记 稳态问题的求解步骤 稳态离散相问题的设定、求解的一般过程如下： 1 求解连续相流场 2 创建离散相喷射源(射流源) 3 求解耦合流动（如果希望计算的话） 4 用PLOT 或REPORT 图形界面来跟踪离散相 非稳态问题的求解步骤 非稳态离散相问题的设定、求解的一般过程如下： 1 创建离散相喷射入口 2 初始化流场 3 设定求解的时间步长和时间步数。在每个时间步，颗粒的位置将得到更新。如果 求解问题是非耦合流动，那么，颗粒的位置在每个时间步计算完成之后得到更新的； 如果是耦合流动，那么，颗粒的位置在每个时间步内的相间耦合迭代计算过程中都 会得到更新。 FLUENT 提供的离散相模型选择如下： 对稳态与非稳态流动，可以应用拉氏公式考虑离散相的惯性、曳力、重力 预报连续相中，由于湍流涡旋的作用而对颗粒造成的影响 离散相的加热/冷却 液滴的蒸发与沸腾 颗粒燃烧模型，包括挥发份析出以及焦炭燃烧模型（因而可以模拟煤粉燃烧） 连续相与离散相间的耦合 液滴的迸裂与合并 应用这些模型，FLUENT 可以模拟各种涉及离散相的问题，诸如：颗粒分离与分级、 喷雾干燥、气溶胶扩散过程、液体中气泡的搅浑、液体燃料的燃烧以及煤粉燃烧。 19.2-19.5 介绍离散相计算中所用到的物理方程；设定、求解和后处理在19.6-19.13 中介绍。 19.1.2 湍流中的颗粒 随机轨道模型或颗粒群模型（19.2.2）可考虑颗粒湍流扩散的影响。在随机轨道 模型中，通过应用随机方法（19.2.2）来考虑瞬时湍流速度对颗粒轨道的影响。 而颗粒群模型则是跟踪由统计平均决定的一个“平均”轨道（19.2.2）。颗粒群中 的颗粒浓度分布假设服从高斯概率分布函数（PDF）。两种模型中，颗粒对连续相 湍流的生成与耗散均没有直接影响。 19.1.3 应用范围 颗粒体积分数的适用范围： FLUENT 中的离散相模型假定第二相（分散相）非常稀薄，因而颗粒-颗粒之间 的相互作用、颗粒体积分数对连续相的影响均未加以考虑。这种假定意味着分散相 的体积分数必然很低，一般说来要小于10-12%。但颗粒质量承载率可以大于10-12%， 即用户可以模拟分散相质量流率等/大于连续相的流动。 模拟连续相中悬浮颗粒的限制 稳态拉氏离散相模型适用于具有确切定义的入口与出口边界条件问题，不适用 于模拟在连续相中无限期悬浮的颗粒流问题，这类问题经常出现在处理封闭体系中 的悬浮颗粒过程中，包括：搅拌釜、混合器、流化床。但是，非稳态颗粒离散相模 型可以处理此类问题。参阅第十八、二十章来确定具体多相流问题中的适用模型。 在FLUENT 的其它模型中应用离散相模型的限制。 一旦使用了离散相模型，下面的模型将不能使用： 1.选择了离散相模型后，不能再使用周期性边界条件（无论是质量流率还是压差 边界条件） 2.可调整时间步长方法不能与离散相模型同时使用 3.预混燃烧模型中只能使用非反应颗粒模型 4.同时选择了多参考坐标系与离散相颗粒模型时，在缺省情况下，颗粒轨道的显 示失却了其原有意义；同样，相间耦合计算是没有意义的。 在 FLUENT 的缺省模式下，重力加速度等于零。 如果要考虑重力的影响，必须在Operating Conditions 面板中设定重力加速度的大小 和方向。 19.2.1 颗粒运动方程 颗粒的力平衡 FLUENT 中通过积分拉氏坐标系下的颗粒作用力微分方程来求解离散相颗粒（液滴 或气泡）的轨道。颗粒的作用力平衡方程（颗粒惯性＝作用在颗粒上的各种力）在 笛卡尔坐标系下的形式（x方向）为： 其中， u为流体相速度， up 为颗粒速度，为流体动力粘度， r为流体密度，ρp为 颗粒密度（骨架密度），dp为颗粒直径，Re为相对雷诺数（颗粒雷诺数）。 对于亚观尺度（直径＝1-10微米）的颗粒，Stokes曳力公式[ 170]是适用的。这种情 况下， F D定义为： 上式中的系数Cc为Stokes曳力公式的Cunningham修正（考虑稀薄气体力学的颗粒壁 面速度滑移的修正），其计算公式为： 布朗力 对于亚观粒子（直径＝1～10微米），附加作用力可包括布朗力。布朗力的分量 可由高斯白噪声过程来模拟，其谱密度Sn，ij, 由文献[ 135]给出： 为考虑布朗力的影响，必须要激活能量方程选项。只有选择了非湍流模型才能 激活布朗力选项。 P611 热泳力（热致迁移力或辐射力） 从1到10微米 对于悬浮在具有温度梯度的气体流场中的颗粒，受到一个与温度梯度相反的作 用力。这种现象被称为热泳。颗粒平衡方程19.2-1 中的其它作用力Fx可包含这种热 泳力： 其中DT , p为热泳力系数。可以定义为常数、多项表达式或用户定义函数。也可以采 用Talbot [ 237]得到的表达式： 其中： Kn = Knudsen数