流体总应力可分解为各向同性压力应力和偏黏性应力。按照连续介质力学中常用的符号约定,柯西应力张量可写为 (\boldsymbol{\sigma} = -p\mathbf{I} + \boldsymbol{\tau}),其中 (p) 为静压或热力学压力,(\boldsymbol{\tau}) 为黏性应力张量。这种分解使中试装置模型能够区分压力梯度力和黏性应力:前者负责输运和加速流体,后者表示由变形引起的动量传递。
关键在于,压力在各个方向上作用相同,而黏性应力取决于速度梯度。二者的散度共同出现在动量平衡方程中,并决定压降、壁面摩擦、混合转矩和流量分布。
应力分解的含义
总应力包括压力和黏性效应
作用于流体微元上的总应力由柯西应力张量表示:
[ \boldsymbol{\sigma} = -p\mathbf{I}+\boldsymbol{\tau} ]
其中,(\mathbf{I}) 为单位张量,(-p\mathbf{I}) 为各向同性压力贡献,(\boldsymbol{\tau}) 为黏性或偏应力贡献。
负号遵循通常的约定,即正压力表示压缩应力,而拉伸法向应力为正。
压力是各向同性部分
压力在 (x)、(y) 和 (z) 方向上产生相同的法向应力:
[ \sigma_{xx}^{(p)}=\sigma_{yy}^{(p)}=\sigma_{zz}^{(p)}=-p ]
由于各向同性压力场没有优先方向,因此不会产生剪切应力。
对于一般应力张量,热力学压力与平均法向应力的关系为:
[ p=-\frac{1}{3} \left( \sigma_{xx}+\sigma_{yy}+\sigma_{zz} \right) ]
因此,黏性应力就是剩余的偏应力部分:
[ \boldsymbol{\tau}=\boldsymbol{\sigma}+p\mathbf{I} ]
因此,黏性应力张量的迹为零:
[ \tau_{xx}+\tau_{yy}+\tau_{zz}=0 ]
这一零迹条件使黏性张量成为偏应力贡献,而不是体积应力贡献。
分解如何进入动量方程
从一般动量平衡开始
线动量的局部平衡方程为:
[ \rho\frac{D\mathbf{v}}{Dt}
\mathbf{f} + \nabla\cdot\boldsymbol{\sigma} ]
其中,(\rho) 为密度,(\mathbf{v}) 为速度,(D/Dt) 为物质导数,(\mathbf{f}) 表示单位体积上的体力,例如重力。
代入 (\boldsymbol{\sigma}=-p\mathbf{I}+\boldsymbol{\tau}) 可得:
[ \rho\frac{D\mathbf{v}}{Dt}
\mathbf{f}
\nabla p + \nabla\cdot\boldsymbol{\tau} ]
这就是柯西动量方程的压力-黏性形式。
压力梯度驱动整体输运
项
[ -\nabla p ]
表示由静压空间变化产生的力。例如,在中试规模管道中,施加的压降会推动流体加速通过管网。
在重力作用下的静止液体中,压力随标高变化。若 (z) 轴向上测量:
[ \nabla p=\rho\mathbf{g} ]
其中 (\mathbf{g}) 指向下方。当流体静止时,该静压场与重力相平衡。
黏性应力阻碍变形
项
[ \nabla\cdot\boldsymbol{\tau} ]
表示由黏性应力空间变化产生的净黏性力。它负责壁面摩擦、速度分布的形成以及相邻流体层之间的动量扩散等效应。
流体静止时可以存在压力,但黏性应力通常需要速度梯度或变形才能产生;复杂材料中的特殊本构效应除外。
黏性应力与速度的关系
牛顿流体关系
对于各向同性牛顿流体,黏性应力与变形速率张量呈线性关系:
[ \boldsymbol{\tau}
2\mu\mathbf{D} + \lambda(\nabla\cdot\mathbf{v})\mathbf{I} ]
其中
[ \mathbf{D}
\frac{1}{2} \left[ \nabla\mathbf{v} + (\nabla\mathbf{v})^T \right] ]
这里,(\mu) 为动力黏度,(\lambda) 为与体积变形相关的第二黏度系数。
对于不可压缩流体,
[ \nabla\cdot\mathbf{v}=0 ]
因此本构关系可简化为:
[ \boldsymbol{\tau}
\mu \left[ \nabla\mathbf{v} + (\nabla\mathbf{v})^T \right] ]
在恒定黏度和密度等假设下,动量方程随即成为熟知的牛顿流体形式。
剪切应力源于速度梯度
对于沿 (x) 方向流动且速度随 (y) 方向变化的简单流动,主要的黏性剪切应力为:
[ \tau_{xy}
\mu\frac{\partial u}{\partial y} ]
这一关系直接适用于中试装置中的管流、搅拌罐、边界层和流道。
产生壁面剪切的同一速度梯度也会导致压力损失和泵功率需求。
也可能出现黏性法向应力
黏性效应并不局限于剪切分量。在加速流、收缩流、扩张流或拉伸流中,(\boldsymbol{\tau}) 的对角分量可能具有重要作用。
对于不可压缩牛顿流动,这些黏性法向应力与法向速度梯度相关,同时保持偏应力张量的零迹条件。
在中试装置设备中的分解应用
管网和流道
在管道中,压力梯度提供主要驱动力,而黏性应力将动量传递给管壁。二者形成的平衡决定速度分布和压降。
对于中试装置中的湍流,模型应力同时包括分子黏性和速度波动的影响。雷诺应力通常表示为:
[ \tau_{xy}^{t}
-\rho\overline{u'v'} ]
其中,(u') 和 (v') 为波动速度分量。
混合罐和反应器
在搅拌容器中,压力影响整体力场和静压分布。黏性应力控制局部能量耗散、剪切作用以及叶轮所需的转矩。
对于非牛顿反应器内容物,局部黏度可能随变形速率变化。应力分解仍然有效,但 (\boldsymbol{\tau}) 的本构模型会变为非线性模型。
分离塔和重力设备
在沉降柱、倾析器和浮选设备中,静压随深度变化,并有助于确定相界面和液位行为。黏性应力控制沉降阻力、界面运动和近壁流动。
因此,压力场和黏性场具有不同的建模用途:压力确定法向力分布,而黏性应力描述与变形相关的动量传递。
换热器和热流设备
在换热器流道中,壁面剪切和压降与速度梯度场相关。同一流动结构还会影响热边界层和对流传热。
在高雷诺数下,湍流动量输运可能显著超过分子黏性输运,因此湍流模型必须考虑额外的表观应力。
理解其中的权衡
压力不等同于总法向应力
一个常见的建模错误是将所有法向应力都视为静压。在可压缩流、拉伸流、快速加速流或非牛顿流中,黏性法向应力可能十分重要。
压力是应力状态的各向同性部分;其余法向应力属于偏应力张量。
湍流应力需要建模,而不是仅使用分子黏度
在湍流中,时间平均动量方程包含由速度波动产生的雷诺应力。这些应力不同于分子黏性应力,但湍流模型通常会通过有效黏度或涡黏度来表示其综合影响。
在湍流中试设备中仅使用分子黏度,可能会低估壁面摩擦、压降、混合需求和泵功率。
必须保持符号约定一致
一些参考资料使用 (+p\mathbf{I}) 表示以压缩压力为正的总应力张量,而另一些资料使用 (-p\mathbf{I})。两种约定都可以使用,但动量方程和法向应力定义必须始终采用同一约定。
标准工程形式为:
[ \boldsymbol{\sigma}=-p\mathbf{I}+\boldsymbol{\tau} ]
不可压缩性不会消除所有黏性应力
对于不可压缩流体,(\nabla\cdot\mathbf{v}=0) 会从牛顿流体本构方程中去除体积膨胀项,但不会消除由其他速度梯度引起的剪切应力或偏法向应力。
非牛顿流体需要合适的本构定律
对于聚合物溶液、浓缩悬浮液和发酵液,黏度可能取决于应变速率的大小。广义模型可写为:
[ \boldsymbol{\tau}
\mu(\Gamma) \left[ \nabla\mathbf{v} + (\nabla\mathbf{v})^T \right] ]
其中,(\Gamma) 表征局部变形速率。
压力-黏性分解仍然适用,但必须根据流体实测的流变行为计算黏性应力。
如何将其应用于您的项目
正确的分解方式取决于所使用的设备、流动状态和流体模型。
- 如果您的重点是管道压降:分别计算压力梯度贡献、壁面黏性应力和体黏性应力;当流动为湍流时,还应计入湍流雷诺应力。
- 如果您的重点是混合或反应器设计:使用黏性应力场评估叶轮转矩、局部能量耗散和剪切作用,同时保留压力以进行整体力平衡和静压平衡分析。
- 如果您的重点是分离设备:利用静压建立随深度变化的法向载荷,并使用黏性应力对沉降、界面运动和阻力进行建模。
- 如果您的重点是非牛顿加工:保持压力分解不变,但选择能够将偏应力与实测速率相关黏度联系起来的本构定律。
- 如果您的重点是瞬态或湍流操作:区分分子黏性应力与建模得到的湍流应力,避免低估压降和动量输运。
将各向同性压力与依赖变形的应力分离,可以为中试装置模型提供清晰的物理基础,从而预测流动、摩擦、混合、分离和放大行为。
总结表:
| 应力分量 | 物理含义 | 主要建模作用 | 中试装置应用 |
|---|---|---|---|
| 静压 | 在所有方向上均等作用的各向同性法向应力 | 驱动流动并建立静压载荷 | 管道、流道、塔、倾析器 |
| 分子黏性应力 | 由速度梯度引起的偏应力 | 对壁面摩擦、动量扩散和变形进行建模 | 泵、换热器、搅拌罐 |
| 湍流雷诺应力 | 由速度波动引起的动量输运 | 表示湍流流动中的额外表观应力 | 高流量管道、混合、放大 |
| 非牛顿应力 | 基于流体流变特性、随速率变化的偏应力 | 预测复杂流体中的剪切、转矩和流动行为 | 聚合物溶液、悬浮液、发酵液 |
需要有关压降、混合、分离或传热设计的可靠中试装置数据?联系 LABPARK,讨论您的需求。LABPARK 为化学工程、生物过程与生物技术以及环境与水处理应用提供教育和职业培训用单元操作中试装置。我们的系统帮助大学、研究机构和企业评估压力与黏性效应、培训操作人员、验证过程模型,并借助实用且面向应用的设备支持工艺放大。