两相流：颗粒群受到的曳力

最新推荐文章于 2024-07-10 11:04:17 发布

CFD_Tyro

最新推荐文章于 2024-07-10 11:04:17 发布

阅读量1.6k

点赞数 1

分类专栏：理论基础文章标签：算法

原文链接：https://www.jianshu.com/p/b6efc44797b9

版权

理论基础专栏收录该内容

11 篇文章 5 订阅

订阅专栏

此文转载自本人在简书平台上发布的文章

颗粒和流体之间的作用力是两相流领域最重要的本构关系。任何科学问题的研究都遵循先易后难的原则，对颗粒流体间作用力，同样是先从单颗粒、低雷诺数的情况研究，这篇文章介绍了处在爬流状态下单颗粒的受力。本文介绍颗粒群与流体间的作用力，尤其是颗粒群所受曳力。

一、颗粒群曳力关联式研究思路

符号说明： $u$ 无穷远处来流速度； $u_0$ 床层内流体真实速度； $U$ 表观速度，即床层进口流体速度，也是化工原理中的“空塔速度”， $U=(1-\phi)u_0$ ；颗粒群床层固含率 $\phi$ ；粒径 $d$

单颗粒曳力

$\phi \rightarrow 0$ ， $Re\rightarrow 0$ ，Stokes定律： $C_D=24/Re$

$F=3\pi \mu d u=\frac{1}{2}\rho_g u^2 \frac{\pi}{4}d^2 C_D$

$\phi \rightarrow 0$ ， $R e < 1$ ，Oseen公式， $C_D=\frac{24}{Re}(1+\frac{3}{16}Re)$
$\phi \rightarrow 0$ ， $500<Re<2\times 10^5$ ，Newton公式， $C_D=0.44$
$\phi \rightarrow 0$ ，标准阻力公式， $C_D=\frac{24}{Re}(1+\frac{1}{6}Re^{2/3})$

颗粒群曳力

流体流过静止的颗粒床层，床层固含率为 $\phi=NV_p/V$ ，底部入口流速为 $U$ ， $U=(1-\phi)u_0$

颗粒受到两种力：由流体运动产生的曳力 $F_d$ ，压力梯度力 $F_p=-V_p\nabla P$ ，因此颗粒受到的合力为 $F_{gs}=F_d+F_p$

流体流过床层的压降损失大小应等于颗粒对其的作用力：

$(P_{bottom}-P_{top})S_{bed}=V\nabla P=NF_{gs}$

因此， $F_d=(1-\phi)F_{gs}$ ，定义颗粒群中单颗粒曳力与Stokes曳力之比为无量纲曳力，

$F(\phi,Re)=F_d/3\pi \mu d U , \quad F(0,0)=1, \quad Re=\frac{dU\rho_g}{\mu}$

工程中常用床层压降来关联曳力，转化成床层压降为

$F(\phi,Re)=-\frac{1-\phi}{\phi} \frac{\nabla P d^2}{18\mu U}$

目前有了单颗粒无量纲曳力的表达式，研究者进一步从以下两种途径来从简单情况推向复杂情况

方式一：从低雷诺数出发

$F(\phi, Re)=F(\phi,0)+\alpha (\phi)Re \tag{1}$
通过第二项来添加 $R e$ 的影响。
研究者提出了以上两个方程的形式：

$F(\phi,0)=\frac{a\phi}{18(1-\phi)^2},\quad \alpha(\phi)=\frac{b}{18(1-\phi)^2} \tag{2}$

Ergun总结了大量实验，得出上式中 $a = 150, b = 1.75$ ，其他研究者也提出了不同的 $ab$ 值，例如著名的Koch-Hill-Ladd曳力。Van der Hoef给出了目前最精确的 $F(\phi,0)$ 表达式：

$F(\phi,0)=\frac{180\phi}{18(1-\phi)^2}+(1-\phi)^2 (1+1.5\sqrt{\phi}) \tag{3}$

方式二：从低固含率出发

$F(\phi, Re)=F(0,Re)(1-\phi)^{-\beta} \tag{4}$
第一项 $F (0, R e)$ 在前文《单颗粒曳力》已介绍过。通过第二项来添加 $\phi$ 的影响，其确定方式常通过不同 $\phi$ 的颗粒群沉降实验测量，测量沉降速度 $v_t(0)/v_t(\phi)$ 的比值，然后将其关联为 $\phi$ 的函数：