2. 不同流态下摩擦系数选取与阻力对比
同一根管、三种流态假设：摩擦系数怎么选
本例取一段瓦斯抽采直管作为固定对象：内径 D = 0.30 m，长度 L = 200 m，管壁绝对粗糙度 ε = 0.15 mm，抽采气体按常温近似取密度 ρ = 1.10 kg/m³、动力黏度 μ = 1.75×10⁻⁵ Pa·s。任务是：在同一流量下分别按层流、过渡流、湍流三种假设取摩擦系数，算出三组沿程阻力，再用雷诺数判断实际应采用的模型。
达西—魏斯巴赫公式为 h_f = λ·(L/D)·(v²/2g)，工程上更常用压降形式 Δp = λ·(L/D)·(ρv²/2)。摩擦系数 λ 的取法取决于流态：层流用 λ = 64/Re；湍流用柯列勃洛克（Colebrook）隐式式或莫迪图；过渡区 λ 随 Re 与 ε/D 剧烈变化，不能简单外推。
三种流态下的计算过程
先固定一个工况流量 Q = 0.070 m³/s（约 252 m³/h），计算流速与雷诺数：
v = Q / (πD²/4) = 0.070 / 0.07069 ≈ 0.990 m/s
Re = ρvD/μ = 1.10 × 0.990 × 0.30 / 1.75e-5 ≈ 1.87×10⁴
相对粗糙度 ε/D = 0.15 / 300 = 5.0×10⁻⁴Re ≈ 1.87×10⁴ 已明显大于 4000，属于湍流。但为了完成“三种假设对比”，仍分别按三种模型取 λ：
- 层流假设：λ = 64/Re = 64/1.87×10⁴ ≈ 0.00342（该值在物理上不成立，仅作对比基准）。
- 过渡流假设：取过渡区常用经验值 λ ≈ 0.035（Re 在 2000~4000 附近波动时的保守估计）。
- 湍流（Colebrook）：解
1/√λ = -2log₁₀(ε/(3.7D) + 2.51/(Re√λ))，迭代得 λ ≈ 0.0268。
代入压降公式 Δp = λ·(L/D)·(ρv²/2)，其中 ρv²/2 = 1.10×0.990²/2 ≈ 0.539 Pa，L/D = 200/0.30 ≈ 666.7：
| 流态假设 | λ 取值 | Δp 计算 | Δp 结果 |
|---|---|---|---|
| 层流 | 0.00342 | 0.00342×666.7×0.539 | ≈ 1.23 Pa |
| 过渡流 | 0.035 | 0.035×666.7×0.539 | ≈ 12.6 Pa |
| 湍流 | 0.0268 | 0.0268×666.7×0.539 | ≈ 9.63 Pa |
三种结果相差近一个数量级：层流假设把阻力严重低估，过渡流假设又偏保守。实际 Re ≈ 1.87×10⁴ 落在湍流区，应取湍流 λ ≈ 0.0268，对应 Δp ≈ 9.6 Pa（200 m 管段）。
流态判定与模型选择逻辑
实际工程中应先算 Re，再决定 λ 的取法，而不是先假定流态。判定阈值：Re < 2000 层流；2000 ≤ Re ≤ 4000 过渡流（结果不确定，需按保守值或实测校核）；Re > 4000 湍流。湍流区再根据 ε/D 与 Re 用 Colebrook 式或莫迪图取 λ。下图给出本例从流量到流态判定、再到三种 λ 取值与阻力对比的完整关系。
结果判读与工程取舍
三种假设的阻力差异说明：摩擦系数不是可以随手取的经验常数，它由流态和相对粗糙度共同决定。本例中：
- 层流公式给出的 1.23 Pa 明显偏低，因为层流时流速剖面为抛物线、黏性剪切主导，而本例惯性力已占主导，公式前提不成立。
- 过渡流取 0.035 是“不知道流态时”的保守估计，用于设计余量可以，但用于精确阻力预测会高估约 30%。
- 湍流 Colebrook 解出的 0.0268 与莫迪图上 Re≈1.9×10⁴、ε/D=5×10⁻⁴ 的位置一致，是本工况应采用的值。
若实际抽采流量在运行中大幅波动（例如从 0.02 m³/s 升到 0.15 m³/s），Re 会跨越 2000~4000 的过渡区，此时应逐工况重算 Re 与 λ，而不是固定一个 λ 值套用全流量范围。对管网阻力计算而言，先判流态、再取 λ、最后代公式是必须保持的顺序。