《两相流与热物理-刘婕》第三章第七节_2_-1.pptVIP

PDF模拟应用实例--同轴突扩通道内喷雾蒸发 轴向平均速度 轴向脉动速度 能正确预报流体和颗粒速度脉动的相对关系 在颗粒浓度低的回流区内得到和主流区一样好的结果 容易处理宽筛分和颗粒的复杂经历 平均直径 § 7、11格子波尔兹曼 LBM 模型及算法 华中科技大学煤燃烧国家重点实验室对格子波尔兹曼 LB 模型进行了研究，在非线性渗流LB模型、燃烧LB模型、微细尺度流动的建模准则等研究具有特色，将格子波尔兹曼方法应用于湍流、多相流、渗流、燃烧、流固耦合中，获得好的结果 [23]。 若假定气相湍流场是局部均匀和各向同性时，则 、 和 可由下式确定: （7.35） 如果不用上述假定条件，可用Boussinesq假定求得 7.36 式中 和 、 为数量级为1的常数。 颗粒的随机轨迹则由下列公式确定： 7.37 为颗粒位置半径 在进行足够多次随机取样计算后，我们可以得到颗粒的统计平均速度和统计平均轨迹。 这样随机轨迹模型运用时需注意的关键问题是怎样确定积分时间 ，根据本方法中取随机数的原则，应取： 7.38 式中 是某随机涡旋生存时间，而 是颗粒穿过某随机涡旋的穿越时间[14]。 随机涡旋的生存时间应为： 7.39 其中 为随机涡的尺寸（即Kolmogorov涡旋尺寸）Hinze曾给出了计算公式： 7.40 而 ，将其代入上式得： 7.41 而 ，故有 7.42 颗粒穿过某随机涡旋的时间 可由下列方法得到，根据颗粒运动微分方程式，在忽略重力和假定颗粒特征响应时间 在颗粒穿过某涡旋中不变，则可得下列方程： 7.43 ＝常数 进一步假定颗粒穿过随机涡旋时，随机涡旋速度不变，上式整理后得： 积分上式 可得 或 （7.44） 颗粒在穿越随机涡旋中的轨迹方程为： 7.45 代入（7、44）积分上式 整理后就可得： 7.46 当 时， 无解，此时 应取 代入计算。 除上述的这种随机处理颗粒湍流扩散的方法外，还有不少其它的颗粒湍流扩散随机模拟方法。 例如，文献[15]曾应用Monte-Carlo方法处理颗粒群的湍流扩散，求解了气-固两相管道流动问题， 文献[16]则用随机函数取样方法求解了煤粉的两相自由射流。 7、6、3 脉动频谱随机轨道模型（FSRT模型） 上节中所述的描述颗粒湍流扩散的各种随机模型都是以气相湍流脉动速度的平均平方根值为基础，给予不同的随机数以达到模拟颗粒的湍流扩散特性。 但是，根据第六节章中所讨论过的关于湍流运动的基本概念可知，湍流运动实际上是由具有各种不同周期、振幅和方向的三元脉动随机地组合在一起的结果，每个湍流涡旋（或湍流微团）有着不同的脉动频率、频谱，振幅和方向等，故用一个平均平方根速度很难反映出这种湍流涡旋的本质。 另外，如第六节中所分析的那样，如果我们以周期性脉动气流近似的代替湍流，研究颗粒在周期性脉动气流中运动的规律可知：在不同脉动频率下，同一尺寸的颗粒运动轨迹随着脉动频率的增加而上下脉动次数增加，但脉动幅度减小，若脉动相对振幅增加，则颗粒脉动幅度也随着增加。在一定的脉动频率和脉动振幅的气流中，小尺寸的颗粒（将随气流脉动飘扬，但随着颗粒直径的增大，颗粒的脉动将受到惯性力的作用而衰减。当颗粒尺寸较大时（，基本上颗粒将不随气流脉动或脉动很小而很快衰减。 由这些定性的分析可见，不同尺寸的颗粒在湍流脉动气流中扩散运动的情况是不同的。 同时，湍流的结构（如脉动频率、频谱、振幅和方向等）对颗粒的湍流扩散将有着很大的影响。 因此简单地用颗粒的浓度梯度来描述颗粒的扩散（如漂移速度或漂移力修正方法）或用气相场脉动速度的平均平方根值随机地考虑颗粒在复杂的湍流运动中的扩散情况都是不够完善的，和实际的情况相差较大。 为了充分地考虑不同颗粒尺寸下所产生的不同湍流扩散效果和更接近地描述颗粒群在复杂的气流湍流运动作用下的情况，在PSIC方法的基础上，参考文献[17]提出了一种新的描述颗粒-气体多相流动的数值计算模型。 该模型用双方程来求解气相湍流场，并用随机的富氏级数来模拟气流的湍流脉动速度，这种随机的模拟方法不仅考虑到气流湍流脉动频谱和频率的影响，而且还考虑到方向的随机性，故这种方法称作脉动频谱随机轨迹模型（Fluctuation-Spectrum-Random-Trajectory Model简称FSRT模型）。 用拉格朗日方法求解不同尺寸的颗粒沿轨迹的运动情况，并用体积平均方法求出颗粒相的浓度场和速度场，同时，该模型也计及了两相之间质量、动量和能量的耦合。 该模型的气相运动方程组及颗粒相源项的计算方法和PSIC方法一样，见前所述。所不同的是对颗粒相运动方程式进行了新的改进和处理。颗粒的Lagrangi