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

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

针对平行板(上板速度 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 取正值——即与上板运动方向相反的逆向压力梯度——泊肃叶分量推动方向与库埃特拖动相反。下板附近这个反向推力胜过上板的拖动,速度变为负值(逆流)。模拟器用红色箭头和红色填充区显示。这在轴承和齿轮的间隙流动中确实可见,潤滑膜内的油会部分反向流动,称为"逆流袋"。点击 Play 让 dp/dx 在 −2000 到 +2000 间扫描,可以一次看清楚:纯库埃特直线→泊肃叶贡献使其鼓凸→出现逆流的过渡。
🙋
壁面剪切应力显示 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 以上转捩为湍流——和你在圆管中熟悉的 2300 不一样。默认 Re=2000 正好处于这个边界附近。把 μ 增大 10 倍(轻质油 0.01 Pa·s)Re 降到 200(确实是层流),把间隙 h 增大 10 倍(20 mm)Re 跳到 20000(完全湍流),此时层流解析解就只能作为参考值。

常见问题

按物理驱动机制区分。库埃特流由"边界运动(壁面速度)"驱动,适用于轴承油膜、齿轮间隙流、旋转粘度计等"运动壁与静止壁夹住粘性流体"的场景。泊肃叶流由"压力差"驱动,适用于管道流量计算、点滴、微流体芯片等"从上游推送出来的流动"。实际中两者经常并存,例如发动机轴承同时存在"轴旋转引起的库埃特驱动"和"外部泵供给潤滑油的泊肃叶驱动",按叠加解分析。本工具显示的正是这一合成解。
在纯库埃特流(dp/dx = 0)中,平行板间的流体既无体积力也无压力梯度,运动方程退化为极简的 $\mu\,d^2u/dy^2 = 0$(二阶导数=0)。在边界条件 u(0)=0、u(h)=U 下,二阶导数为 0 → 一阶导数(剪切率)为常数 → 速度是 y 的线性函数,得到 $u(y) = U y/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 在试样内全程均匀,从而精确测量非牛顿流体(触变性流体、宾汉流体、幂律流体等)的剪切率依赖黏度 $\mu(\dot\gamma)$。化妆品、食品、涂料、血液等行业广泛使用。

实际应用

潤滑工程与摩擦学:滑动轴承、滚动轴承、齿轮间隙的潤滑油流动正是库埃特流+泊肃叶流的合成。流体潤滑理论的核心方程——雷诺方程,就是从本工具显示的一般解 $u(y) = U y/h - (1/(2\mu))(dp/dx)\,y(h-y)$ 出发推导的。摩擦损失 W = τ·U·面积、发热量、最小油膜厚度、轴承承载能力等都由该式计算。从发动机轴承、智能手机精密电机到高速列车曲轴的油膜,本式是所有潤滑膜回转机械分析的中心。

流变学与黏度测量:库埃特型圆筒流变仪和圆锥-平板流变仪都基于库埃特流原理工作。通过让 du/dy 在样品内全程均匀,既可精确测量牛顿流体的黏度,也可测量非牛顿流体(触变性流体、宾汉流体、幂律流体等)的剪切率依赖黏度函数 $\mu(\dot\gamma)$。化妆品、食品、涂料、印刷油墨、血液、聚合物溶液等工业品质控制中不可或缺。

微流体与芯片实验室(LOC):微流体芯片中间隙 h ≈ 10〜100 μm 的薄通道内,电渗流在壁面附近生成薄薄的库埃特状层,泵驱动的泊肃叶流作用于主流;二者并存。其合成剖面用于设计细胞分选、蛋白质分离、PCR 装置内的流路。本工具的"逆流出现条件"也用于刻意构造逆流域以促进混合的微通道设计。

地球动力学与冰川流动:冰川流动、地幔对流、地壳板块运动在长时间尺度下都可视为粘性流。冰川的边界条件不同(底部固定基岩、顶面自由),但底部附近的速度剖面与库埃特剖面定性相似。黏度 μ ≈ 10¹⁵ Pa·s 极大且流速极慢(每年数米),但本工具同样的运动方程在地球物理中亦在起作用。

常见误解与注意事项

最常见的误解是"库埃特流总是线性剖面"。剖面线性只在 dp/dx = 0 的纯库埃特流时成立。实际的轴承、齿轮、流变仪中几乎总有压力梯度:例如轴承中需要在油膜内产生压力分布以承载负荷,这就贡献了泊肃叶分量。在本模拟器中拖动 dp/dx 滑块,可以体验剖面从直线鼓凸为抛物线、再出现逆流的过程。标准的轴承设计使用包含该压力分布的雷诺方程求解。

第二个易错点是"以为平行板的临界 Re 与圆管同为 2300"。平面剪切驱动流(纯库埃特)的著名结果是"对所有 Re 线性稳定",但实验中 Re ≈ 360〜1500 即可转捩为湍流(亚临界转捩)。本工具始终显示层流解析解,因此 Re 较高时与实际流动会偏离。同时要记住不同流动类型的临界值不同:圆管约 2300、平板上的边界层约 5×10⁵ 等。

最后请注意"将单位宽度流量 Q'(m²/s)与普通体积流量 Q(m³/s)混淆"。平行板流被处理为流路在纸面深度方向无限延展的二维问题,所以流量为"深度方向单位宽度",单位为 m²/s。乘以实际通道宽度 W [m] 即可得到三维体积流量 Q = Q'·W [m³/s]。本工具同时显示 Q'(单位宽度流量)与 V_avg(平均流速),按需要选用。CAE 分析中混淆 2D 假设与 3D 换算会引起 W 倍的流量误差,请务必注意。