C1-04 制冷剂路径压降与均流设计
课程代码 C1-04 · 板块 C 乘员舱空调与 HVAC / C1 制冷系统与循环 时长 约 5 小时（7 讲 + 1 次仿真作业） 适合对象 换热器 / 仿真 / 结构工程师；负责换热器芯体与集流管设计、系统仿真诊断 前置 C1-02 四大件选型匹配；A1-02 对流换热 h 估算；A2-01 管网水力学（并联铁律与流量分配，第 2 讲 2.1 直接承接、不重推）；A2-03 两相流基础（一根管里的两相流物理，本课只讲几十根管之间怎么分） 本讲义定位 讲师授课蓝本 / 学员自学讲义，是 C1-04 大纲的完整展开版
引言：换热器有没有被"用满"
换热器选型算出来的换热面积，装到车上却达不到设计制冷量——这类问题里，最常见、也最难查的一个根因不是芯体做小了、不是风量不够，而是制冷剂在换热器内部根本没被均匀地分配到每一根扁管里。有的管分得多、有的管分得少，甚至有的管进的偏气、有的管进的偏液。结果是：一部分换热面在满负荷干活，另一部分在空转或被淹没。整器的名义换热面积没变，能真正参与换热的"有效面积"却被打了折。
这就是偏流（maldistribution），是换热器性能达不到设计值的头号隐形杀手。之所以"隐形"，是因为拆开来看每个零件都正常：扁管没堵、集流管没裂、制冷剂没少，可整器就是不给力；而且你用单管算法去核算怎么算都达标，因为单管算法根本没把"分配"这件事算进去。
这门课的主线是三件事：压降怎么组成、为什么两相不能用单相公式糊弄；流体在集流管和扁管之间怎么分配、为什么两相分配天生易畸变；以及怎么用均流手段和仿真诊断把偏流揪出来、按住。核心判断始终是一句话：换热器的性能，一半取决于芯体选型，另一半取决于它有没有被均匀地用满。
第 1 讲 两相流压降：分配这件事的计价单位
这门课的射程是换热器整器内部的冷媒分配——几十根并联扁管之间，谁分到多少流量、分到什么气液比例。压降之所以要放在第一讲讲，不是因为本课要教怎么算压降，而是因为各扁管分到多少流量，是由它两端的压差与它自身的阻力特性共同定的：压降在本课是分配的计价单位，不是目的。
所以这一讲与 A2-03 两相流基础：沸腾、冷凝与流型 的分工要先划死：A2-03 两相流基础：沸腾、冷凝与流型 讲一根管里的两相流（干度与空泡率、流型图、沸腾与冷凝、两相压降关联式的族谱与适用范围），本课讲几十根管之间怎么分。本讲会用到两相压降的量，但一律只讲「这些量在分配问题里怎么用」，⛔ 不重推两相流本身的物理。不划这条线会有两种坏法：一是把 A2-03 两相流基础：沸腾、冷凝与流型 整讲重讲一遍，读者学两遍同样的关联式而始终没学到分配；二是反过来，以为「会算单管压降就会算换热器」——而单管算法里根本没有「分配」这个变量。
还有三件事不在本课射程内，但必须在这里显式占一行，⛔ 不许因为不讲就整类不出现：
- 并联通道的动态失稳（Ledinegg 流量漂移／压降型振荡／密度波）——本课讲的全是静态（时不变）偏流，摆动型的机理与稳定判据去 A2-03 两相流基础：沸腾、冷凝与流型 第 6 讲，本课⛔ 不给任何稳定判据的限值；
- 干涸与 CHF 的判据本身——本课只用「干涸之后换热系数骤降」这个后果，判据去 A2-03 两相流基础：沸腾、冷凝与流型；
- 冷冻油对分配与压降的影响——去 C1-03 充注量优化与系统油循环管理，本课只写去向。
最后是两个贯穿全课的书写口径，从这一讲起生效：本课凡说「流量」一律指质量流量 ṁ（kg/s），⛔ 不是体积流量（两相体积流量沿程剧变，两者完全不是一个数）；换热量一律写作「换热量」或 Q̇（W），⛔ 正文与图上禁止裸写 Q 表示换热量——第 2 讲要用的不均匀度符号 γ 里那个下标量，指的是各流路的质量流量，两者只有分开写才不会串。
1.1 压降账是四项不是三项，第四项才是均流手段的支出方
常见的三项式写法是
Δp_total = Δp_fric + Δp_acc + Δp_grav
它是对的，但它描述的是沿程分布项之和。换热器里还有一整类集中发生的压降，沿「摩擦／加速／重力」这条线索往下数是数不到它的，因为它根本不在这条线上。本课把压降账记成四项，第四项是局部（集中）阻力 Δp_local——⚠ 这一项是本课按全集自检补入的，⛔ 不是任何标准或既有条款里既有的写法。
它至少发生在五个具体位置，读者可以拿自己的换热器逐个对照数：
- 扁管进口的突缩（集流管截面大、扁管 port 截面小）；
- 入口分配孔板／节流件；
- 扁管出口的突扩；
- 多流程的折返弯；
- 集流管进出口接管处的转向与分／合流。
漏掉这一项不是精度问题，是记账口径问题。本课最核心的那件均流硬件——扁管入口分配孔板——本身就是一个局部阻力件：它的收益（分配变均）是显性的，而它的支出（多出来的那份压降）在三项式的账上根本不出现。一笔只看得见收益、看不见支出的交易，必然被过度使用，或被完全误解。而这件硬件同时出现在「压降账」与「均流手段」两张表上——同一件硬件两个归因，这条会在第 4 讲再敲一次。
四项各自的性质并不一样，先按「要不要分段积分」把它们分开：
摩擦项 Δp_fric——流体与壁面之间的剪切耗散。两相流里气相占据流道中心高速流动、把液相挤成贴壁薄膜，气液界面粗糙、相互拖曳，剪切耗散远比单相剧烈。它随局部干度变，必须沿程分段积分；它的量化方式（两相摩擦倍率）是 1.3 的内容。
加速项 Δp_acc——相变让比容沿程剧变，同样的质量流量要求流速改变，这份动量变化必须由压力梯度提供。它有一条容易被忽略的性质：它是端点量。均相口径下
Δp_acc = G²·(1/ρ_out − 1/ρ_in)
分相口径下它是进出口两个状态的 [x²/(ρ_v·α) + (1−x)²/(ρ_l·(1−α))] 之差（α 为局部空泡率）。两种口径下它都只含端点状态，与中间那条 x(z) 走的是什么路径无关。本课用三条端点相同、中间形状不同的干度路径（线性／前慢后快 x∝z²／前快后慢 x∝√z）代入均相式做过数值验证，三条路径给出完全相同的 Δp_acc（该验证只用于说明「与路径无关」这一件事，其带量纲的结果含教学假设的输入，⛔ 本课不印那个数、也不得被取用为任何压降参考值）。
⇒ 于是「哪些项必须分段」这个问题就答死了：摩擦项与重力项必须分段，加速项不必。⚠ 前提是流通截面 A 与质量流速 G 沿程不变；多流程换 pass 时 G 变了，每个 pass 各自贡献一项、pass 过渡处另有一项（接第 4 讲）。⚠ 而正因为它不用分段，它极容易在「按段列表」的算法里被整段整段地跳过，最后一项都没算——加速项是最常被整项漏掉的一项，而漏掉整项才是真正的低估来源（1.4 再收一次口）。方向可以给死：蒸发段是压降（比容增大、流体被加速），冷凝段是压力回升（比容减小、动量被回收，即负压降）；⛔ 不要把冷凝段的回升当成算错了而人为取绝对值，那会把总账算大。这个符号在第 7 讲换向时会翻过来。
重力项 Δp_grav——形如 ρ_mix·g·Δz·sinθ，其中 ρ_mix 由局部空泡率定，而空泡率随干度沿程变 ⇒ 与摩擦项一样是沿程积分量，⛔ 不能用纯液或纯气密度代替混合密度。⚠ 更要紧的是水平情形：沿程分量小不等于重力不起作用——它转化为集流管横截面内的上下气液分层，直接改写各扁管分到的气液比例。竖置时它改的是压降账，水平时它改的是相分配；安装姿态在这里第一次成为设计变量，第 3 讲会把它讲透。⛔ 别把「水平放置重力项可略」这句从压降账搬到相分配账上——那是把一个结论搬出了它自己的射程。
局部项 Δp_local——由各处的局部损失系数与当地动压定，可略与否是算出来的结论，不是默认值：有孔板／节流件时它不可略；纯光管、无孔板的短流程里它可能可略。⛔ 不要把局部损失系数 ζ（量纲一，部分手册写作 K_L）从单相手册直接搬来用在两相段——两相的突缩、突扩另有修正，各自的形式与射程去 A2-01 管网水力学：沿程/局部压降与流量分配 与 A2-03 两相流基础：沸腾、冷凝与流型 查，本课不给任何 ζ 值。⚠ 这里必须先把一处撞名钉死：本课第 2 讲 2.3 起会反复出现支路特性式 Δp_i ＝ K·ṁ_i^n，那个 K 是「支路特性系数」，量纲为 Pa/(kg/s)^n，与这里量纲一的 ζ ⛔ 不是同一个量（两者只是习惯上都被写成 K）。「本课不给值」这句话在这两处的理由也不同：ζ 是有手册可查而两相段不适用，支路特性系数 K 是本来就只能由你自己那条支路的 Δp–ṁ 特性曲线定的平台相关量。
1.2 均相还是分相由流型与质量流速定，选错是量级问题不是精度问题
两相流的精确描述极其复杂，工程上靠两类简化：
- 均相模型（homogeneous）：假设气液充分混合、以同一速度流动（无滑移），把两相当成一种取平均密度与平均黏度的伪单相流体。它适合高质量流速、雾状／环状流、微通道——管径小、表面张力强，气液难以各走各的，天然接近这个假设。
- 分相模型（separated）：承认气液各占一部分流道、各有各的速度（滑移比大于 1），分别列动量方程再耦合。Lockhart-Martinelli 族属于这一类，适合低质量流速、气液明显分离（分层流、大管径）。
本课要给出的不是「两个模型的介绍」，而是选错的代价有多大。用同一组物性、同一个干度，两条常用关联式给出的两相摩擦倍率 φ_lo²（定义见 1.3）在中高干度段可以差一倍以上：
| 干度 x | 均相 φ_lo² | 分相 φ_lo² |
|---|---|---|
| 0.10 | 6.26 | 15.37 |
| 0.30 | 13.91 | 33.96 |
| 0.50 | 20.33 | 44.54 |
| 0.70 | 26.12 | 47.67 |
| 0.90 | 31.52 | 42.31 |
（物性基准态＝R134a 饱和 5 ℃：p_sat = 349.7 kPa、ρ_l = 1278.1 kg/m³、ρ_v = 17.131 kg/m³ ⇒ ρ_l/ρ_v = 74.61；μ_l = 2.501e−4 Pa·s、μ_v = 1.091e−5 Pa·s ⇒ μ_l/μ_v = 22.92。⚠ 这个饱和温度是教学假设值，选它只为把低压侧的对照拉开，⛔ 不对应任何平台的实际蒸发温度，⛔ 禁止被当作可取用值；分相曲线的极大值出现在 x ≈ 0.681、φ_lo² ≈ 47.70。）
这个量级差意味着一件事：「先随便选一个模型、不确定度靠安全系数补」的做法不成立。⚠ 而这张对照表不是「哪个更对」的判据——哪个更对由流型定，流型判据本课不判，去 A2-03 两相流基础：沸腾、冷凝与流型。
三个易错点：① 微通道套分相、粗管路套均相，正好选反；② 在同一门课、同一份计算书里中途换模型——同一处修正的两种参数化只能取一，⛔ 不得中途换用、⛔ 不得叠乘；③ 把两个模型的结果取平均——它们是两套前提，平均出来的数没有任何前提。
1.3 两相摩擦倍率的折算基准必须写死：全课只用 φ_lo²
两相摩擦压降的工程写法是把它折算回一个熟悉的单相基准再乘一个倍率：
Δp_fric = f·(L/D)·(ρu²/2)·φ_lo²
φ_lo² ≡（实际两相摩擦压降）/（假设总质量流量全部以液相流过同一流道时的摩擦压降），下标 lo ＝ liquid only，量纲一。
⚠ 这里有一个通则级别的坑：另一个常见基准是 φ_l²，它的分母换成「只有液相那一部分质量流量单独流过同一流道」。两者只差一个下标、量纲相同、名字都叫「两相摩擦倍率」，而换算关系里含 (1−x) 的高次幂，在高干度段差出很多倍。拿一个基准的数代进另一个基准的式子，得到的是量级级别的错。⇒ 本课全程只用 φ_lo²（全液基准），⛔ 正文、图、图注与口头讲述里一律禁止出现不带下标的 φ²，⛔ 禁止与 φ_l² 混用。从文献抄一条倍率关联式时，第一件事是核它的基准，而不是核它的系数。
本课全程固定用下面这两条关联式（⛔ 不中途换用）：
- 均相（McAdams 混合黏度 + Blasius，f ∝ Re^−1/4）：φ_lo² = [1 + x(ρ_l/ρ_v − 1)]·[1 + x(μ_l/μ_v − 1)]^(−1/4)
- 分相（Lockhart-Martinelli 湍-湍 tt、C = 20，写成 Chisholm 的 φ_lo² 形式）：Γ² = (ρ_l/ρ_v)·(μ_v/μ_l)^(1/4)、B = C/Γ、φ_lo² = 1 + (Γ² − 1)·[B·x^0.875·(1−x)^0.875 + x^1.75]
逐步手算复核（用 1.2 的物性，x = 0.10）：均相 —— 1 + 0.10×(74.61 − 1) = 8.361；1 + 0.10×(22.92 − 1) = 3.192；3.192^(−1/4) = 0.7481；8.361 × 0.7481 = 6.26。分相 —— Γ² = 74.61 × (1/22.92)^(1/4) = 74.61 × 0.4570 = 34.10，Γ = 5.839，B = 20/5.839 = 3.425；方括号 = 3.425 × 0.1^0.875 × 0.9^0.875 + 0.1^1.75 = 3.425 × 0.13335 × 0.91193 + 0.01778 = 0.4343；φ_lo² = 1 + 33.10 × 0.4343 = 15.37。两个数与 1.2 表里的一致。
φ_lo² 在两相段一般远大于 1，这给出一条红线：⛔ 绝不能用单相公式算完再乘一个安全系数就交差——差的是量级，不是系数。⚠ 但反过来也要注意：「φ_lo² 常比单相大数倍」里的「单相」指的是全液基准，⛔ 不是「同样气液两相但按纯气算」。⛔ 具体设计取值须按实际工质、质量流速与流型另查关联式，并核该关联式自己的适用范围。
1.4 一刀切的误差方向由 φ_lo²(x) 的凹凸决定，⛔ 不是靠直觉判的
干度 x ＝ m_vapor / m_total（气相质量占总质量的比例）在换热器里不是常数：蒸发器进口是节流后闪发出来的低干度两相，沿程吸热、干度一路升高；冷凝器进口是过热气，冷凝段干度从 1 降到 0，出口是过冷液。⚠ 而且 x(z) 本身是非线性的——各段热流密度不同，干涸之后换热系数骤降，x(z) 在干涸点前后是两条不同的曲线（干涸判据本身去 A2-03 两相流基础：沸腾、冷凝与流型）。⇒ 「x 沿程均匀分布」这个前提本来就不成立，摩擦项必须沿程分段。
于是就有了那个几乎人人都会问的问题：如果偷懒用一个「平均干度」一刀切，算出来是偏大还是偏小？
流传最广的直觉是：一刀切抹平了高干度段的陡增，所以算出来偏小，再乘个安全系数才稳妥。⚠ 这个方向是反的。
理由不在直觉里，在曲线的形状里。设沿程积分的平均值为 ⟨φ_lo²⟩，用平均干度直接代入得到的值为 φ_lo²(x̄)。按 Jensen 不等式，曲线凹（二阶导为负）⇒ φ_lo²(x̄) > ⟨φ_lo²⟩，即一刀切是高估。而本课这两条关联式的实算结果是：在 x∈(0,1) 上二阶导为负的采样点占 100%（均相的二阶差分最大值为 −3.3、分相为 −183.7——即最接近零的那个采样点也仍是负的），两条曲线都严格凹。
直觉之所以会反，是因为读者记住的是「φ_lo² 随干度大幅上升」，于是脑子里画的是一条凸（越来越陡地上翘）的曲线；而图 2 里这两条曲线的实际形状是低干度段陡升、中高干度段变缓甚至回落——陡的那一段在左边，不在右边。
实算的高估量级（正号一律表示高估）：
| 干度区间 | 均相 | 分相 |
|---|---|---|
| x∈[0.05, 0.95] | +2.9% | +17.1% |
| x∈[0.02, 0.98] | +3.5% | +20.2% |
⚠ 这里最要紧的一句不是那个方向，而是判据本身：⛔ 不得把它写成「一刀切一定高估」这样的全称——本结论只对上面这两条关联式、这个教学假设饱和态与这两个干度区间成立；换一条关联式（含流型分区或 Froude 数的形式），曲线形状可能不同，方向必须对你自己所用的关联式现算。⛔ 这两个百分数也只覆盖摩擦这一项，⛔ 不得被当作整体压降的通用修正系数取用。
⇒ 由此把「两相压降算不准」的成因分成两层，两层的处置完全不同：
- 第一层：整项缺失——漏掉加速项、漏掉局部项、φ_lo² 直接取 1（用单相公式）。这三种是量级级别的错。
- 第二层：积分方式——一刀切 vs 分段。本课实算是百分之几到二十几，且方向要现算。
不分层就会出现最坏的一种做法：知道「两相压降容易算不准」，于是笼统地乘一个安全系数。⚠ 安全系数补得了第二层、补不了第一层（少了一整项，乘多少都补不出来）；更坏的是，它在补第一层之前，先把第二层已经存在的高估再放大一次——压降预算被放得过大，据此把流道做大，就正好接上本课反钉子那条链（第 2 讲会说明为什么把流道做大反而会让分配更不均）。⇒ 本课的立场是：用把项算全来代替安全系数，⛔ 本课不给任何安全系数取值。
干度沿程变化还有第二重意义，它把本课的前后两半接了起来：干度不只是压降账里的一个自变量，它还是分配畸变的驱动源。并联的每一根扁管都有自己的一条 x(z) 曲线——只要它们的进口干度或热负荷不同，这些曲线就不同；而曲线不同 ⇒ 各管沿程的 φ_lo² 积分与加速项不同 ⇒ 各管的阻力不同 ⇒ 反过来改写流量分配。
⇒ 到了并联分配问题里，x 从自变量变成了耦合变量：它既被流量决定（同样热负荷下，流量少的管干度升得更快），又反过来决定流量。⚠ 所以并联两相通道的分配问题在数学上是耦合的，⛔ 不能先解流量再算干度，必须迭代；用一条平均的 x(z) 代表所有管，等于把耦合项直接消掉，算出来的分配永远比实际均匀。
关于这个耦合的方向，本讲只给两句、⛔ 不多给：① 在沸腾侧，干度更高的那根管其沿程 φ_lo² 积分与加速项一般更大、阻力更大、于是分得更少——⛔ 但这不是一条可以外推到全干度区间的单调关系（图 2 里分相曲线在 x ≈ 0.681 之后回落），高干度段的方向要按你自己那条曲线现判；② 在冷凝侧这条耦合的符号是相反的，第 5 讲会把两侧的符号并起来讲，第 7 讲会说明换向之后同一根扁管上这个符号会翻转。⚠ 而「这个耦合会不会把支路特性推成非单调、乃至失稳」是另一条判据链，⛔ 不能从本条读出来——它去 A2-03 两相流基础：沸腾、冷凝与流型 第 6 讲；本课第 2 讲只把「支路特性 Δp–ṁ 是否单调递增」当成手算判据的一个硬前提来用。
1.5 定 G 有三条约束不是两条，而第三条是否决线
选管径、定流程数，本质上都是在定管内质量流速 G ＝ ṁ/A_c（kg/(m²·s)）。工程上通常只在两条约束之间找 G：
- 换热（性能账）：G 高则管内换热系数高；
- 压降（能耗账）：G 高则压降大（摩擦项含 G² 依赖，加速项也含 G²）。
冷媒侧还有第三条常被漏掉的约束——噪声（感知账）：两相段流速过高时，节流后的闪发、气液界面扰动、气泡的生成与崩溃直接产生可闻流噪，经管路、支架与空调箱壳体传进座舱。
三条约束的性质不同，这才是关键，也是本课在这一节要立的判断（⚠ 这条判断是本课归纳的，⛔ 不是任何标准既有）：压降是一项可以做加减法的代价——它换来的换热收益可以和它相抵；而噪声是一条否决线——越线就是整改，没有「多一点噪声换多一点换热」这种交易。尤其蒸发器离乘员耳朵最近，整车静置或低速时又没有掩蔽噪声，感知账常常比能耗账先爆。
⇒ 由此得到求解顺序：先用噪声划出可行域，再在域内对换热与压降做加减法；⛔ 不要把噪声当成一项可以被换热收益抵消的代价，也⛔ 不要等到整车 NVH 阶段才发现越线——那时芯体早已冻结。
三条纪律，与这一节配套：
- ⛔ 本课不给 G 的工程取值区间——它由换热需求、压降预算与整车 NVH 目标三条约束共同定，须由本项目的系统仿真与台架确定。
- ⛔ 本课不给流噪随流速上升的指数，只承接一条方向性判断：流致噪声的声功率随流速高次幂上升、比压降账更陡。⚠ Lighthill 的 u⁶／u⁸ 标度是单相气动声源的结论，两相闪发、气液界面扰动与气泡崩溃是另一类声源，跨射程引用不成立；⛔ 既然不给单点值，就不得写「噪声比压降高一到两个数量级」这类倍数比较，也⛔ 不得拿一个指数去反算「G 抬多少、噪声涨多少」。
- ⛔ 所选两相关联式自带 G 的适用范围，代进去之前先核关联式自己的射程，⛔ 不得外插；关联式的 G 适用范围去 A2-03 两相流基础：沸腾、冷凝与流型，本课不转述。
这一条会在第 6 讲的压降—换热—噪声三方权衡里再收一次口——到那里还要加上一层：孔板这类加阻力的均流件本身也是流噪源，所以「买均匀」这笔交易付的从来不只是压降。
本讲小结：压降在本课是分配的计价单位不是目的，一根管里的两相流物理承 A2-03 两相流基础：沸腾、冷凝与流型、本课不重推；压降账是四项不是三项——摩擦与重力必须分段积分，加速项是只由端点决定的端点量（不必分段，却最常被整项漏掉，冷凝时是压力回升），第四项局部阻力是本课按全集自检补入的，而入口分配孔板正是它与「均流手段」共用的那件硬件；均相与分相由流型与质量流速定，同一干度上可差一倍以上，⛔ 不得中途换用、不得取平均；两相摩擦倍率全课只用全液基准 φ_lo²，⛔ 禁止裸写不带下标的形式、禁止与 φ_l² 混用；用平均干度一刀切的误差方向由 φ_lo²(x) 的凹凸决定，本课这两条关联式都严格凹 ⇒ 是高估（+2.9% / +17.1%），⛔ 但不得写成全称，换关联式必须现算——真正的低估来源是整项缺失，那是安全系数补不了的；干度沿程变化同时是分配畸变的驱动源，并联两相通道的分配在数学上是耦合的、必须迭代，这是第 2、3 讲的物理起点；定 G 有换热、压降、噪声三条约束，前两者是可加减的代价、噪声是否决线（本课归纳），⛔ 本课不给 G 区间、不给流噪指数。
后面还有 6 讲正文 · 关键公式 · 案例拆解 · 常见误区 · 动手做