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

幂律流体模拟器 — 非牛顿流体的奥斯特瓦尔德-德瓦勒模型

基于奥斯特瓦尔德-德瓦勒幂律 $\tau = K\,\dot{\gamma}^{n}$ 实时计算剪切应力 $\tau$、表观粘度 $\mu_{\text{app}}$ 以及圆管内表观雷诺数。调节稠度系数 $K$、幂指数 $n$、剪切速率 $\dot{\gamma}$ 与流体密度 $\rho$，即可对比剪切变稀（假塑性）、牛顿与剪切增稠（膨胀）流体的流动曲线和圆管内速度剖面形状。

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

表观雷诺数采用固定的直径 $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
—
流体类型
—
表观雷诺数（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 时趋于尖锐。

理论与主要公式

奥斯特瓦尔德-德瓦勒幂律本构方程（$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\,\dot{\gamma}^{\,n-1}$$

圆管层流速度剖面（$R$＝内半径、$u_{\max}$＝中心线速度）：

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

圆管表观雷诺数（直径 $D$、平均速度 $U$、流体密度 $\rho$）：

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

$n < 1$＝假塑性（剪切变稀）／$n = 1$＝牛顿／$n > 1$＝剪切增稠（膨胀）。在本工具的默认值（$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$。

什么是幂律流体（奥斯特瓦尔德-德瓦勒模型）

🙋
流变学教材里反复出现「幂律流体」这个词，到底是什么意思呢？
🎓
简单来说，就是剪切应力 $\tau$ 与剪切速率 $\dot{\gamma}$ 不再是线性关系，而是按 $\tau = K\,\dot{\gamma}^{n}$ 的幂次关系联结的流体。$K$ 是稠度系数，决定粘度量级；$n$ 是幂指数，决定流体性质。$n = 1$ 时退化为牛顿流体（水、空气）；$n < 1$ 为假塑性（剪切变稀），$n > 1$ 为膨胀流体（剪切增稠）。在本工具的默认设置（$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。
🙋
「剪切变稀」具体有什么好处？日常生活中有例子吗？
🎓
油漆、血液、番茄酱、洗发水、酸奶都是剪切变稀流体。慢慢提起刷子时油漆不滴落（低 $\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，中心被「压平」，靠近壁面处梯度陡升，呈塞流（plug-like）形态。$n > 1$ 时指数小于 2，中心更尖，壁面附近剪切减弱。下方画布可直接观察，这一性质对聚合物挤出与动脉血流分析至关重要。

常见问题

因为它在两端都会失效。低 $\dot{\gamma}$ 时 $\mu_{\text{app}} = K \dot{\gamma}^{n-1}$ 趋于无穷（$n < 1$），但真实高分子溶液在足够低剪切速率下会稳定在零剪切粘度 $\mu_0$ 平台。高 $\dot{\gamma}$ 时表观粘度趋近零，无法体现实际存在的无限剪切极限 $\mu_\infty$。在目标工艺的剪切速率窗口（典型 1〜1000 1/s）内拟合幂律是合理的，但若需要跨越两个平台，则要使用 Carreau 模型（$\mu = \mu_\infty + (\mu_0 - \mu_\infty)[1 + (\lambda\dot{\gamma})^2]^{(n-1)/2}$）或 Cross 模型。
本工具显示的表观雷诺数用当地壁面剪切速率下的 $\mu_{\text{app}}$ 进行非牛顿修正，对层流估算足够，但无法准确预测层湍流过渡点。Metzner-Reed 广义雷诺数 $Re_{\text{gen}} = \rho U^{2-n} D^n / [K\,((3n+1)/(4n))^n\,8^{n-1}]$ 由幂律流体的 Hagen-Poiseuille 层流解推导，可像牛顿流体一样以 2100 作为过渡判据。工业 CFD 求解器（Fluent、OpenFOAM）会在内部计算广义雷诺数，本工具显示的值仅作初步评估。
温度升高加剧热运动，破解高分子链的缠结，因此 $K$ 通常遵循 Arrhenius 关系 $K(T) = K_0 \exp(E_a / RT)$，活化能 $E_a$ 典型值为 20〜50 kJ/mol，升温 10 K 时粘度通常下降 20〜30%。幂指数 $n$ 对温度的敏感性较弱（变化幅度约 ±0.05），在常规加工温度范围内可近似为常数。注塑成型仿真中通常配合 $K(T)$ 与独立的冷却模型，以预测流动前沿的凝固。
奥斯特瓦尔德-德瓦勒是两参数模型，不含屈服应力。对于必须先克服屈服应力 $\tau_y$ 才能流动的流体，需要使用 Bingham（$\tau = \tau_y + \mu_p \dot{\gamma}$，屈服后线性）或 Herschel-Bulkley（$\tau = \tau_y + K \dot{\gamma}^n$，屈服后幂律）。新拌混凝土、牙膏、钻井泥浆都具有屈服应力，仅用幂律模型会在低剪切区显著偏离。只有当 $\tau_y$ 可以忽略或工艺始终处于高剪切区时，幂律模型才足够。

