3. 内流道几何特征与流动阻力关联

内流道的压降并不是一个只由流量决定的单一数值,而是几何形状与流动状态共同作用的结果。工程上常用达西—魏斯巴赫公式作为分析骨架:

Δp = (f · L / D) · (ρ · v² / 2) + ΣK · (ρ · v² / 2)

其中:
  f   —— 沿程阻力系数(与Re、相对粗糙度ε/D有关)
  L   —— 直管段长度
  D   —— 水力直径
  ρ   —— 流体密度
  v   —— 截面平均流速
  ΣK  —— 各局部阻力系数之和

这个公式把压降拆成两部分:沿程摩擦损失和局部形状损失。理解几何特征对压降的影响,本质上就是理解几何如何改变 f、L/D、v 和 ΣK 这四个量。下面按几何因素逐一展开。

管径:四次方级别的敏感参数

管径的影响远比其他几何参数剧烈。在质量流量 ṁ 恒定时,v = ṁ /(ρ·A),而 A = πD²/4,因此动压项 ρv²/2 与 D⁻⁴ 成正比;同时 L/D 又与 D⁻¹ 成正比。两者叠加后,湍流区压降近似与 D⁻⁵ 相关。也就是说,管径缩小 10%,压降可能上升约 60%~70%。这是内流道设计中“宁粗勿细”的根本原因。

但管径并非越大越好。管径增大后流速下降,若流体含颗粒或易结垢,低流速会导致沉积;同时大管径带来体积、重量和成本上升。工程上通常先由允许压降反推最小可行管径,再校核流速是否落在合理区间(液体一般 1~3 m/s,气体 10~20 m/s,含固介质适当降低)。

长度:最直观的线性关系

沿程损失与长度成正比,这一点没有非线性陷阱。真正需要注意的是“有效长度”的统计口径:弯头、阀门、进出口等局部件在工程估算中常折算成当量长度 L_eq = K·D/f,与直管长度相加得到总有效长度。若建模时只画了直管而漏掉这些当量长度,压降会被系统性低估。

弯曲:曲率半径决定局部损失

弯头的局部阻力系数 K 主要受两个参数控制:弯曲半径与管径之比 R/D,以及弯转角。以 90° 弯头为例,R/D = 1(急弯)时 K 约为 0.9,R/D = 2 时降至约 0.3,R/D = 5 时可低至 0.15 左右。也就是说,把急弯改成缓弯,单个弯头的损失可以降到原来的三分之一以下。

弯头还会在出口下游形成二次流和分离区,扰动可持续 10~20 倍管径的距离。如果两个弯头相距很近(小于 5D),它们的损失不能简单相加,而应视为组合件整体评估。CFD 分析中,这类“弯头—弯头”“弯头—突扩”耦合区域往往是压降预测偏差最大的位置。

突扩与突缩:不可逆损失的重灾区

突扩的局部阻力系数可用 Borda-Carnot 关系估算:K_exp = (1 - A₁/A₂)²。当面积比很小时,K 趋近于 1,意味着几乎全部动压被耗散。突缩则分两部分:缩流断面处的加速损失和再附着的混合损失,工程上常用 K_con ≈ 0.5·(1 - A₂/A₁) 估算。

关键认识是:突扩突缩的损失主要来自流动分离和涡旋耗散,属于不可逆损失,无法通过下游扩压回收。因此在内流道设计中,应尽量用渐扩/渐缩段替代突变,渐扩角控制在 7° 以内、渐缩角控制在 15°~30° 以内,可显著降低 K 值。

截面形状:水力直径是统一尺度

非圆截面用水力直径 D_h = 4A/P(A 为流通面积,P 为湿周)来统一描述。相同流通面积下,形状越“扁”,湿周越大,D_h 越小,沿程阻力越高。例如正方形截面 D_h = a,而同等面积的扁矩形(宽高比 10:1)D_h 只有约 0.55a,压降明显更高。

此外,尖角截面(如矩形、三角形)的角隅区存在低速涡,会额外增加损失;圆角或椭圆截面则更接近圆管特性。层流条件下,不同截面的 f·Re 常数也不同:圆管为 64,正方形约 57,扁矩形可低至 40 以下,说明形状对层流阻力的影响不可忽略。

