TMS BOOK · ACADEMY 讲义

电池包液冷流道 CFD 与流量均匀性

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

J3-02 电池包液冷流道 CFD 与流量均匀性

课程代码 J3-02 · 板块 J 仿真分析 CAE / 换热器/部件级 CFD 时长 约 4.0 小时(8 讲) 前置 J3-01 换热器芯体 CFD 与多孔介质模型;A2-01 管网水力学与流量分配;D1-02 电芯温度不均与温差控制目标 本讲义定位 讲师授课蓝本 / 学员自学讲义,是 J3-02 大纲的完整展开版


按任务导航

你手上的事 直接看 还要配
定节流孔该加在哪一端 第 1 讲:画两条集管静压沿程曲线判抢流端 第 6 讲(手段按「改压差链哪一段」选)· 第 3 讲(入口给错,这条判据是哑的)
选建模层级并列全边界清单(开算前) 第 2 讲:七种层级、逐支路 Re、共轭传热域 第 3 讲(十二类边界+跨模型交接口径)
写一条流量均匀度的验收判据 第 4 讲:γ 的三条口径、它为什么只是过程指标 第 5 讲(温差目标→流量离散→γ 的反解链)
审一份「已调匀」的仿真报告 第 7 讲:验证与确认两侧的证据清单、解析解验设置 第 4 讲(γ 与压降的取值口径)· 第 5 讲(温差归因两级)
判额定点调匀的方案在低温与老化后还成不成立 第 8 讲:工况漂移四域、阻力律指数变了方案跟着变 第 6 讲(节流孔权重随工况变)· 第 7 讲(敏感性核查)

本课不覆盖:流道拓扑与结构优化(→ J3-03 冷板/流道拓扑优化)· 冷板部件级压降与温差设计(→ D2-01 电池液冷冷板设计:流道、压降与温差均匀性)· 芯体多孔介质等效与 j/f 标定(→ J3-01 换热器芯体 CFD 与多孔介质模型)· 系统管网与泵选型(→ A2-01 管网水力学:沿程/局部压降与流量分配 J1-06 电池热管理一维建模与快充温升仿真)· 电芯产热模型与温差控制目标(→ D1-01 动力电池产热机理与产热速率建模 D1-02 电芯温度不均与温差控制目标)· 接触热阻与 TIM(→ D2-05 冷板与电芯的接触热阻与 TIM 应用)· 包壳外侧空气边界(→ J2-03 整车外流场与气动—热耦合仿真)· 分配台架与试验方法(→ L2-03 冷却系统台架与流量分配试验)· 两相/沸腾/直冷(本课只算单相液冷)· 分区阀控制律(本课只给接口)· 求解器面板与选项(→ 求解器现行版本手册)· 整包热扩散与热安全试验(→ 归 GB 38031 的射程,以现行目录为准;⛔ 它不为流量分配与电芯温差提供对标判据);以及任何一个具体的目标值、参量阈值与本底限值——γ 限值、电芯温差目标与最热点温升上限、产热值、压降与泵功预算、节流孔孔径、包壳外侧 h_ext、R_TIM、网格量与 GCI 目标、低温工况温度点,一律由本项目自己的设计验证计划与所引标准现行版本给出,⛔ 本课一个数都不给(课内带数算例全为教学假设值,⛔ 禁止照抄)。

引言:一条经常对的口诀,比一条经常错的贵得多

一个十来片冷板并联的电池包,进出口开在异侧(Z 型)。台架标定下来,靠近进液口那几片下面的电芯明显更热。会上很快就有了结论,而且几乎是异口同声的一句:「Z 型远端流量多」——远端抢流、近端吃亏,那就在远端支路上加节流孔,把它多吃的那份逼回来。

这句话就是本课要拆掉的第一件东西。它的麻烦不在于它错,而在于它经常对:在常见的短粗集管构型上,它给出的方向确实是对的。于是它先赢得信任,再在某一个平台上安静地反过来。一条经常错的口诀会很快被人丢掉;一条经常对的口诀会被人一直用到它出事的那一天,而且出事的时候没有人怀疑它。

它为什么会被读反,其实是可以指出来的:这条经验的完整形式带着一个前提——集管动压回升主导时远端可用压差大、流量偏多;沿程摩擦主导时近端流量偏多。人记住的是「Z 型 / U 型」这两个字母,忘掉的是前面那半句「谁在主导」。而字母是几何图纸上一眼就能看见的,主导项不画曲线看不见。

一个不会报错的错误动作

拿到这个 Z 型包,不画进出集管的静压沿程曲线,直接判「远端抢流」,于是把节流孔加在远端。

如果这个包实际是沿程摩擦主导、近端在抢流,这一刀就正好加在流量本来就少的那一端。被抬阻的那支流量进一步下降,其余支路按并联关系分走这份流量 ⇒ 分配的离散度不降反升;同时系统总压降上升、泵功跟着上升——三个指标一起朝坏的方向走

而最贵的不是这一步,是下一步:「我是按 Z 型的规律加的节流孔」这句话会让人判定方向是对的、只是量不够,于是继续收紧孔径,朝着错误的方向迭代。⚠ 这类错在 CFD 上不会以不收敛的形式暴露:残差照掉、云图照出、总压降还是个看起来合理的数。它没有报错的位置。

判方向的四步,今天就能对手上的包用

先看你的算例本身:计算域含不含集管?入口是不是单一总进口给总流量?——若各支路入口各给了一份流量(或用周期段/单板模型代替全包),下面三步全是哑的:分配已经被当成输入给掉了,你会得到一张完美平齐的柱状图,而它什么也没证明(第 2、3 讲)。

沿集管轴向逐站取截面平均静压,画出进(分流)集管的 p_div(x) 与出(合流)集管的 p_com(x) 两条曲线。⚠ 方向全课钉死:x=0 与支路序号 i=1 都取在进液口一侧,「近端」= i 小的一侧。

看两条曲线的差 Δp_b(x) 沿 x 是增还是减:增 ⇒ 远端可用压差大 ⇒ 远端抢流;减 ⇒ 近端抢流。方向由这一步给出,布置形式给不出——U 型与 Z 型只决定这两条曲线怎么叠(异侧进出时两条曲线的沿程变化相加,同侧进出时部分抵消),⛔ 不决定叠出来的朝向。

