A3-05 湿工况换热器计算:析湿、湿翅片效率与焓差驱动的 ε-NTU
课程代码 A3-05 · 板块 A 基础理论与共性技术 / A3 换热器与工质基础 时长 约 3.5 小时(5 讲 + 1 次湿工况正算与台架对账实操 45 分钟) 适合对象 仿真 / 结构工程师,尤其是做部件级 1D 与芯体模型标定的人(必修) 前置 A2-04 湿空气物理与焓湿图:空调负荷计算基石;A3-01 换热器传热基础:ε-NTU 法与 LMTD 法;A3-02 翅片管/扁管/微通道换热形式与传热强化 本讲义定位 讲师授课蓝本 / 学员自学讲义,是 A3-05 大纲的完整展开版
引言:翅片一开始淌水,模型就有两处不对,不是一处
一台蒸发器的 1D 模型,在常温中湿工况下与焓差台架对得很好;一进高温高湿 pull-down,它同时冒出两个偏差——模型报的出风干球温度偏低(温降被算大了),而总制冷量又偏低。很多人第一反应是"标定没调准",于是回头去改空气侧换热系数、去乘一个湿工况系数。改不动,因为错的不是系数,是口径。这门课要敲的第一根钉子就是:析湿一次改两处——驱动势由温差换成焓差,翅片效率由干态的 η_f,dry 掉到湿态的 η_f,wet;只改一处就是半个模型。凡是"湿工况就是干工况乘个系数"的做法,补的都只是驱动势那一半,总表面效率 η₀ 的折损照样漏掉,剩下的偏差仍然是系统性的、方向固定的,重复试验平均不掉。
讲到这里,最容易顺口冒出来的一句是:"析湿系数 ξ 恒大于 1,湿表面比同温差的干表面多带走好几倍热,所以拿干工况算法去算湿盘管,是偏保守、偏安全的。"这句话必须当场拆掉,它是本课最容易被读反的一条。读反之后会做出两个具体的错误动作。其一,把"用干工况 η₀ 算湿盘管"当成留了余量,于是不改口径、只加一个安全系数——而实际方向是高估:干态 η₀ 高估了有效面积,于是高估热导、高估出风温降。ξ 大于 1 多出来的那部分是潜热,它去凝水,不降干球温度;所以"湿态换热更强"推不出"出风更冷",更推不出"干算法偏安全"。后果落在出风温度、除湿量、防雾这些恰恰是验收指标的量上,全部偏乐观。其二,把总制冷量的偏差与温降的偏差当成同一个方向去调参——总热是低估(潜热不进账),温降是高估,方向相反;任何"整体乘一个系数"的标定,都必然把其中一项调得更糟。正确的读法只有一句:湿工况总热更大,但多出来的那部分热不降温;干工况口径带来的是一低一高、方向相反的两个系统性偏差,不是一个可以靠安全系数吸收的单向裕度。
第二条同样容易读反的是:"Le≈1 只是个方便的近似,差个 10% 无所谓。"它不是精度问题,而是整个焓势法能不能成立的定义性前提——把显热与潜热两条通道合并成一条焓差通道,靠的正是这一步代换,前提没了,合并本身就不成立。而且刘易斯因子的真值小于 1(常温常压下约 0.90~0.93,取值随所选扩散系数关联式而变,引用时必须声明来源),这意味着按 Le=1 算是低估传质通道——真实的传质系数比 h_c/c_p,m 大约 8%~11%,方向同样和"近似总是偏保守"的直觉相反。它的正确处置不是弃用焓势法,而是加一项修正,第 1 讲给出判据与量级。
第二根钉子是口径纪律,它的好处是量纲当场可判:温差驱动的 UA 量纲是 W/K,焓差驱动的 U_h·A_a 量纲是 kg/s(U_h 的单位是 kg/(m²·s),不是 W/(m²·K));温度效能 ε 与焓效能 ε_h 也不是同一个效能,本课全篇分用、不写裸的效能符号。同样地,台架上那条 Q = ṁ_a·(h_in − h_out) 是测量口径——进出口一测就得总热,但它不回答"这台芯体为什么只能换这么多";焓势法是正算口径,由几何与表面温度把热量算出来。两者形式相似而性质不同。⛔ 本课不会出现"拿干态效能反算湿态 UA""湿态 UA 与干态 UA 直接比大小"这两种写法:它们连量纲都不同,是量级级别的错,不是精度问题;交付表里这两个量必须分列两行、各带单位。
最后交代射程。本课是跨产品通用的方法课,只给算法与判据,不碰结构、不碰回路、不碰控制。蒸发器的结构形式、冷媒分配均流与出风温度均匀性在 G2-03 蒸发器设计与出风温度均匀性;箱内布置、凝水管理与"湿态压降上升、选风机必须按湿态"那本账在 C2-05 箱内蒸发器布置与冷凝水管理,排水防霉与自清洁在 C2-12 冷凝水排放、防霉与蒸发器自清洁;热泵室外机的双工况角色互换在 G2-09 直接式热泵室外换热器(冷凝/蒸发双向工况)部件设计;除湿再热的回路构型与解耦控制在 C3-06 热泵除湿再热模式:冷媒侧构型、中压与双 EXV 解耦控制。焓湿图判读、过程线与显热比拆分是前置课 A2-04 湿空气物理与焓湿图:空调负荷计算基石 的射程,本课假定你已能在图上定点,一句不重讲;台架怎么搭、数据怎么判有效在 L1-01 换热器性能试验(风洞焓差法)(部件级芯体风洞焓差法)与 L2-01 焓差实验室与空调系统性能试验(整机/系统级),本课只消费它们产出的数据。射程的下边界是 0℃:表面温度跌破冰点、转为凝华结霜之后,析湿系数、饱和焓斜率与饱和焓—温关系这三套量全部失效——那不是精度下降,是换了机理,交 A3-04 污垢、结霜与传热衰减机理。
五讲的走法与两根钉子一一对上:第 1 讲讲清湿表面上到底发生了什么(两条通道占用同一块面积、析湿系数的含义与量级、刘易斯关系的前提与失效);第 2 讲把驱动势换成焓差,建立饱和空气等效焓与焓势法主式,并声明全课的压力口径;第 3 讲把翅片效率换成湿态、重算总表面效率,说清"用干态值是系统性高估而非随机误差";第 4 讲把这两处一起送进焓差驱动的 ε-NTU,并处理部分湿盘管的分区与迭代;第 5 讲与 ADP-BF、焓差台架对账,把偏差逐类归因,再用多个高湿工况点反标定。课末的 45 分钟动手做把这条链走一遍。请记住这门课要交付的东西:不是一组标定完就当常数用的参数,而是一条换了风量、换了蒸发温度、换了进风含湿量还能重算的正算通路。
第 1 讲 湿表面上发生了什么:显热与潜热在同一片翅片上耦合
干工况的换热器计算里,空气侧只做一件事——把温度降下来。而表面温度一旦低于进风露点,水蒸气就开始在翅片上凝结,空气侧同时在做两件事:降温,脱水。本讲不给算法,只回答"这一层上到底发生了什么",因为后面四讲每一条式子的合法性都压在本讲的三个判断上——两条通道共用同一块面积、析湿系数把潜热折算成显热的一个倍数、刘易斯关系是把两条通道并成一条的定义性前提。这三条里任何一条被当成"差不多就行",后面算出来的 UA、总表面效率与效能都是建在错地基上的。
1.1 两条通道,一块面积——显热与潜热不是两台并联的换热器
是什么。湿表面上同时进行两件事:传热,驱动势是温差 (T_a − T_w);传质,驱动势是含湿量差 (d_a − d_s,w)。这里 T_a、d_a 是主流空气的干球温度与含湿量,T_w 是表面(凝水膜面)温度,d_s,w 是与 T_w 相对应的饱和含湿量。凝结在膜面上放出的汽化潜热并不是"另一条路"上的热——它同样要穿过同一层空气侧边界层、同一片翅片、同一段管壁才能交给冷媒。
本课符号约定:含湿量统一记 d,入公式一律用 kg/kg,报数用 g/kg;汽化潜热统一记 h_fg(部分中文资料记作 r);空气侧对流换热系数记 h_c,单位 W/(m²·K);传质系数记 h_D;湿空气比热记 c_p,m,按每 kg 干空气计,单位 kJ/(kg·K)。空气侧总面积统一记 A_a,全课的参照面积口径都取它。
为什么这件事必须先钉死。因为两条通道共用同一串热阻,所以不存在"显热面积"与"潜热面积"之分。把显热与潜热当成两台并联的换热器、各自配一块面积再相加,会把所需总面积算成实际的两倍——这不是精度问题,是结构性的错。反过来,也不能把凝结那部分热当成"额外负荷"再去追加面积:它没有多占一寸翅片。
工程量级。取一个高温高湿的教学算例:进风 38℃ / 60%RH(h_a = 103.90 kJ/kg 干空气、d_a = 25.55 g/kg),表面温度取 T_w = 12℃(教学假设值,不对应任何平台,见下)。此工况下总热是显热的 2.55 倍——也就是说显热只占总热的 39.3%、潜热占 60.7%,而潜热通道一寸额外面积都没有占。(39.3% / 60.7% 由未舍入的倍数复算给出,与"2.55"这个两位显示值直接相除会差 0.1 个百分点,属显示舍入。)湿空气焓与饱和空气焓一律取 ASHRAE 口径:干空气 0℃、液态水 0℃ 为焓零点,压力 101.325 kPa,每 kg 干空气;本讲全部物性值由 CoolProp 8.0.0 的 HAPropsSI 逐点算出。
易错点。两个方向相反、后果都很实:① 先用析水量乘汽化潜热算出潜热,再把它当"额外负荷"去加一块面积——面积被算多;② 只用 c_p·ΔT 记账,潜热完全不进总量——总热被算少。两种错都源自同一个没想清楚的问题:这两条通道到底占不占同一块面积。
1.2 析湿系数 ξ 是这两条通道的总账;ξ < 1 不是"小于 1 的系数",是判据在报警
是什么。析湿系数定义为总热与显热之比:
ξ = Q_总 / Q_显 = 1 + h_fg·(d_a − d_s,w) / [c_p,m·(T_a − T_w)]
把 1.1 里的传质系数按刘易斯关系代换掉(那一步见 1.3),它有一个更常用的等价形式:
ξ = (h_a − h_s,w) / [c_p,m·(T_a − T_w)]
其中 h_s,w 是"贴着湿表面、温度为 T_w 的那层饱和空气"的焓,即 h_s(T_w)。它不是表面上真实存在的某股空气,只是一个由 T_w 算出来的换算量——它的构造与它为什么不可测,是第 2 讲的内容,本讲只借用记号。
为什么要有这么一个系数。因为它把潜热通道折算成显热通道的一个倍数:只要拿到 ξ,干工况那一套式子换一个系数就能复用,不必另起一套潜热的方程。第 3 讲的湿翅片效率、第 4 讲的湿态效能,走的都是这条路。
工程量级(全部输入已标注,⛔ 不存在可照搬的通用值)。四个逐点复算的算例,压力 101.325 kPa、ASHRAE 焓基准,表面温度均为教学假设值:
| 进风状态 | 进风露点 | 表面温度 T_w(教学假设值) | ξ |
|---|---|---|---|
| 27℃ / 60%RH(h = 61.54 kJ/kg、d = 13.483 g/kg) | 18.58℃ | 10℃ | 1.84 |
| 38℃ / 60%RH(h = 103.90 kJ/kg、d = 25.55 g/kg) | 28.88℃ | 12℃ | 2.55 |
| 35℃ / 40%RH(h = 71.64 kJ/kg) | 19.39℃ | 10℃ | 1.64 |
| 25℃ / 40%RH(h = 45.31 kJ/kg、d = 7.915 g/kg) | 10.48℃ | 12℃ | 0.838 ← 越界 |
方向由此可读:进风含湿量升高、表面温度降低,ξ 增大。同一台芯体换个工况,ξ 能差一倍以上(1.64 → 2.55)。⛔ 因此本课不给 ξ 的通用值——它取决于进风状态与表面温度,必须按自己这台芯体的工况点算,或由焓差台架数据反标定。
易错点一:把 ξ 当成可以照搬的常数。上表四行用的是同一条式子,结果跨了三倍多。任何"湿工况就按 ×× 倍算"的做法,都是把一个随工况走的量当成了器件属性。
易错点二(本课归纳的判据):算出 ξ < 1 时,正确解读是"式子被用在了它的适用域之外",不是"得到了一个小于 1 的系数"。看上表最后一行:25℃ / 40%RH 的露点只有 10.48℃,而设的表面温度是 12℃——表面根本是干的,压根没有析湿,同一条式子当然会算出小于 1 的结果。此时若照样把 ξ = 0.838 代进第 3 讲的翅片参数 m_wet = √(2·ξ·h_c/(k·δ)),就会得出"析湿反而提高了翅片效率"这种明显反物理的假结论。ξ < 1 是干湿分区判据在报警(分区怎么做在第 4 讲),不是一个数值结果。
易错点三,也是本课最容易被读反的一条。读到"ξ 恒大于 1、湿表面比同温差的干表面多带走好几倍热",顺手就会推出一句错话:"那用干工况算法算湿盘管就是少算了,偏保守、偏安全。"这句话推不出来,而且方向恰好相反:
- ξ > 1 多出来的那部分热是潜热,它去凝水,不降干球温度。所以"湿态换热更强"推不出"出风更冷";
- 用干工况的总表面效率去算湿盘管,是高估了有效面积,因而高估温降(这一重高估的机理与量级在第 3 讲);
- 而总热若仍按 c_p·ΔT 记账,潜热整条不进账,是低估——本节那个 38℃ / 60%RH 算例里,被漏掉的正是总热的 60.7%。
⇒ 干工况口径带来的是一低一高、方向相反的两个系统性偏差,不是一个可以靠安全系数吸收的单向裕度。读反之后会做的那个具体错误动作是:不改口径,只加一个安全系数——而任何"整体乘一个系数"的标定,必然把总热与温降这两项里的一项调得更糟。这条指纹会在典型案例里再出现一次。
1.3 把两条通道并成一条,靠的是刘易斯关系——它是定义性前提,不是精度问题
是什么。刘易斯数是热扩散率与质扩散率之比,Le = α / D_AB;工程上真正用到的是刘易斯因子 Le_f ≡ Le^(2/3),它出现在传热与传质系数的关联式里:
h_c / (h_D · c_p,m) = Le^(2/3) = Le_f
当 Le_f → 1 时这条关系退化成 h_D = h_c / c_p,m,把 1.1 那两支通道相加,含湿量差与温差就合并成一个焓差,得到第 2 讲要建立的那条主式 q″ = (h_c/c_p,m)·(h_a − h_s,w)。
为什么它是"定义性前提"而不是"一个近似"。焓势法不是换个变量重写一遍——它是靠这条关系把传质系数 h_D 消掉。没有这一步代换,(h_a − h_s,w) 这个合成驱动势根本凑不出来,两条通道就只能各带各的系数分开算。所以 Le 的地位与"某个系数取 0.9 还是 0.95"完全不同:它决定的是后面四讲的整套式子成不成立,而不是它们准不准。⇒ "Le≈1 只是个方便的近似,差个 10% 无所谓"这句话是错的——差的那 10% 不落在结果的精度上,而落在"能不能写出这条式子"上。
工程量级。空气—水蒸气体系在常温常压下 Le ≈ 0.85~0.90,对应 Le_f ≈ 0.90~0.93。⚠ 这两个区间随所选的 D_AB 关联式而变,取值必须声明来源(α 由 k/(ρ·c_p) 算;D_AB 取自 ASHRAE Handbook—Fundamentals 的传质章或所用关联式,版本以现行有效版为准)。
易错点一:以为这个近似"偏保守"。方向恰好相反。因为 Le_f < 1,真实的传质系数是 h_D = h_c / (c_p,m · Le_f),比 h_c / c_p,m 大约 8%~11%——按 Le = 1 算是低估潜热通道,不是留了余量。"近似总是偏保守"这条直觉在这里不成立。
易错点二:把 Le 与 Le_f 混用。本课区间下 Le 是 0.85~0.90 而 Le_f 是 0.90~0.93,两者只差几个百分点——公式里写错哪一个,误差都会被当成"物性不确定度"吞掉,查不出来。⛔ 因此本课正文一律分列 Le 与 Le_f,不写不带下标的歧义形式。
回看 1.1 那张剖面图右侧的合并框:那一格写的就是这一步代换。前提没了,合并这一步就不成立——而不是"合并出来的结果差一点"。
1.4 Le_f 偏离 1 加一项修正就够了;真正让焓势法失效的是"表面状态变了"的三类场景
是什么。Le_f ≠ 1 时,主式带一个修正项:
q″ = (h_c/c_p,m)·[ (h_a − h_s,w) + (1/Le_f − 1)·h_fg·(d_a − d_s,w) ]
Le_f = 1 时后项归零,退回 1.3 那条主式。
为什么这样处置。有一种严格写法流传很广——"Le≈1 这个前提不成立,后面的式子一条都不能用"。这条太严了,会把可以用的工况一起挡掉。工程上 Le_f 稳定在 0.9 附近,加一项修正即可继续用;真正必须换方程的,是表面状态本身变了的场景,不是 Le_f 略偏离 1。这条口径是本课归纳的处置办法,不是某部标准的既有条款。
工程量级。取 Le_f = 0.91,把修正项与驱动焓差相比 —— (1/Le_f − 1)·h_fg·(d_a − d_s,w) ÷ (h_a − h_s,w),h_fg 取 10~12℃ 的 2470~2477 kJ/kg:
- 27℃ / 60%RH、T_w = 10℃(教学假设值):修正项占驱动焓差的 4.4%;
- 38℃ / 60%RH、T_w = 12℃(教学假设值):占 5.9%。
这个量级与湿工况 j 因子取值的不确定度同一档。⇒ 本课的口径是二选一,没有第三条路:要么把这一项加上,要么在交付物里显式声明"忽略了 Le_f 修正"。⛔ 不许默默按 Le_f = 1 算却不写。以上两个百分比是算例,随工况变,不是通用修正量。
真正让焓势法失效的三类场景。它们的共同特征是"表面不再是那个 0℃ 以上的连续凝露面":
- 表面温度接近或高于当地露点,盘管趋于全干。此时该回到温差驱动的常规 ε-NTU(干工况本体在 A3-01 换热器传热基础:ε-NTU 法与 LMTD 法)。⚠ 判据里的"当地露点"指该处空气状态对应的露点,它沿程会随析湿而下降,不是进风露点;这也正是 1.2 里 ξ < 1 报警的那种情况。
- 极低含湿量的干冷工况。潜热通道本身微不足道,合成驱动势里的相对误差被放大,此时双通道分算或干脆按全干处理更稳。
- 表面温度跌破 0℃,凝露转为凝华结霜。这一类整体交 A3-04 污垢、结霜与传热衰减机理,本课止于 0℃ 以上凝露。⚠ 结霜之后不是"精度下降"而是换了机理:析湿系数 ξ、饱和空气焓 h_s(T) 及其斜率这三套量全部失效,霜层还会自己长厚、自己改变通道截面。
易错点。最常见的一种是把"Le≈1 不成立"与"表面不再是湿的"混为一谈——前者加一项修正就够,后者必须换方程。判断顺序也不能颠倒:先判是不是出了本课射程(结霜),再判要不要用焓势法(全干/极低含湿量),最后才判修正项怎么取。
1.5 焓湿图那一侧的东西,本课一句不重讲
是什么。本课与前置课 A2-04 湿空气物理与焓湿图:空调负荷计算基石 都会出现含湿量差、装置露点、旁通因子这些词,但分工是清楚的:焓湿图判读、六族线、过程线、显热比拆分、湿源盘点、析水量估算——全部是 A2-04 湿空气物理与焓湿图:空调负荷计算基石 的射程;本课只讲"空气与湿表面之间那一层"的耦合机理,以及由芯体几何正算出总热与出风状态的那条通路。
为什么要在开讲就划这条线。两门课的词表重叠得厉害,不划线就会在第 1 讲复述一遍焓湿图"以便读者衔接",把课时耗在前置内容上,而真正缺的正算通路仍然没讲。本课假定你已经能在焓湿图上定点、能读出一个状态的焓、含湿量与露点——本讲用到的进风状态三元组(焓 / 含湿量 / 露点)就是这么来的,本课不再重讲怎么定点,只讲这些量进到表面耦合里之后发生了什么。
顺带把另外两条边界一并说清,免得后面反复解释:
- 凝水在流动侧的账不在本课。凝水滞留会抬高空气侧压降、进而影响风量,这条账(含"选风机必须按湿态"这条结论)在 C2-05 箱内蒸发器布置与冷凝水管理;本课只算换热侧,一句不重讲。
- 由实测反推的集总参数不是本课的产出。装置露点与旁通因子必须由焓差台架实测或经标定的仿真给出(这是 G2-03 蒸发器设计与出风温度均匀性 已经钉死的红线,不存在可以照搬的通用数值);本课交付的是由几何与表面温度正算的那条通路,两者是互校关系,第 5 讲会把它们对起来。
这一讲留下的三件事。第一,驱动势要从温差换成焓差,那个 h_s,w 到底怎么构造、它的口径是什么——第 2 讲。第二,ξ 一旦大于 1,翅片上的等效换热系数就被放大,翅片效率必然往下掉——第 3 讲。第三,这两处改完之后一起进效能法,才能正算出湿工况的出风干球温度与含湿量——第 4 讲。只改其中一处,就只是半个模型。
后面还有 4 讲正文 · 关键公式 · 案例拆解 · 常见误区 · 动手做