3. 一维热流体建模的物理基础
一维热流体建模的本质，是把三维连续介质中的守恒规律沿流动主方向做降维处理：保留沿程方向的空间变化，把垂直于流动方向的截面效应通过经验关系或集中参数折算进来。理解质量、能量、动量三大守恒方程在一维框架下的简化形式，以及何时该用集中参数、何时必须用分布参数，是后续搭建管路、换热器、泵阀与冷板模型的前提。
一、控制体与符号约定
先统一参照约定，全文与图示保持一致。取流动方向为 x 正方向，控制体为管段内介于 x 与 x+dx 之间的微元，截面积 A、水力直径 Dh、湿周 P。质量流量 ṁ 沿 x 正向为正，换热器热侧放热记为负、冷侧吸热记为正。工质密度 ρ、比焓 h、流速 u、压力 p、壁温 Tw、流体温度 Tf。所有方程均以“进入控制体为正、离开为负”的净通量形式书写。
二、质量守恒：一维连续性方程
对可压缩一维流动，微元控制体的质量守恒为：
∂(ρA)/∂t + ∂(ṁ)/∂x = 0
其中 ṁ = ρ·u·A在热管理系统中，绝大多数液体回路（乙二醇水溶液、冷却油、制冷剂液相段）可视为不可压缩，ρ 沿程近似不变，于是简化为 ∂ṁ/∂x = 0，即稳态下沿程质量流量恒定。这一条看似简单，却是建模时最实用的约束：管路串联段流量相等，并联支路按节点压力分配。对制冷剂两相段或气态回路，密度随压力显著变化，必须保留 ∂(ρA)/∂t 项，此时质量守恒与动量、能量方程强耦合，通常需要联立求解。
工程上还有一个常被忽略的细节：当管径沿程变化（如扩口、集流管）时，即使不可压缩，流速 u 也会随 A 变化，u₁A₁ = u₂A₂。这直接影响后续对流换热系数与压降的计算，建模时应在几何参数表中显式记录每个节点的 A 与 Dh。
三、能量守恒：从焓方程到换热器简化式
一维能量守恒的通用形式为：
∂(ρA·h)/∂t + ∂(ṁ·h)/∂x = q' - ṁ·g·dz/dx + Φ
q' = P·α·(T_w - T_f) （单位长度换热量）其中 α 为对流换热系数，Φ 为黏性耗散（液体低速时通常可忽略）。稳态、忽略位能与耗散后，得到管路与换热器中最常用的形式：
ṁ·(h_out - h_in) = q_total
q_total = ∫ P·α·(T_w - T_f) dx对单相液体，dh ≈ cp·dT，于是 ṁ·c_p·(T_out - T_in) = q_total，这就是换热器热平衡计算的核心。以某液冷板为例：冷侧流量 ṁ = 0.15 kg/s，cp = 3800 J/(kg·K)，进口 25 ℃、出口 31 ℃，则吸热量 q = 0.15×3800×6 ≈ 3420 W。若热侧为 45 ℃ 进水、40 ℃ 出水、流量 0.12 kg/s、cp = 4180 J/(kg·K)，则放热量约 2508 W。两者不相等说明存在向环境的散热或测量误差，建模时应把外壳散热作为独立热阻支路补入，而不是强行让两侧相等。
换热器建模的另一个关键选择是“ε-NTU”还是“分布参数”。当两侧温差不剧烈、单相且流道结构简单时，用 ε-NTU 集中参数即可；当出现相变、大温差或需要捕捉沿程温度分布时，必须把换热器离散为若干沿程单元，逐段用上面的焓方程推进。
四、动量守恒：压降的组成与简化
一维动量方程可写为：
∂(ṁ)/∂t + ∂(ṁ²/(ρA))/∂x = -A·∂p/∂x - τ_w·P - ρ·g·A·sinθ稳态、不可压缩、水平管路下，简化为压降平衡：
Δp = Δp_friction + Δp_local + Δp_acceleration
Δp_friction = f·(L/D_h)·(ρu²/2)
Δp_local = ΣK·(ρu²/2)
Δp_acceleration = (ρu²)_out - (ρu²)_in摩擦因子 f 按流态选取：层流 f = 64/Re，湍流用 Blasius 或 Colebrook。局部阻力系数 K 需按实际管件查表，弯头约 0.3～0.9、突扩约 1.0、突缩约 0.5。加速压降在单相液体中通常很小，但在两相或高流速气体回路中不可忽略。建模时建议把每个管段的 L、Dh、粗糙度、K 汇总成参数表，避免在方程里硬编码。
五、集中参数与分布参数的选择依据
选择依据可归纳为三条判据：一是毕渥数 Bi = αLc/λ，Bi < 0.1 时固体内部温度均匀，可用集中参数；二是沿程温升与总温差的比值，若单段温升小于总温差的 5%～10%，集中参数误差可接受；三是时间尺度，若流体穿越时间远小于热惯性时间，动态过程可用集中参数近似。反之，出现相变、大温差、强瞬态或需要定位热点时，必须用分布参数。
六、建模落地建议
实际搭建一维模型时，建议按以下顺序推进：先建立几何与物性参数表，记录每段的 A、Dh、L、粗糙度、K 与 cp、ρ、λ；再按“质量→能量→动量”的顺序逐段写出稳态方程，确认流量分配与温度推进方向；最后根据 Bi 与温升判据决定离散段数。对液冷板这类对象，通常沿程离散 10～20 段即可捕捉温度分布，段数过多会显著增加求解时间而收益有限。对制冷剂蒸发器，则必须按干度分区离散，每段内用平均干度计算 α 与压降。