库埃特流 模拟器 返回
粘性流体模拟器

库埃特流 模拟器 — 平行平板间的粘性流

分析平行平板(上板速度 U、下板静止)间的定常层流,通过改变上板速度、隙间、粘性和压力梯度,实时可视化剪应力和单位宽度流量,以及库埃特+泊肃叶的合成速度分布。

参数设置
上板速度 U
m/s
隙间 h
mm
粘性系数 μ
Pa·s
压力梯度 dp/dx
Pa/m

假定流体为水(ρ = 1000 kg/m³)。当 dp/dx = 0 时为纯库埃特流,当 U = 0 时为纯泊肃叶流。

计算结果
平均壁面剪应力
单位宽度流量
平均流速
雷诺数
平板间速度分布

上=动板(速度 U)/下=静止板/箭头=速度矢量(库埃特直线+泊肃叶抛物线的合成)/反流用红色表示

正规化分布 u/U vs y/h

横轴=y/h [0, 1]/纵轴=u/U/蓝实线=合成分布/灰虚线=纯库埃特直线(参考)

理论与主要公式

平行平板间(上板速度 U、下板静止、隙间 h)的定常层流可从纳维叶-斯托克斯方程获得解析解,为库埃特直线和泊肃叶抛物线的叠加。包含压力梯度 dp/dx 的一般解为:

$$u(y) = U\,\frac{y}{h} - \frac{1}{2\mu}\,\frac{dp}{dx}\,y(h-y)$$

上下板的壁面剪应力(牛顿粘性律 $\tau = \mu\,du/dy$):

$$\tau_{\mathrm{top}} = \frac{\mu U}{h} - \frac{h}{2}\frac{dp}{dx},\qquad \tau_{\mathrm{bot}} = \frac{\mu U}{h} + \frac{h}{2}\frac{dp}{dx}$$

单位宽度体积流量和平均流速:

$$Q' = \int_{0}^{h} u\,dy = \frac{Uh}{2} - \frac{h^{3}}{12\mu}\frac{dp}{dx},\qquad V_{\mathrm{avg}} = \frac{Q'}{h}$$

雷诺数(以隙间 $h$ 为特征长度,以 $U$ 为特征速度):

$$Re = \frac{\rho\,U\,h}{\mu}$$

其中 $U$ 为上板速度 [m/s]、$h$ 为隙间 [m]、$\mu$ 为粘性系数 [Pa·s]、$dp/dx$ 为流动方向压力梯度 [Pa/m]、$\rho$ 为密度 [kg/m³]、$y$ 为下板处的高度。

库埃特流模拟器简介

🙋
我第一次听说"库埃特流"。这和哈根-泊肃叶流有什么区别吗?
🎓
泊肃叶流是"用压力推动"的流,而库埃特流是"用板子拉动"的流。下板固定,上板以速度 U 运动,中间的流体因粘性被拖动,形成层状滑动。如果压力梯度为零,速度就是从下(u=0)到上(u=U)的漂亮直线——这就是纯库埃特流。公式就一行:$u(y) = U y/h$。
🙋
啊,不是抛物线而是直线?我记得圆管的时候是抛物线……
🎓
对,因为驱动机制不同,形状也不同。库埃特流的驱动机制只是板子"移动"这个边界条件,内部没有体积力,所以运动方程就是 $\mu\,d^2u/dy^2 = 0$——结果是线性。加上压力梯度后,右边变成常数 $\mu\,d^2u/dy^2 = dp/dx$,就出现抛物线解。模拟器显示的正是这两者的合成解:$u(y) = U y/h - (1/(2\mu))(dp/dx)\,y(h-y)$。把 dp/dx 拉到 0 就变成线性,把 U 拉到 0 就变成纯泊肃叶流。
🙋
压力梯度滑块出现的"反流"是什么意思?
🎓
当压力梯度设成 dp/dx > 0,也就是"与上板拉动方向相反"时,泊肃叶成分会反向推流体。下板附近,这个反向推力会超过上板的拉力,导致速度变成负数(反向流动)。这在轴承和齿轮间隙中确实会发生——润滑油会部分地反向流动。用 Play 按钮让 dp/dx 从 −2000 扫到 +2000 Pa/m,你就能看到:纯库埃特直线→被泊肃叶抛物线弯曲→反流出现,全过程一目了然。
🙋
说壁面剪应力是 0.5 Pa,这个大吗?小吗?
🎓
用默认值(U=1 m/s、h=2 mm、水)算的话,τ = μU/h = 0.001·1/0.002 = 0.5 Pa,这是个小数值。但要是把隙间 h 改成 0.02 mm(润滑膜的典型值),τ 就跳到 50 Pa。轴承的摩擦损失、油膜温度上升,计算发热功率 W = τ·U·面积 时就用这个公式。你手机里的精密轴承、汽车曲轴上的油膜,本质上都是这个库埃特流公式在支配。
🙋
显示 Re = 2000,这个作为层流的分析解够吗?
🎓
很好的问题。平行平板的剪应力驱动流在 Re ≈ 1500 以上可能过渡到湍流,默认值的 Re=2000 正好在边界附近。把粘度乘以 10(比如轻油的 0.01 Pa·s),Re 就跌到 200,肯定是层流。反过来把隙间乘以 10(20 mm),Re 就蹿到 20000,进入湍流区。模拟器显示的是层流分析解,Re 高的时候要慎重对待。