节流孔加在 ③ 判出来的那一端——方位是结果,⛔ 不是判据。

还有一个连 CFD 都不用跑的先手估算:进集管里,支路取流带来的动压回升项与沿程摩擦项都正比于 ½ρu_man²(u_man = 集管的截面平均流速),所以两者之比与总流量一阶无关,只由集管的长径比和取流回升的强弱定。⇒ 集管越细长,摩擦越占优(近端抢流倾向);越短粗,回升越占优(远端抢流倾向)。 第 1 讲把这条比值写成一个两行可复算的式子,并给两组只差集管长径比的教学算例——同样是 Z 型,两组判出来的抢流端相反

把这件事推到底:本课的两条主线

主线一:并联分配是「集管压降」与「支路压降」之间的份额之争。 写成一个可判真假的量:Λ ≡ 集管侧压力沿程变化的幅度 / 单支路压降(★ 本课归纳的判据量,⛔ 不是任何标准或手册既有的定义;分子分母各由什么构成,第 1 讲与第 4 讲逐项拆开)。Λ→0 时各支路两端的压差趋于相同、分配趋于均匀;Λ 越大越不均。⇒ 从此每一个均流方案都要回答同一句话:你压低的是 Λ 的分子哪一项,还是抬高了分母? 变截面集管与加大集管截面压的是分子,节流孔与改流道截面抬的是分母,分区供液从两条路径一起压分子(⛔ 分母不动),导流片治的则是分子里的一个局部畸变项。答不出这一句的方案,本课认为它还没有被理解,只是被试出来的。⚠ Λ 只判量级与方向,⛔ 不得拿它反算各支路流量。

主线二:流量不均本身不是缺陷,只有被翻译成电芯温差的那一部分才是。 流量不均匀度 γ(口径见第 4 讲)是过程指标——同样的 γ,产热分布不同、电芯在包内的位置不同,最热电芯的温度可以差很多。⇒ 验收判据必须落在电芯温差与最热点上,⛔ 不落在 γ 上。而这次翻译有两级,两级的口径⛔ 不能混着说,否则会得出「流量减半、温升翻倍」这一类系统性高估(据此定的流量余量偏大、泵功跟着偏大);第 5 讲给出正确的放大倍数、它的代数式与它失效的两个边界。

⚠ 顺带把一件读者一定会去找的事说明白:流量不均匀度的限值与电芯温差目标都没有公开统一的标准,多由主机厂/电池厂的设计验证计划内控规定 ⇒ ⛔ 本课不给 γ 的任何限值,也⛔ 不许从别的项目抄一个来当验收线。本课给的是替代物:从本项目的温差目标往回推到允许的流量离散、再推到允许的 γ 的那条反解链(第 5 讲)。

本课交出什么

输入是一份已经有几何的整包液冷方案(集管走向、冷板数量与流道形式、进出口位置);输出是三份带口径的东西:① 各支路流量分布与不均匀度 γ、γ_R;② 最热电芯的位置与全包电芯温差 ΔT_cell;③ 系统压降 Δp_sys 与水力功率 P_hyd 的拆账(三者取的是哪一层、哪一段、哪一种功率,口径分别在第 4、5 讲钉死)。三份合起来才构成一次可复用的结论,第 1 讲逐份说明缺了哪一份下游具体做不了什么。

⇒ 全课的落点可以用一句话概括:把那条口诀,换成 Δp_b(x) 这条曲线的斜率符号。 口诀给你一个答案,曲线给你一个判据——判据在你换一个包、换一个工况、换一个供应商的时候仍然成立,答案不会。

第 1 讲 方向看两条静压曲线的差,量级看阻力比 Λ

这一讲不讲网格、不讲边界条件、不讲求解器面板——那是第 2、3 讲的事。它只做一件事:把「几十片冷板并联,流量为什么分不匀、往哪一端偏」这个问题,从一句可以背下来的口诀,换成一条你能在自己的算例上画出来、而且能判真假的判据。这一步做不到,后面七讲都是在一个错的方向上做精细化:网格再细、模型再贵,节流孔只要加错那一端,流量不均匀度、系统压降、泵功会一起朝坏的方向走,而这类错在 CFD 上不会以不收敛的形式暴露——残差照掉、云图照出、总压降还是个看起来合理的数。

这一讲要建立三个判断:① 分配由每一支路自己的两端压差决定,而那个压差沿集管是变的 ⇒ 各支路的阻力特性完全一样,流量也可以差很多,不均的根在集管、⛔ 不在支路;② 偏向的方向由进、出两条集管的静压曲线之差 Δp_b(x) 沿程是增是减决定,布置形式(U 型/Z 型)只决定这两条曲线怎么叠,⛔ 不决定叠出来的方向;③ 不均的量级由阻力比 Λ = 集管侧压力沿程变化幅度 ÷ 单支路压降 决定 ⇒ 第 6 讲要排的每一个均流手段,从此都能被问同一句话:你压低的是 Λ 的分子哪一项,还是抬高了分母?答不出这句话的方案,本课认为它还没被理解,只是被试出来的。

本讲开头先钉三条约定(全课通用,⛔ 中途不换写法) · 方向:沿集管的轴向坐标 x(m)取 x=0 在进液口一侧,支路序号 i=1..N 也从离进液口最近的那一支起数;「近端」= i 小的一侧、「远端」= i 大的一侧。⛔ 全课的正文、图与图注一律照此——方向一反,这一讲的结论会整条跟着反。 · 流量与热:支路流量一律写质量流量 ṁ_i(kg/s,i 为支路序号,N 为并联支路数);Q 在全课只表热量,⛔ 一律不用 Q 表流量;确需体积流量时写 V̇ 并注明温度(密度随温度变)。 · 摩擦因子:全课压降式一律用 Darcy 摩擦因子 f_D:Δp = f_D·(L/D_h)·½ρu²,⛔ 全课不出现不标类型的裸 f。⚠ 本讲会用到 f_D、水力直径 D_h、流量不均匀度 γ、系统压降 Δp_sys、水力功率 P_hyd 这几个量,但只用它们的定性含义;它们的取值口径(标准差取哪一种、分母是什么、取值区间从哪到哪、与 Fanning 因子怎么换算)在第 4 讲统一钉死,⛔ 本讲不另立一套。

