› 幂律流体 模拟器 返回
非牛顿流体模拟器

幂律流体 模拟器 — 非牛顿流体的 Ostwald-de Waele 模型

基于 Ostwald-de Waele 幂律 $\tau = K\,\dot{\gamma}^{n}$，实时计算剪切应力 $\tau$、表观粘度 $\mu_{\text{app}}$、圆管内见观 Reynolds 数。通过滑块改变稠度系数 $K$、幂指数 $n$、剪切速率 $\dot{\gamma}$、流体密度 $\rho$，对比可视化 shear-thinning（假塑性）／Newtonian／shear-thickening（膨胀性）的流动曲线和速度分布。

参数设置
稠度系数 K
Pa·sⁿ
幂指数 n
剪切速率 γ̇
1/s
流体密度 ρ
kg/m³

圆管内见观 Re 假设直径 $D = 10$ mm、平均速度 $U = 1$ m/s。$\mu_{\text{app}} = K\,\dot{\gamma}^{\,n-1}$、$Re_{\text{app}} = \rho U D / \mu_{\text{app}}$。

计算结果
—
剪切应力 τ
—
表观粘度 μ_app
—
流体类型
—
见观 Re（D=10mm, U=1m/s）
流动曲线 τ vs γ̇（log-log）

横轴 = $\log_{10}\dot{\gamma}$（1/s）／纵轴 = $\log_{10}\tau$（Pa）／多个幂指数 $n$（0.3/0.5/0.7/1.0/1.3/1.7）的 $\tau(\dot{\gamma})$ 曲线叠加／黄● = 当前 $(\dot{\gamma},\tau)$。斜率等于 $n$，截距等于 $\log K$。

圆管内速度分布 u(r)

轴对称圆管横截面（半径 $R$）／横轴 = 径向位置 $r/R$／纵轴 = 正规化速度 $u/u_{\max}$／$u(r) = u_{\max}\left[1 - (r/R)^{(n+1)/n}\right]$／n=1 时为抛物线，n < 1 时变平，n > 1 时变尖。

理论和主要公式

Ostwald-de Waele 幂律（$K$ 是稠度系数 Pa·sⁿ，$n$ 是幂指数，$\dot{\gamma}$ 是剪切速率 1/s）：

$$\tau = K\,\dot{\gamma}^{\,n}$$

表观粘度为 $\mu_{\text{app}} = \tau/\dot{\gamma}$，所以对于牛顿流体（n=1），$\mu_{\text{app}} = K$ 为常数；其他情况下依赖于剪切速率：

$$\mu_{\text{app}} = K\,\dot{\gamma}^{\,n-1}$$

圆管层流的速度分布（$R$ = 内半径，$u_{\max}$ = 中心线速度）：

$$u(r) = u_{\max}\left[1 - \left(\dfrac{r}{R}\right)^{(n+1)/n}\right]$$

圆管内见观 Reynolds 数（直径 $D$、平均速度 $U$、流体密度 $\rho$）：

$$Re_{\text{app}} = \dfrac{\rho\,U\,D}{\mu_{\text{app}}}$$

$n < 1$ = shear-thinning（假塑性）／$n = 1$ = Newtonian／$n > 1$ = shear-thickening（膨胀性）。本工具默认值（$K = 1$ Pa·sⁿ、$n = 0.5$、$\dot{\gamma} = 10$ 1/s）时，$\tau = \sqrt{10} \approx 3.16$ Pa、$\mu_{\text{app}} = 10^{-0.5} \approx 0.316$ Pa·s、$Re_{\text{app}} \approx 31.6$。

幂律流体（Ostwald-de Waele 模型）概述

