J1-03 冷却液回路水力与热仿真建模
课程代码 J1-03 · 板块 J 仿真分析 CAE / 一维系统仿真 时长 约 3.5 小时(6 讲) 前置 A2-01 管网水力学;A2-02 泵曲线匹配;J1-01 1D 建模范式;A2-05 空气侧阻力链与进风量 本讲义定位 讲师授课蓝本 / 学员自学讲义,是 J1-03 大纲的完整展开版
引言:模型不会报错,它会告诉你「这个现象不存在」
冷却液回路是整车热管理里唯一一张把几乎所有热源和几乎所有散热面连在一起的网:电池、电驱、电控挂在水侧,低温散热器(LTR)通向空气,Chiller 通向制冷剂侧,中间靠一台或几台电子水泵、一段段管路和一组多通阀,串并联成若干个可以互相切换的模式。给这张网做一维仿真,表面上是一件相当程式化的活——把泵曲线放进去、把每段管路的阻力算出来、把换热器降阶成一个 UA、连线、按下运行;工具会收敛,会吐给你一张流量分配表和一组节点温度,曲线通常还挺好看。
难就难在这里:它几乎总会给你一个结果,而那个结果对不对,求解器不负责。本课要建立的,正是替求解器负这个责的三种判断力——这张网的方程是怎么闭合的(未知量、方程数、参考点,缺一样就有一类结论永远不会出现);水力与热是怎么互相咬住的(哪些反馈可以断、断了以后偏差朝哪个方向);以及一份 1D 结论的射程到哪里为止(它能回答哪一类问题、不能回答哪一类,必须写进交付说明而不是留在建模者脑子里)。
钉子①:第一位的缺陷来源不是参数取得准不准,而是自由度建没建
本课从头到尾要敲的第一根钉子是这一句:
冷却液回路 1D 模型的第一位缺陷来源不是「参数取得准不准」,而是「自由度建没建」——一个没有被建进模型的现象,模型不会报错,它会如实地报告「这个现象不存在」。
这句话之所以要当钉子,是因为它和大多数人对「仿真误差」的直觉正好错位。参数不准是一个精度问题:它会随着你标定、加密、收紧容差而改善,而且通常有一个可感知的量级。自由度缺失是一个结构问题:它不随任何参数变化。把网格加密十倍、把容差收紧三个数量级,那个现象照样是零;而求解器全程一个字都不会说,因为在它看来,模型里根本没有承载那个现象的状态量。把 8 片并联冷板挂在同两个节点上,不均匀度就恒为零;把回路建成刚性不可压,压力波的峰值就恒为零——不是「小」,是「零」。
由这根钉子直接推出本课的三条工作纪律。
推论 A(怎么操作):每讲一个建模对象,先答一句「不建它,模型会漏报哪一整类结论」。本课把这一问归纳成十条自由度的全集(水力、定压、充液含气、热容、输运、换热、空气侧、拓扑切换、物性、控制),第 1 讲立起来,后面每讲兑现一到两条,逐条给出「不建 ⇒ 漏报什么」。⚠ 这份全集是本课自己归纳的,不是任何一份标准或既有规范里现成的条目;它的价值不在于条数,而在于它是按「建模决策」而不是按「硬件盒子」切的——沿着「泵—管—换热器」这条最顺手的硬件线索去数,会整整齐齐地漏掉定压、含气、拓扑切换、控制这四整类,因为它们不对应任何一个硬件框。而恰恰是它们决定了气蚀判据能不能算、水锤能不能出现、切换瞬态能不能看见、以及你看到的结果属于哪一类。
推论 B(怎么读结果):模型跑通、收敛、曲线好看,都不构成「这个风险不存在」的证据——它只说明「你建进去的那几个自由度上没有问题」。同理,「和试验对上了」也不能反证模型结构完整:对上的只是你测了的那几个量(第 6 讲会把这一层拆开讲)。⇒ 因此每一份 1D 结论都必须附一句射程声明:本模型能回答哪一类问题、不能回答哪一类。这一条在第 6 讲会落成一份可以照着填的交付口径,⛔ 不是一句自觉。
推论 C(和前一门课划清界):J1-01 1D 热管理仿真导论(KULI/AMESim/GT/Flowmaster) 里也有一根关于自由度的钉子,但那是另一件事,本课不重讲、也⛔ 不得把两者写成同一条。J1-01 1D 热管理仿真导论(KULI/AMESim/GT/Flowmaster) 讲的是 1D 这个方法本身没有的自由度——横截面上永远只有一个值,你再加控制体也长不出第二个值来。本课讲的是你在建模时自己选择不建的那些自由度:定压点、气腔、金属热容、控制体链、离散阀位、空气侧闭环——它们只要你愿意建就有,代价是建模工作量和收敛难度,不是方法的边界。两者的处置方向完全相反:前者要换方法或改分工(交给三维或交给试验),后者要改模型(就在你这张网里补上)。判不清这一点的人,会把本来该自己补的东西推给 CFD,或者反过来指望在 1D 里看见它根本装不下的东西。
钉子②:水力与热是双向咬合的,「先流量后温度」是迭代顺序不是传递方向
第二根钉子把本课的两半内容连起来:
水力与热在冷却液回路里是双向咬合的,「先解流量、再解温度」是一个迭代顺序,⛔ 不是一次性的单向传递。
至少三条反馈把它闭成环:UA(ṁ) 让温度反过来依赖流量分配;μ(T) 让流量分配反过来依赖温度(低温段尤其剧烈,Re 会跨过层流/湍流分界,连阻力指数都跟着变);第三条更隐蔽——空气侧的风量与 air-on 也在这个环里,它们不是外部给定的边界,而是要和水侧联立解出来的(第 3 讲)。
⇒ 由此得到一条贯穿全课的判据:模型里每断开一条反馈(物性按定温取、UA 按额定流量取、air-on 按环温取),就产生一处方向可判的系统性偏差。请注意本课的立场:断开是允许的,工程上常常必须断——但偏差的方向、以及它成立的条件,必须写进交付说明,⛔ 不许用「影响不大」四个字带过。判不出方向的解耦,等于把一个未知量伪装成了零。
由钉子② 还要顺出本课在能量侧的统一口径,它是本课最容易被含糊过去的一处:稳态、热源功率给定、且未触发温度保护时,Σ排热 ≡ Σ产热——这是能量守恒的直接后果,没有讨论余地。⇒ 加流量、加 UA、改拓扑,改变的是达成这个平衡所需要的温度水平(各节点温度、回路温差),⛔ 不是稳态排热量本身。而「散热能力」是另一个量——它指的是给定温度上限下能排掉的最大功率,这个量确实随流量与 UA 增大。两句话都对,混起来就出错:「加流量提高散热能力」对,「加流量提高稳态散热量」错。第 1、2 讲会用同一组算例把这两句话分别兑现一次。
反钉子:等压降不是求解器附送的定律,是你要自己建进去的约束
本课最容易被读反的是这一条:
「并联支路两端压降相等」,是一条你必须自己建进模型的约束,且只在特定拓扑下严格成立——⛔ 不是一条把支路连上去就自动生效、可以直接拿来当结论用的定律。
它之所以几乎必然被读反,是因为「流量分配」「压降相等」「串联叠加」这几句话各自都对,合起来却给人一个印象:这是物理定律,接上线它就成立了。真相分两层,两层都要点破。
第一层:它是列方程时要加的约束。这张网的方程闭合靠三样东西凑齐——节点流量守恒、每条支路自己的 Δp-Q 特性、以及一个压力参考点。三样缺一样,方程组就不定;而 1D 工具在方程不定时通常不会拒绝求解,它会用别的东西把它补齐——最常见的是你自己在某条支路上放了一个「定流量边界」,或者工具的默认初值把解锁在了某一支上。⇒ 结果是模型照常收敛、照常给出一张漂亮的流量分配表,而那张表和水力无关。自查只有一条:数一数未知量和方程数,并逐条检查有没有哪条支路的流量是「我给的」而不是「解出来的」。
第二层:这条铁律的射程是「各支路真正共用同一对节点」。总管/集流管本身有沿程压降和动压回升时,各支路接入点的可用压差沿流向是变化的——靠近进液口的先天占优、末端天然吃亏。而在 1D 建模里,「共用一对节点」不是实物的性质,是你连线时做的选择:把电池包里 k 片并联冷板画成挂在同两个节点上的 k 条支路,求解器当然会还给你 k 个完全相等的流量。⇒ ★ 他以为模型在「验证」分配,实际上模型在「假设」分配。这就是读反的那一下。
读反之后会做错的那个具体动作:面对一个包内温差超标的问题,1D 模型给出 k 片冷板流量严格相等、出水温几乎相同 ⇒ 工程师据此判定「流量分配没问题」,把温差归因到电芯产热一致性或冷板流道结构,于是去改电芯分选、去改冷板内部流道。而真因(集流管造成的近端多、远端少)在这个模型里从头到尾不可能出现——那个自由度在连线的那一刻就被取消了。⇒ 改进方向整个指反,一轮样件白做。正解是:要在 1D 里看得见分配不均,必须把总管/集流管本身建成带节点的支路(让各支路接在不同节点上),或者干脆把这一层交给部件级三维计算去定(前指 J3-02 电池包液冷流道 CFD 与流量均匀性)。
⚠ 反向同样要防:⛔ 不得由此推出「等压降是错的」。总管压降相对支路压降可以忽略时,它是很好的近似;而且它永远是列方程的正确做法。错的只是把它当成「分配没问题」的证据。判据是:先算总管压降占支路压降的比例——比例小到什么程度才可忽略,取决于你要分辨多大的不均,这是本课要你自己算的一步,⛔ 本课不给一个通用门槛(它是平台相关量)。
全课共用一个教学算例,每个数你都能自己复算
为了让上面这些判断都能被验证而不是被相信,本课全程共用一组教学算例(称作 T2):一组冷却液物性、一套回路几何与容积、两条支路加一段串联段的阻力、一条泵曲线、LTR 两侧的 hA、两档空气侧流量与环温、两个热源功率、一个控制体数。⚠⚠ 这一组全部是教学假设值,不对应任何平台、任何车型、任何在售产品,⛔ 禁止照抄取用;它存在的唯一目的,是让本课出现的每一个数都能被你自己拿计算器复算出来、并且彼此严格自洽。真实工程取值必须由你所在平台的实测或仿真确定。
它顺带把钉子① 变成可量化的:同一个 T2 回路,漏建金属热容时热时间常数从 102.19 s 掉到 65.50 s;把 air-on 拍成环温时 LTR 进水温度低算 6.63 K(回路各节点同步低算同一个 6.63 K,其中回路最高点由 74.87 ℃ 低算成 68.24 ℃)。这两个偏差都不是「误差」——它们是两个自由度没建的必然后果,标定再多轮也标不回来。
本课的射程,以及不归本课管的那些事
★ 一门课把边界写清楚,比多讲两条更重要。本课只负责「把冷却液回路建成一个能解、能读、能交付的水力—热耦合一维模型」,下面这些相邻的事各有归属,本课只写去向、不重讲:
- 通用数值内核(微分代数方程与一致初值、代数环、容差取法、事件定位、守恒残差自检、刚性与积分器选型)归 J7-08 一维系统模型的数值内核与「跑飞」排查:DAE/代数环、步长与容差、一致初值与守恒自检;本课只讲它在水力域长什么样、从模型侧怎么规避。
- 相关性验证的流程、指标与阈值归 J7-01 仿真—试验相关性验证(Correlation);本课只给「按层标定、标定与验证不用同一组数据」这条纪律。
- 阀件与集成流道的阻力曲线怎么生产(多通阀、阀岛、集成模块壳体内流道的 ζ/Kv-Cv 曲线族)归 J3-05 阀件与集成流道内流 CFD:ζ/Kv-Cv 特性曲线族的生成与交付;本课只定接口口径——要什么形状的曲线族、按什么参数成族。
- 集流管尺度的流量均匀性归 J3-02 电池包液冷流道 CFD 与流量均匀性;本课只负责判断「在 1D 里要不要把总管建成带节点的支路」。
- 拓扑与模式矩阵的设计归 B3-01 冷却液回路拓扑:串/并联、多回路耦合与模式切换,多通阀的密封、内漏与阀位重叠几何归 G4-03 多通阀(四通/八通)结构与密封设计,集成模块的边界划分归 G8-01 集成热管理模块(TMM/阀岛/超级水壶)架构;本课只负责把它们翻译成算例,并校核每一个模式下网络还闭不闭合。
- 空气侧的阻力链、格栅有效通流面积、层叠顺序决策归 A2-05 空气侧压降、风量与风阻的工程处理/G7-03 前端进气量、格栅开度与散热能力匹配/G7-02 多换热器层叠布置与空气侧串扰,机舱三维反演归 J2-02 机舱热管理与前端模块 CFD;本课只负责气—液联立的算法与收敛判据。
- 制冷剂侧的两相离散与物性归 J1-02 制冷循环一维建模与制冷剂物性,Chiller 的选型与承压归 G3-01 Chiller(板式换热器)设计与选型,电芯级热模型归 J1-06 电池热管理一维建模与快充温升仿真;本课只定接口。
- 控制策略本体(泵转速、阀开度、风扇的标定与逻辑)归 K 板块的相应课程;本课只到「它是模型的一个自由度,闭不闭环决定了你的结论属于哪一类」,⛔ 不给策略、不给标定值。
- 往下的交接见 J1-04 整车热管理系统级联合仿真(整车联合仿真接口)。
关于标准,本课取一条克制的写法:涉及的三类文件——冷却液产品本身的国家标准、计算传热的验证与确认框架、以及冷却系统性能试验体系——本课一律只写去向、不写条款内容,因为本课并未逐条核对它们的原文;版次与年份以现行目录为准,引用前须自行核对,⛔ 本讲义不把任何一个年份写死当作依据。另有一条本课自己作出的射程判定(不是这些文件里的原话):验证与确认这一对概念里,「验证」那一半的做法归 J7-08 一维系统模型的数值内核与「跑飞」排查:DAE/代数环、步长与容差、一致初值与守恒自检、「确认」那一半归 J7-01 仿真—试验相关性验证(Correlation),本课只借用它「两者不可互相替代」这一条结论;冷却液的国家标准管的是冷却液这个产品,⛔ 管不到回路仿真的方法与允差;试验体系管的是台架怎么做,⛔ 不给建模判据。具体件号按项目确认,⛔ 本讲义不编。
六讲怎么走
- 第 1 讲 水力网络:把方程数清楚,把参考点定下来——支路特性从哪来(本课归纳的七类来源)、泵曲线与系统阻力曲线求交、并联分配与配平、定压点与含气这两个最容易漏的自由度、低温这个第二工况族,以及水力域专属的七条不收敛成因。
- 第 2 讲 换热部件的水侧降阶——ε-NTU 在回路里到底扮演什么角色、UA 是工况量不是部件属性、水侧加流量的收益上限、C_min 会翻边、Chiller 的特例与它的边界。
- 第 3 讲 空气侧不是给定边界——阻力链与风扇工作点解出风量、层叠串扰把上游的 air-on 传给下游、气—液闭环怎么迭代、以及收敛判据该取在哪两个量上。
- 第 4 讲 多回路拓扑与模式切换——串联与并联是两种不同的方程结构、阀在模型里的六种表达方式、拓扑重构与离散事件、止回阀这类由模型自己解出状态的元件、每个模式都要单独校核闭合性。
- 第 5 讲 时间尺度——回路热时间常数的分母到底该取什么(本课在这里对一个流传很广的写法作一次订正)、输运时滞是纯延迟而不是一阶惯性、管段必须离散成控制体链、要不要算水锤的判据、四个时间尺度放在同一根轴上看。
- 第 6 讲 标定与交接——按层标定(先水力后热,⛔ 不许用一个综合量同时标两层)、证据分三层共八类、「模型温度与实测对不上」的八类排查、两个错互相抵消的模型为什么最危险、以及一份 1D 交付里必须写清的射程声明。
课程末尾是一次双回路水力实操:建一个「电子水泵—散热器—电池冷板—电驱—三通阀」的双回路模型,解并联分配与稳态温度,切阀看重分布,再阶跃热源量出时间常数。四步分别对应本课的四条目标,每一步都有一个可以判真假的交付物——⚠ 其中量时间常数那一步,同时也是验证本课那条订正是否成立的证据。
第 1 讲 水力网络建模:流量是解出来的,不是给出来的
这一讲要建的是冷却液回路的水力层——泵、管段、阀件、并联支路,以及它们背后那张节点—支路网络。它看起来是全课最“基础”的一讲,实际上是本课钉子最集中的一讲:1D 冷却液回路模型的第一位缺陷来源不是「参数取得准不准」,而是「自由度建没建」——一个没有被建进模型的现象,模型不会报错,它会如实地报告「这个现象不存在」。水力层正是这句话最容易兑现的地方,因为水力网络的所有约束都靠你自己连线、自己列方程,而工具几乎从不拒绝求解一个不完整的方程组。
本讲的推进路线是:先把这句钉子说透并给出它的操作化清单(1.1),然后把网络方程的闭合条件摊到桌面上、让读者能自己数(1.2);接着交代支路特性从哪来(1.3),把工作点与并联分配真正解一遍(1.4);再用一组三方案的对照把“配平”这件事讲透(1.5~1.8)——这四节里藏着本课最容易被读反的一条;随后是泵的建法(1.9)、两条常被当成“部件”而其实是自由度的东西(1.10 定压点、1.11 含气),最后是低温这个第二工况族(1.12)与水力域专属的不收敛与错解(1.13)。
本讲全程使用一个教学算例,下文称它为本课教学算例。⚠ 它的每一个输入都是教学假设值,存在的唯一理由是让本课出现的每一个数都能被你在一张纸上复算出来;它不对应任何平台,⛔ 禁止照抄取用,真实取值须由本项目按 BOM、管路走向、部件数据表与实测/系统仿真确定。
1.1 没有被建进去的现象,不会以「报错」的形式暴露
是什么。1D 网络模型里,每一个可观测量都对应一个状态变量。一个物理现象要能在结果里出现,模型里必须先有承载它的状态量:气蚀要有气腔体积,压力波要有流体可压缩性,支路间的流量不均要有集流管上的中间节点,切换过渡态要有离散阀位。这些状态量不存在时,对应的现象不是“算得不准”,而是在这份模型里根本不可能出现。
为什么这不是「精度不够」而是结构性缺失。精度问题会随参数调整而改善:网格加密、容差收紧、物性取得更细,误差都会往下走。结构性缺失不随任何参数变化——把控制体加密十倍、把容差收紧三个数量级,那个现象照样恒为零,而求解器全程不报一个字。这就是两者最要命的区别:精度问题会用“还差多少”的形式告诉你它在,结构缺失只会用“没有这个问题”的形式让你放心。
操作化:每讲一个建模对象,先答一句「不建它,模型会漏报哪一整类结论」。把这一问在冷却液回路上问一遍,得到的是十条必须显式决定“建还是不建”的自由度——⚠ 这十条是本课按「先画全集、再逐项核射程」的做法归纳出来的,⛔ 它不是任何标准、手册或既有分类里现成的清单,引用时须注明出处是本课:
① 水力自由度:节点压力 + 支路流量。不建(给定流量边界)⇒ 抢流量与配平永远不出现,模型不能回答「改了会怎样」。 ② 定压自由度:全网一个压力参考点(膨胀壶/去气罐接点)。不建 ⇒ 不可压网络的压力场只确定到相差一个常数 ⇒ 一切依赖绝对压力的判据(气蚀、局部沸腾、压力盖开启、承压校核)全部失效,而模型照常给出流量与温度。 ③ 充液/含气自由度:液位、气腔体积、流体可压缩性。不建(刚性不可压)⇒ 压力波、水锤、气堵、吸空在模型里没有承载它们的状态量,永远不会出现。 ④ 热容自由度:冷却液 + 金属件(管路、壳体、芯体、冷板基板)分别建。不建金属 ⇒ 时间常数偏小、瞬态偏快。 ⑤ 输运自由度:管段离散成控制体链。管段当零容积 ⇒ 纯延迟消失,温度传播变成瞬时,控制回路的相位裕度算错。 ⑥ 换热自由度:UA/ε 随两侧流量与温度变。写成常数 ⇒ 低流量段的换热被系统性高估。 ⑦ 空气侧自由度:风量与进风温度由阻力链、风扇工作点与层叠串扰解出。拍成常数 ⇒ 结论偏乐观且方向一致。 ⑧ 拓扑切换自由度:阀位是离散状态、切换是事件、切换会改变网络的连接关系。写成一个连续可变阻力 ⇒ 过渡态(瞬时短路/断流)与切换冲击永远不出现。 ⑨ 物性自由度:ρ/c_p/μ/k 随温度(与乙二醇体积比)变,尤其 μ(T)。定温取值 ⇒ 低温工况下阻力与分配全错。 ⑩ 控制自由度:泵转速、阀开度、风扇由控制器按温度闭环给。开环拍常数 ⇒ 你看到的是「物理上能不能」,而不是「这台车实际会怎样」;两者的结论可以相反。⚠ 控制策略本体归 K 板块,本课只到「它是模型的一个自由度,闭不闭环决定了结论的性质」,⛔ 不给策略与标定值。
⛔ 这十条里有四条是靠“看硬件”永远想不到的。只沿「泵—管—换热器」这条最熟悉的硬件线索枚举,必漏 ②③⑧⑩ 四整类——它们不对应任何一个硬件框,全是建模决策。而 ② 决定气蚀判据能不能算、③ 决定水锤与气堵能不能出现、⑧ 决定切换瞬态、⑩ 决定你看到的是哪一类结果。⚠ 漏掉的方式不是“少了一条”,是那一整类根本不出现在选项里,所以清单看起来始终是完整的。
工程量级。本课后面几讲会用四处实例把这句钉子量化(数值均由本课教学算例的教学假设值算出,⛔ 不对应任何平台):把 8 片并联冷板挂在同两个节点上,分配不均匀度恒为 0;把回路建成刚性不可压,压力峰值的上界 1.07 MPa 在结果里恒为 0;漏建金属热容,回路热时间常数从 102.19 s 掉到 65.50 s;把进风温度拍成环温 40 ℃,算出来的回路温度整体低了 6.63 K。四处的共同点是:模型全都收敛,曲线全都好看。
与前置课的分界。J1-01 1D 热管理仿真导论(KULI/AMESim/GT/Flowmaster) 讲的是 1D 这个方法本身没有的自由度——横截面上永远只有一个值,加控制体也长不出来;本课讲的是 你在建模时自己选择不建的那些自由度,后者只要你愿意建就有,代价是建模工作量与收敛难度,而不是方法边界。⚠ 两者的处置完全相反:前者要换方法或换分工,后者要改模型。
易错点。
- 把「模型收敛且结果好看」当成「这个风险不存在」的证据。它只说明你建进去的那些自由度上没有问题。
- 用「和试验对上了」来反证模型结构完整。对上的只是你测了的那几个量——这正是第 6 讲按层标定要处理的事。
- 交付结论时不写射程声明。不写,读者就默认这份模型什么都能答。⇒ 每一份 1D 结论都必须附一句:本模型能回答哪一类问题、不能回答哪一类。
1.2 未知量与方程数必须能数出来:三样凑不齐,方程组不定
是什么。不可压回路的水力问题是一组代数方程(不含时间导数)。把回路抽象成节点—支路网络之后,未知量是 N_节点 个内部节点的压力 + M 条支路的流量,闭合它需要三样东西凑齐:
① 每个节点上的流量守恒 Σ Q = 0(N_节点 条); ② 每条支路自己的 Δp-Q 特性(M 条,它从哪来见 1.3); ③ 一个绝对压力参考点(1 个,消掉压力场的任意常数,见 1.10)。
为什么必须自己数。非因果的物理网络工具会替你自动装配方程组,你看不到方程本身。于是“少了一条约束”这件事不会以「你少写了一条方程」的形式报出来,它只表现为:收敛变慢、对初值敏感、或者换一个初值得到另一个解。⚠ 这三种表现都很容易被归因成“数值问题”,然后你去调容差——而调容差治不好结构性问题(1.13 会再讲一次)。
工程量级。本课教学算例是能说明问题的最小拓扑:泵与主管段(含低温散热器)为串联段,一个分流节点与一个汇合节点之间挂两条并联支路(电池冷板支路、电驱支路)。于是未知量 4 个(2 个节点压力 + 2 条支路流量),方程 4 条(2 条节点守恒 + 2 条支路特性),再加 1 个压力参考点。这个规模可以在一张纸上列完并手解——1.4 就把它解出来。
易错点。
- 在某条支路上加「定流量边界」,好让它拿到“该拿的”流量。那是把一个未知量当成了输入:方程数超了,而模型从此再也不能回答「改这条支路的阻力会怎样」。⚠ 这是本讲反钉子第一层最常见的具体形态。
- 把泵同时建成「定压升」和「定流量」。同样是过约束。
- 忘了定压点。全场压力浮动一个常数,流量与温度照样解得出来,而气蚀判据算不出来。
1.3 支路 Δp-Q 特性有七种来源,而手册只覆盖其中一类
是什么。1.2 的第 ② 样东西——每条支路的 Δp-Q 特性——是整张网络里唯一“要你去找数据”的部分。按「这条曲线从哪来」而不是按硬件分,它有七类来源(⚠ 这一分类同样是本课归纳的,⛔ 不是手册既有):
① 可解析式:直管段沿程,Darcy-Weisbach 加 f(Re, 相对粗糙度)(承 A2-01 管网水力学:沿程/局部压降与流量分配,本课⛔ 不重推)。⚠ 这里刻意写「相对粗糙度」四个字而不用 ε——本课的 ε 只表示换热器效能,⛔ 全课不用同一个字母表示粗糙度(见〈全课符号与口径表〉)。 ② 手册 K 值:标准弯头/三通/突扩突缩/常规阀。⚠ 三通的分流与合流 K 随支路流量比变化甚至变号,⛔ 不得用单一常数(A2-01 管网水力学:沿程/局部压降与流量分配 已写死)。 ③ 供应商随件交付的 Δp-Q 曲线:换热器芯体、冷板、成品泵、成品阀。用了它就不许再叠一遍它内部的局部损失。 ④ 部件级 CFD 生产的 ζ / Kv-Cv 曲线族:多通阀、阀岛、热管理集成模块壳体内的集成流道——手册一个字都没有。 ⑤ 台架实测反演:样件、非标件,以及要把内漏含进去的场合。 ⑥ 整机水力平衡反标:装车后测总流量与若干处压差,反推整段未知阻力(第 6 讲的水力层标定)。 ⑦ 主动元件的特性:泵曲线族(转速/PWM)、比例阀开度—Kv 曲线。⚠ 它们不是阻力,是把一条外来曲线塞进网络,因此有“曲线外区段”这个专属陷阱(1.13 第 ② 条)。
⛔ 只沿「手册查 K 值」这条线索枚举,必漏 ④⑤⑥——而现代集成模块(阀岛、热管理集成模块)的主要压降恰恰全落在 ④ 里。集成度越高,回路里“查得到”的那部分压降占比越小;而模型对这些件通常只能先塞一个猜的 K 值进去,然后一切结论都建在那个猜的数上。
本课的射程。①②③⑦ 在本课射程内并给用法;④ 本课只给接口口径,生产方法归 J3-05 阀件与集成流道内流 CFD:ζ/Kv-Cv 特性曲线族的生成与交付;⑤⑥ 只给“拿它标哪一层”,试验方法本身不在本课。
工程量级。⛔ 集成流道的 ζ 与 Kv-Cv 属平台相关量,本课不给数。本课只给接口口径:一维侧要的是按 (阀位 × 端口对 × 流量) 成族的 Δp-Q(或 ζ-Re)曲线,并且必须随曲线族交付五项——参考截面、Δp 类型(全压还是静压,口径承 A2-06 流体力学基础:连续性、伯努利与静压/动压/全压口径)、介质与比重折算、取压位置、阀位定义;缺任何一项,一维侧都没法消费这条曲线。
易错点。
- 拿一个「同口径球阀」的手册 K 值代多通阀——量级都不对。
- 拿到一条曲线就当它覆盖全阀位。曲线族的每一支对应一个阀位与一个端口对,缺哪一支哪个模式就没法算。
- 把含内漏与不含内漏的两份曲线混用。前者是“在计算域里留了微间隙”算出来的,后者不是,⛔ 两者不可互换。
1.4 工作点是解出来的,并联分配也是——而等压降是你要加进去的约束
是什么(工作点)。泵工作点是泵曲线 Δp_pump(Q) 与系统阻力曲线 Δp_sys(Q) 的交点,由不得你挑(承 A2-02 泵与风机的工作点、特性曲线与匹配,曲线读法与相似定律本课⛔ 不重讲)。在 1D 里它不是“查图”而是“解方程”:系统阻力曲线本身随阀位、温度(μ(T))、结垢与主动进气格栅状态变,工作点会沿泵曲线滑动。⇒ 模型里它必须每一步现解,⛔ 不能在建模时固化成一个流量数。
工程量级(基线水力解)。本课教学算例取泵曲线 Δp_pump = 110 − 0.1·Q²(kPa,Q 以 L/min 计)、串联段 R_s = 0.075、并联段等效 R_par = 0.100(kPa/(L/min)²),于是系统曲线 Δp_sys = 0.175·Q²。令两者相等:
110 − 0.1·Q² = 0.175·Q² ⇒ 0.275·Q² = 110 ⇒ Q = 20.00 L/min、Δp = 70.0 kPa
主管内径 25 mm ⇒ 流速 u = 0.679 m/s;取 ρ = 1050 kg/m³、μ = 1.2×10⁻³ Pa·s(约 60 ℃ 量级的 50 vol% 乙二醇水溶液,⚠ 教学假设值,实际须查所用牌号的物性表,查用口径承 A1-07 热管理常用工质的热物性与查用)⇒ Re = 1.49×10⁴,在湍流区,故取阻力指数 n = 2。
是什么(并联分配)。并联的两条支路挂在同一对节点上,节点压力唯一 ⇒ 两条支路的 Δp 必然相等。⚠ 这是拓扑连接的后果:在网络方程里,它体现为“两条支路的特性方程共用同两个节点压力未知量”。
★ 反钉子第一层:它是你列方程时要加的约束,⛔ 不是求解器附送的定律。很多人把「并联支路两端压降相等」当成一条可以直接取用的结论——“把支路连上去它就成立,所以几条并联支路会自动按 √R 分好”。真相是:网络方程的闭合靠 1.2 那三样东西凑齐,而 1D 工具在方程不定时通常不会拒绝求解,它会用别的东西把方程数补齐——最常见的是你自己在某条支路上加了一个定流量边界,或者工具的默认初值把解锁在了某一个分支上。⇒ 结果是模型照常收敛、照常出一张漂亮的流量分配表,而那张表与水力无关。
判据只有一条,而且很便宜:数一数未知量与方程数(1.2 那三样),并检查有没有哪条支路的流量是“我给的”而不是“解出来的”。
工程量级(基线分配)。基线并联段 Δp = 40.0 kPa,两条支路的阻力系数为 R₁ = 40.0/144 = 0.278(电池冷板支路)、R₂ = 0.625(电驱支路),于是
Q₁ = √(40.0/0.278) = 12.00 L/min Q₂ = √(40.0/0.625) = 8.00 L/min
两者相加恰为 20.00 L/min——这是你手算复核的第一个钩子。
易错点。
- 用「各支路阻力比」直接分流量。Δp ∝ R·Q² 时正确关系是 Q_i ∝ 1/√R_i,⛔ 不是 Q_i ∝ 1/R_i;用错会把分配的差距放大成平方。
- 两条支路指数 n 不同时仍套 1/√R。一条已进层流(n = 1)、另一条还在湍流(n ≈ 2)时,那条公式的前提就不成立了(1.12)。
- 只求出总流量就以为水力解完了。分配还没解。
1.5 配平的两条路:同一个分配比例,结果完全不同
是什么。把 1.4 的水力解接到热上,只差一步除法:并联支路进口同温(都来自低温散热器出口),出口温度 = 进口 + Q̇_i/(ṁ_i·c_p)。取 c_p = 3400 J/(kg·K),两条支路的热容流率是
ṁ·c_p = 714.0 W/K(电池支路,12.00 L/min) 与 476.0 W/K(电驱支路,8.00 L/min)
热源取电池 3.0 kW、电驱 5.0 kW(⚠ 教学假设值),于是
ΔT_电池 = 3000/714.0 = 4.20 K ΔT_电驱 = 5000/476.0 = 10.50 K 全回路 ΔT = 8000/1190.0 = 6.72 K
基线的温度场是:低温散热器出水 57.74 ℃ ⇒ 电池支路出口 61.94 ℃、电驱支路出口 68.24 ℃(回路最高点);两股混合后 (12.00×61.94 + 8.00×68.24)/20.00 = 64.46 ℃,与低温散热器进水一致——这是第二个自查钩子,混合温度对不上就说明你哪一步的流量或温升算错了。
要解决的问题。回路最高点 68.24 ℃ 落在电驱支路出口。想把它降下来,最直觉的动作是“配平”——把分配从 12.00/8.00 调到热负荷比 3.0∶5.0(即 37.5%/62.5%,此时两支路温升相等)。达到同一个分配比例有两条路:
- 方案 A:在电池支路加节流孔,把 R₁ 从 0.278 加到 1.736;
- 方案 B:把电驱支路的流道阻力降下来,把 R₂ 从 0.625 降到 0.100。
为什么两者结果不同。加阻力让系统阻力曲线变陡,工作点沿泵曲线左移,总流量下降;减阻力让曲线变缓,工作点右移,总流量上升。⇒ 分配比例一模一样,而回路拿到的总流量完全不同。
工程量级(三方案对照,全部由教学假设值算出,⛔ 不对应任何平台)。
| 项 | 基线 | 方案 A(加孔板) | 方案 B(降阻力) |
|---|---|---|---|
| 总流量 Q_total | 20.00 L/min | 16.20 L/min(−19%) | 22.67 L/min(+13%) |
| 分配(电池/电驱) | 12.00/8.00 | 6.07/10.12 | 8.50/14.17 |
| 两支路温升 | 4.20/10.50 K | 均 8.30 K | 均 5.93 K |
| 回路最高点 | 68.24 ℃ | 65.71 ℃ | 63.85 ℃ |
| 泵压升 | 70.0 kPa | 83.8 kPa | 58.6 kPa |
| 水力功率 Δp·Q | 23.33 W | 22.61 W | 22.15 W |
| 稳态排热 | 8.0 kW | 8.0 kW | 8.0 kW |
⇒ 方案 B 在每一列上都优于方案 A:最高点低 1.86 K、总流量高 40%、泵压升还低 25.2 kPa——而两者的分配比例一模一样。这就是本课要拿量化证据支撑的那条结论:配平手段的排序(先改流道、再调阀、最后才加孔板,承 A2-01 管网水力学:沿程/局部压降与流量分配)不是偏好,是这张表。
⚠ 两条必须随这张表一起写出来的限定。① 三个方案的泵转速未变,走的是同一条泵曲线,水力功率都落在 22~23 W ⇒ ⛔ 不得写成「方案 B 省泵功」;B 的优势是同样的泵功买到更多流量与更低的温度。② 三个功率数如此接近,是本算例那条泵曲线的产物(三个工作点都落在它的水力功率峰附近),⛔ 不是通则——换一条更平或更陡的泵曲线,结论会变。
能量侧的口径(本讲要立、后面几讲反复用)。表的最后一行三个 8.0 kW 一分不差,这不是巧合而是能量守恒的直接后果:稳态、热源功率给定、且未触发温度保护时,Σ排热 ≡ Σ产热。⇒ 流量、UA、拓扑改变的是达成这个平衡所需的温度水平(各节点温度与回路温差),⛔ 不是稳态排热量本身。⚠ 与它成对的另一句同样成立:给定温度上限下的最大可排热功率(也就是“散热能力”)确实随流量与 UA 增大。两句话都对,混起来就出错——「加流量提高散热能力」对,「加流量提高稳态散热量」错。⇒ 交付时 ⛔ 不得据这张表报「改了之后散热量提高了」。
1.6 「抢到流量」不等于「更凉」
是什么。1.5 的那组数里藏着一条很容易滑过去的判据。基线上,电池支路拿了 60% 的流量(12.00/20.00),却只带 37.5% 的热(3.0/8.0 kW);电驱支路拿 40% 的流量,却带 62.5% 的热。两条支路进口同温 57.74 ℃,出口分别是 61.94 与 68.24 ℃,相差 6.30 K。
为什么要在水力这一讲就点破。读者从「阻力小的支路抢流量」这句正确的话,很容易滑到「抢到流量的支路就凉」——而决定支路温升的是 Q̇/(ṁ·c_p) 这个比值,⛔ 不是流量本身。两条支路的流量比与热负荷比不匹配时,多拿流量的那条反而可能是最凉的,少拿的那条最热。基线正是这个形态。
⇒ 这一条直接决定诊断顺序:看到某条支路水温偏高,先算它的热负荷占比与流量占比,再决定要不要动流量。⚠ 若不先算这一步就照直觉去动阀,很容易把方向做反——比如在“水温高”的那条支路上加节流孔,它流量更少、温升更大,问题被做反了一轮。
易错点。
- 看到某支路水温高就去查它的流量,而没有先算它的热负荷占比。
- 把「流量配到相等」当成配平目标(1.8 正面点破)。
- 忘了「两条支路进口同温」这个前提。串联或有旁通时它就不成立了(第 4 讲)。
1.7 ★ 反钉子第二层:等压降铁律的射程是「各支路真正共用同一对节点」
是什么。1.4 立的那条约束有一个前提:两条支路挂在同一对节点上。而总管/集流管本身有沿程摩擦压降与动压回升时,各支路接入点的可用压差沿流向变化——靠近进液口的先天占优、末端天然吃亏。⇒ 等压降只是“总管压降相对支路压降可忽略”这一退化情形下的近似(承 A2-01 管网水力学:沿程/局部压降与流量分配 已写死的边界,本课⛔ 不重推)。
为什么它在 1D 建模里比在手算里更危险。手算时,“共用一对节点”是一个你知道自己做了的近似;建模时,它是连线的默认结果——把电池包 k 片并联冷板画成挂在同两个节点上的 k 条支路,求解器当然会给出 k 条完全相等的流量,而没有任何一步会提醒你做过这个假设。
⇒ ★ 他以为模型在「验证」分配,实际模型在「假设」分配。这就是读反的那一下。
读者会做错的那个具体动作。一个包内温差超标的问题,1D 模型给出 k 片冷板流量严格相等、出水温几乎相同 ⇒ 工程师据此判「流量分配没问题」,把温差归因到电芯产热一致性或冷板结构,转去改电芯分选与冷板流道;而真因(集流管造成的近端多、远端少)在这个模型里从头到尾不可能出现——那个自由度在连线的一瞬间就被取消了。⇒ 改进方向整个指反,一轮样件白做。
正解。要在 1D 里看到分配不均,必须把总管/集流管本身建成带节点的支路链(各冷板支路接在沿流向不同的节点对上),或者直接用部件级手段定这一层(前指 J3-02 电池包液冷流道 CFD 与流量均匀性)。
⚠ 反向也要防:⛔ 不得由此推出「等压降铁律是错的」。总管压降相对支路压降可忽略时,它是很好的近似;而且它永远是列方程的正确做法。错的只有一件事——把它当成“分配没问题”的证据。
工程量级。判据是「总管压降占支路压降的比例」,而这个比例要多小才可忽略,取决于你要分辨多大的不均 ⇒ ⛔ 属平台相关量,本课不给通用门槛,只给算法:把总管建成带节点的支路,比较建与不建两版的分配差,差小于你要分辨的量级才允许退化。
易错点。
- 据一张「分配完美」的 1D 结果判「流量没问题」——上面那个错误动作的入口。
- 反过来因为怕漏而把每一处并联都建成集流管网络。建模量与收敛难度会失控;判据仍是“你要分辨多大的不均”。
- 把 U 型与 Z 型集管的分配偏向按布置形式直接套结论。它由集管压力沿程曲线的主导项决定(J3-02 电池包液冷流道 CFD 与流量均匀性 已写死)。
1.8 配平的目标函数是各支路各自的限值裕度
是什么。1.5 给了配平的量化后果,1.6、1.7 给了两条容易读反的前提。剩下最后一问:配平到底要配到什么?本课的答案是——目标函数是各支路各自的限值裕度,⛔ 不是「流量相等」,也⛔ 不是「温升相等」(⚠ 这条判据是本课归纳的,须在引用时注明出处)。
为什么“温升相等”不是目标。方案 A 确实做到了两支路温升相等(均 8.30 K),代价是电池支路出口从 61.94 ℃ 抬到 65.71 ℃。若电池的温度限值比电驱严,这次配平就把“一个超温”换成了“两个都逼近限值”——指标看起来更漂亮了,裕度反而更薄。
为什么“流量相等”更不是目标。1.6 已经给出了理由:温升由 Q̇/(ṁ·c_p) 决定,把流量配到相等,等于按一个与热负荷无关的量去配。
决策顺序(可执行)。① 先算这条支路的热负荷占比与流量占比——一张纸就能算,且它决定后面所有动作的方向;② 再按各支路各自的限值裕度定目标,⛔ 不拿“相等”当目标;③ 最后才选手段,三条并列:改流道/降阻力(工作点右移、总流量升)、调比例阀开度、加节流孔(工作点左移、总流量降)——三条对总流量、最高点温度、泵压升的方向互不相同,1.5 那张表就是它们的量化对照。
⚠ 最后一条纪律:配平是一次零和的重分配,外加一笔水力损失,⛔ 不凭空创造裕度。给一条支路加了流量,另一条就少了;⇒ 每次配平之后必须回头复核另一条支路的裕度有没有被吃掉。
工程量级。⛔ 各支路的温度限值、裕度目标与分配不均匀度的目标值不给数——它们由电池、电驱、座舱各自的限值与本项目的裕度策略给定,属平台相关量;均匀度目标的回标口径前指 J3-02 电池包液冷流道 CFD 与流量均匀性。本课只负责把“目标函数是什么”与“决策顺序”定下来。
1.9 电子水泵是一族曲线加一条控制律,不是一个可调的流量输入
是什么。电子水泵给出的是一族随转速/PWM 参数化的 Δp-Q 曲线。在模型里,泵是“一条随控制量变形的曲线”,工作点仍然是求交解出来的;按需流量则是“这族曲线 + 一条控制律”。相似定律(Q ∝ N、Δp ∝ N²、P ∝ N³)用于在族内换算,其边界(Re 跨度大、效率岛随转速漂移时会偏)承 A2-02 泵与风机的工作点、特性曲线与匹配,本课⛔ 不重讲。
为什么要在建模层专门强调。把泵建成“给定流量源”是最常见的偷懒法,它同时毁掉两件事:① 支路分配对总流量的响应(1.5 那三个方案的差别全部消失,因为总流量被你钉死了);② 阀切换后工作点沿泵曲线的自然滑动(第 4 讲切串联时,总流量本该大幅下降)。⚠ 这正是 1.1 第 ① 条自由度“不建 ⇒ 模型不能回答「改了会怎样」”的具体形态。
与“加泵速”的关系。「靠加泵速解决水温高」有没有用、有多大用,取决于控制热阻在哪一侧——这一条第 2 讲会给出量化的收益上限,本讲只到“泵是一族曲线,加转速改变的是工作点位置”。
工程量级。⛔ 具体的转速—流量标定表属平台相关量,本课不给数。本课教学算例只演示“同一条泵曲线上工作点随系统阻力滑动”(1.5 的三个方案,转速全程未变)。
易错点。
- 用相似定律外推到远离标定转速的区段而不做台架校核。
- 把 PWM 占空比线性映射成流量。
- 在低转速段用高转速标定的效率反算泵功。
1.10 定压点是模型的一个自由度,不是一个部件
是什么。不可压网络的方程里只出现压差,因此压力场有一个任意常数——解出来的每个节点压力都可以整体加上同一个数而方程仍然成立。膨胀壶/去气罐的接点就是把这个常数钉死的地方。⇒ 它在模型里的身份不是“一个部件”,而是 1.2 里的第 ③ 样闭合条件。
不建它会漏报哪一整类结论。不设参考点,流量与温度照样解得出来(它们只依赖压差),而一切依赖绝对压力的判据——气蚀、局部沸腾、压力盖开启、承压校核——全部失效。⚠ 注意这句话的形态:不是“算得不准”,是“这一类判据在这份模型里根本没有依据”。
为什么接点位置在模型里是结构性的。把参考点从泵入口挪到泵出口,全场压力就整体平移一个泵压升;而汽蚀余量只看泵入口那一个绝压。⇒ 同一份模型、同一组流量与温度,接点画错位置就能把气蚀结论翻过来。
工程量级。⛔ 压力盖开启压力、膨胀壶安装高度、接点到泵入口那段的沿程与局部损失都属平台相关量,本课不给数。本课只给结构判据两条:① 参考点必须画在实车膨胀壶真正的接点上;② 接点到泵入口那段主流管路的损失要老老实实建出来——它是从泵入口绝压里扣掉的(G10-02 膨胀水壶/去气罐(degas)设计与排气 已写死「接点位置的选择本质是拿 Δp_f 换汽蚀余量」)。裕量判据与失效机理归 G5-04 泵的气蚀、空气分离与可靠性,去气罐容积、液位与压力盖选型归 G10-02 膨胀水壶/去气罐(degas)设计与排气,⛔ 本课不给限值与裕度门槛。
易错点。
- 全场按表压建模、参考点随手放在“回路里某个点”。
- 同时设了两个定压点。过约束(1.13 第 ④ 条)。
- 把定压点当成“给系统加压的手段”。它只钉住一个数,⛔ 不产生流量。
1.11 不可压刚性回路里没有「气」这个自由度
是什么。常规冷却液回路模型把液体当不可压、把气腔当不存在。于是节点压力可以被解成低于饱和蒸气压、甚至解成负绝压,而方程组仍然完全自洽、迭代照常收敛。⇒ 气蚀、气堵、吸空在这份模型里没有承载它们的状态量,永远不会出现。
为什么它是本课钉子最锋利的一个实例。气蚀在实车上是一个响亮的失效:噪声、流量塌陷、叶轮点蚀,谁都跑不掉。而在模型里它是完全静默的——它不表现为报错,表现为“没有这个问题”。⚠ 这两种表现之间没有任何过渡地带,也没有任何一个求解器会替你把它们区分开。
正解:补自由度 + 后处理自检,两件事都要做。要在模型里看到这一类现象,必须先补上气腔/可压缩性这个自由度(1.1 第 ③ 条);而无论建不建,都必须做后处理自检——把每个节点、每一步的绝对压力与该点温度下的饱和蒸气压 p_v(T) 相减,出现非正裕度就报出来。⚠ 判据本身不含任何门槛数:它只回答“有没有越过 p_v”,⛔ 不回答“该留多少裕度”(那属平台相关量,机理与裕量判据归 G5-04 泵的气蚀、空气分离与可靠性)。
工程量级。⛔ 饱和蒸气压 p_v(T) 随冷却液牌号与乙二醇体积比变化,本课不给数,只给来源:所用牌号的物性数据表。⚠ 冷却液本身的性能要求与试验方法按 GB 29743 系列现行版本核——本课未取原文,只写去向,版次与年份以现行目录为准、引用前须核。★ 另需说明一条射程(这是本课自己作出的判定,须注明出处):GB 29743 管的是冷却液这个产品,⛔ 管不到回路仿真方法与允差;⇒ 它能给你物性边界,⛔ 给不了任何建模判据。
易错点。
- 指望求解器报警。它不会,因为方程组是自洽的。
- 只在稳态查一次。最危险的是泵启动、阀切换、下坡回充这类瞬态时刻。
- 把「模型没算出负压」当成「不会气蚀」。模型里根本没有最低液位与气腔——第 ③ 条自由度没建,这个结论就没有依据。
1.12 低温是水力模型的第二个工况族,不是同一个工况换个温度
是什么。Darcy 摩擦系数在湍流段近似 f ∝ Re^(−0.25)(Blasius 型)⇒ Δp ∝ Q^1.75~2;层流段 f = 64/Re ⇒ Δp ∝ Q^1(非圆截面按水力直径 d_h = 4A/P,⚠ 扁流道沿用 64 会把阻力低估到约 2/3,A2-01 管网水力学:沿程/局部压降与流量分配 已写死,本课⛔ 不重推)。⇒ 温度变化不只是让某个系数变大变小,它会换掉阻力律的指数。
为什么它影响的是分配而不只是总压降。并联分配比例由两条支路的指数与系数共同决定。若一条支路已进层流、另一条还在湍流,比例会随总流量与温度双向漂移——⛔ 1.4 那条 Q_i ∝ 1/√R_i 的前提(两条支路同指数、n = 2)已经不成立了。⇒ 你在常温下配好的分配,在低温下不是“稍微差一点”,而是另一组分配。
工程量级。本课教学算例常温工况:主管 u = 0.679 m/s、d = 25 mm、μ = 1.2×10⁻³ Pa·s(约 60 ℃ 量级)⇒ Re = 1.49×10⁴,湍流。低温(−20 ℃ 量级)50 vol% 乙二醇水溶液的动力黏度高约一个数量级——⛔ 只给方向与量级,具体值须查所用牌号的物性表——同一条管路的 Re 落到 10³ 量级,已在 2300 以下的层流区。
⚠ 层流/过渡/湍流的分界要按区间用:2300~4000 之间是典型灰区,A2-01 管网水力学:沿程/局部压降与流量分配 已写死“灰区里关联式都不可靠”。⇒ ⛔ 不得把 2300 当成一个可以据以切换关联式的确定门槛;模型在灰区内的结果只作趋势,且必须在交付说明里声明这一点。
易错点。
- 全程用一组定温物性。1.1 第 ⑨ 条自由度没建的具体形态。
- 用常温标定的泵曲线算低温工况(A2-02 泵与风机的工作点、特性曲线与匹配 已写死)。
- 在 2300~4000 的灰区里取任一关联式当准值。那一段的结果只能当趋势。
1.13 「不收敛」与「收敛到一个错解」是两件事
是什么。水力域专属的不收敛与错解成因有七条(⚠ 本课按网络方程的闭合条件逐项推得,须注明出处是本课):
① 并联支路的等压降约束漏设(方程欠定); ② 泵曲线被外推到曲线之外(Q < 0 或 Q > Q_shutoff 段,工具会线性外推出一条无物理意义的曲线,还会给出“解”); ③ 阀全关后支路成孤立节点(该节点的压力不再被任何方程决定); ④ 定压点缺失、或不小心设了两个(前者压力场不定,后者过约束); ⑤ 零流量与流向反转处的奇异(Δp ∝ Q² 在原点导数为零,牛顿迭代的雅可比在那里退化,稳态求解器常卡在“几乎不流”的支路上); ⑥ 状态相关阀(止回阀)造成的解不唯一或在两状态间抖振; ⑦ 支路给定流量边界与泵曲线同时存在(过约束)。
⛔ 只沿“并联约束、泵曲线外推、阀全关”这三条最常被提到的线索枚举,会漏 ④⑤⑥⑦——④⑦ 是“约束数不对”的另外两种形态,⑤ 是最常见的原点奇异,⑥ 是模型内生的离散。
★ 更要紧的是那条分岔:「不收敛」是响亮的,「收敛到一个错解」是安静的。①④⑦ 三条既可能让求解器失败,也可能让它照常收敛——欠定时被默认初值或某条定流量边界补齐,过约束时被工具用最小二乘或优先级规则化解。⇒ 右边这一支必须主动去查,⛔ 不能等工具提醒。两支的处置方向也相反:左支多为数值与结构奇异,模型侧有规避法;右支是约束数不对,靠调容差治不好。
工程量级。⛔ 无门槛数。给三条可执行的模型侧规避:
① 给每条支路一个极小的最小流通面积(避免全关成孤立节点),或者在阀全关时把该支路整体从网络里摘掉; ② 在 Q → 0 附近给支路特性一段线性化(把 Δp = R·Q|Q| 在原点附近换成 Δp = R·ε·Q)——⚠ 这是数值处理,它会改变极小流量段的解,必须在交付说明里写明换过; ③ 泵曲线在 Q < 0 与 Q > Q_shutoff 段显式给出外推规则,⛔ 不让工具默认线性外推。
射程。通用数值机制——DAE 与指标、一致初值、rtol/atol 取法、事件定位、守恒残差自检——是所有 1D 模型共有的,归 J7-08 一维系统模型的数值内核与「跑飞」排查:DAE/代数环、步长与容差、一致初值与守恒自检,本课⛔ 不重讲;本课只负责“它在水力域长什么样、从模型侧怎么规避”。⚠ 两条判据链互不推导:⛔ 不得用数值侧的结论去判水力结构,也⛔ 不得反过来。
易错点。
- 把「收敛了」当成「解对了」。过约束与欠定都可能给出一个收敛的错解。
- 靠调容差治结构性问题。
- 把工具报的「非线性方程组求解失败」直接当成「模型物理有错」。它常常只是初值问题,诊断入口在 J7-08 一维系统模型的数值内核与「跑飞」排查:DAE/代数环、步长与容差、一致初值与守恒自检。
本讲小结:水力层要立的第一件事,是把流量解出来而不是给出来。① 网络方程的闭合靠三样——节点守恒、每条支路的 Δp-Q 特性、一个绝对压力参考点;三样缺一样方程组就不定,而工具通常不会拒绝求解 ⇒ 判据是“数一数,并检查有没有哪条支路的流量是我给的”。② 支路特性有七类来源,手册 K 值只覆盖其中一类,集成流道那一类手册一个字都没有,本课只定接口口径、⛔ 不给数。③ 工作点是泵曲线与系统曲线的交点,本课教学算例基线 Q = 20.00 L/min、Δp = 70.0 kPa、分配 12.00/8.00 L/min,全部可在一张纸上复核。④ ★ 反钉子两层:「并联支路 Δp 相等」是你要加进方程组的约束而不是求解器附送的定律(第一层);它的射程是“各支路真正共用同一对节点”,把 k 片冷板挂在同两个节点上就等于在模型里取消了流量不均这个自由度(第二层)⇒ 他以为模型在验证分配,实际模型在假设分配;⚠ 而反过来 ⛔ 不得推出“等压降铁律是错的”,它永远是列方程的正确做法。⑤ 配平的目标函数是各支路各自的限值裕度,⛔ 不是流量相等、⛔ 不是温升相等;同一个分配比例,用加阻力达到与用减阻力达到,结果全面不同(最高点 65.71 ℃ vs 63.85 ℃),⚠ 而三方案水力功率都在 22~23 W ⇒ ⛔ 不得写“省泵功”。⑥ 三方案稳态排热恒为 8.0 kW ⇒ 流量与拓扑改变的是达成平衡所需的温度水平,⛔ 不是稳态排热量本身。⑦ 「抢到流量」不等于「更凉」——温升由 Q̇/(ṁ·c_p) 决定,60% 的流量只带 37.5% 的热。⑧ 定压点与含气是两条常被当成“部件”的自由度:前者不建则一切绝压判据失效,后者不建则气蚀在模型里恒不出现,而两种情形下模型都照常给出漂亮的流量与温度。⑨ 低温是第二个工况族,μ(T) 会换掉阻力律的指数、让分配比例本身漂移。⑩ 水力域专属的不收敛有七条,而最贵的那一类不是“不收敛”,是“收敛到一个错解”。⇒ 水力解出来之后,下一个自由度是换热:UA 到底是部件属性还是工况量、加流量的收益上限在哪——那是第 2 讲。
后面还有 5 讲正文 · 关键公式 · 案例拆解 · 常见误区 · 动手做