1.1 本课的坐标系:输入是一份已有几何的方案,输出是三份带口径的交付物

是什么。本课的输入是一份已经有几何的整包液冷方案——集管走向与截面、冷板数量与流道形式、进出口位置都已经画出来了。输出是三份可以被下游直接取用的交付物:

  1. 各支路流量分布(柱状图)与流量不均匀度 γ、γ_R(标准差归一与极差归一两个数,口径在第 4 讲钉死);
  2. 最热电芯的位置与全包电芯温差 ΔT_cell(=同一时刻全包电芯最高温 − 最低温,这一层的取法在第 5 讲钉死);
  3. 系统压降 Δp_sys 与水力功率 P_hyd 的账——拆到集管、支路、局部损失各占多少。

为什么是三份而不是一份。这三样合起来才构成一次可复用的结论,少一份下游就有一步做不了:只交第 ① 份,γ 是过程指标,同样的 γ 在不同产热分布、不同电芯位置下对应的最热电芯温度可以差很多 ⇒ 判不了「够不够好」;只交第 ② 份,你知道热点在哪,却不知道该改哪一段(集管?支路?某个局部损失?);只交第 ③ 份,回答不了「均匀性是不是用泵功买来的」。⇒ 验收判据必须落在电芯温差与最热点上,⛔ 不落在 γ 上——这是全课第二根钉子,第 4、5、6、8 讲各会敲一次。

图1 本课三份交付物与各自回答的问题。三列关系图。左列三个方框是本课的三份交付物:①各支路流量分布柱状图加流量不均匀度 γ 与 γ_R;②最热电芯位置加全包电芯温差 ΔT_cell;③系统压降 Δp_sys 与水力功率 P_hyd 的拆账。中列写每一份回答的问题:分得匀不匀/热点在哪、差多少/匀是不是用泵功买来的。右列写只交这一份时下游做不了的事,逐条一句:只交①,γ 是过程指标,判不了够不够好;只交②,知道热点在哪却不知道该改哪一段;只交③,回答不了均匀性的代价从哪来。①指向②的箭头上标「两级翻译:Q_gen 除以 ṁ 乘 c_p 的项,与 Q_gen 乘 R_tot 的项」,②指向③的箭头上标「同一比较条件」。右下角一个禁止框写:三份缺一不可,不得只报其中一份作为验收。图上不出现任何数值。本图为交付物之间的关系示意,不代表任何平台的实际交付物清单,不得据图读取任何数值或指标限值。
图1 三份交付物与各自回答的问题。从图上要读出三个判断:①交付物是三份不是一份——左列三个框各配了一个不同的问题;②①→②那支箭头上写着「两级翻译」,说明流量分布本身不是验收依据,中间还隔着两级口径不同的换算(第 5 讲);③右列每一行都写明了「缺了它,下游具体做不了哪一步」。三列的对应关系与「缺了做不了什么」为本课归纳。本图为交付物之间的关系示意,⛔ 不得据图读取任何数值或指标限值。

工程量级。本课的限值类数字一律是平台相关量、⛔ 一个都不给,但每一个都带着去处:γ 的限值无公开统一标准,多由主机厂/电池厂的设计验证计划内控规定 ⇒ 不写这句,读者会去找一个并不存在的标准;本课给的替代物是从温差目标反解 γ 的链条(第 5 讲)。电芯温差目标与最热点温升上限 ⇒ D1-02 电芯温度不均与温差控制目标 与本项目的工况定义;电芯产热值 ⇒ D1-01 动力电池产热机理与产热速率建模 / D1-02 电芯温度不均与温差控制目标;系统压降与泵功预算 ⇒ J1-06 电池热管理一维建模与快充温升仿真 / A2-01 管网水力学:沿程/局部压降与流量分配。⇒ 而三份交付物的结构、以及它们之间的判据关系与平台无关,可以直接照用。

今天就能用的一份产物规格(把上面这段变成可核对的三行,交付前逐行对):

  • 分配柱状图:横轴 i=1..N 且 i=1 标「近端」,纵轴归一化(或写明量纲由本项目给);必须同时标出 N,以及这次比较是在同总流量/同泵工作点/同泵功中的哪一个条件下做的(为什么必须声明,见第 3 讲);γ 与 γ_R 两个数并列、各带名字,⛔ 不写一句笼统的「不均匀度是多少」。
  • 温度云图:标出最热电芯的编号与它所属的支路,并写明报的温差是「全包电芯最高−最低」这一层,⛔ 不与「支路间平均温差」「单颗电芯内部温差」混报——这三个数在同一个包上可以差好几倍,而它们在中文里都叫「温差」。
  • 压降-泵功对比表:γ(或 γ_d)、最热电芯温差、P_hyd 三量并列,并注明比较条件;只报相对变化,⛔ 不把某一次的绝对幅度写成经验值。

易错点。

  • 把本课当成某个求解器的操作手册去读。具体面板与选项归求解器文档,本课交付的是判据——判据换一个求解器仍然成立,面板不成立。
  • 以为「算出流量分布」就完成了。没有被翻译成温差与最热点的分配结果,不能作为验收依据。

1.2 支配方程只有一条:每一支路都在它自己的那个压差下取流量

是什么。对每一条支路都有

Δp_b,i = p_div(x_i) − p_com(x_i),i = 1..N,且 Σ ṁ_i = ṁ_total

其中 p_div(x) 是进(分流)集管在位置 x 的静压(Pa),p_com(x) 是出(合流)集管在同一位置的静压(Pa),Δp_b,i 是第 i 条支路两端的压差(Pa)。每一支路自身的 Δp_b,i(ṁ_i) 是它的阻力特性 ⇒ 这是 N 个非线性方程加一个总流量约束,解出来的就是分配。

为什么把它放在最前面。为了立刻取消一个很自然的错误直觉——「各支路做得一模一样,所以流量一样」。各支路的阻力特性可以完全相同,而驱动它们的压差并不相同,因为 p_div(x) 与 p_com(x) 都沿 x 变。⇒ 分配不均的根在集管,⛔ 不在支路。