表面粗糙度:从水力光滑到完全粗糙

粗糙度通过相对粗糙度 ε/D 影响摩擦系数。在层流区,粗糙峰被黏性底层覆盖,粗糙度几乎不起作用;进入湍流后,随着 Re 增大,粗糙度影响逐渐显现。当 Re·√(f/8)·(ε/D) > 70 时进入完全粗糙区,此时 f 只取决于 ε/D,与 Re 无关。

工程上常见的等效砂粒粗糙度参考值:拉制钢管约 0.0015~0.002 mm,商用钢管约 0.045 mm,铸铁约 0.26 mm,混凝土约 0.3~3 mm。对于小管径流道,即使绝对粗糙度不大,ε/D 也可能相当可观。例如 D = 5 mm 的商用钢管,ε/D ≈ 0.009,已接近中等粗糙水平。

在 CFD 建模中,若采用壁面函数法,粗糙度通过等效砂粒高度参数输入;若采用解析到壁面的网格,则需保证第一层网格 y⁺ ≈ 1,否则粗糙度信息无法被正确捕捉。

几何因素影响规律总览

内流道几何特征 → 压降影响规律 管径 D Δp ∝ D⁻⁵(湍流,定质量流量) D↓10% → Δp↑约60%~70% 敏感度最高,优先控制 长度 L Δp ∝ L(线性) 需计入局部件当量长度 L_eq 漏算当量长度 → 系统性低估 弯曲 R/D 90°弯头:R/D=1→K≈0.9 R/D=2→K≈0.3,R/D=5→K≈0.15 缓弯可降至急弯的1/3以下 突扩 / 突缩 突扩:K=(1−A₁/A₂)² 突缩:K≈0.5·(1−A₂/A₁) 不可逆损失,渐扩角≤7° 截面形状 D_h = 4A/P,湿周越大阻力越高 层流 f·Re:圆管64,方管57 尖角涡 → 额外损失 表面粗糙度 ε/D 层流:几乎无影响 湍流:Re↑ → 影响↑ → 完全粗糙区 小管径 ε/D 放大效应显著 压降敏感度排序(工程经验) 管径 D > 突扩突缩 > 弯曲 R/D > 粗糙度 ε/D > 截面形状 > 长度 L(线性可控) CFD 分析中的几何处理要点 管径:确认水力直径定义,非圆截面用 D_h = 4A/P 长度:统计所有直管段,局部件折算当量长度 L_eq = K·D/f 弯曲:相邻弯头间距 < 5D 时按组合件整体评估,不单独叠加 粗糙度:壁面函数法输入等效砂粒高度;解析到壁面需 y⁺ ≈ 1

工程经验与建模建议

综合以上规律,内流道压降分析可以遵循一条清晰的优先级:先控管径,再控突变,后控弯曲,最后校核粗糙度。管径决定了压降的量级,突扩突缩和急弯决定了局部损失的集中程度,粗糙度则在小管径或高 Re 工况下成为不可忽略的修正项。

在 CFD 建模中,几何清理阶段就应识别出所有 R/D < 2 的急弯、面积突变比超过 2:1 的突扩突缩、以及相邻距离小于 5D 的连续局部件。这些位置应加密网格并作为压降监测的重点区域。对于长直管段,若沿程摩擦占主导,可采用一维管网计算快速估算,再用三维 CFD 校核局部件附近的流动细节,这样能在保证精度的同时显著降低计算成本。

需要特别注意的是:当多个局部件距离很近时,它们的阻力系数不能简单相加。例如弯头紧接突扩,下游的分离区会被上游二次流强化,实际 K 值可能比单独件之和高出 20%~50%。这类耦合效应只能通过 CFD 或专门实验数据来捕捉。

最后,几何特征对压降的影响并非孤立作用。管径变化会同时改变 Re 和 ε/D,弯曲会改变下游速度分布进而影响后续直管段的摩擦系数,粗糙度在低 Re 下几乎无效但在高 Re 下可能使 f 翻倍。因此,实际分析中应把几何参数作为一组耦合变量来考虑,而不是逐个单独调整。建立参数化模型、批量扫描关键几何变量、观察压降响应曲线,是形成工程直觉和优化设计的有效路径。