粘性应力张量连接了中试工厂仪器的测量结果与流体动量传输的方式。 在牛顿流体中,粘性应力是根据局部速度梯度和动力粘度计算得出的。在化学工程中试工厂中,该公式用于解释壁面剪切、速度分布、压降、混合行为、换热器流动以及放大决策。
应力张量将流体变形转化为机械应力。对于不可压缩牛顿流体,压力代表总应力的各向同性部分,而粘性张量描述由速度梯度引起的剪切和变形。
粘性应力张量的含义
流体点的应力
流体单元会受到来自周围流体的表面力。这些力随表面的方向而变化,因此它们由二阶张量而非单个标量来表示。
总柯西应力张量通常分解为
[ \boldsymbol{\sigma}=-p\mathbf{I}+\boldsymbol{\tau}, ]
其中 (p) 是热力学压力,(\mathbf{I}) 是单位张量,(\boldsymbol{\tau}) 是粘性应力张量。
上述符号约定将正压力视为压缩。其他文献可能使用不同的符号,因此在实验或模型中必须保持约定一致。
张量的对称性
对于没有分布体偶力的常规流体,应力张量是对称的:
[ \tau_{xy}=\tau_{yx},\qquad \tau_{xz}=\tau_{zx},\qquad \tau_{yz}=\tau_{zy}. ]
尽管 (3 \times 3) 张量包含九个条目,但这种对称性留下了六个独立的应力分量:三个正应力和三个剪应力。
这一特性在中试工厂计算中非常重要,因为作用在表面上的应力和相应的力矩平衡必须在物理上兼容。特殊的极性流体或受显著体扭矩影响的流体需要更通用的理论,但这超出了大多数常规单元操作实验的范围。
速度梯度如何产生粘性应力
牛顿本构方程
对于各向同性牛顿流体,粘性应力与变形速率张量线性相关:
[ \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) 是第二粘度系数,(\mathbf{v}) 是速度场。张量 (\mathbf{D}) 捕获了与流体变形相关的速度梯度部分。刚体旋转改变了流体单元的方向,但其本身并不会产生牛顿粘性应力。
不可压缩液体操作
对于不可压缩流动,
[ \nabla\cdot\mathbf{v}=0. ]
因此本构方程简化为
[ \boldsymbol{\tau}=2\mu\mathbf{D}. ]
对于 (x) 方向的速度随 (y) 变化的简单剪切流,
[ v_x=v_x(y), ]
相关的剪应力变为
[ \tau_{xy}
\mu\frac{dv_x}{dy}. ]
这是在层流管道流动演示、壁面摩擦装置、粘度计实验和简化混合分析中最直接应用的公式。
正应力和剪应力分量
非对角项(如 (\tau_{xy}))描述剪应力。它们出现在相邻流体层以不同速度运动时。
对角项描述粘性正应力。在不可压缩、稳态、恒定粘度的简单剪切流中,这些项通常消失或变得可以忽略不计,但在拉伸流、收敛通道、射流和复杂的混合器几何结构中,它们可能很重要。
张量如何进入动量平衡
局部动量方程
对无穷小流体单元应用动量平衡得到
[ \rho\frac{D\mathbf{v}}{Dt}
\rho\mathbf{g} + \nabla\cdot\boldsymbol{\sigma}, ]
其中 (\rho) 是密度,(\mathbf{g}) 是重力加速度,(D\mathbf{v}/Dt) 是物质加速度。
使用应力分解得到
[ \rho\frac{D\mathbf{v}}{Dt}
\rho\mathbf{g}
\nabla p + \nabla\cdot\boldsymbol{\tau}. ]
压力梯度代表驱动或抵抗流体整体运动的各向同性力。粘性张量的散度代表由粘性应力空间变化引起的单位体积净粘性力。
与纳维-斯托克斯方程的联系
对于恒定粘度的牛顿流体,当流动不可压缩时,粘性项简化为
[ \nabla\cdot\boldsymbol{\tau}
\mu\nabla^2\mathbf{v} ]
动量方程变为
[ \rho\frac{D\mathbf{v}}{Dt}
\rho\mathbf{g}
\nabla p + \mu\nabla^2\mathbf{v}. ]
该方程是解释中试规模设备中速度和压力测量结果的数学基础。
表面牵引力和张量运算
作用在单位法向量 (\mathbf{n}) 表面上的牵引力由应力张量获得:
[ \mathbf{t}=\boldsymbol{\sigma}\cdot\mathbf{n}. ]
结果是一个包含该表面上法向和切向力分量的向量。在计算管道壁、填料表面、叶轮叶片、换热器壁和其他设备边界上的力时,此运算至关重要。
由于张量-向量乘法取决于运算顺序,因此必须始终使用所采用的索引和矩阵约定。对于常规牛顿流体的对称应力张量,相关的物理关系与力矩平衡一致。
在中试工厂单元操作中的应用
管道流动和流体摩擦实验
在充分发展的圆管中,轴向压力梯度与剪应力的径向变化相平衡。壁面剪应力与测量管道长度上的压降直接相关:
[ \tau_w
-\frac{R}{2}\frac{dp}{dx}, ]
对于半径为 (R) 的管道,服从稳态、充分发展流动的通常假设。
在层流状态下,将此平衡与牛顿本构关系结合会产生抛物线速度分布。因此,可以将中试工厂的流量和压降测量值与预测的壁面应力、摩擦系数和泵需求进行比较。
泵和管道网络
张量公式确定了动量如何通过固体壁面附近和管件中的速度梯度耗散。这些局部损失构成了泵必须提供的压力头。在中试规模下,由此产生的压降数据有助于确定泵、管径、阀门布置或运行流量是否适用于更大的工艺。当比较不同粘度的流体时,应力模型特别有用。
换热器
在换热器通道中,粘性应力决定了速度分布和水力阻力。速度场随后影响停留时间、边界层发展和对流换热系数。仅靠应力模型无法预测传热。它提供了必须与能量方程耦合的动量场,并在相关时考虑温度相关的粘度和密度。
混合罐和反应器
在搅拌罐中,叶轮产生强烈的速度空间变化。粘性应力张量描述了局部变形率,而由此产生的流场影响混合、循环、悬浮和能量耗散。对于高粘度流体,层流或过渡流的应力预测可以与混合器扭矩和功率密切相关。对于湍流操作,平均流应力还必须包括速度波动的影响。
填料塔和分离设备
在填料床和塔中,速度梯度出现在填料表面周围和狭窄的流动通道内。相关的应力有助于摩擦压降并影响液体分布。当局部不均匀、沟流或壁面效应很重要时,张量框架最有用。仅基于平均流量的相关性可能足以进行常规设计,但它们无法揭示局部应力集中。
从层流到湍流
雷诺应力
在湍流中,瞬时速度被分解为平均部分和波动部分:
[ u_i=\overline{u_i}+u_i'. ]
对动量方程进行平均会引入雷诺应力,例如
[ \tau_{xy}^{R}
-\rho,\overline{u'v'}. ]
这些不是分子粘性应力。它们代表由湍流速度波动引起的动量传输,在求解平均流时必须对其进行建模或测量。
表观动量传输
与仅使用分子粘度相比,雷诺应力增加了有效动量传输。这就是为什么通常不能在不进行修改的情况下,通过使用层流本构方程来预测湍流压降、壁面摩擦和混合性能的原因。
在工程模型中,湍流闭合、涡粘模型或实验拟合的摩擦相关性被用来表示这种额外的传输。
近壁面行为
不透水壁面附近的粘性底层由分子粘度主导,而远离壁面时湍流传输变得越来越重要。数值场方法可以直接解析该区域,或使用壁面函数来连接近壁面和外流区域。
中试工厂可以使用测得的压降、速度分布或局部探针来评估所选的湍流模型是否准确地代表了设备。
理解权衡
张量模型是局部的,仪器通常是平均的
应力张量是逐点定义的,但中试工厂仪器通常报告平均量,如总流量、进出口压力、扭矩或整体温度。测得的压降可以验证应力场的综合结果,但它不能唯一确定每个局部应力或速度值。
牛顿假设可能失效
关系式 (\boldsymbol{\tau}=2\mu\mathbf{D}) 假设牛顿流体的粘度与变形历史和剪切速率无关。聚合物溶液、浆料、乳液、糊状物和某些生物流体可能是剪切变稀、剪切变稠、粘弹性或屈服应力流体。对于这些材料,必须用适当的流变模型代替本构方程。使用恒定粘度可能会对压降、混合器功率和放大行为产生误导性预测。
不可压缩性有限制
当密度变化很小时,条件 (\nabla\cdot\mathbf{v}=0) 通常适用于液体单元操作。它可能不适用于高速气体、强压力变化、可压缩流动或具有显著热膨胀的工艺。当可压缩性影响实验时,正确的公式必须考虑密度变化。
放大并不能保留所有量
仅匹配几何尺寸并不能保证动态相似的应力场。雷诺数、流动状态、粘度、粗糙度、叶轮条件、停留时间和相关的无量纲组都可能影响结果。当中试实验的运行条件和测量输出被选择用于支持明确定义的放大标准时,其价值最高。
如何将其应用于您的项目
当数学公式与可测量的操作变量和特定的设计问题联系起来时,它就变得很有用。
- 如果您的主要重点是压降预测: 测量流量、流体性质、管道或通道几何形状以及压力损失,然后使用粘性应力平衡来评估壁面剪切和摩擦行为。
- 如果您的主要重点是层流速度分布: 将牛顿本构关系应用于测得的速度梯度,并将预测的分布与局部速度数据进行比较。
- 如果您的主要重点是湍流管道或通道流动: 通过适当的湍流模型或经过验证的相关性包含雷诺应力,而不是将分子粘度视为唯一的动量传输机制。
- 如果您的主要重点是混合性能: 将局部变形率和能量耗散与叶轮转速、扭矩、粘度和浓度或示踪剂测量值联系起来。
- 如果您的主要重点是换热器设计: 使用应力公式获得速度场,然后将其与能量方程耦合以预测对流换热。
- 如果您的主要重点是放大: 保留控制相关应力、压降、混合或传热行为的无量纲组和流动状态。
中试工厂将粘性应力张量从抽象的连续介质力学方程转化为测量、验证和放大流体传输的实用工具。
总结表:
| 应用 | 应力张量有助于确定什么 | 关键测量值或输入 |
|---|---|---|
| 管道流动 | 壁面剪应力、速度分布、摩擦和压降 | 流量、压力损失、半径、粘度 |
| 泵和管道 | 水力损失和泵需求 | 管道布局、管件、流量、流体性质 |
| 换热器 | 速度分布、水力阻力和传热流动条件 | 通道几何形状、流量、温度、粘度 |
| 混合罐和反应器 | 变形率、扭矩、功率、循环和混合行为 | 叶轮转速、扭矩、粘度、示踪剂数据 |
| 填料塔 | 摩擦压降、液体分布和沟流 | 填料几何形状、流量、流体性质 |
| 湍流系统 | 来自雷诺应力的额外动量传输 | 雷诺数、速度数据、湍流模型 |
| 放大研究 | 动态相似性和相关流动行为的保留 | 无量纲组、操作条件、几何形状 |
将流体力学理论与中试工厂结果联系起来
使用 LABPARK 教育和职业单元操作中试工厂,将粘性应力张量计算转化为可靠的工程数据。我们的系统帮助大学、研究机构和企业在可测量的操作条件下研究压降、壁面剪切、混合、传热、湍流、流变学和放大。LABPARK 为化学工程、生物工艺和生物技术以及环境和水处理应用提供中试工厂解决方案。使用实际实验来验证模型、比较流体行为、培训学生和操作员,并做出更好的设备和工艺决策。
联系 LABPARK 讨论您的实验目标,并确定适合您教学、研究、工艺开发或工业培训需求的中试工厂。