⚠ 这条平衡还常被写成「进集管压降 − 出集管压降 = 支路压降」的形式。字面读下去像是两条集管的总压降相减,于是读者会以为各支路的 Δp_b 是同一个数——那正是本讲第一步要取消的东西。它实际说的是同一个轴向位置上两条集管的静压之差:p_div(x_i) − p_com(x_i)。⇒ 判断一份报告有没有踩这个坑,只要问一句:你的 Δp_b 是一个数,还是 N 个数

图2 整包并联流道的压力链拓扑。平面拓扑示意:左端进液法兰接进(分流)集管,集管沿 x 方向向右,其上等间距挂 N 条支路(画出 5 条并用省略号表示其余),各支路下端汇入出(合流)集管,右端接出液法兰。图下方画 x 轴,x 等于 0 与 i 等于 1 标在进液口一侧并写「近端」,i 等于 N 的一侧写「远端」。两条集管上分别标 p_div(x) 与 p_com(x)。第 i 支两端用一对箭头标出 Δp_b,i 等于 p_div(x_i) 减 p_com(x_i),五对箭头长度各不相同;每支中间标 ṁ_i。右侧竖排两个公式框:Δp_b,i 等于 p_div(x_i) 减 p_com(x_i)(i 从 1 到 N),以及 ṁ_i 求和等于 ṁ_total。框下一行禁止语:Δp_b 不是全局常数,把它当常数等于假设集管无损失。图上不出现任何尺寸或数值。本图为拓扑示意,箭头长度不代表任何真实压差比例,不得据图量取尺寸或压差。
图2 并联流道的压力链与支配方程。从图上要读出三个判断:①五对箭头长度不同且各带 Δp_b,i 标签 ⇒ 各支路的阻力特性可以完全一样,而驱动它们的压差并不一样;②约束只有一条 Σṁ_i=ṁ_total(写在公式框里)⇒ 分配就是「N 条阻力曲线各自在自己的压差下取值」;③x=0、i=1、「近端」三个标签画在同一侧——全课的方向约定就钉在这张图上。本图为非比例拓扑示意,⛔ 不得据图量取任何尺寸或压差。

工程量级。各支路 Δp_b,i 之间的差异,往往只占系统压降 Δp_sys 的百分之几到十几;而这点压差差被支路的阻力律放大或压缩:层流段 Δp ∝ ṁ¹ ⇒ 压差差 10% 就是流量差 10%;湍流段 Δp ∝ ṁ^1.8 ⇒ 压差差 10% 只换来 1.1^(1/1.8) − 1 ≈ 5% 的流量差。(★ 这条量级关系是本课归纳的;指数 1.8 是 1.75~1.8 区间内的取值——光滑管湍流的经典关联式给 f_D ≈ 0.316·Re^(−0.25),即 Δp ∝ ṁ^1.75,来源等级为行业手册/教科书、本课未逐页复核 ⇒ 全课统一取 1.8 做量级演算,⛔ 它不是任何流道的实测指数、⛔ 不作设计输入。)⇒ 湍流支路天然比层流支路好调匀,同一份不均的压差,在层流支路上会原封不动地变成流量不均。这也预告了第 8 讲:低温高黏会把支路从过渡带推进层流带,阻力律的指数从约 1.8 变成 1,均流方案是按原来那种形状调的

从进液法兰到出液法兰,压力依次经过哪几段——影响并联流量分配的因素全集。这个全集我是沿几何路径数出来的:从进液法兰出发,压力依次经过分流集管、支路、合流集管。⚠ 但不在这条几何线索上、却同属这个集合的还有两类:全局耦合类(泵工作点、热-流反馈——它们不在管路的某一段上,而是把整条链接成了环)与单相稳态 CFD 算不出、但台架上真实存在的类(气泡)。只沿几何线索数,这两类会成建制地不出现,而清单看起来是完整的。

  • 进(分流)集管侧:①沿程摩擦 ②支路取流造成的动量回升(静压回升)③分流三通的局部损失 ④集管截面沿程变化(渐缩/渐扩)⑤进口射流冲击与集管弯头
  • 支路侧:⑥入口收缩与出口扩张的局部损失 ⑦支路沿程摩擦(含入口段修正)⑧支路内转向(蛇形流道的弯头)⑨节流孔/调阻件 ⑩加工偏差、钎焊余料与堵塞造成的支路间阻力散布
  • 出(合流)集管侧:⑪合流三通的局部损失 ⑫速度沿程递增造成的动量损失 ⑬沿程摩擦 ⑭出口截面变化
  • 全局耦合:⑮泵曲线与系统曲线的交点(⛔ 总流量不是给定量)⑯温度→黏度→阻力的反馈 ⑰高差与重力项
  • 台架有而本课单相稳态 CFD 算不出:⑱气泡滞留与排气不净

⇒ 射程:①~⑭ 在本课全解(第 4 讲把 ①~⑨、⑪~⑭ 拆成压降账;⚠ ⑩ 不进那张账——它不是一条支路上某一项的大小,而是同一项在支路之间的差异,要用最坏组合或统计抽样另做一套,落在第 8 讲);⑮⑯⑰ 本课只给判据与接口(⑮ 去向 J1-06 电池热管理一维建模与快充温升仿真A2-01 管网水力学:沿程/局部压降与流量分配,⑯ 在第 8 讲判方向);⑱ 只给去向 L2-03 冷却系统台架与流量分配试验

易错点。

  • 把 Δp_b 当成一个全局常数。那等于假设集管无损失——而集管有没有损失,正是这门课要算的东西。
  • 忘了总流量约束本身也不是给定的。ṁ_total 由泵曲线与系统曲线的交点决定,均流改造会让这个交点移动(第 3 讲)。

今天就能用的一步:拿到任何一版分配结果,先把 Δp_b,i 与 ṁ_i 两列并排列出来,再按这条支路的阻力律指数把压差差折成流量差(层流 1:1、湍流约 2:1)。两列对不上,说明后处理取的不是同一位置的静压——先查量具,⛔ 别先怀疑模型。

1.3 进(分流)集管的静压会沿程上升:两项反向博弈,谁主导由一个比值定

是什么。进集管沿 x 不断向支路分流,集管流速 u_man(x)(m/s,u_man,0 指进集管入口处的值)沿程递减。一维动量平衡里因此出现两个方向相反的项:沿程摩擦让静压降,速度递减带来的动量回升让静压升。⇒ p_div(x) 到底是升还是降,取决于谁大。

