|
|
一分四卜型岔管:过岔分处总水头异常增大的解决办法
(雷诺应力模型,总水头 H = 位置水头 + 总压/ρg,支管进口到分岔处总水头沿程减小,过了分岔处变大)
一、先定性:总水头变大是"物理上不可能"的
总水头(机械能)H = z + P_total/(ρg),无外部能量输入时沿流线只能单调减小(粘性耗散)。所以过了分岔处总水头变大,100% 是计算/数值问题,不是真实物理。不用怀疑物理模型,直接查数值。
二、最可能的原因(按概率排)
1. 总压的平均方式错了(最常见)
Fluent 里"面平均"有几种,用错就出问题:
算术平均:P_tot_avg = (ΣP_tot)/N —— 在分岔处容易这样算,会高估
质量加权平均:P_tot_mw = Σ(mᵢ·P_totᵢ)/Σmᵢ —— 正确(能量应按质量通量加权)
分岔处速度分布极不均(有的支管快、有的回流),用算术平均算总压,动压项 ½ρv² 会被高估,总水头就"变大了"。
✅ 解法:总水头用质量加权平均算,别用算术平均。
2. 动压里"v"的取法
P_total = P + ½ρ|v|²,|v| 应该用单元真实速度。若你用截面平均速度 |v̄|² 代替 |v|²,在分离/回流区会严重低估或高估动压。分岔处回流区尤其明显。
✅ 解法:P_total 直接用 Fluent 的 total-pressure 变量(内部已正确处理),别自己拼 P + ½ρv²。
3. 回流区污染采样
分岔处有回流、涡旋,在这些区域取总压样本,局部负速度/异常压力会干扰。
✅ 解法:排除回流区(v·n<0 的区域)再统计。
4. RSM 在强分离处的数值不稳定
雷诺应力模型对强曲率/分离敏感,分岔处(强弯曲、回流)容易:压力-应变项预测不准、收敛慢、总压振荡。
✅ 解法:
分岔处网格加密(壁面 y+ 达标)
残差收到 10⁻⁵ 以下,并监测总压是否稳定(不只是看残差)
若 RSM 不稳,先用 k-ω SST 算一遍对比,确认是不是模型引入的
三、具体解决步骤(照做)
1.确认总水头是质量加权平均,不是算术平均
2.P_total 用 Fluent 自带 total-pressure,不自己拼
3.分岔处网格加密,排除回流区统计
4.残差收到 10⁻⁵,监测分岔口总压曲线是否收敛平稳
5.若还异常:k-ω SST 对比 + 检查边界条件(支管出口是否设对)
四、一句话总结
分岔处总水头变大 = 数值问题。先查总压是不是用了质量加权平均、是不是自己拼动压,再查网格和 RSM 收敛。物理上总水头不可能增大。 |
|