超常颗粒两相流动—林建忠.ppt

6.2 不同长径比颗粒的沉降 长径比 为 5， 固粒轴 线的稳 定取向 是水平 方向 长径比 为 7.5。 长径比 为 2.8 的颗粒 转向速 度最快 6.3 颗粒沉降的实验结果 6.4 沉降速度与倾斜角的关系 6.5 不同长径比颗粒终了雷诺数 6.6 不同长径比颗粒的漂移速度 6.7 颗粒旋转到水平位置所需时间 6.8 两个颗粒不同初始位置的沉降 6.9 两个颗粒沉降实验装置 6.10 两个颗粒沉降实验结果 6.11 两个颗粒沉降结论 初始水平平行，颗粒将逐渐保持平行 地分开，最后达到一稳定状态。初始垂 直平行，后者将落在前者尾流中，相互 作用导致沉降加快。初始交叉放置，后 面颗粒在尾流作用下旋转并倾斜，一旦 赶上前面颗粒，两颗粒相对滑动。颗粒 作用对沉降的影响依赖于作用的模式。 7. 颗粒对于湍流场的影响 圆柱状颗粒对湍流场的影响比圆球颗 粒的情形大，因为相同体积下，前者影 响范围大。研究颗粒对湍流场影响，要 用双向耦合，耦合不仅体现在平均量， 而且要体现在脉动量，即要对雷诺平均 运动方程、湍动能、耗散率方程修正。 纤维悬浮湍流瞬时运动基本方程 a颗粒取向张量，ε是应变率张量，μf是 表观粘度，与颗粒浓度和长径比有关。 瞬时量表示为平均加脉动代入上式 7.1 湍流瞬时运动基本方程 纤维悬浮湍流瞬时运动基本方程 7.2 湍流平均运动基本方程 和 为零。 纤维瞬时取向概率分布函数方程 Ψ(p) 是颗粒瞬时取向概率分布函数， p 是颗粒取向矢量， Ψ(p) dp是位于p与 dp之间颗粒的概率分布。 旋转角速度，ω涡张量，λ长径比相关 7.3 颗粒取向概率分布函数方程 纤维平均取向概率分布函数方程 代入以下瞬时取向方程后取平均 线位移与角位移扩散系数及最终方程 7.4线、角位移扩散系数及最终方程 雷诺应力表示式 需要确定含颗粒的 k, ε方程 7.5 雷诺应力表示式 含纤维作用项的湍动能方程 7.6 含颗粒作用项的湍动能方程 含纤维作用项的湍流耗散率方程 7.7 含颗粒作用项的耗散率方程 最终方程组 7.8 最终方程组 7.9 含圆柱状颗粒湍射流场 这是一种很普遍的流场，颗粒对于射 流场的平均速度分布、湍动能和湍耗散 率都会产生影响。 ○: C=0.002% □:C=0 : C=0.01% 颗粒浓度对平均速度的影响 ?: C=0.006% Re=7000, λ=15 ○: C=0.002% □:C=0 : C=0.01% 颗粒浓度对湍动能的影响 ?: C=0.006% Re=7000, λ=15 3.7 阻力、升力、力矩的拟合式 将圆柱置于风洞中，改变圆柱与来流 夹角和来流速度，可得实验数据拟合式： (1) 阻力 (2) 升力 (3) 力矩 4. 颗粒取向分布函数方程及求解 含圆柱状颗粒两相流的基本方程为： ε是应变率张量，μf 是表观粘度，与颗粒 密度和长径比有关， a是颗粒取向张量 p是颗粒取向矢量，Ψ(p)取向分布函数 纤维瞬时取向概率分布函数方程 Ψ(p) dp是取向角位于p与dp之间的颗 粒的概率分布： 是旋转角速度，ω是涡张量，λ是与长 径比相关的参数。 4.1 颗粒取向分布函数方程 4.2 联立方程 两相流方程： 颗粒取向张量方程和分布函数方程： 4.3 方程的求解 以上方程可耦合求解，也可迭代求解， 求解后可得到颗粒的取向分布和两相流 的速度分布，但不能得到单颗粒的轨迹 和空间分布。若要得到颗粒运动信息， 则要由牛顿第二定律建立颗粒的动力学 方程求解，或采用格子玻尔兹曼等方法 进行直接数值模拟。 5. 圆柱颗粒在流场中的取向分布 与圆球颗粒不同，圆柱状颗粒取向分布 是这类多相流研究的重点之一，因为取向 分布决定多相流的整体特性，影响产品性 能。流体运动特性影响颗粒取向分布，而 颗粒取向分布也会影响流体流动特性。所 以不同流动特性会导致不同颗粒取向分布。 5.1 混合层中的取向分布 U1 U2 流体方程和谱方法求解 颗粒方程 根据细长体理论的颗粒点力和速度关系 颗粒受力和力矩方程 混合层卷起成涡 Stokes=0.01的颗粒运动及取向 St=颗粒松弛时间/流场特征时间 Stokes=0.1的颗粒运动及取向 Stokes=1.0的颗粒运动及取向 计算与实验结果对比 颗粒偏离中心的程度 颗粒轴线与水平的夹角