为什么它能写成一个与流量无关的比值。两项都正比于 ½ρu_man² ⇒ 它们的比值不含流量的一阶依赖。在「等间距等流量取流、u_man 沿 x 线性递减」这一教学假设下逐项积分:摩擦项带一个 ∫₀¹(1−ξ)²dξ = 1/3,得 (f_D·L_man/(3·D_h,man))·½ρu_man,0²;取流回升项得 K_r·½ρu_man,0²。其中 L_man 是集管从第 1 支到第 N 支的有效长度(m),D_h,man 是集管的水力直径(m),K_r 是集管动压回升系数——支路取流处静压回升占动压变化的比例(无量纲,★ 本课记法)。两者相除:

回升项 ÷ 摩擦项 ≈ 3·K_r·D_h,man/(f_D·L_man)

⇒ 系数 3 是这条积分给出的代数结果,⛔ 不是经验系数、⛔ 不是修正因子。⇒ 集管越细长(f_D·L_man/D_h,man 大)摩擦越占优,越短粗回升越占优。

把它算成一个手算题(⚠ 以下两组全部是教学假设值,⛔ 不对应任何平台、⛔ 不得取用)。两组共同取 K_r = 0.4、集管 f_D = 0.03(⚠ 这是一个便于手算的整值,落在光滑管湍流 f_D 的量级内,⛔ 不是由任何关联式在某个 Re 下复算出来的值、⛔ 不得当作集管摩擦因子的取用值):

  • 短粗集管:D_h,man = 20 mm、L_man = 0.5 m ⇒ 3×0.4×0.020/(0.03×0.5) = 1.6 ⇒ 取流回升占优 ⇒ 远端抢流倾向
  • 细长集管:D_h,man = 10 mm、L_man = 1.0 m ⇒ 3×0.4×0.010/(0.03×1.0) = 0.4 ⇒ 沿程摩擦占优 ⇒ 近端抢流倾向

⇒ 两组的抢流端相反,而两组可以是同一种布置形式。这就是 1.4 那条反钉子的可手算形式。

图3 进(分流)集管静压沿程的两项博弈。单幅曲线图。横轴 X 等于 x 除以 L_man,范围 0 到 1,X 等于 0 在进液口一侧。纵轴为归一化静压,即 p_div 减 p_out 再除以 ½ρu_man,0 的平方,无量纲。三条实线对应比值 3·K_r·D_h,man 除以 f_D·L_man 分别等于 1.6、1.0、0.4,各带文字标签:回升主导时静压沿程上升;两项相当时近乎平坦;摩擦主导时静压沿程下降。其中 1.6 与 0.4 两条就是正文 1.3 的短粗集管与细长集管两组。另有两条细虚线画出其中一组的两项分量:摩擦项为负的 A 除以 3 再乘以 1 减去 1 减 X 的三次方,回升项为 K_r 乘以 1 减去 1 减 X 的平方,各带名称与符号,其中 A 每次出现都带着它的定义 A 等于 f_D 乘 L_man 除以 D_h,man,即集管的摩擦参数。图内信息框列出三组的 A 与 K_r 取值及其比值,并加粗写明这些是教学假设值、不对应任何平台。曲线由本课自设的一维集管教学模型现算,不代表任何真实集管,不得据图取压力值、不得反算流量不均匀度或任何支路流量。
图3 进集管静压的两项博弈。从图上要读出三个判断:①进集管的静压可以沿程上升,与「流速降下来了、压力当然也降」的直觉相反(上升那条曲线标着「回升主导」);②两条虚线说明升与降是两项相反的贡献,谁主导由 3·K_r·D_h,man/(f_D·L_man) 定,三个比值直接印在图上;③三条曲线用的是同一个 u_man,0(同一总流量)却给出三种走向 ⇒ 方向与总流量一阶无关。图上 A 与 K_r 为教学假设值,⛔ 不对应任何平台、⛔ 不得取用;曲线是教学模型的解,⛔ 不得据图取压力值、⛔ 不得反算 γ 或任何支路流量。

K_r 怎么拿到。本课不给 K_r 的数——它随支路开孔形式、开孔比、三通几何与取流比而变,属构型相关量,须由本项目标定;⛔ 禁止照抄别的项目或上一代的取值。替代物是两条标定路径:① 后处理反解——在全包算例里沿集管逐站取截面平均静压与动压,用相邻两站之间的静压回升量 ÷ 动压变化量即得该站的 K_r(本课给 K_r 的定义就是这个比值);② 局部算例——单独建一段含一到两个取流三通的集管算例,扫开孔比与取流比,得到 K_r 随工况的变化。⚠ 两条路径都必须写明取的是哪个截面、哪种平均方式(面积平均静压还是质量流量加权总压)——换一种平均方式就是另一个 K_r。

射程(写在这里,因为它决定了这个比值能不能进报告)。该比值只用于判「谁主导」,⛔ 不得反算 γ、⛔ 不得反算各支路流量。取流不等间距、或集管本身变截面时,上面那条积分要重新做,⛔ 不许原式照搬。

易错点。

  • 以为「集管里流速降下来了、压力当然也降」。那是把动量回升整项漏掉了——而它恰恰是远端抢流的物理来源。
  • 以为「把集管做粗就能扭转方向」。加粗让两项同时按 ½ρu_man² 缩小 ⇒ Δp_man 下降、幅度改善;但比值 3·K_r·D_h,man/(f_D·L_man) 里 D_h,man 就在分子上——摩擦项还要再按 1/D_h,man 多缩一档 ⇒ 比值只朝回升占优那一侧移动一点,⛔ 不会把抢流端换到另一头(除非这个包原本就落在临界附近)。⇒ 改的是幅度、⛔ 不是方向。第 6 讲 6.3 会再敲一次这条。

今天就能用的四步:量出你这个包的 D_h,man 与 L_man → 取集管 f_D 的量级 → 按上面两条路径之一标出 K_r → 代进 3·K_r·D_h,man/(f_D·L_man)。大于 1 ⇒ 回升占优、远端抢流倾向;小于 1 ⇒ 摩擦占优、近端抢流倾向。这一步只花几分钟,而它决定了 1.4 那一刀该往哪一端下。

1.4 出(合流)集管一定下降;U 型与 Z 型只决定两条曲线怎么叠

