TMS BOOK · ACADEMY 讲义

人字通道内的冷媒侧两相:流型、关联式与 CFD 能回答什么

大纲的完整展开版——讲师授课蓝本 / 学员自学材料

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 下就制造湍流」。⇒ 单位体积上的剪切远大于圆管。

人字波纹通道与圆管的几何对比示意图。画面分左右两栏,中间一条竖直分隔线。左栏标题「圆管」:画一根水平的长直管,管壁为两条平行水平线,管内均匀,沿流向画三个等大的虚线圆表示三个截面,三个截面下方各标注「A_c 恒定」;管内画一条水平的液膜带贴在下壁、上部留白表示气核,标注「相分布沿流向连续 ⇒ 流型可定义」;左栏底部一行小字写「截面恒定 · 无内部障碍 · 扰动来自壁面剪切」。右栏标题「人字波纹通道」:画一条水平流道,上下边界为两条相位相反的正弦波曲线,两条曲线周期性相交,把流道分成一串首尾相接的透镜形腔体,相交处即通道最窄处;在其中三个相交处画实心橙色圆点,框下一行小字写「红点=接触点(钎焊点),流被劈开又合并;A_c 沿流向最小与最大交替」;两个腔体中部各画一小段淡蓝色短条代表被反复撕开重铺的液膜;框下另两行写「相分布特征长度 ≈ 波纹节距」与「⇒『这一段是什么流型』缺少稳定的观察尺度」;右栏底部方框内三行写「截面沿流向周期性收缩扩张」「接触点反复劈分与合并」「扰动来自几何强制混合」,方框下一行小字写「当量直径毫米量级,而扰动远强于同径光滑圆管」。画面最下方横跨两栏一条横向结论条,写着「三条几何差别各取消圆管的一条前提:截面恒定 / 相分布连续 / 扰动机制」。本图为结构示意,⛔ 不得据图取值。
图1 圆管与人字波纹通道的三条几何差别:截面是否沿程恒定、有无接触点反复劈分、扰动来自何处。这张图要读出的结论是:把圆管结论搬到板换上失效的不是「精度」,而是流型这个概念赖以成立的前提——相分布的连续性与稳定的观察尺度在这里都不成立。本图为结构示意,⛔ 不得据图取值。

⚠ 关于当量直径,这里必须插一句会直接影响你算出来的数的话。板换的 d_h 在文献里有两套常见定义:一套取2 倍板间距(把通道近似成两块平行板之间的缝隙),另一套取 4×流通面积/湿周(严格的水力当量直径定义,会把波纹展开带来的湿周增加算进去)。两者在同一副板上可以差到接近一倍,因为波纹把湿周放大了大约 φ 倍而流通面积几乎没变。⇒ Re 与 G 的计算结果、以及关联式射程里那一栏「适用 d_h 范围」,都会跟着差。⛔ 引用任何一条板换关联式之前,必须先在原文里确认它的 d_h 是哪一套定义;确认不了,这条关联式就不该用作设计依据。

d_h 的量级(行业典型区间):车用与商用钎焊板换的通道当量直径典型落在毫米量级,⚠ 具体值随板型、板间距与所采用的定义而变,⛔ 本课不给具体数——它由你手上那副板的几何直接算得,且必须与关联式的定义口径一致。

1.2 惯性、表面张力与剪切的相对权重变了

两相流型由几个力的竞争决定:惯性(把相打散、把液膜甩上壁)、重力(让密度大的相沉下去)、表面张力(让界面收成最小面积、维持液桥与液塞)、以及相间剪切(气核拖动液膜)。圆管流型图(无论以 G 与 x 为坐标,还是以某种无量纲群为坐标)本质上就是这几个力的比值在某个几何上的分区结果。

人字通道把这几个比值同时改了,而且不是同比例改的:

  • d_h 变小 ⇒ 表面张力的相对权重上升(这是所有小通道两相研究的共识方向:通道越小,重力越退居次位,表面张力越显眼);
  • 扰动变强 ⇒ 惯性与剪切的权重上升;
  • 接触点强制分流 ⇒ 引入了一个圆管里根本不存在的机制:相分布被几何周期性地重置。

⛔⛔ 请注意最后一条的性质与前两条不同。前两条是「同一个坐标系里点挪了位置」,理论上还有救——你可以找一张覆盖更小 d_h、更高扰动的流型图。第三条是坐标系里没有这条轴。无论你怎么调 G 与 x,圆管流型图上都没有一个位置对应「相分布每隔一个波纹节距被重置一次」。⇒ 这就是为什么本课的判据是「只能当定性参考」而不是「要做修正」——修正的前提是差别可以映射到坐标轴上,而这里的主要差别映射不上去。

