4. 基于允许压降反算最大可行流量

为什么“反算”比“正算”更贴近选型现场

正算压降是“给定流量求 Δp”，而工程上更常见的约束是反过来：泵或系统只允许换热器消耗某个压降上限 Δp_allow，问“最多能走多大流量”。由于湍流区压降近似与流量平方成正比（Δp ≈ K·Q²），层流区近似与流量一次方成正比，因此不能简单线性外推，必须按流态分段反算，再取管程、壳程两个上限中的较小值作为整台设备的可行流量。

本节采用分段幂律模型：Δp = a·Qⁿ，其中 n 由流态决定（层流 n≈1，湍流 n≈2，过渡区取 1.75 作工程近似）。反算公式为 Q_max = (Δp_allow / a)^1/n。关键点是：先假设流态→反算流量→用反算流量复核雷诺数→若流态假设不成立则换段重算，形成一次收敛迭代。

教学示例：一台固定几何的换热器

设管程与壳程的压降系数已由前序计算或厂家数据拟合得到（示例假设值，非实测）：

系统允许换热器总压降上限 Δp_allow = 60 kPa（管程、壳程各自独立受此上限约束，取严者）。

反算过程与流态复核

管程：Q_t = (60 / 0.85)^(1/1.85) = (70.588)^0.5405 ≈ 9.86 m³/h
壳程：Q_s = (60 / 0.42)^(1/2.0)  = (142.857)^0.5  ≈ 11.95 m³/h

取小值：Q_max = min(9.86, 11.95) = 9.86 m³/h（管程控制）

流态复核（示例几何：管程流通面积 A_t = 0.0012 m²，ρ=990 kg/m³，μ=0.0008 Pa·s）：
v = Q/(3600·A_t) = 9.86/(3600×0.0012) ≈ 2.28 m/s
Re = ρ·v·d/μ，取 d=0.02 m → Re ≈ 990×2.28×0.02/0.0008 ≈ 5.6×10⁴
Re > 4000，湍流假设成立，n=1.85 可用，反算结果有效。

若复核发现 Re 落入层流区（<2300），则需把 n 改为 1 重新反算，流量会显著增大，此时往往由壳程转为控制侧——这正是必须做流态复核的原因。

可直接套用的反算小工具（Python）

把上面的分段幂律与流态复核写成函数，输入拟合系数与允许压降即可得到控制侧与最大流量，便于批量核算不同工况。

def q_max_from_dp(dp_allow, a, n):
    """Δp = a·Q^n 反算 Q_max"""
    return (dp_allow / a) ** (1.0 / n)

def re_check(Q_m3h, A_m2, d_m, rho=990.0, mu=8e-4):
    v = Q_m3h / 3600.0 / A_m2
    return rho * v * d_m / mu

# 示例参数（假设值）
dp_allow = 60.0
tube = dict(a=0.85, n=1.85, A=0.0012, d=0.02)
shell = dict(a=0.42, n=2.00, A=0.0045, d=0.025)

Qt = q_max_from_dp(dp_allow, tube["a"], tube["n"])
Qs = q_max_from_dp(dp_allow, shell["a"], shell["n"])
Qmax = min(Qt, Qs)
ctrl = "管程" if Qt < Qs else "壳程"

print(f"管程 Q_t = {Qt:.2f} m³/h, Re = {re_check(Qt, tube['A'], tube['d']):.2e}")
print(f"壳程 Q_s = {Qs:.2f} m³/h, Re = {re_check(Qs, shell['A'], shell['d']):.2e}")
print(f"控制侧 = {ctrl}, Q_max = {Qmax:.2f} m³/h")

运行输出（示例预期）：管程 Q_t≈9.86 m³/h、Re≈5.6×10⁴；壳程 Q_s≈11.95 m³/h、Re 处于湍流；控制侧为管程，Q_max≈9.86 m³/h。若某侧 Re 落入层流，应把该侧 n 改为 1 后重算并再次复核，直至流态与 n 自洽。