是什么。出(合流)集管里,流量沿流动方向不断汇入,速度递增 ⇒ 动量损失与沿程摩擦同号,静压单调下降,越靠近出口越低。⇒ 进集管两项相反(可能升也可能降),出集管两项同向(一定降)——两条曲线的性格不同,这才是 U 型与 Z 型差别的物理根源,而不是「Z 型天生不均」。

顺带一个可用的推论:出集管的压降通常比进集管更可预测(没有正负号之争)⇒ 拿一维估算去校核出集管那一段,比校核进集管更容易对得上。⚠ 但⛔ 不要把这句话读成「出集管可以用一维关联式顶掉」:出集管里全是合流三通,而三通与突扩正是局部损失系数最不可靠的地方 ⇒ 第 2 讲那一刀是两条集管都留给三维、把一维用在规整的支路流道上(2.2;三维/一维的界面画在三通处,见 2.7 验收第 2 项)。⇒ 这里说的「更可预测」只影响你对结果的信心,⛔ 不改变那一刀切在哪儿。

两种布置怎么叠。

  • Z 型(异侧进出):两条集管的流动方向相同 ⇒ p_div 与 p_com 的沿程变化在做差时相加,Δp_b(x) 沿 x 单调放大;
  • U 型(同侧进出):两条集管流向相反 ⇒ 两者的沿程变化在做差时部分抵消

⇒ 这解释了为什么在动压回升主导时 Z 型远端可用压差大、流量偏多,而在沿程摩擦主导时近端流量偏多;也解释了为什么 U 型在很多常规构型上看起来更均——它靠的是抵消,而抵消的程度依赖两条曲线的相对幅度 ⇒ 抵消可以不完全,也可以过头。⛔ 本课不给「U 型比 Z 型均多少」的任何数——那完全取决于 K_r 与 f_D·L_man/D_h,man。

★ 这一节要当场算破的那句话:「Z 型远端流量多、U 型近端流量多——记住这条口诀就能判方向。」

为什么它必然被读反(诱因就在这条结论平常的写法里,⛔ 不是本课立个靶子来打)。这条结论在本课程体系里的完整措辞,是把主导项布置形式用括号绑在一起说的:「集管动压回升主导时(Z 型异侧进出更典型)远端可用压差大、流量偏多;沿程摩擦主导时(U 型同侧进出更典型)近端流量偏多」。读者记住的是括号那两个字母(Z/U),忘掉的是括号那个前提(哪一项在主导)。⚠ 更糟的是:这条口诀在常见的短粗集管构型上碰巧是对的(1.3 的第一组:比值 1.6,回升占优,Z 型确实远端抢流)⇒ 它不会立刻出错,只在集管细长或取流回升偏弱的构型上翻过来(1.3 的第二组:比值 0.4)。一条经常对的口诀,比一条经常错的口诀危险得多——它会先赢得信任,再在某一个平台上安静地反过来。

读者会做错的那个具体动作。拿到一个 Z 型集管的包,不画 p_div(x) 与 p_com(x) 两条沿程曲线,直接判「远端抢流」,于是在远端支路加节流孔(或把集管远端做渐缩)。若这个包实际是沿程摩擦主导、近端抢流,这一刀正好加在流量本来就少的那一端

代价(结构性,⛔ 不给数)。被抬阻的那一支流量进一步下降,其余支路按并联关系分走这份流量 ⇒ γ 不降反升;同时系统总压降上升、泵功上升 ⇒ 三个指标同时朝坏的方向走。而最贵的不是这一步,是下一步:「我是按 Z 型的规律加的节流孔」这句话会让人判定方向对、只是量不够,于是继续收紧孔径——朝错误的方向迭代。⚠ 这类错在 CFD 上不会以不收敛的形式暴露:残差照掉、云图照出、总压降还是个看起来合理的数。

正确的判据(★ 本课归纳,⛔ 不是任何标准或手册既有的条款)。

  1. 判方向只看一件事——把 p_div(x) 与 p_com(x) 两条沿程静压曲线画出来,看它们的差 Δp_b(x) 沿 x 是增还是减。布置形式只决定这两条曲线怎么叠,⛔ 不决定叠出来的方向。
  2. 主导项之比可以先估:按 1.3 的 3·K_r·D_h,man/(f_D·L_man) 算一个数,判出谁主导——它与总流量一阶无关(f_D 随 Re 只弱变),所以在总流量还没定死的阶段就能用。⛔ 该比值只用于判「谁主导」,⛔ 不得反算 γ、⛔ 不得反算各支路流量。
  3. 不许拿别的项目的结论直接套:同样是 Z 型,集管长径比与取流回升强弱不同,方向可以相反。
  4. 节流孔一律加在抢流的那一端(由上面第 1 条判出来的那一端),⛔ 不是「远端」也⛔ 不是「近端」——方位是结果,⛔ 不是判据
图4 布置形式与主导项的四种组合。二乘二四联图。行为布置形式:上行 Z 型异侧进出,下行 U 型同侧进出。列为集管参数:左列回升主导,取正文 1.3 的短粗集管组(比值 1.6);右列摩擦主导,取细长集管组(比值 0.4)。每格纵轴为归一化支路压差,即 Δp_b(X) 除以 ½ρu_man,0 的平方;横轴 X 等于 x 除以 L_man。格内一条实线画 Δp_b(X),两端各打点并标「近端 i 等于 1」与「远端 i 等于 N」,Δp_b 最大的那一端另加箭头标「抢流端,节流孔加在这里」。下行左格另加一条虚线,标「合流侧动量系数与分流侧不等时」,其最大值落在另一端。四格标题各写该格的参数组与一句结论,四格共用同一纵轴范围。图底一行写:布置形式只决定两条集管曲线怎么叠,不决定叠出来的方向。四组参数为教学假设值、不对应任何平台,不得据图取值,不得把任何分界当判定阈值。
图4 同一种布置形式、两种主导项,抢流端相反。从图上要读出三个判断:①同为 Z 型,左右两格的「抢流端」箭头指向不同,右格的曲线还在中段最低 ⇒ 「Z 型远端多」只是回升主导那一格的形态,「沿 x 单调放大」在左上格成立,右上格给出的是它的射程;②下行的虚线在同一格里把最大值挪到了另一端 ⇒ U 型的方向也不是布置形式给的;③四格的判据是同一条——看 Δp_b(X) 沿 X 是增还是减。四组参数为教学假设值、⛔ 不对应任何平台,⛔ 不得据图取值、⛔ 不得把任何分界当判定阈值。