🙋
流变学教科书中总是出现的「幂律流体」，到底是什么流体？
🎓
简单来说，就是剪切应力 $\tau$ 和剪切速率 $\dot{\gamma}$ 的关系不是直线，而是通过 $\tau = K\,\dot{\gamma}^{n}$ 这样的幂律关联的流体。$K$ 是稠度系数，决定粘度的大小；$n$ 是幂指数，决定流体的性质。当 $n$ 恰好等于 1 时就是牛顿流体（水、空气）。$n < 1$ 时是假塑性（shear-thinning），$n > 1$ 时是膨胀性（shear-thickening）。在本工具的默认值（$K = 1$ Pa·sⁿ、$n = 0.5$、$\dot{\gamma} = 10$ 1/s）下，$\tau = \sqrt{10} \approx 3.16$ Pa，表观粘度 $\mu_{\text{app}} \approx 0.316$ Pa·s，在「计算结果」中可以看到。
🙋
「shear-thinning」有什么实际用处吗？有哪些常见的流体？
🎓
油漆、血液、番茄酱、洗发水、酸奶都是假塑性流体。用刷子涂油漆时，缓慢提起时油漆不会滴下来（低 $\dot{\gamma}$ 下粘度高），用力刷动时则光滑地铺开（高 $\dot{\gamma}$ 下粘度低）。没有这种性质，油漆就涂不了，血液也不能顺畅流动。试试用滑块把 $n$ 调到 0.3，看看流动曲线的斜率变得缓和（也就是说，即使剪切速率增加，应力也上升不快，表观粘度下降），这就是假塑性的表现。
🙋
那 $n > 1$ 的膨胀性流体呢？
🎓
典型代表是淀粉水溶液（oobleck）。慢慢把手指放进去会沉下去，但狠狠一击就会像固体一样反弹。这是因为当 $\dot{\gamma}$ 增大时，胶体粒子形成「水合簇」，表观粘度急剧上升。试试用滑块把 $n$ 调到 1.7，会看到流动曲线陡峭地上升，这说明剪切速率越高，粘度越大。在防弹背心的液体装甲材料研究中也有应用。
🙋
为什么速度分布的形状会随 n 大幅变化？
🎓
圆管层流的解为 $u(r) = u_{\max}\left[1 - (r/R)^{(n+1)/n}\right]$。当 $n = 1$ 时指数等于 2，也就是抛物线分布（Hagen-Poiseuille 流）。当 $n < 1$ 时指数大于 2，中心部分变平，壁近处陡峭变化，呈现「塞子流动」的样子。当 $n > 1$ 时指数小于 2，中心变尖，壁近处的剪切变弱。看下面的画布，改变 $n$ 可以直观看到这些变化，这对理解挤压成形和血管内血流至关重要。

常见问题

因为 Ostwald-de Waele 在两端都会物理破裂。在低 $\dot{\gamma}$ 下，$\mu_{\text{app}} = K \dot{\gamma}^{n-1}$ 会趋于无穷（当 $n < 1$ 时），无法再现现实高分子溶液中存在的零剪切粘度 $\mu_0$。在高 $\dot{\gamma}$ 下则相反，$\mu_{\text{app}}$ 接近 0，也无法再现极限粘度 $\mu_\infty$。实际工艺的剪切速率范围（通常 1～1000 1/s）内才是近似有效。要想在广泛范围内覆盖，需要用 Carreau 模型（$\mu = \mu_\infty + (\mu_0 - \mu_\infty)[1 + (\lambda\dot{\gamma})^2]^{(n-1)/2}$）或 Cross 模型。
本工具显示的见观 Re 是基于壁面剪切速率 $\dot{\gamma}_w$ 处 $\mu_{\text{app}}$ 的简单定义，对于层流区的粗估计足够，但无法准确预测乱流过渡点。Metzner-Reed 的广义 Reynolds 数 $Re_{\text{gen}} = \rho U^{2-n} D^n / [K\,((3n+1)/(4n))^n\,8^{n-1}]$ 是从幂律流体的 Hagen-Poiseuille 解推导出来的，过渡点可以像牛顿流体那样用 2100 左右的值。在实际 CFD 工程中，Fluent 或 OpenFOAM 内部计算的是广义 Re，所以本工具的值只能作为初期检验用。
温度上升激活热运动，高分子链的缠结得以解除，因此 $K$ 大致按 Arrhenius 型 $K(T) = K_0 \exp(E_a / RT)$ 下降。活性化能 $E_a$ 通常在 20～50 kJ/mol，温度上升 10°C 粘度一般下降 2～3 成。相比之下，幂指数 $n$ 的温度依赖性较弱（±0.05 左右），在室温到成形温度的范围内基本可视为常数。在射出成形模拟中，通常结合温度相关的 $K(T)$ 和独立的冷却模型，预测流动前沿的凝固。
Ostwald-de Waele 是不含屈服应力的 2 参数模型。屈服应力 $\tau_y$（应力超过这个值才能流动）的流体需要用 Bingham（$\tau = \tau_y + \mu_p \dot{\gamma}$，线性）或 Herschel-Bulkley（$\tau = \tau_y + K \dot{\gamma}^n$，幂律）来描述。新鲜混凝土、牙膏、钻井液等都有屈服应力，用 Ostwald-de Waele 就不适用，特别是底端行为会严重偏离。只有当屈服应力可以忽略，或者 $\dot{\gamma}$ 始终很高时，Ostwald-de Waele 才足够。

