2. 流体力学基础与压降机理
内流道压降分析的本质，是把三维黏性流动的控制方程在管道几何与工况约束下逐步简化，最终得到可计算、可校核的压降表达式。理解这一简化链条，才能在CFD建模时正确选择物理模型、设置边界条件，并判断仿真结果是否合理。
一、三大守恒方程及其在内流道中的简化
流体流动受质量、动量、能量三大守恒定律支配。在CFD中，它们分别对应连续性方程、动量方程（Navier-Stokes方程）和能量方程。
连续性方程表达质量守恒。三维瞬态形式为 ∂ρ/∂t + ∇·(ρu) = 0。对于内流道中常见的定常、不可压缩流动（液体或低速气体，马赫数低于0.3），密度为常数，方程简化为 ∇·u = 0，即速度场散度为零。工程上更常用的是其一维积分形式：ṁ = ρAV = 常数。这意味着管道截面积缩小处流速必然升高，截面积扩大处流速降低——这是理解局部损失和变径段流动的起点。
动量方程表达牛顿第二定律。不可压缩牛顿流体的三维形式为 ρ(∂u/∂t + u·∇u) = −∇p + μ∇²u + f。方程左侧是惯性项，右侧依次为压力梯度力、黏性力和体积力。对于充分发展的管内流动，惯性项中对流加速度为零，方程退化为压力梯度与黏性力之间的平衡：dp/dx = μ∇²u。这一平衡直接决定了沿程损失的物理来源——压力能持续转化为克服黏性摩擦的耗散。
能量方程表达能量守恒。在内流道中，若忽略传热和做功，且流体不可压缩，能量方程可简化为机械能守恒的伯努利形式：p₁/ρg + V₁²/2g + z₁ = p₂/ρg + V₂²/2g + hL。其中 hL 为水头损失，代表机械能不可逆地转化为内能。这正是压降分析的核心：总压降 Δp = ρg·hL，由沿程损失和局部损失两部分组成。
二、沿程损失与局部损失的产生机理
沿程损失源于流体与管壁之间以及流体内部各层之间的黏性摩擦。在充分发展的管流中，壁面处流速为零（无滑移条件），管轴心处流速最大，形成速度梯度。黏性切应力 τ = μ(du/dy) 持续做功，将机械能转化为热能。沿程损失用达西-魏斯巴赫公式计算：h_f = λ(L/d)(V²/2g)，其中 λ 为沿程阻力系数，L 为管长，d 为管径，V 为截面平均流速。λ 的取值取决于流态：层流时 λ = 64/Re，仅与雷诺数有关；湍流时 λ 由雷诺数和相对粗糙度 ε/d 共同决定，常用Colebrook方程或Moody图求解。
局部损失产生于流道几何突变处——弯头、渐缩/渐扩段、阀门、三通、孔板等。其机理是：流体经过这些位置时，速度大小或方向发生剧烈变化，产生边界层分离和漩涡区。漩涡消耗机械能，且这种耗散不可逆。局部损失用 h_j = ζ(V²/2g) 计算，ζ 为局部阻力系数，通常由实验测定。例如，90°标准弯头 ζ ≈ 0.9，全开闸阀 ζ ≈ 0.2，突扩管 ζ = (1−A₁/A₂)²。在CFD中，局部损失区域往往需要更密的网格来捕捉分离和回流。
三、层流与湍流的判别
流态判别依据雷诺数 Re = ρVd/μ = Vd/ν。对于圆管内的充分发展流动，工程上采用以下判据：
- 层流：Re < 2300。流体分层有序流动，速度剖面为抛物线形，沿程阻力系数 λ = 64/Re。层流中黏性力主导，动量方程中的对流项可忽略，CFD可采用层流模型直接求解。
- 过渡区：2300 ≤ Re ≤ 4000。流动不稳定，可能间歇出现湍流斑。CFD中此区域最难模拟，通常需要瞬态计算或采用低雷诺数湍流模型。
- 湍流：Re > 4000。流体出现随机脉动，速度剖面趋于平坦（1/7次方律），沿程阻力系数显著增大且与粗糙度相关。CFD需选用合适的湍流模型（如k-ε、k-ω SST、Spalart-Allmaras等）。
在内流道CFD分析中，判别流态是选择物理模型的第一步。若误将湍流按层流计算，压降会被严重低估；反之则可能高估压降并浪费计算资源。实际工程中，可先根据设计流量和管道尺寸估算Re，再决定建模策略。
四、从方程到CFD：简化链条的工程意义
将三大方程在内流道中简化的过程，本质上是在回答三个问题：流体如何保持质量守恒（连续性）、压力如何驱动流动并克服阻力（动量）、机械能如何耗散为热（能量）。CFD求解器正是对这三个方程的离散化数值求解。理解简化形式，有助于判断：
- 何时可以用一维工程公式快速估算，何时必须做三维CFD；
- 边界条件应设速度入口还是压力入口，取决于已知量和方程形式；
- 网格加密应优先放在速度梯度大、漩涡强的区域（壁面附近、弯头下游、阀门后）；
- 结果验证时，沿程损失应与达西公式对比，局部损失应与ζ系数对比，偏差过大说明网格或模型有问题。