常见问题

按物理驱动机制来区分。库埃特流以"边界运动(壁面速度)"为驱动源,用于轴承油膜、齿轮间隙流、旋转粘度计等"动面与静面夹间的粘性流"场景。泊肃叶流以"压力差"为驱动源,用于管道流量计算、输液、微流体芯片等"从上游被推动的流"场景。实际应用中两者往往共存,例如引擎轴承就同时有"转轴库埃特驱动"和"外部泵泊肃叶驱动"的叠加。本模拟器正是展示这种合成解。
纯库埃特流(dp/dx = 0)时,平板间流体不受体积力和压力梯度,运动方程简化为 $\mu\,d^2u/dy^2 = 0$。这个二阶导数=0 的方程在边界条件 u(0)=0、u(h)=U 下的解就是 u(y) = U y/h,是线性的。物理上,粘性将剪应力从上板均匀传递到下板,流体各层以恒定的速度梯度 du/dy = U/h 滑动,像"叠叠乐推牌"那样。剪应力在深度方向保持不变,这是设计公式的基础。
加入 dp/dx ≠ 0 后,运动方程变成 $\mu\,d^2u/dy^2 = dp/dx$,抛物线成分叠加到库埃特直线上。一般解为 $u(y) = U y/h - (1/(2\mu))(dp/dx)\,y(h-y)$,上下板的壁面剪应力变成非对称:τ_top = μU/h − (h/2)·dp/dx、τ_bot = μU/h + (h/2)·dp/dx。当 dp/dx 足够大且为正(逆向)时,下板附近会出现反流。这与轴承空蚀现象和润滑膜不稳定性有关。雷诺方程(流体润滑理论的基础方程)正是以库埃特+泊肃叶的合成解为出发点推导的。
圆筒式旋转粘度计(库埃特型流变仪)正是应用库埃特流的典型。内圆筒(或外圆筒)以恒定角速度旋转,在内外圆筒间的薄隙间放入试样,通过测量旋转扭矩来测粘度。当隙间 h 远小于圆筒半径时,可将圆筒面视为平行平板近似,用本模拟器的公式 τ = μU/h 反算粘度 μ。圆锥-平板型流变仪也用同样原理,通过在全范围内保证剪速率 du/dy = U/h 的均匀性,精密测量非牛顿流体的粘度(剪速率依赖性)。广泛应用于化妆品、食品、涂料、血液等产业品质管理。

现实应用

润滑工程·摩擦学:滑动轴承、滚珠轴承、齿轮间隙中的润滑油流动,正是作为库埃特+泊肃叶合成流来处理。雷诺方程是流体润滑理论的基础方程,以本模拟器显示的一般解 $u(y) = U y/h - (1/(2\mu))(dp/dx)\,y(h-y)$ 为出发点推导。摩擦损失 W = τ·U·面积、发热量、油膜最小厚度、轴承承载能力等,都从这个方程计算。从引擎轴承、智能手机精密电机到新干线曲轴,所有有润滑膜的旋转机械的分析中心都是这个公式。