实际应用

高分子和塑料成形：聚乙烯、聚丙烯熔体通常表现为 $K = 10^3 \sim 10^5$ Pa·sⁿ、$n = 0.3 \sim 0.6$ 的强假塑性。在射出成形时，金型浇口处的剪切速率可达 $10^3 \sim 10^4$ 1/s，表观粘度下降 100 倍，从而保证流动性。在本工具中输入 $K = 100$、$n = 0.4$、$\dot{\gamma} = 1000$ 1/s，可以直观看到表观粘度的数量级下降。流动分析软件（Moldflow、Moldex3D）内部就是用这个模型和参数。

食品和化妆品工业：番茄酱 $n \approx 0.3 \sim 0.5$、酸奶 $n \approx 0.5 \sim 0.7$、洗发水 $n \approx 0.5 \sim 0.8$。从容器中的可挤出性、用刷子的涂布感、口腔内的口感（食品感官评价）都由幂指数决定。产品设计时需要用流变仪测量 $K, n$，通过增粘剂、乳化剂的调整达到目标口感。在本工具中对比 $n = 0.3$ 和 $n = 0.8$ 时的速度分布，平坦化和尖锐化的差异很容易看出。

血液流变学（生物力学）：血液因为红细胞的凝聚和变形，呈现弱假塑性，典型 $K \approx 0.01 \sim 0.1$ Pa·sⁿ、$n \approx 0.7 \sim 0.9$。在大动脉的高剪切域粘度接近水，在毛细血管的低剪切域粘度增加。动脉瘤、血管狭窄的血流 CFD 标准上就用这个模型，通过壁面剪切应力分布预测血栓形成风险。本工具的圆管分布功能可以对比 $n = 0.8$ 和 $n = 1.0$（牛顿近似）的差异。

钻井液和地热流体：石油和地热井的钻探中使用的钻井液是膨润土悬浊液，典型 $n \approx 0.4 \sim 0.7$。将钻屑（cuttings）从井底运到地表的能力（孔洞清理能力）依赖于流动前沿的速度分布，$n$ 越小越平坦，输运效率越高。另一方面，泵的驱动功率由高 $\dot{\gamma}$ 域的表观粘度决定，$K$ 和 $n$ 的平衡成为现场优化的关键指标。

常见误区和注意事项

最常见的误区是：「$n < 1$ 说明是假塑性，那在任何低剪切速率下粘度都会下降」。实际上真实的高分子溶液有「零剪切粘度平台」，在 $\dot{\gamma} < 0.01$ 1/s 左右粘度会保持在常数 $\mu_0$。Ostwald-de Waele 在数学上当 $\dot{\gamma} \to 0$ 时 $\mu_{\text{app}} \to \infty$，所以用它在低剪切速率下做 CFD 分析会导致求解器发散。OpenFOAM 的 powerLaw 类中需要用 nuMin / nuMax 来裁剪，本工具之所以设定 $\dot{\gamma} = 0.1$ 1/s 的下限也是这个原因。只在实际工艺的 $\dot{\gamma}$ 范围内才有效。