易错点。

  • 把出集管当成进集管的镜像来处理。两者的动量项符号不同,它不是镜像。
  • 出口边界紧贴最后一支的合流处 ⇒ 合流损失被出口边界污染,而最末几支恰恰是分配判断最关心的(第 3 讲给可判真假的自检法)。

1.5 不均的量级由阻力比 Λ 定;量级走到极端,就是两种短路形态

是什么。定义 Λ ≡ Δp_man / Δp_b(无量纲),其中 Δp_man = Δp_b(x) 沿 x 的极差(Pa,即集管侧压力沿程变化的幅度),分母是单支路压降。Λ→0 时各支路两端压差趋于相同 ⇒ 分配趋于均匀;Λ 越大越不均。★ 这个判据量是本课归纳的,⛔ 不是任何标准或手册既有的定义。

为什么值得给它一个名字。它把「集管形式、集管尺寸、支路阻力、支路数、节流孔」这些看似互不相干的变量收进同一个判据里 ⇒ 从此每一个均流方案都能被问同一句话:你压低的是 Λ 的分子哪一项,还是抬高了分母?(第 6 讲会照这句话把所有手段排一遍:变截面集管与加大集管截面压的是分子,节流孔与改流道截面抬的是分母,分区供液从两条路径一起压分子(⛔ 分母不动)。)

射程。Λ 与 γ 的关系是单调的定性关系,⛔ 本课不给 Λ→γ 的定量换算式——那需要具体的支路阻力律与集管几何,且在层流段与湍流段形状不同;Λ 的量级本身是平台相关量、⛔ 不给数(替代物就是 1.3 的比值加这里的定义式,两者都能在你自己的算例上算出来)。

两种极端形态。Λ 大到一定程度,分配会落进两种极端,而它们的补救方向相反

  • 近端短路:靠进液口的几支吃掉大部分流量,远端支路接近滞流;
  • 远端短路:远端几支抢流,近端偏少。

两者在分配柱状图上一个左高右低、一个左低右高;而更可靠的识别特征是集管压力沿程曲线的斜率符号——柱状图只给结果,曲线给的是成因

⚠ 极端形态下,滞流那一支的电芯温升由 Q_gen/(ṁ·c_p) 项主导且 ṁ 很小 ⇒ 温升会非线性放大:φ(冷却液温升在电芯相对进液温升中的占比,口径与推导在第 5 讲)在滞流支路上 →1,放大倍数逼近 2、也就是逼近线性反比 ⇒ 它是全课温升被放大得最厉害的那一支。⚠ 但⛔ 别把「放大得最厉害」读成「1+φ 在这里不好使」——恰恰相反:滞流支路 Re 已经很低、入口段很短,1+φ 的前提在那里反而更接近成立(第 5 讲 5.3);1+φ 真正失效的是另外两条边界(入口段主导、支路落在过渡/湍流带),两条都让实际放大倍数略大于 1+φ,同在 5.3。⛔ 具体数值是平台相关量、不给。

图5 近端短路与远端短路两种极端形态。两列对照。左列为近端短路,右列为远端短路。每列上行是 Δp_b(X) 沿程曲线,带斜率方向标注与「抢流端」箭头;下行是同一算例解出的各支路流量柱状图,支路数 N 等于 10,横轴 i 从 1 到 10 且 i 等于 1 标「近端」,纵轴为 ṁ_i 除以支路平均流量。最大与最小两支柱顶标出数值,柱状图右上角写出该例的 γ 与 γ_R 两个数,各带名字并写出 N。两列下方各一行「补救方向」,箭头指向抢流端并写「节流孔加在这一端」,两支箭头方向相反。图底一行:只看柱状图判不出成因,柱状图给结果,曲线给成因。本图为教学算例的示意,两组分配由教学模型现算、参数为教学假设值,不对应任何平台,不得据图取用图上任何 γ、γ_R 或流量比。
图5 两种极端形态,以及为什么补救方向相反。从图上要读出三个判断:①两列的柱状图左右镜像,而两支补救箭头方向相反 ⇒ 判错形态,这一刀就加在流量本来就少的那一端;②γ 与 γ_R 两个数并列写出且互不相等(各带名字与 N)⇒ 报不均匀度必须说清是哪一个;③曲线的斜率方向才是识别特征,柱状图只是它的结果。本图为教学算例的示意,两组分配由教学模型现算、参数为教学假设值、⛔ 不对应任何平台,⛔ 不得据图取用任何 γ、γ_R 或流量比。

易错点。

  • 把 Λ 当成能直接算出 γ 的公式去用。它只判量级与方向。
  • 只盯分母(加节流孔)而忘了分子。抬支路阻力会同时抬高 Δp_sys 与泵功,而改集管往往更便宜(第 6 讲的权衡)。
  • 只看柱状图就动手,把两种形态的补救方向搞反。
  • 看到「近端电芯更热」就直接断定该在近端做文章。现象与成因之间还隔着一次判据:近端更热确实对应近端流量少,但流量少是因为远端在抢 ⇒ 该动的是远端。本课那个 12 片冷板 Z 型案例的反直觉之处正在这里。

今天就能用的一张识别表(三行,从现象走到动作):柱状图形态(左高右低/左低右高)→ Δp_b(x) 曲线的斜率符号(增/减)→ 抢流端(远端/近端)→ 节流孔与渐缩加在抢流端。⚠ 前两行必须填,只填第一行就动手,正是上面那条反钉子的入口。

1.6 支路数 N 越多越难调匀——增加 N 把 Λ 的两头同时往坏里推

是什么。设集管截面与总流量不变,把并联支路数从 N 增到 2N(节距不变 ⇒ L_man 翻倍):

  • 分子侧:集管更长 ⇒ 沿程摩擦项 ∝ f_D·L_man/D_h,man ×2;而取流回升项与 L_man 无关 ⇒ ×1 ⇒ Δp_man 增大;
  • 分母侧:每支路分到的流量减半 ⇒ 单支路压降 Δp_b 下降:层流支路 0.5^1 = 0.5 倍,湍流支路 0.5^1.8 ≈ 0.29 倍。