实际应用

聚合物与塑料加工：聚乙烯、聚丙烯熔体的典型参数为 $K = 10^3 \sim 10^5$ Pa·sⁿ、$n = 0.3 \sim 0.6$，呈强烈剪切变稀。在注塑模具浇口处，局部剪切速率可达 $10^3 \sim 10^4$ 1/s，表观粘度下降约两个数量级，正是这种特性使得充模成为可能。在本工具中输入 $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 与 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$。钻屑能否被携带至地面（hole-cleaning capability）取决于流动前沿速度剖面：$n$ 越小，中心越平坦，输送效率越高。但泵送功率由高剪切区的表观粘度决定，因此 $K$ 与 $n$ 的平衡是井场的实际优化课题。

常见误解与注意事项

最常见的误解是「$n < 1$ 就意味着任意低剪切速率下都剪切变稀」。真实高分子溶液存在零剪切粘度平台，当 $\dot{\gamma} < 0.01$ 1/s 时粘度稳定在 $\mu_0$。奥斯特瓦尔德-德瓦勒在 $\dot{\gamma} \to 0$ 时数学上给出 $\mu_{\text{app}} \to \infty$，CFD 求解器在低剪切区会发散。OpenFOAM 的 powerLaw 类需要 `nuMin` / `nuMax` 进行裁剪，本工具将最小剪切速率设为 0.1 1/s 也是同样的考虑。请只在目标工艺的 $\dot{\gamma}$ 范围内使用本模型。

第二个误解是「只要测得幂指数 $n$，就能完全描述流体行为」。$n$ 是局部斜率，剪切速率窗口不同时 $n$ 也会不同。例如某聚合物溶液在 1〜10 1/s 范围内 $n = 0.5$，而在 100〜1000 1/s 范围内 $n = 0.7$。正确做法是在目标工艺的 $\dot{\gamma}$ 范围内进行局部拟合，把「目录上的 $n$」直接外推到全范围会导致浇口和流道处的预测严重偏差。覆盖多个数量级的工艺通常需要 Carreau 模型。

最后是「用奥斯特瓦尔德-德瓦勒模型的表观雷诺数 2100 判断层湍流过渡」的误解。幂律流体的过渡雷诺数依赖于 $n$，按 Metzner-Reed 计算，$n = 0.5$ 时过渡值约 2400，$n = 1.5$ 时约 1800。本工具显示的表观 Re（使用当地 $\mu_{\text{app}}$）仅用于教学，实际工程过渡判断应使用 Metzner-Reed 公式或 Ryan-Johnson 修正。点击播放按钮扫描 $n$，可直观看到即使 $K$、$\dot{\gamma}$、$\rho$ 不变，表观 Re 也会跨数量级变化。

相关工具

宾汉塑性流体模拟器 — 具有屈服应力的非牛顿流体 宾汉塑性流体 τ = τ_y + μ_p γ̇ 的屈服应力、塑性粘度、表观粘度、宾汉数与 BP 雷诺数实时计算工具，可视化牙膏、油漆、新鲜混凝土、钻井泥浆等屈服应力流体的流动曲线与圆... 游泳阻力与功率模拟器 计算游泳者所受流体阻力和所需推进功率。调节速度、体型和水温，分析不同泳姿的能量消耗，适用于运动科学和流体力学教育。 射流泵（喷射器）模拟器 CAE技术文章《射流泵（喷射器）模拟器》介绍了一款用于模拟和分析喷射器工作过程的仿真工具。该模拟器支持实时可视化工作原理，用户可通过调节动力流体压力、流量和吸入条件，即时计算并获取关... 滑动轴承模拟器 本文介绍基于短轴承近似雷诺方程的滑动轴承流体润滑模拟器，可实时计算油膜压力分布、最小油膜厚度及Sommerfeld数，为轴承设计与性能评估提供高效CAE分析工具。 伯努利定理仿真器 使用伯努利定理仿真器，在线调节管道入口、喉部及出口直径与流量，实时计算并可视化三截面的流速与压力分布。通过管道截面图与压力柱状图，直观展示文丘里效应，深入理解流体在变径管道中的能量转... Navier-Stokes流动可視化 Navier-Stokes方程流体模拟器：在线实时仿真盖驱动方腔流（Re=100/400/1000）及绕流。通过色阶与流线可视化2D不可压缩流动，是学习计算流体力学（CFD）与流体动...