第二个常见错误是：「测到了幂指数 $n$ 就能完全决定流体行为」。实际上 $n$ 是局部斜率，改变剪切速率范围，$n$ 也会变。某聚合物溶液在 1～10 1/s 可能 $n = 0.5$，但在 100～1000 1/s 就变成 $n = 0.7$。Ostwald-de Waele 应该在目标工艺的剪切速率域内做局部拟合，不能拿「样本表中的 $n$」直接用在全过程。要涵盖多个 $\dot{\gamma}$ 域就需要 Carreau 模型做标准。

最后一个常见误解是：「用见观 Re 的 2100 来判断乱流过渡」。幂律流体的乱流过渡点随 $n$ 变化，在 Metzner-Reed 广义 Re 中 $n = 0.5$ 约 2400，$n = 1.5$ 约 1800。本工具的见观 Re（用局部 $\mu_{\text{app}$）只是为教学目的呈现，实际过渡判定需要用 Metzner-Reed 公式或 Ryan-Johnson 修正。试试扫描 $n$（按播放键）时，即使 $K$、$\dot{\gamma}$、$\rho$ 不变，见观 Re 也会整个数量级摇摆，就会直观理解这一点。

关联工具

Bingham 塑性流体 模拟器 — 具有屈服应力的非牛顿流体 Bingham 塑性流体 τ = τ_y + μ_p γ̇ 的屈服应力、塑性粘度、表观粘度、Bingham 数、BP Reynolds 数的实时计算。牙膏、油漆、新鲜混凝土、钻井液的... 游泳阻力和推进力 模拟器 游泳中的流体阻力（F=½ρCdAv²）和所需推进力的实时计算。速度、体型、泳法的不同阻力系数可视化。 喷射泵（喷射器）模拟器 喷射泵（喷射器）的流量、压力、引射比的简单计算模拟器。只需调整驱动流体的条件，动作就会实时可视化，混合效率和压缩比瞬间... 滑动轴承 模拟器 滑动轴承的流体润滑模拟在线运行。Sommerfeld 数或最小油膜厚度、承载能力、摩擦系数的实时计算。油膜压力分布可视化，支持轴承设计的高效评检... 伯努利定理 在变截面管路流动中直观理解伯努利定理的模拟器。实时动画可视化压力和速度变化，学习能量守恒原则。流体力学学习和教育中...

使用指南

  1. 输入流体常数 K（Pa·sⁿ）和流动指数 n（无量纲）。番茄酱约为 K=100、n=0.33；增稠剂混聚合物溶液约为 K=1.2、n=1.5
  2. 设定剪切速率 γ̇（s⁻¹）和流体密度 ρ（kg/m³）。口腔内流动约 γ̇=100～1000 s⁻¹，血液 ρ=1060 kg/m³ 可作参考
  3. 实时更新显示遵循幂律模型 τ=Kγ̇ⁿ 的表观粘度 μₐₚₚ、剪切应力 τ、流动曲线，以及 D=10mm 管路 Re 见观值

具体计算例

发油（K=50 Pa·sⁿ、n=0.5）在 γ̇=500 s⁻¹ 下剪切时：τ=50×500⁰·⁵=1118 Pa，μₐₚₚ=1118÷500=2.24 Pa·s。同流体密度 ρ=900 kg/m³、特征速度 U=1 m/s、D=10mm 圆管流时，见观 Re=(900×1×0.01)÷2.24≈4.0（层流）；而牛顿流体（n=1、μ=2.24 Pa·s）的 Re≈400，差异巨大

工程应用注意

  1. n<1（假塑性）的流体在高速泵送和金型填充时粘度大幅下降。医用甘油制剂（n=0.8）用注射针推动时，低 γ̇ 粘度高、驱动扭矩增加，须提前确认
  2. n>1（膨胀性）的番茄酱或玉米淀粉悬浊液（n=1.2～1.8），快速加速时粘度反而上升，容器急倾时流动反而困难，应该缓缓倾斜为正解
  3. 见观 Re 重算时，不仅液体比重会影响，n 值对有效粘度的改变也直接影响乱流过渡（Re≈2300 标准不适用）。n<0.5 的食品反复加热膏状物即使见观 Re<1，局部加速也会诱发乱流混合，配管设计前需确认全 γ̇ 范围的流动曲线