非牛顿流体 — 故障排除指南

非牛顿流体 — 故障排除指南

非牛顿流体计算的典型故障

1. 粘度发散导致计算停止

原因：在幂律模型中，当 $\dot{\gamma} \to 0$ 时 $\eta \to \infty$（$n < 1$ 的情况）。在流动停滞区域发生。

对策：

• 设置粘度上下限：Fluent中指定 Minimum Viscosity 和 Maximum Viscosity。通常设置 $\eta_{\text{min}} = \eta_\infty$、$\eta_{\text{max}} = \eta_0$
• 改用Carreau模型：零剪切粘度 $\eta_0$ 会自动设置上限
• 添加背景粘度：在幂律粘度上添加微小的常数粘度

2. 屈服应力模型的收敛困难

原因：屈服面位置（$\tau = \tau_y$）在各迭代间变化，"固体⇔液体"的转换不稳定

对策：

• 分阶段增加Papanastasiou正则化参数 $m$：从 $m = 10 \to 100 \to 1000$ 逐步提高
• Bi-viscosity模型：假设屈服前具有非常高的粘度（$10^3 \sim 10^5\,\text{Pa}\cdot\text{s}$）
• 降低松弛系数：设置粘度的 under-relaxation 为0.3-0.5

3. 剪切速率的计算精度不足

原因：网格粗糙，壁面的剪切速率 $\dot{\gamma}$ 计算不准确

对策：

• 充分细化壁面第一层（管内流中 $R/\Delta r > 20$ 为目标）
• 使用二阶或更高阶梯度方案
• 在壁面布置垂直于壁面的棱柱层

4. 温度依赖粘度的热失控

原因：高粘度流体（高分子熔体等）的高速剪切中粘性散逸占主导地位

对策：

• 激活能量方程的Viscoeus Dissipation项
• 充分减小时间步长
• 确认粘度的温度依赖系数（Arrhenius型：$\eta = \eta_0 \exp[E_a/(RT)]$）

调试流程图

1. 首先用牛顿流体（$\eta = \eta_0$）获得定常解

2. 应用幂律，从 $n = 0.9$ 开始逐步降至目标值

3. 若存在屈服应力，从 $\tau_y = 0$ 逐步增加

4. 各阶段使用前一步骤的解作为初始值

5. 温度依赖性在最后添加