📚 课程目录第2章管路流动基础与压降机理

2. 管路流动基础与压降机理

在 Fluent 中做管路流阻与压降仿真，最容易犯的错误不是网格画得不好，而是把物理机理搞错了：把该用湍流模型算的工况按层流算，或者把沿程阻力与局部阻力的量级判断反了，导致压降结果偏差几倍甚至一个数量级。本节把管路压降的物理基础梳理清楚，为后续建模、边界条件设置和结果判读提供依据。

一、压降的两个来源：沿程阻力与局部阻力

流体在管路中流动，总压降由两部分构成：沿程阻力（摩擦阻力）和局部阻力（形体阻力）。

沿程阻力来自流体与管壁之间的黏性摩擦，以及流体内部的黏性剪切。它沿整段直管均匀分布，与管长成正比。管壁越粗糙、流速越高、黏性越大，沿程损失越大。在长直管段占主导的管路（如长输管线）中，沿程阻力往往是压降的主要贡献。

局部阻力来自管路中的弯头、三通、变径、阀门、进出口等局部构件。这些构件迫使流体改变方向或速度大小，产生旋涡、流动分离和二次流，把机械能不可逆地耗散为热能。局部阻力集中在很短的距离内，与管长无关，只与构件形式和流速有关。在设备密集、弯头阀门多的管路（如泵站、换热器进出口管路）中，局部阻力可能超过沿程阻力。

工程判断经验：直管段总长与管径之比 L/D 很大（如大于 1000）时，沿程阻力通常占主导；L/D 较小而弯头、阀门密集时，必须认真统计局部阻力，不能忽略。

二、层流与湍流的判别：雷诺数

黏性流动有两种基本流态：层流和湍流。层流时流体质点分层有序运动，只有分子尺度的动量交换；湍流时流体质点做随机脉动，存在涡团尺度的动量交换，摩擦和掺混都显著增强。

判别流态的准则数是雷诺数 Re，定义为惯性力与黏性力之比：

Re = ρ · v · D / μ = v · D / ν

其中：
  ρ  流体密度          [kg/m³]
  v  管道截面平均流速   [m/s]
  D  管道内径          [m]
  μ  动力黏度          [Pa·s]
  ν  运动黏度 = μ/ρ    [m²/s]

对于圆管内的充分发展流动，工程上采用如下判别：

Re < 2300：层流，流动稳定，扰动会被黏性耗散掉。
2300 < Re < 4000：过渡区，流动状态不稳定，可能间歇出现湍流斑，仿真结果对扰动和网格敏感。
Re > 4000：湍流，流动充分发展，需采用湍流模型。

举例：20 ℃ 水在 D = 50 mm 圆管中流动，ν ≈ 1.0×10⁻⁶ m²/s。若 v = 0.04 m/s，则 Re = 0.04×0.05/1.0×10⁻⁶ = 2000，属层流；若 v = 0.1 m/s，Re = 5000，已进入湍流。可见同一根管子，流速变化 2.5 倍就可能跨越流态边界。在 Fluent 中，层流用 Laminar 模型，湍流需选 k-ε、k-ω SST 等模型，选错会直接导致压降和速度剖面失真。

三、达西-魏斯巴赫公式

沿程压降用达西-魏斯巴赫（Darcy-Weisbach）公式计算：

Δp_f = f · (L / D) · (ρ · v² / 2)

其中：
  Δp_f  沿程压降        [Pa]
  f     达西摩擦系数     [-]
  L     管段长度        [m]
  D     管道内径        [m]
  ρ     流体密度        [kg/m³]
  v     平均流速        [m/s]
  ρv²/2 动压            [Pa]

该公式把压降表达为动压乘以一个由几何和摩擦决定的系数 f·L/D。层流时 f 有解析解：f = 64/Re，与粗糙度无关。湍流时 f 同时依赖 Re 和相对粗糙度 ε/D，没有简单解析式，需要借助穆迪图或经验公式（如 Colebrook 公式、Swamee-Jain 显式近似）求解。

局部压降用阻力系数法：

Δp_m = K · (ρ · v² / 2)

其中 K 为局部阻力系数，由构件形式决定，例如：
  90° 标准弯头        K ≈ 0.75
  90° 长半径弯头      K ≈ 0.45
  突然扩大（A1→A2）   K ≈ (1 - A1/A2)²
  闸阀全开            K ≈ 0.2
  截止阀全开          K ≈ 6～10

总压降为各段沿程压降与各局部压降之和。注意：局部阻力系数通常以局部构件下游或上游的动压为基准，不同手册的基准截面可能不同，引用时必须核对，否则 K 值会差出数倍。

四、穆迪图与阻力系数体系

穆迪图（Moody Chart）把 f 表示为 Re 和 ε/D 的函数，是湍流沿程阻力的核心工具图。图中可读出几个关键规律：

层流区（Re < 2300）：所有粗糙度曲线重合为 f = 64/Re。
过渡区：曲线分散，f 随 Re 变化剧烈，工程上常取偏保守值。
湍流光滑管区：f 随 Re 增大而减小，粗糙度影响不明显。
完全粗糙区（阻力平方区）：f 只取决于 ε/D，与 Re 无关，压降与 v² 成正比。

在 Fluent 中，壁面粗糙度通过壁面边界条件的 Roughness Height（ε）和 Roughness Constant 设置，等效于穆迪图中的 ε/D。若仿真对象是旧钢管、铸铁管，ε 可取 0.2～0.5 mm；新不锈钢管可取 0.01～0.05 mm。粗糙度设置不当，是压降仿真与实测偏差的常见原因之一。

五、从机理到 Fluent 设置的对应关系

把上述机理映射到 Fluent 操作，可以形成一条清晰的设置链：

算 Re 定流态：根据设计流量、管径和介质物性算出 Re，判断层流还是湍流，据此选择 Viscous 模型。过渡区工况建议做层流与湍流两套对比，观察压降差异。
沿程阻力靠壁面：直管段的沿程阻力由壁面剪切决定，需正确设置壁面粗糙度 ε 和粗糙度常数，并保证近壁网格 y+ 与所选湍流模型匹配（如 k-ε 增强壁面处理要求 y+ 在 1 量级或 30 以上）。
局部阻力靠几何：弯头、阀门、变径的局部损失由几何解析度决定。弯头处网格过粗会抹平分离涡，导致 K 值偏小、压降偏低。建议在弯头内侧和阀门缩口处加密网格。
结果判读：仿真得到的进出口总压差即总压降。可与达西-魏斯巴赫公式加局部阻力系数法的手算值对比，若偏差在 10%～15% 以内，通常认为设置合理；若偏差过大，优先检查流态选择、粗糙度设置和局部几何是否被简化。

常见误区：把整段管路都按层流算，或忽略阀门、弯头的局部阻力，都会使压降被严重低估。另一个误区是直接套用手册 K 值而不核对基准截面——同一阀门，以进口动压为基准和以出口动压为基准，K 值可能相差一倍以上。

掌握本节内容后，面对一条具体管路，应当能够先估算 Re 判断流态，再用达西-魏斯巴赫公式和局部阻力系数法手算一个压降量级，最后在 Fluent 中按对应模型和边界条件复现。手算值不是用来替代仿真，而是用来判断仿真结果是否可信的参照系。