2. 流体力学基础回顾

管路系统的一维流动仿真，本质上是在“质量守恒、动量守恒、能量守恒”三条物理定律的约束下，求解沿程压力、流速与温度的变化。本节把三维的流体力学方程沿管路轴向做截面平均，得到工程上可直接使用的简化形式，并给出层流与湍流的判别方法，为后续建模与参数设置打基础。

一、控制体与一维化假设

取管路中长度为 dx 的一段作为控制体，截面为 A(x)。一维化假设包括：同一截面上压力、温度、密度均匀；流速沿轴向变化，截面平均速度为 v；忽略径向速度分量。符号约定统一为：沿流动方向为正方向，x 从入口指向出口，v 与 x 同向时为正，压力 p 为绝对压力，重力势能以水平管为基准。

二、连续性方程：质量守恒

对可压缩流体，连续性方程的微分形式为：

∂ρ/∂t + ∂(ρvA)/∂x = 0

稳态流动时 ∂ρ/∂t = 0，得到 ρvA = 常数，即质量流量 ṁ 沿程不变。若流体不可压缩（ρ 为常数），则简化为：

A₁v₁ = A₂v₂ = Q = 常数

这说明管径缩小则流速升高，管径扩大则流速降低。例如内径 50 mm 的管段流速为 2 m/s，若变径到内径 25 mm，则流速升至 8 m/s。仿真中，连续性方程通常作为流量守恒约束，用于检查各支路流量分配是否合理。

三、动量方程：压力与阻力的关系

沿轴向的一维动量方程可写为：

∂(ρv)/∂t + ∂(ρv²)/∂x = -∂p/∂x - ρg·sinθ - (f/2D)ρv|v|

稳态、水平管（θ = 0）时，惯性项与压力梯度、摩擦项平衡。工程上更常用的是积分形式，即伯努利方程加损失项：

p₁/ρg + v₁²/2g + z₁ = p₂/ρg + v₂²/2g + z₂ + h_f + h_m

其中 h_f 为沿程损失，h_m 为局部损失。沿程损失用达西–魏斯巴赫公式计算：

h_f = f · (L/D) · v²/(2g)

摩擦系数 f 的取值直接依赖流态：层流时 f = 64/Re，湍流时需用 Colebrook 公式或 Moody 图。这正是下一节判别层流与湍流的工程意义所在。

四、能量方程：温度与传热

一维稳态能量方程（忽略动能与势能变化对温度的贡献）可简化为：

ṁ·cp·(T₂ - T₁) = q'·L

其中 q' 为单位管长的热流密度（W/m），cp 为定压比热。若管壁与流体换热，q' = h·πD·(T_wall - T_fluid)。对于绝热管段，q' = 0，则 T₂ = T₁。仿真中，能量方程用于求解沿程温度分布，并与动量方程耦合——因为温度变化会改变黏度 μ，进而影响 Re 和摩擦系数 f。

五、层流与湍流的判别

流态由雷诺数 Re 判别：

Re = ρvD/μ = vD/ν

对于圆管内的流动，工程判别标准为：

以 20 ℃ 水为例，ρ = 998 kg/m³，μ = 1.002×10⁻³ Pa·s。若管内径 D = 0.05 m，流速 v = 0.05 m/s，则 Re = 998×0.05×0.05 / 1.002×10⁻³ ≈ 2490，处于过渡区，仿真中应避免使用单一层流或湍流模型，建议采用低雷诺数湍流模型或直接数值模拟。若流速提高到 0.5 m/s，Re ≈ 24900，属于充分湍流，可用标准 k-ε 模型。

六、三方程在管路仿真中的协同

连续性方程确定流量与流速分布，动量方程确定压力沿程变化，能量方程确定温度分布。三者通过物性参数（密度 ρ、黏度 μ、比热 cp）相互耦合。仿真求解时，通常先假设初始温度场，计算 Re 和 f，求解动量方程得到压力与流量，再代入能量方程更新温度，迭代至收敛。

实际仿真中，建议先根据设计流量和管径估算 Re，判断流态后再选择对应的摩擦系数模型。若 Re 落在过渡区，应避免使用单一层流或湍流模型，可考虑低雷诺数湍流模型或直接数值模拟。对于非圆截面管道，用水力直径 D_h = 4A/P 代替 D 计算 Re。三方程通过物性参数耦合，迭代求解时通常先假设温度场，计算 Re 和 f，求解动量方程得到压力与流量，再代入能量方程更新温度，重复至收敛。