A2-03 两相流基础:沸腾、冷凝与流型
课程代码 A2-03 · 板块 A 基础理论与共性技术 / A2 流体、两相流与湿空气 时长 约 4 小时(5 讲 + 1 次两相压降与流型算例) 适合对象 仿真 / 系统工程师;做蒸发器、冷凝器、直冷板与微通道设计的人(必修) 前置 A1-04 蒸气压缩循环与 log p-h 图判读;A2-01 管网水力学(单相压降与并联分配);A1-07 热管理常用工质的热物性与查用 本讲义定位 讲师授课蓝本 / 学员自学讲义,是 A2-03 大纲的完整展开版
引言:贴壁还是分层,决定这段管子在不在干活
冷媒进了蒸发器,从一头的过冷液走到另一头的过热气,中间大半段是气液两相。这段管子里,冷媒是液是气、气相贴壁还是分层沉底,直接决定它是在换热还是白占空间。同样一根管、同样的进出口温度,内部流型不同,壁温可以差一大截——而只测进出口,你根本看不出来。
这门课把最基本也最要命的几件事讲清:用干度与空泡率描准任一截面的两相状态、在流型图上定流型;看懂核态沸腾—CHF—干涸的过渡,解释某段壁温为什么突跳;把冷凝换热里膜状冷凝、剪切主导与不凝气体讲透;再用均相与分相两种模型估两相压降、说清量级差从哪来。最后全部落到整车最现实的战场——制冷剂直冷电池板与微通道并联蒸发。一句话贯穿全课:两相流最难的从来不是传热,是均流。关联式再准,只要几条并联支路分配不均、某几路提前干涸,整块板的温差就崩了。
第 1 讲 两相流的描述量
单相流里一个流速、一个温度、一个压力就把截面说清了。两相流不行——同一截面上两种相、两种速度、两种密度,得先有专门的描述量,后面所有关联式才有入口参数。
1.1 干度 x:热力学干度与真实干度的分道
干度 x = ṁ_g / (ṁ_g + ṁ_l),是气相质量流量占总流量的份额,两相状态的第一描述量。冷媒在蒸发器里一路吸热,x 从进口接近 0 逐段涨到出口,所有两相换热与压降关联式,第一个要的输入就是当地 x。
这里埋着最常见的一个坑:热力学干度和真实干度不是一回事。热力学干度是拿当地焓做能量平衡、假设气液处于热力学平衡态反算的;真实干度是截面上实际的气相质量份额。两者对不上的根子在热力学非平衡:进口过冷段主流温度还低于饱和温度,按焓平衡算出的热力学干度是零甚至负值,可只要壁面过热到一定程度就早已起泡(过冷沸腾),近壁一层有真实气相。你若信了那个"零"、当成单相液体去算,就漏掉了起泡带来的换热强化和额外压降。出口干涸之后主流还有没蒸完的液滴、气相已过热,热力学干度又可能算出大于 1。落到操作上:进口过冷段、出口过热段要结合壁面过热度判断有没有起泡,别只信焓平衡算出的干度。
1.2 空泡率 α 与滑移比 S:为什么 x 很小 α 已经很大
空泡率 α = A_g / (A_g + A_l),截面上气相占的面积份额。它和干度的根本区别是:x 论质量,α 论面积(体积),两者隔着一个巨大的密度比。气相密度远小于液相:R134a、R1234yf 这类常规制冷剂在 0~10℃ 蒸发段,ρ_l/ρ_g 大致在 50~90 的量级——同样一份质量,气相要占的体积大几十倍。但这个比值随工质与压力变化很大,别当常数背:CO₂ 在 5℃ 亚临界下 ρ_l/ρ_g 只有 6 的量级(物性查用承接 A1-07),CO₂ 系统里「x 很小 α 就已很大」的悬殊会明显收窄。写成式子就是 α = 1/[1 + ((1−x)/x)·(ρ_g/ρ_l)·S],其中 S=u_g/u_l 是滑移比。⚠ 这里有个方向容易记反:撑开 α 的是密度比,滑移 S>1 是把 α 往回拉的——气相跑得快,意味着同样的质量流只需要更小的过流面积,所以 α 反而更小;均相(S=1)才是 α 的上限。本节下方那张图取的正是 S=1,故为上限曲线;滑移 S>1 时实际 α 略低于图示,但「x 小、α 已经很大」的结论不变。拿 R134a 5℃ 蒸发算一遍(ρ_l/ρ_g ≈ 75):S=1 时 x=0.1 得 α≈0.89,按 Zivi 关系取 S=(ρ_l/ρ_g)^(1/3)≈4.2 也还有 α≈0.67——所以哪怕 x 只有 0.1,α 仍在 0.7~0.9 的量级。
这个反直觉的事实直接决定两相密度算得对不对:两相段平均密度由 α 决定,不是由 x 决定。谁图省事拿一个小 x 当"九成液一成气"去估密度,算出来会严重偏高,压降、加速项跟着离谱。记牢:干度小不等于气少占地,凡涉及截面体积份额、平均密度,一律用 α,别用 x 顶替。
1.3 质量流密度 G:比"流速"更通用的自变量
质量流密度 G = ṁ / A_c,单位流道截面积上的质量流量,几乎所有两相关联式的核心自变量。为什么不用"流速"?因为两相里"流速"含糊——是气相、液相还是混合速度?而且 x 一路涨、混合物密度一路降,同一 G 下表观速度一路加快。G 就没这个歧义:质量流量除以流道面积,沿程守恒、与当地相态无关,是干净可传递的自变量。
这对选型有直接价值:把管径换掉——比如常规圆管换成更细的微通道扁管——别急着套原关联式。A_c 一变,同样质量流量下 G 会成倍跳,而每个关联式都在某个 G 范围内标定。换管径第一件事,是核 G 有没有跑出关联式的适用范围;跑出去了就得换标定范围合适的关联式或重新实验。这是仿真最容易忽略又最容易翻车的一步。
1.4 表面张力与 Bond / Weber 数:微通道里重力让位
到了微通道,一件在常规管里可忽略的事突然成主角:表面张力。用两个无量纲数把握——Bond 数衡量重力与表面张力之比(正比于特征尺寸平方),Weber 数衡量惯性与表面张力之比。常规大管尺寸大、Bond 数大,重力压过表面张力,气液在重力下分层;一旦流道缩到微通道量级,Bond 数骤降,重力让位于表面张力,表面张力把液相箍在流道周壁上,重力再没能力把气液拉分层。
这就是为什么微通道流型和常规管不同、不能照搬常规管流型图——常规管低流量下必现的分层流,在微通道里可能根本不出现。这是第 2 讲流型图"适用边界"的物理根源,也是微通道成为直冷板主流选择的原因之一:它天然规避了分层这个杀手(代价是均流更难,见第 5 讲)。
本讲小结:x 论质量、α 论面积,被密度比撑开、被滑移拉回(S=1 的均相值是 α 的上限),x 小 α 已经很大;G 比流速通用,换管径先核 G 出没出界;微通道里 Bond 数小、表面张力压过重力,流型规律与常规管不同。
后面还有 4 讲正文 · 关键公式 · 案例拆解 · 常见误区 · 动手做