圆管流型图不可直接移植到板换的示意图。画面主体是一张二维流型图坐标系,横轴标注「干度 x(0 到 1)」,纵轴标注「质量流速 G」。坐标区内用三条弯曲的实线分界线划出四个区域,四个区域从左下到右上依次用浅色底纹标注「分层 / 波状」「弹状 / 塞状」「环状」「雾状 / 弥散」,坐标区左上角有一个小标签写「标定载体:光滑圆管」。图中画一条粗虚线折线代表某台板换沿程的工况轨迹,从左侧中部起点标注「进口 x_1」向右上延伸到右侧标注「出口 x_2」,轨迹整体穿过弹状区与环状区;轨迹上叠加一个半透明的灰色阴影带,阴影带明显比轨迹线宽,标注「板换上这条线读出来的分区没有依据」。坐标区右侧引出三个说明框,自上而下:第一框写「d_h 更小 ⇒ 表面张力权重上升 ⇒ 分界线整体移动(可映射)」,第二框写「扰动更强 ⇒ 惯性与剪切权重上升 ⇒ 分界线整体移动(可映射)」,第三框加粗边框写「接触点周期性重置相分布 ⇒ 这张图上没有对应的坐标轴(⛔ 不可映射)」,第三框用一条带叉号的箭头指向坐标区。画面底部一条横向结论条写「前两条差别可以靠换一张标定更近的图来缓解;第三条差别在这个坐标系里根本不出现 ⇒ 只能当定性参考,⛔ 不能作为设计判据」。本图为定性示意,⛔ 不得据图取值。
图2 圆管流型图上的三类差别:前两类只是把分界线整体挪动(还能靠换图缓解),第三类在这张图的坐标系里没有对应的轴。这张图要读出的结论是:板换流型判断的困难不是精度不够,是坐标系不全 ⇒ 不存在「加一个修正系数」这条出路。本图为定性示意,⛔ 不得据图取值。

1.3 强制混合:高换热系数的结构原因,也是流型难移植的原因

这里有一个值得停一下的观察:板换两相换热系数高,和板换两相流型难判,是同一件事的两面。

板式换热器之所以在两相侧表现好,很大程度上就是因为波纹几何把气液强制混合了:液相不容易在某一侧长时间堆积成厚膜,气核也不容易长时间占住通道中心。液膜被反复撕开重铺 ⇒ 平均膜厚小 ⇒ 导热热阻小 ⇒ h_tp 高。这是好事。

而同一个机制,让「这一段是什么流型」失去了稳定答案。⇒ 你不能一边享受强制混合带来的换热收益,一边指望用一个建立在「相分布稳定」之上的工具去描述它。

⚠ 这条观察在工程上有一个直接推论,值得写进你的判断清单:当有人告诉你某台板换「两相侧换热系数比手册值低很多」时,先别急着怀疑关联式,先问强制混合是不是被削弱了。削弱强制混合的常见原因有三类:① 实际质量流速远低于设计值(低负荷或分配不均导致某些通道流量极低);② 板型换成了低 β 的软板而关联式仍按原板型;③ 通道内有异物或钎料堆积,把局部流道堵成了近似光滑的窄缝。这三类都不是关联式的问题,而是你以为的几何和实际的几何不一样了。

1.4 重力方向:竖直安装的上升流与下降流不是同一件事

板式换热器在整车与储能系统里通常竖直安装——端口在上下两端、流道沿重力方向。而绝大多数公开的圆管流型图,要么是水平管的(这一类最多),要么是垂直管的,而且垂直管的图往往只针对上升流。⇒ 光是「用哪张图」这个问题,在板换上就已经不唯一。

竖直流道里,上升流与下降流的行为差别很实在:

  • 上升流:重力与主流方向相反,液相被拖着往上走,需要足够的气相剪切才能把液膜带上去。低质量流速时液相回落,会出现局部的液体积聚甚至流动不稳定。重力位压降 Δp_grav 为正(消耗压头)。
  • 下降流:重力与主流同向,液相自然下落,液膜倾向于变薄且更快。重力位压降 Δp_grav 为负(提供压头)。低质量流速下的降膜形态与上升流完全不同。

⇒ 三条工程含义:

  1. 同一台板换,多流程排布下会同时存在上升段与下降段(这正是 G3-05 多流程与过冷式板换:在一台换热器里做分段 分段设计要面对的事),⇒ 各段的两相行为、压降符号都不同,一维分段计算必须按段给出流向;
  2. Δp_grav 在两相段不可忽略,且它需要空泡率 α 才能算(重力位压降由两相平均密度决定,而平均密度由 α 决定,不是由 x 决定——这是 A2-03 两相流基础:沸腾、冷凝与流型 §1.2 反复强调的那条);
  3. 安装姿态变更是一次真正的设计变更。把一台原本竖直安装的板换改成水平安装,两相侧的行为会变,而它不会触发任何一张变更单——包络对得上、接口对得上、单相侧性能几乎不变。⛔ 这是本课列进常见误区的一条。

1.5 判据:圆管流型图在这里只能当定性参考

把这一讲收成一条可执行的判据:

判据 1:把圆管流型图直接画到板换工况上,得到的分区没有依据。它可以作为定性参考(用来提示「这段大概处在低干度还是高干度、大概偏向液主导还是气主导」),⛔ 不能作为设计判据——不能据它判定某段是否会干涸、不能据它选关联式、更不能据它下「结构没问题」的结论。

⚠ 那么该用什么?三条替代路径,按可信度从高到低:

  1. 用本板型的台架数据直接给出行为边界:在你关心的工况范围上扫 h_tp 与 Δp 对 x、G 的响应,把拐点位置量出来。这是唯一一条能给出本板型结论的路。代价是台架时间。
  2. 用以板换为载体标定的公开关联式与流型判据,并严格核对四项射程(第 3 讲)。代价是你必须找得到覆盖你板型的那一条,而这一条常常不存在。
  3. 用圆管结论做定性提示,并在交付物上显式声明它是定性的。这条路唯一的价值是让你不至于完全没有方向感;⛔ 它的结论不得进入设计冻结。

后面还有 5 讲正文 · 关键公式 · 案例拆解 · 常见误区 · 动手做

会员专属

后续为会员深水区内容——四库数据与深度拆解。

查看会员方案