流变学·粘度测量:库埃特型流变仪(圆筒粘度计)和圆锥-平板型流变仪都基于库埃特流原理工作。通过在样品全域保持均匀的剪速率 du/dy = U/h,精密测量牛顿流体粘度,也能测量非牛顿流体(触变流体、宾汉流体、幂律流体等)的粘度函数 $\mu(\dot\gamma)$。应用于化妆品、食品、涂料、印刷油墨、血液、高分子溶液等的产业品质控制,是不可或缺的工具。

微流体·芯片实验室(LOC):微流体芯片中,隙间 h ≈ 10~100 μm 的薄通道内,电渗流(在壁面附近生成的库埃特成分)和泵驱动流(泊肃叶成分)同时存在。两者的合成分布用于细胞分选、蛋白质分离、PCR 芯片内流路设计等。本模拟器显示的"反流条件",也被有意利用在微流体中促进混合。

地壳变动·冰川流动:冰川流动、地幔对流、地壳板块运动在长时间尺度上可视为粘性流体。冰川以基岩为底(固定)、大气为顶(应力为零),边界条件不同,但近基底的速度分布与库埃特流类似。冰川粘度约 10¹⁵ Pa·s,极其巨大,流速也极缓(年数米),但用同样的运动方程解析。

常见误解与注意

最常见的误解是"库埃特流总是线性分布"。纯库埃特流(dp/dx = 0)确实线性,但实际轴承、齿轮、流变仪内必然存在压力梯度。例如轴承要承受负荷就会产生油膜内部压力分布,这就是泊肃叶成分叠加上去。本模拟器的 dp/dx 滑块可以体验:从线性→被抛物线弯曲→反流出现的全过程。雷诺方程求解轴承时必须包含压力分布,这是标准做法。

第二常见的误解是"平板的临界雷诺数和圆管一样是 2300"。平行平板剪应力驱动流的线性稳定分析表明"各雷诺数都线性稳定"这个著名结果,但实验显示 Re ≈ 360~1500 就开始湍流过渡(次临界过渡)。本模拟器只显示层流分析解,Re 高时与实际流动可能有偏差。值得注意的是不同流型有不同临界值:圆管 2300、平板上界面层 Re ≈ 5×10⁵,要按流型查表。

第三常见的误解是"单位宽度流量 Q'(m²/s)和体积流量 Q(m³/s)混淆"。平板间流把 2D 问题处理,流量的单位是"单位宽度下的值",即 m²/s。乘以实际流路宽度 W [m] 才能换算成 3D 体积流量 Q = Q'·W [m³/s]。本模拟器同时显示 Q'(单位宽度流量)和 V_avg(平均流速),按需应用。CAE 分析中如果混淆 2D 假设和 3D 换算,就会在流量计算中引入 W 倍的误差,要特别注意。

使用指南

  1. 上板速度 (slUVal) 在 0~2 m/s 范围内设置,下板保持静止
  2. 平板间隙 (slHVal) 在 0.1~10 mm 范围内指定,决定通道高度
  3. 流体粘度 (slMuVal) 以 mPa·s 单位输入,水约 1 mPa·s,油约 100 mPa·s
  4. 压力梯度 (slDpDxVal) 以 Pa/m 单位设置,控制外部驱动力有无
  5. 模拟运行后,实时显示速度分布、剪应力分布,可直观观察流场特征

具体计算示例

SAE 10W-30 发动机油(μ=100 mPa·s)充填在平板间隙 H=0.5 mm 的圆周轴承内,上板以 U=1.5 m/s 速度驱动,纯库埃特流(压力梯度 dp/dx=0 Pa/m)场景下,壁面剪应力 τ_w=μ(U/H)=0.1×(1.5/0.0005)=300 Pa=0.3 kPa(μ=100 mPa·s=0.1 Pa·s)。单位宽度流量 q=U×H/2=1.5×0.0005/2=0.000375 m³/(m·s),雷诺数 Re=ρUD/μ=(860×1.5×0.001)/0.1≒12.9,属层流区。

工程实务注意

  1. 轴承设计中压力梯度通常为正值,同时存在泊肃叶复合流,流量、温升等性能会改变,需与承载条件联动模拟
  2. 粘度温度依赖性:油温每升高 10℃,粘度约下降 5~15%,需配合定常热分析重新计算
  3. 非牛顿流体(高分子添加油等)中,剪应力与剪速率非线性关系 τ=f(du/dy),本模拟器不适用
  4. 微尺度(<100 μm)需进行滑移边界修正,检查努森数(Knudsen number)