⇒ Λ = Δp_man/Δp_b 两头都往上走:分子按摩擦项主导取 ×2,粗算得 2/0.5 = 约 4 倍(层流支路)与 2/0.29 ≈ 约 7 倍(湍流支路)。⚠ 若这个包是回升项占优,分子的放大会小于 2(回升项那一段不随 L_man 变长)⇒ 结论的方向不变、幅度变小——⇒ 支路数这条影响主要走的是摩擦项那条路。

图6 支路数加倍如何把阻力比 Λ 的两头同时推高。左右两栏拆账条图。左栏是 Λ 的分子,即集管侧压力沿程变化幅度 Δp_man;右栏是分母,即单支路压降 Δp_b。每栏两根条,分别对应 N 与 2N,前提是节距不变因而集管有效长度翻倍、集管截面与总流量不变。左栏两根条按两项分色堆叠,分别是沿程摩擦项与动量回升项,并各标名称与倍数:摩擦项乘 2、回升项乘 1。右栏两根条标乘 0.5(层流支路,压降正比于流量的一次方)与乘 0.29(湍流支路,压降正比于流量的 1.8 次方),两种阻力律各画一根。两栏中间画 Λ 等于分子除以分母的箭头,下方写出结果:Λ 放大约 4 倍(层流支路)与约 7 倍(湍流支路)。图底一行写清三条前提:集管截面不变、总流量不变、节距不变。倍数由图上写明的三条前提直接推出,是教学量级关系,不代表任何平台的实际 Λ,不得据图读取 Λ、Δp_man 或 Δp_b 的绝对值。
图6 支路数加倍,Λ 的两头同时朝坏的方向走。从图上要读出三个判断:①Λ 的两头同时变坏(左栏条变长、右栏条变短);②分子里的回升项那一段没有变长 ⇒ 支路数的影响主要走摩擦项这条路,⇒ 想抵消它就要动集管长度,⛔ 动截面不够;③层流与湍流两个结果并列(约 4 与约 7)⇒ 支路阻力律的指数会改变这件事的量级。倍数由图上写明的三条前提直接推出,是教学量级关系,⛔ 不代表任何平台的实际 Λ,⛔ 不得据图读取 Λ、Δp_man 或 Δp_b 的绝对值。

为什么这条值得单独讲。它给出了分区供液为什么有效的物理理由:把 N 支拆成 n_z 组各自带短集管的小并联,每组的 L_man 与组内支路数都变小 ⇒ Λ 的分子从两条路径一起降。⚠ 但⛔ 别把它读成上面那笔账的反向:分母不动——总流量与支路总数都没变,每一支仍然取到 ṁ_total/N,同一条流道、同一个流量 ⇒ Δp_b 不变(第 6 讲 6.5 展开,那里还要处理它带来的三项代价)。

工程量级。整包液冷的并联冷板数常在十的量级(典型)。⚠ 本课程体系案例里的 12 片与动手做的 8~12 支路场景设定值,⛔ 不是推荐值、⛔ 不是典型带的端点——本课不由它们推出任何结论,换成别的支路数结论照样成立。⛔ 具体 N 是构型相关量(由模组布置与冷板分片决定),本课只给它对 Λ 的作用方向。⚠ 另外:报 γ 必须同时报 N——不同 N 之间的 γ 可比性有限。

易错点。

  • 以为「支路多了平均下来更均匀」。统计直觉在这里是反的:这不是抽样问题,是压差沿程累积问题。
  • 拆分区时只拆支路不拆集管(还是一根长集管上挂几组)⇒ L_man 没变,Λ 的分子没降。

今天就能用:上面那笔账的价值不在「4 倍」「7 倍」这两个数,在算法本身——换成你自己的支路数比、你自己的节距变化、你自己那条支路阻力律的指数,代进「分子倍数 ÷ 分母倍数」,答案自己就出来了。⚠ 三条前提(集管截面不变、总流量不变、节距不变)任何一条不成立,这笔账要重列。

1.7 什么时候可以不做 CFD:先估 Λ,再定建模层级

是什么。若支路阻力远大于集管阻力(Λ ≪ 1),各支路两端压差近似相同 ⇒ 分配近似只由各支路自身阻力决定,用一维并联网络(支路阻力律 + 集管沿程/局部损失关联式)就能给出足够准的分配,⛔ 不必上三维。

为什么这条判据重要。它把 CFD 从「必做的仪式」变成有条件的选择,也把网格预算从「越细越好」变成「花在能改变结论的地方」(第 2 讲)。⇒ 先估 Λ,再决定建模层级——这正是第 2 讲那张建模层级表(图7)下面那条决策链的第一步。

判据结构(⛔ 本课不给 Λ 的判别阈值)。「≪ 1」到底是多少,取决于你要分辨多大的流量差 ⇒ 本课给的是把两个量放在同一把尺上比的做法:

  1. 写下你要分辨的支路间流量差(例如:小于百分之几的差异对本项目的温差目标没有意义);
  2. 用 1.2 的折算把它换成允许的支路间压差差:层流段压差差多少就是流量差多少,湍流段压差差 10% 只换来约 5% 的流量差 ⇒ 湍流支路允许的压差差更大;
  3. 把「Λ 引起的支路间压差差」与上一步得到的允许值比——前者远小于后者 ⇒ 一维并联网络够用;
  4. 否则按「我要回答的问题需要解出哪一层」去选建模层级(第 2 讲的七种层级与它们各自抹掉了什么)。

这条判据的失效边界(必须一起写出来,否则它会被当成万能)。集管有弯头、突扩、非对称进口,或支路阻力律跨越层流与过渡带时,一维关联式最不可靠——⇒ 那正是必须上 3D 的场合,而不是「Λ 小所以可以省」。

易错点。

  • 把这条读成「Λ 小就不用算了」。不均可以小到不影响分配,却仍然影响最热点——产热在支路间的不均会与流量不均叠加或抵消,最热电芯未必落在流量最少那一支上(第 5 讲)。⇒ Λ 小免掉的是分配这一问,⛔ 不是温度那一问。
  • 忽略失效边界,把一维网络用到集管带弯头、突扩或非对称进口的构型上——那几处恰恰是局部损失系数最不可靠的地方。

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

会员专属

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

查看会员方案