G3-06 人字通道内的冷媒侧两相:流型、关联式与 CFD 能回答什么
| 项 | 内容 |
|---|---|
| 课程代码 | G3-06 |
| 所属板块 | G 总成级产品专题 · G3 Chiller 与板式换热器 |
| 难度 | L3 专家 |
| 时长 | 3.5 学时(6 讲 + 1 次两相计算与验证方案实操) |
| 适合对象 | 做板式换热器两相侧性能预测与异常定位的部件工程师(主要受众);承担板换两相 CFD 任务、需要先判断「该不该做、能回答什么」的仿真工程师;设计两相台架试验点与测点布置的试验工程师 |
| 前置 | A2-03 两相流基础:沸腾、冷凝与流型 两相流基础:沸腾、冷凝与流型 · G3-02 板式换热器钎焊结构与流道优化 板式换热器钎焊结构与流道优化 · G3-04 板式换热器正向设计指导:从工况谱到板型、板片数与流程的定量闭环 板式换热器正向设计指导。体系外前置需自备:两相流的基本参数定义(干度、空泡率、质量流速)与制冷剂物性库的使用 |
| 本讲义定位 | 这是一门判断力课,不是公式课。它不教你算某一条两相关联式,它教你判断:在人字波纹通道这个具体对象上,你手里那条关联式还成不成立、你要不要做 CFD、做出来的结果凭什么可信、以及你打算用什么试验去验证它。凡是能在别的课里学到的机理与方法本体,本课一律引用不重讲 |
引言:三门课把地基打好了,中间这块没人踩
全课程体系里,两相内容有两块地基。
第一块是 A2-03 两相流基础:沸腾、冷凝与流型。它把两相流的描述量、沸腾与冷凝的机理、流型与流型图讲透了:干度 x 是什么、空泡率 α 为什么远大于 x、两相压降为什么必须拆成三项、干涸是怎么发生的。这些是本课的全部入口参数。但 A2-03 的流型判据与流型图,几乎全都建立在圆管上——那是两相流研究的标准试验件,也是绝大多数公开流型图的载体。
第二块是 C1-04 制冷剂路径压降与均流设计。它讲并联管路之间的相分配:一根集管把两相分给十几根并联扁管,谁进气谁进液,怎么用 VOF 去看。那是一个通道与通道之间的问题。
而人字波纹通道,既不是圆管,也不是并联扁管。
它的截面沿流向周期性地收缩又扩张;相邻两张板片的波纹交叉,在接触点上把流道反复劈开又合并;当量直径只有毫米量级,而扰动强度远高于同当量直径的光滑圆管。你要研究的现象发生在一条通道内部,而这条通道的几何在整个体系里没有第二个对象与它相似。
于是就出现了这样一种常见局面:一位工程师读完了 A2-03,手上有一张漂亮的圆管流型图;他打开供应商给的板换资料,找到质量流速与干度,把点标到图上,读出「环状流」三个字,然后基于这三个字去解释一个高负荷工况下的过热度失控。这一整套动作里,每一步的操作都是对的,唯独第一步的前提不成立。而它不会报错——流型图不会因为你换了几何就拒绝作图。
本课要补的就是这一层。六讲的走向是:
- 第 1 讲说清几何到底改变了什么,以及为什么圆管流型图在这里只能当定性参考;
- 第 2 讲把冷媒侧沿程拆成相区,说明各区的主导机理与两端各自的失效形态;
- 第 3 讲处理关联式:形式族、四项射程、以及超出射程之后该怎么办;
- 第 4 讲给三条路线(一维分段 / 关联式加修正 / CFD)的分工判据;
- 第 5 讲讲板换两相 CFD 的三类困难,以及它们怎么表现为「算出来但不可信」;
- 第 6 讲讲验证:为什么整器指标不足以验证内部两相行为。
⚠ 本课有一条贯穿全篇的自限,现在就说清楚:本课不给任何两相关联式的系数,也不给流型转变的具体判据数值。原因不是藏私,是那些数都带着射程,而射程一旦脱离原文就不可核。本课给的是形式、去向与判断方法——形式和去向不会因为你换了板型而失效,系数会。
这门课不讲什么
| 内容 | 去哪门课 | 边界说明 |
|---|---|---|
| 两相流基础、圆管流型图、沸腾与冷凝机理 | A2-03 两相流基础:沸腾、冷凝与流型 | 本课把它的结论当输入,⛔ 不重讲。你必须已经能区分热力学干度与真实干度、能说清 α 与 x 的差别 |
| 两相仿真的方法本体:VOF 与欧拉双流体怎么选、相变源项模型、时间步与数值稳定性 | J3-04 两相流、沸腾与冷凝仿真 | 本课只讲「在板换这个对象上该不该做、做了能回答什么」,⛔ 不讲怎么做 |
| 并联管路之间的相分配及其 VOF 用法 | C1-04 制冷剂路径压降与均流设计 | 那是通道之间的问题,本课是通道内部的问题 |
| 端口到流道的分配不均与分配器结构 | G3-02 板式换热器钎焊结构与流道优化 第 3 讲 | 结构根因在那里,本课引用 |
| 过热度控制、MSS 与均流诊断的控制侧 | G3-03 水—冷媒换热的过热度与均流控制 | 控制回路怎么调不在本课,本课只在验证一节引用它的诊断手段 |
| 正向设计链条、单相关联式的射程管理 | G3-04 板式换热器正向设计指导:从工况谱到板型、板片数与流程的定量闭环 | 本课把射程管理的框架搬到两相侧,⛔ 不重讲框架本身 |
| 本体分段与过冷段的结构实现 | G3-05 多流程与过冷式板换:在一台换热器里做分段 | 本课在第 2 讲说明「为什么要分段」,怎么分是那门课 |
| 换热器芯体级 CFD 与多孔介质等效 | J3-01 换热器芯体 CFD 与多孔介质模型 | 本课只讲两相,且只讲通道级 |
本课的符号约定(含与兄弟课的口径对齐)
⛔⛔ 这一节请务必读完再往下。本课要跨三门课取用量,其中有两个符号在兄弟课里被赋予过不同含义,直接互代会出系统性偏差。
直接沿用兄弟课口径的符号:
| 符号 | 含义与口径 | 来源 |
|---|---|---|
| x | 干度,气相质量流量占总流量的份额,x = ṁ_g /(ṁ_g + ṁ_l) | 沿用 A2-03 两相流基础:沸腾、冷凝与流型 §1.1。⚠ 须分清热力学干度与真实干度,见下文 |
| α | 空泡率,截面上气相占的面积份额,α = A_g /(A_g + A_l) | 沿用 A2-03 两相流基础:沸腾、冷凝与流型 §1.2。⛔ 本课不用 α 表示任何换热系数 |
| S | 滑移比,S = u_g / u_l | 沿用 A2-03 两相流基础:沸腾、冷凝与流型 §1.2 |
| G | 质量流速,kg/(m²·s)。本课以单个板间通道的最小流通截面为基准 | 沿用 G3-04 板式换热器正向设计指导:从工况谱到板型、板片数与流程的定量闭环。⚠ 与 A2-03 的通用式 G = ṁ/A_c 一致,差别只在 A_c 取哪一个,见下文 |
| β | 波纹角,自流动方向量起。大 β = 硬板 H,小 β = 软板 L | 沿用 G3-02 板式换热器钎焊结构与流道优化 §1.1 |
| φ | 面积放大系数(波纹展开面积 ÷ 投影面积) | 沿用 G3-02 板式换热器钎焊结构与流道优化 |
| h | 表面换热系数,W/(m²·K) | 沿用 G3-01 Chiller(板式换热器)设计与选型 / G3-02 板式换热器钎焊结构与流道优化 / G3-04 板式换热器正向设计指导:从工况谱到板型、板片数与流程的定量闭环。⚠ 与 G3-03 水—冷媒换热的过热度与均流控制 的用法冲突,见下文 |
| U | 总传热系数,W/(m²·K) | 沿用 G3-01 Chiller(板式换热器)设计与选型 |
| n_pass | 流程数(单侧) | 沿用 G3-04 板式换热器正向设计指导:从工况谱到板型、板片数与流程的定量闭环 |
| Δp_fric / Δp_acc / Δp_grav | 两相压降的摩擦项 / 加速项 / 重力位项 | 沿用 A2-03 两相流基础:沸腾、冷凝与流型 的下标写法 |
| Φ_l² | 两相摩擦倍率,Φ_l² = (dp/dz)_tp /(dp/dz)_l | 沿用 A2-03 两相流基础:沸腾、冷凝与流型 |
本课首次定义或首次写足口径的符号:
| 符号 | 含义 | 口径(⛔ 这一栏才是重点) |
|---|---|---|
| i | 比焓,J/kg | 本课用 i,⛔ 不用 h——理由见下面那条跨课提醒 |
| h_tp | 两相侧当地换热系数,W/(m²·K) | 「当地」=某一个干度切片上的值,⛔ 不是两相段平均值 |
| h̄_tp | 两相段平均换热系数,W/(m²·K) | 按面积加权的平均。它是一个交付量,不是一个物理量——第 2 讲会说明它掩盖了什么 |
| d_h | 通道当量直径,m | ⛔ 板换的 d_h 有两套常见定义(按 2×板间距,或按 4×流通面积/湿周),两者在同一副板上可以差到接近一倍。凡引用一条关联式,必须先确认它用的是哪一套 |
| x_th | 热力学干度 | 拿当地焓做能量平衡、假设气液处于热力学平衡态反算的 x |
| x_1、x_2 | 两相段的进口干度与出口干度 | 用于界定「两相段」的两个端点 |
| x_dry | 干涸起始干度 | 蒸发侧 h_tp 由升转降的那个拐点对应的干度。它是平台相关量,取值见第 2 讲 |
⚠ 三处必须说清的跨课口径提醒——这三处如果不说,读者互代之后不会报错,只会得到一个错的数。
提醒一:h 在 G3-03 水—冷媒换热的过热度与均流控制 里被用作比焓,在本课里是换热系数。G3-03 水—冷媒换热的过热度与均流控制 §5.3 用 x = (h − h_f)/(h_g − h_f) 计算支路干度,那里的 h 是比焓;而 G3-01 Chiller(板式换热器)设计与选型、G3-02 板式换热器钎焊结构与流道优化、G3-04 板式换热器正向设计指导:从工况谱到板型、板片数与流程的定量闭环 三门课都用 h 表示表面换热系数。本课属于后者那条链(本课通篇在讨论换热系数随干度的变化),所以 h = 换热系数,而比焓在本课写作 i。⇒ 于是 G3-03 水—冷媒换热的过热度与均流控制 那条式子在本课的写法是 x_th = (i − i_f)/(i_g − i_f)。⛔ 读 G3-03 时不要把那里的 h 当成本课的 h,两者量纲都不同(J/kg 对 W/(m²·K)),一旦代错量纲会当场对不上,这是唯一一处「代错了会立刻暴露」的地方——反倒是最安全的一处。
提醒二:G 的分母必须写明取哪一个截面。A2-03 两相流基础:沸腾、冷凝与流型 给的是通用定义 G = ṁ/A_c,它没有、也不需要指定 A_c;G3-04 板式换热器正向设计指导:从工况谱到板型、板片数与流程的定量闭环 在板换语境下把 A_c 钉死为单个板间通道的最小流通截面,本课沿用这一条。⛔⛔ 危险在于板换有三个都能被叫做「流通面积」的东西:端口面积、单通道最小截面、以及全部并联通道的截面之和。若拿端口面积或整器截面之和当 A_c,算出来的 G 与本课口径可能差一个数量级(差的倍数 = 并联通道数 ÷ n_pass 这一类因子)。而两相关联式的射程恰恰是按 G 来划的 ⇒ 分母取错,你会以为工况落在射程内,其实早已出界。⇒ 凡在交付物上写 G,必须同时写明分母是哪一个截面。
提醒三:x 有两种,本课默认指真实干度,凡用焓平衡算出来的一律写 x_th。A2-03 两相流基础:沸腾、冷凝与流型 §1.1 已经把这件事讲透:进口过冷段主流还没到饱和温度、按焓平衡算出的 x_th 是零甚至负值,而壁面早已起泡;出口干涸之后主流还有液滴而气相已过热,x_th 又可能算出大于 1。板换里这两处非平衡区比圆管更难忽略,原因在第 1 讲——波纹通道的强制混合会让近壁与主流之间的交换比圆管快,非平衡的表现形式与圆管不同。⇒ 本课凡写 x 而未加下标处,指真实干度;凡是从焓平衡来的,一律写 x_th。
第 1 讲 波纹通道不是圆管:几何改变了什么
这一讲只做一件事:把「板换两相不能照搬圆管结论」这句话,从一句口号变成一条可以说出理由、也可以说出射程的判断。
1.1 三个几何特征,每一个都取消了圆管的一条前提
先把人字波纹通道的几何摊开。两张相邻板片,波纹方向相对于主流各偏一个角度 β,扣在一起,就形成了一条被反复劈开又合并的流道网络。它与圆管的差别可以归成三条,每一条都取消了圆管两相研究的一条隐含前提。
特征一:截面沿流向周期性收缩与扩张。圆管的截面沿程恒定,所以「当地质量流速 G」这个量在圆管里是明确的——同一个截面上处处一样,沿程也不变。人字通道里,截面在一个波纹节距之内就走完一次收缩—扩张,⇒ 「当地速度」这个量在一个波纹节距内变化的幅度,可能与整段流道上干度变化带来的速度变化同一量级。你写下的那个 G,是按某个约定截面(本课约定为最小流通截面)算出来的名义值,⛔ 不是通道里任何一处的真实速度。
特征二:接触点把流道反复劈开又合并。两张板的波纹交叉处互相顶住(钎焊板换在这些点上还有钎缝),流体在这里被迫分流、绕行、再汇合。⇒ 一段液膜刚在某个波谷里铺开,走过一个接触点就被撕成两股。这一条取消的前提是相分布的连续性:圆管流型的定义(分层、弹状、环状……)默认相分布在流向上有一定的连续性,才能说「这一段是环状流」。人字通道里,相分布的特征长度与波纹节距同量级 ⇒ 「这一段是什么流型」这句话本身就缺少一个稳定的观察尺度。
特征三:当量直径小而扰动强。板换通道的当量直径处在毫米量级,而波纹带来的扰动强度远高于同当量直径的光滑圆管——这正是 G3-02 板式换热器钎焊结构与流道优化 §1.3 讲的「波纹为什么能在低 Re 下就制造湍流」。⇒ 单位体积上的剪切远大于圆管。
⚠ 关于当量直径,这里必须插一句会直接影响你算出来的数的话。板换的 d_h 在文献里有两套常见定义:一套取2 倍板间距(把通道近似成两块平行板之间的缝隙),另一套取 4×流通面积/湿周(严格的水力当量直径定义,会把波纹展开带来的湿周增加算进去)。两者在同一副板上可以差到接近一倍,因为波纹把湿周放大了大约 φ 倍而流通面积几乎没变。⇒ Re 与 G 的计算结果、以及关联式射程里那一栏「适用 d_h 范围」,都会跟着差。⛔ 引用任何一条板换关联式之前,必须先在原文里确认它的 d_h 是哪一套定义;确认不了,这条关联式就不该用作设计依据。
d_h 的量级(行业典型区间):车用与商用钎焊板换的通道当量直径典型落在毫米量级,⚠ 具体值随板型、板间距与所采用的定义而变,⛔ 本课不给具体数——它由你手上那副板的几何直接算得,且必须与关联式的定义口径一致。
1.2 惯性、表面张力与剪切的相对权重变了
两相流型由几个力的竞争决定:惯性(把相打散、把液膜甩上壁)、重力(让密度大的相沉下去)、表面张力(让界面收成最小面积、维持液桥与液塞)、以及相间剪切(气核拖动液膜)。圆管流型图(无论以 G 与 x 为坐标,还是以某种无量纲群为坐标)本质上就是这几个力的比值在某个几何上的分区结果。
人字通道把这几个比值同时改了,而且不是同比例改的:
- d_h 变小 ⇒ 表面张力的相对权重上升(这是所有小通道两相研究的共识方向:通道越小,重力越退居次位,表面张力越显眼);
- 扰动变强 ⇒ 惯性与剪切的权重上升;
- 接触点强制分流 ⇒ 引入了一个圆管里根本不存在的机制:相分布被几何周期性地重置。
⛔⛔ 请注意最后一条的性质与前两条不同。前两条是「同一个坐标系里点挪了位置」,理论上还有救——你可以找一张覆盖更小 d_h、更高扰动的流型图。第三条是坐标系里没有这条轴。无论你怎么调 G 与 x,圆管流型图上都没有一个位置对应「相分布每隔一个波纹节距被重置一次」。⇒ 这就是为什么本课的判据是「只能当定性参考」而不是「要做修正」——修正的前提是差别可以映射到坐标轴上,而这里的主要差别映射不上去。
1.3 强制混合:高换热系数的结构原因,也是流型难移植的原因
这里有一个值得停一下的观察:板换两相换热系数高,和板换两相流型难判,是同一件事的两面。
板式换热器之所以在两相侧表现好,很大程度上就是因为波纹几何把气液强制混合了:液相不容易在某一侧长时间堆积成厚膜,气核也不容易长时间占住通道中心。液膜被反复撕开重铺 ⇒ 平均膜厚小 ⇒ 导热热阻小 ⇒ h_tp 高。这是好事。
而同一个机制,让「这一段是什么流型」失去了稳定答案。⇒ 你不能一边享受强制混合带来的换热收益,一边指望用一个建立在「相分布稳定」之上的工具去描述它。
⚠ 这条观察在工程上有一个直接推论,值得写进你的判断清单:当有人告诉你某台板换「两相侧换热系数比手册值低很多」时,先别急着怀疑关联式,先问强制混合是不是被削弱了。削弱强制混合的常见原因有三类:① 实际质量流速远低于设计值(低负荷或分配不均导致某些通道流量极低);② 板型换成了低 β 的软板而关联式仍按原板型;③ 通道内有异物或钎料堆积,把局部流道堵成了近似光滑的窄缝。这三类都不是关联式的问题,而是你以为的几何和实际的几何不一样了。
1.4 重力方向:竖直安装的上升流与下降流不是同一件事
板式换热器在整车与储能系统里通常竖直安装——端口在上下两端、流道沿重力方向。而绝大多数公开的圆管流型图,要么是水平管的(这一类最多),要么是垂直管的,而且垂直管的图往往只针对上升流。⇒ 光是「用哪张图」这个问题,在板换上就已经不唯一。
竖直流道里,上升流与下降流的行为差别很实在:
- 上升流:重力与主流方向相反,液相被拖着往上走,需要足够的气相剪切才能把液膜带上去。低质量流速时液相回落,会出现局部的液体积聚甚至流动不稳定。重力位压降 Δp_grav 为正(消耗压头)。
- 下降流:重力与主流同向,液相自然下落,液膜倾向于变薄且更快。重力位压降 Δp_grav 为负(提供压头)。低质量流速下的降膜形态与上升流完全不同。
⇒ 三条工程含义:
- 同一台板换,多流程排布下会同时存在上升段与下降段(这正是 G3-05 多流程与过冷式板换:在一台换热器里做分段 分段设计要面对的事),⇒ 各段的两相行为、压降符号都不同,一维分段计算必须按段给出流向;
- Δp_grav 在两相段不可忽略,且它需要空泡率 α 才能算(重力位压降由两相平均密度决定,而平均密度由 α 决定,不是由 x 决定——这是 A2-03 两相流基础:沸腾、冷凝与流型 §1.2 反复强调的那条);
- 安装姿态变更是一次真正的设计变更。把一台原本竖直安装的板换改成水平安装,两相侧的行为会变,而它不会触发任何一张变更单——包络对得上、接口对得上、单相侧性能几乎不变。⛔ 这是本课列进常见误区的一条。
1.5 判据:圆管流型图在这里只能当定性参考
把这一讲收成一条可执行的判据:
判据 1:把圆管流型图直接画到板换工况上,得到的分区没有依据。它可以作为定性参考(用来提示「这段大概处在低干度还是高干度、大概偏向液主导还是气主导」),⛔ 不能作为设计判据——不能据它判定某段是否会干涸、不能据它选关联式、更不能据它下「结构没问题」的结论。
⚠ 那么该用什么?三条替代路径,按可信度从高到低:
- 用本板型的台架数据直接给出行为边界:在你关心的工况范围上扫 h_tp 与 Δp 对 x、G 的响应,把拐点位置量出来。这是唯一一条能给出本板型结论的路。代价是台架时间。
- 用以板换为载体标定的公开关联式与流型判据,并严格核对四项射程(第 3 讲)。代价是你必须找得到覆盖你板型的那一条,而这一条常常不存在。
- 用圆管结论做定性提示,并在交付物上显式声明它是定性的。这条路唯一的价值是让你不至于完全没有方向感;⛔ 它的结论不得进入设计冻结。
后面还有 5 讲正文 · 关键公式 · 案例拆解 · 常见误区 · 动手做