TMS BOOK · ACADEMY 讲义

1D 热管理仿真导论(KULI/AMESim/GT/Flowmaster)

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

J1-01 1D 热管理仿真导论(KULI/AMESim/GT/Flowmaster)

课程代码 J1-01 · 板块 J 仿真分析 CAE / J1 一维系统仿真 时长 约 3.5 小时(6 讲 + 1 次最小回路实操) 适合对象 CAE / 系统仿真工程师;本课是 J1 模块的起步课,不是零前置的入门课 前置 A1-06 热阻网络法与集总参数建模入门;A2-01 管网水力学:沿程/局部压降与流量分配。正文还会用到 A2-02(泵与风机的工作点)、A3-01(ε-NTU 法)、A1-07(工质热物性与查用),建议一并备好 本讲义定位 讲师授课蓝本 / 学员自学讲义,是 J1-01 大纲的完整展开版


引言:控制体从 20 个加到 2000 个,横截面上仍然只有一个温度

一维系统仿真在整车热管理里的位置,一句话就能说完——它是数字台架:不用实车、几分钟跑一个工况,可以把高温快充、冬季续航、瞬态 pull-down 一次扫遍。真正卡住人的从来不是这句话,而是它后面的三个具体场合。第一个场合:模型搭好了也跑通了,有人来问「电芯上哪一颗最热」「冷板里哪条流道流量最少」,于是开始加控制体、压残差、把时间步往小里调,一轮一轮下去,那个数始终不出现——不是调得不够,是这个模型里根本没有承载那个答案的变量。第二个场合:评审桌上有人问「你这个结果能信到什么程度」,答出来的是求解器的残差降到了多少——而残差回答的是「这组方程解没解对」,它对「这组方程写得对不对、里面的参数是从哪借来的」一个字都没说。第三个场合:结果与试验对不上,第一反应是减时间步、加控制体、把容差压紧,一周之后曲线几乎没动——因为差的那一部分从头到尾记在别的账上。这门课要处理的,就是这三件事。

本课交出四项判断力,每一项都对应一个当场可验的动作。其一,判这个问题该不该用 1D——交付物是一句话写清的响应量,外加「本模型答得出与答不出的问题各列三条」。其二,看懂四款主流工具的建模范式差异并按项目选型——交付物是一份按六个维度写成的提问清单,每一问都要求对方回答「哪个版本、怎么验证的」,⛔ 而不是一张打勾打叉的能力表。其三,搭一个最小可运行的冷却回路模型并跑通稳态——交付物是工作点流量与各节点温度,而且这几个数你自己能用一张纸复算出来。其四,判 1D 结果的可信区间——交付物是一张按敏感度排序的参数表,以及一份「哪些量必须靠试验或 3D 标定、哪些量本模型永远给不了」的清单。

全课六讲,顺序就是一个人真正建一个模型的顺序。第 1 讲把 1D 与 3D 的分界线画清楚,并把「建模层级」从「1D 还是 3D」的两分展开成一条八档的谱;第 2 讲讲建模范式——0D 集总、1D 分布式、准 2D 分层网络三档各自的自由度,因果式信号流与非因果式物理网络的区别,端口上到底守恒的是什么,以及物性用预插值表还是在线调用;第 3 讲做工具选型,把差异拆成六个维度,再讲四款工具各自的血统决定了它们擅长回答什么问题;第 4 讲是求解器与数值基础——稳态水力工作点怎么求、1D 热流网络在数学上到底是什么形态、刚性系统为什么必须用隐式变步长、CFL 在 1D 里的两种身份,以及初值给不好会污染整段瞬态的机理;第 5 讲讲落图,从架构图到拓扑网络的六步,定压点该放哪,以及空气侧边界为什么不能拍成常数;第 6 讲收口到结果可信度——敏感度排序、UA 的热阻拆账、误差指标的分母、参数可辨识性,最后是「结果与试验对不上」的六类原因与排查顺序。其后的典型案例把六讲的判据合成一条链,动手做则用一条最小回路把上面四条目标各验一次。

钉子①(全课主线):选 1D 还是 3D,第一问不是「要多准」,而是「我要的答案是不是一个控制体平均量」。1D 模型把流道沿流向切成一串控制体,每个控制体上写质量、动量、能量三条守恒,未知量是这个控制体的截面平均压力、平均质量流量与平均温度——横截面方向上不存在离散,也就不存在「截面内的分布」这个概念。所以「哪一点最热」「哪条流道流量最少」「壁面膜温多高」「气流在哪里分离」这几类问题,1D 不是算不准,是自由度根本不存在:把控制体从 20 个加到 2000 个,沿流向的分辨率确实提高了两个数量级,而横截面上永远还是一个值。⇒ 由此得到本课反复使用的选型判据:先把要回答的问题写成一个响应量,再问它是不是某个控制体上的平均量。是 ⇒ 1D 是正解,不是将就;不是 ⇒ 加控制体、减时间步、压残差、换更贵的求解器全都无效,只剩三条出口——换方法(3D / 共轭传热)、换分工(1D 出边界条件、3D 看局部)、换问题(先判平均量能不能替代)。这道判据在第 1 讲 1.2 节末画成一张自上而下的判断链,从那里往后,每一讲在「这一层要不要建」上都会回到它。⇒ 由钉子① 还要顺带推出 1D 与 3D 的关系:两者不是精度阶梯,是分工——三方之间真正流动的只有两类东西,3D 与试验往 1D 喂降阶参数(UA/ε、阻力系数、多孔系数、风扇曲线、热容),1D 往 3D 喂边界条件(各支路流量、进出口温度、该看哪个工况的哪个时刻)。⛔ 这根钉子直接派生出本课的一条硬纪律:本课不给任何「1D 比 3D 差几个百分点」式的精度对比——那种说法预设了两者在回答同一个问题,而钉子① 说的恰恰是它们回答的不是同一个问题。

钉子②(与钉子① 配套):一个收敛的 1D 结果,它的可信区间由降阶子模型的参数定,不由求解器定。1D 里的换热器不是被算出来的,是被 UA 或 ε-NTU 代表的;管路不是被算出来的,是被阻力系数代表的;泵与风扇不是被算出来的,是被一条特性曲线代表的。这些参数全都从别处借来——试验、3D、供应商数据表——而借来的东西不会比它的源头更准。⇒ 由此推出结果对不上时的排查顺序:先查拓扑,再查参数,再查边界,再查取值口径,最后才查数值。⛔ 反过来做(先压残差、先减时间步、先加控制体)是本课要拦下的最常见动作,它把时间花在通常最小的那一本账上。⇒ 再推出「该标定什么」的判据:由敏感度排序决定,⛔ 不由「哪个参数我最没把握」决定。本课的教学算例里,四个参数各扰动一次之后,对回路最高温的影响排序是热源功率 > 散热器 UA > 空气侧风量 > 系统阻力系数——那个大家最没把握、最想去翻手册查一查的阻力系数,恰恰是最不要紧的一个。⚠ 这个排序是本算例的结论,⛔ 不是通用排序;具体数值、扰动口径(解耦与耦合两种)以及它们为什么必须并列给出、⛔ 不得只报一个,都在第 6 讲 6.1。

⚠ 两根钉子是配套的,但不得互相推导:钉子① 回答「这个问题该不该用 1D」,判据是响应量的类型;钉子② 回答「用了 1D,结果差在哪」,判据是参数的敏感度与来源。一个模型可以完全满足钉子①(问的确实是平均量)而仍然错得离谱(钉子② 那一侧的参数全是拍的);也可以参数极准而答非所问。两者各由各的判据成立,谁也不是谁的推论。

★ 本课最容易被读反的一条与瞬态步长有关,第 4 讲 4.3~4.6 会把它连同代价一起当场算破,这里先把正读法摆到台面上:热时间常数 τ 是输出侧的时间尺度——它定的是仿真要跑多久、数据存多密、什么时候可以判稳;步长是输入侧的数值参数——它由局部误差容差、当前活跃的最快模态与离散事件定位三者定。两者之间没有换算关系。一个热管理回路里同时活着一族时间尺度:回路集总热容的响应、管路输运的纯延迟、单个控制体的输运时间、控制器的采样周期,以及阀切换与泵启停这类本身没有时间尺度、但步长必须在它上面落一步的离散事件。变步长积分器盯的是这一族里最快的那几个,⛔ 从来不是最慢的那一个。⚠ 更麻烦的是,步长取错在 1D 里不会以报错的形式暴露:隐式求解不发散、不报警、残差照样降到底,你拿到的是一条又光滑又好看、而快过程整段已经没了的曲线。第 4 讲会把这条错误读法派生出的四个具体动作逐个点破,并给出三条必须同时满足的可执行判据,同时把「那就用小步长」这半句话的代价一并算给你看。

全课共用一组教学假设值,本讲义里统一称它 T1。它存在的唯一理由,是让本课出现的每一个结论都能被读者用一张纸复算出来,并且彼此严格自洽——同一个量在正文、图、图注与练习里同值同舍入。⛔ T1 的每一个数都是教学假设值,不对应任何平台,禁止照抄取用;⛔ 它们不得以「典型值」「参考值」的名义被引用,也不得被搬到别的课去。

T1 取值(⛔ 教学假设值,不对应任何平台)
冷却液 50 vol% 乙二醇水溶液、约 60 ℃ 量级:ρ = 1050 kg/m³、c_p = 3400 J/(kg·K)、μ = 1.2×10⁻³ Pa·s
回路 总容积 6 L;金属件等效热容 8000 J/K
管路 主管内径 25 mm;当量总长 6 m;局部损失系数之和 ΣK = 12
泵与系统阻力 泵曲线 Δp = 100 kPa×(1 − Q/40 L/min);系统阻力以「20 L/min 时 60 kPa」定 K
散热器 UA = 350 W/K;空气侧质量流量 0.60 kg/s(干空气 c_p,air = 1006 J/(kg·K))
边界 热源功率 5.0 kW;进风温度 40 ℃
离散 控制体数 N = 20

⚠ 关于这张表要先说三句。第一,冷却液物性的来源级别只到工程手册与物性库,本课未取标准原文,实际取值随乙二醇体积比与工作温度显著变化,须查所用牌号的物性表(查用口径承 A1-07 热管理常用工质的热物性与查用);黏度对温度尤其敏感,换一个工作温度必须重取。第二,金属件等效热容这一项在稳态里完全不出现,只在瞬态里主导——只做过稳态的人最容易把它整项漏掉,而漏掉它系统会显得比真实的「快」。第三,表里没有出现的东西同样是有意的:本课在多处会先画一份「全集」再逐项核射程,一共六份,正文里按序号标出、后文引用时照序号找:全集(一)建模层级八档(第 1 讲 1.1)· 全集(二)1D 元件八类(第 2 讲开篇)· 全集(三)求解的数学形态六种(第 4 讲 4.1)· 全集(四)结果对不上的六类原因(第 6 讲 6.6)· 全集(五)工具选型六维度(第 3 讲 3.1)· 全集(六)验收证据七类(第 6 讲 6.7)。这六份全集都是本课按工程逻辑归纳的,⛔ 不是任何标准或手册的既有清单,也不是行业公认划分,引用时请连同这句一起带走。

本课对数字有三条固定处理:物性与标准值给值并写清来源级别;行业典型区间只给区间并标「典型」,⛔ 端点不得当成可承诺值;平台相关量一律不给数。第三类在本课里有四组,先声明在这里,免得读到正文时以为是漏写:① 1D 与 3D 的单工况机时,以及物性在线调用相对查表的耗时倍数——⛔ 不给数,须在本团队的机器与模型上实测;本课只给结构结论:两者相差若干个数量级,足以让「工况数」从一项约束变成一个自由变量。② 四款工具的能力评分、版本能力清单、性能跑分与许可价格——⛔ 一律不给:来源级别只到供应商数据表且本课未独立核对,何况它们随版本与商务条件变;本课交付的是六维提问清单,⛔ 正文与图上都不会出现任何工具的分数、排名与勾叉表。③ 能量不平衡自检的阈值——⛔ 不给,须由项目定;本课只给「报不平衡量时分母必须写清是全系统总输入功率还是最大单支热流」这条做法。④ 案例里的电芯生热功率、电池入口水温上限、压缩机可用制冷量——⛔ 不给,案例只交付判据结构(逐段温差拆账 → 找最大那一段 → 判瓶颈)。

本课引用三份标准:ASME V&V 20(计算传热的验证与确认)、SAE J2765(车用空调系统台架 COP 测量)、GB/T 18386.1(电动车能耗与续驶里程的工况边界)。⚠ 三份的版次与年份以现行目录为准,引用前须核;⛔ 本课三份都未取原文,条款内容与限值一律不写,只写去向。⚠ 其中 SAE J2765 的射程是移动式空调系统在台架上的 COP 测量规程,⛔ 不得跨射程用到冷却液回路、电池回路与整车工况上去。

最后把本课的边界摆出来,写清楚是为了让你知道哪些东西不在这里找数值内核——微分-代数方程组的指标与代数环、相对/绝对容差的工程取法、一致初值的求法、事件定位、守恒残差自检——全部归 J7-08 一维系统模型的数值内核与「跑飞」排查:DAE/代数环、步长与容差、一致初值与守恒自检,本课只讲到「它们存在,以及它们如何决定步长与初值」这一层,⛔ 不给任何取值建议。仿真与试验相关性验证的完整流程(测点选择、误差指标定义、验收阈值)归 J7-01 仿真—试验相关性验证(Correlation),本课只交付「两侧取值口径必须对齐」这一条纪律和两个可当场判真假的例子。不确定度的量化传播J7-06 仿真不确定度量化(UQ)与稳健/可靠性设计:从确定性单点到带分布的裕度,本课只到敏感度排序这一层。3D 的方法本身属 J2 板块(入口 J2-01 CFD 基础、湍流模型与网格划分(Fluent/STAR-CCM+)),本课只作分工声明。试验侧的台架不确定度与测量方法L7-01 测量技术:热电偶/红外/PIV/流量测量L2-03 冷却系统台架与流量分配试验,本课只声明「它存在,且必须与仿真侧的取值口径对齐」。至于具体回路怎么建,本课一条都不重讲:制冷循环与制冷剂物性在 J1-02 制冷循环一维建模与制冷剂物性,冷却液回路的水力与热在 J1-03 冷却液回路水力与热仿真建模,热泵与模式切换在 J1-05 热泵系统一维建模与模式切换仿真,电池与快充温升在 J1-06 电池热管理一维建模与快充温升仿真,座舱多节点在 J1-08 座舱一维多节点热模型与参数辨识(含 CFD→1D 分区映射),把它们联起来跑在 J1-04 整车热管理系统级联合仿真,模型作为资产怎么管在 J7-07 仿真任务的输入、交付与模型资产管理。本课交付的是范式、判据与流程——把它学扎实,上面这些课你会用同一套判据去读。

第 1 讲 1D 答得出什么、答不出什么:第一问不是「要多准」

这一讲不教任何一款软件的操作。它要把一个几乎每个项目都会被问、而多数人答错方向的问题掰开:这件事该用 1D 还是该上 3D。绝大多数人是按「要多准」来答的——精度要求高就上 3D,精度要求低就用 1D。这条路径从第一步就走岔了:1D 与 3D 的分界根本不在精度,在自由度。它们回答的不是同一类问题,因此也不存在「谁比谁准几个百分点」这回事。本讲先把建模层级的全谱摊开(1.1),再给出唯一的那一问判据(1.2),然后把判「否」之后读者最想做的那个动作——加控制体——当场算破(1.3);接着交代 1D、3D 与试验三方之间到底有哪两类东西在流动(1.4),并由此立起本课的第二根钉子(1.5);最后两节回到正面,说清 1D 真正强的两个主场(1.6、1.7),并把第一根钉子推到验收那一层(1.8)。后面五讲的每一个判断——建什么、用哪款工具、步长怎么定、拓扑怎么落图、结果信到什么程度——都建在这一讲的两根钉子上。

1.1 建模层级不是「1D 或 3D」两档,是一条八档的谱

是什么(全集(一))。把整车热管理里用得到的建模层级摊开,是八档而不是两档。⚠ 这八档是本课按「先画全集、再逐项核射程」的做法归纳出来的,⛔ 它不是任何标准、任何手册、也不是任何既有手册里现成的分类,引用时须注明出处是本课:

稳态手算与表格法——热阻串并联、ε-NTU 手算、Excel 能量平衡(承 A1-06 热阻网络法与集总参数建模入门A3-01 换热器传热基础:ε-NTU 法与 LMTD 法)。它不是「简陋版 1D」,而是在方案筛选与量级判断阶段唯一够快的工具。 ② 0D 集总参数——一个部件一个热容节点,没有空间维度,只有时间维度。 ③ 1D 分布式——沿流向离散成控制体链,横截面取平均。本课的主体。 ④ 准 2D / 分层网络——换热器按流程分区、电池包按模组分支路、座舱按分区节点。形式上仍是 1D 网络,但在拓扑层面表达了空间不均。 ⑤ 3D CFD / 共轭传热——J 板块的 CFD 一支(J2-01 CFD 基础、湍流模型与网格划分(Fluent/STAR-CCM+) 起)。 ⑥ 降阶模型 ROM——从 3D 或试验数据蒸馏出来的代理(响应面、POD、线性状态空间、数据驱动代理)。 ⑦ 纯经验查表模型——直接用台架 map 插值,不写任何守恒方程。 ⑧ 实时 / HIL 模型——建模选择由「每一步必须在固定时间内算完」倒推:定步长、禁迭代、禁事件重启。

每一档的区别,说到底是它的代表变量是什么:0D 是「一个部件一个温度」,1D 是「沿流向一串截面平均量」,准 2D 是「一组分区平均量」,3D 是「一个场」。这一行字比档位名更有用——你要的答案是哪一种形态,就该落在哪一档。

为什么准 2D 值得单独列一档。它是「在 1D 框架里回答空间不均问题」的正解,而不是一个折中方案。它的做法是用拓扑表达空间不均:每一区内部仍然是集总的,但区与区之间有各自独立的流量与温度。于是「我需要知道模组之间差多少」这个问题,就从「需要一个场」降格成了「需要几个平均量」——重新落回下一节那条判据的射程之内。这一档在 1D 工具里的具体建法,座舱那一支前指 J1-08 座舱一维多节点热模型与参数辨识(含 CFD→1D 分区映射),电池包那一支前指 J1-06 电池热管理一维建模与快充温升仿真

工程量级。⛔ 分几区、区界划在哪,是平台相关量——它取决于流道布置与你要回答的那个问题,本课不给数,须由所回答的问题倒推。本课只给判据:分区数由「你要分辨的最小差异单元」定,不由网格习惯定。要分辨模组之间的差异就按模组分,要分辨换热器进出流程的差异就按流程分;⛔ 不存在一个通用的「切几段比较好」。

易错点。

  • 只沿「1D vs 3D」这条最熟悉的线索枚举。这样枚举会整类漏掉 ④⑥⑦⑧ 四档,而漏掉的方式不是「少了一条」,是那一整类根本不出现在选项里——清单看起来是完整的。直接后果:遇到「模组温差」就上 3D,成本涨两个数量级,而结论只需要几个平均值。
  • 把分区数一路加下去当成「细网格」用。加到最后既没有 3D 的物理(横截面上仍然只有一个值),也没有 1D 的速度。
  • 忘了 ⑥⑦ 两档其实已经在你的模型内部。1D 模型里的换热器 map、压缩机 map,本质上就是代理模型与查表;于是一个半查表的模型被当成纯机理模型去外推,而查表模型的边界外行为完全不可信(ROM 的构建方法前指 J7-02 模型降阶(ROM)与控制导向建模,数据驱动那一支前指 J7-05 机器学习/AI 代理模型在热管理中的应用)。

本课的射程。①②③④ 在本课射程内,作为选型谱讲清并给判据;⑤ 只作分工声明,方法全部在 J 板块的 CFD 一支;⑥⑦ 只作类型声明与去向,⛔ 本课不给构建方法与精度指标;⑧ 只声明它的取舍逻辑与 ①~⑦ 方向相反(那边是由算力预算倒推模型,这边是由问题正推模型),⛔ 不给实时化的实现手段,前指 J7-04 整车热管理数字孪生与实时模型

一张定位图,横轴是空间分辨率(从无空间维度到三维场,定性轴无刻度数值),纵轴是可覆盖的工况数乘以可跑的时间跨度(同样是定性轴无刻度数值)。图上八个方框分别标注稳态手算与表格法、0D 集总、1D 分布式、准 2D 分层网络、3D CFD 与共轭传热、降阶模型 ROM、纯经验查表、实时 HIL,每个方框下方一行小字写出它的代表变量,例如 0D 是一个部件一个温度、1D 是沿流向一串截面平均量、准 2D 是分区平均量、3D 是场。两个不同底色的色块把本课射程内的前四档与只作类型声明与去向的后四档分开,后四档各带一个箭头指向课号标签。本图为定性示意,两轴均无刻度,不代表任何平台的实际值,不得据图读取任何数值或倍数
图1 建模层级的八档谱与各自的代表变量。该读出三件事:① 1D 与 3D 之间还夹着一档准 2D 分层网络,它是在 1D 框架里回答空间不均问题的正解,而它恰是只按两分法枚举时必漏的那一类;② 沿谱向右走,换来的是空间分辨率、付出的是工况数与时间跨度——这是一次交换,⛔ 不是一次升级,因此谱上没有「更高级」的一档;③ 后四档不在本课射程内,各自的去向已标在图上。⚠ 八档的划分是本课归纳的,⛔ 不是任何标准或既有手册里现成的分类。本图为定性示意图,两轴均无刻度,方框之间的距离不代表任何倍数关系,⛔ 不得据图读取任何数值或倍数。

1.2 判据只有一问:把问题写成一个响应量,再问它是不是控制体平均量

是什么。1D 模型把流道沿流向切成控制体,在每个控制体上写质量、动量、能量三条守恒。它求解的未知量,是这个控制体的截面平均压力、截面平均质量流量、截面平均温度。横截面方向上不存在离散,因此模型里根本没有「截面内的分布」这个概念——不是分辨得粗,是没有这个变量。

于是本课的第一根钉子可以写成一句话:选 1D 还是 3D,第一问不是「要多准」,而是「我要的答案是不是一个控制体平均量」。

由它推出的判据只有一步:先把要回答的问题写成一个响应量,再问它是不是某个控制体上的平均量。是 ⇒ 1D 是正解,而不是将就;不是 ⇒ 加控制体、加时间步、压残差、换更贵的求解器全都没用,必须换方法或换分工。

为什么这是自由度问题而不是精度问题。「精度不够」意味着这个量算出来了、只是值不对,可以靠更细的离散、更小的容差、更好的格式往真值上逼。而「哪一排扁管最热」这类问题在 1D 里连一个对应的未知量都没有——没有哪个数可以被逼近。把控制体从 20 个加到 2000 个,沿流向的分辨率确实提高了 100 倍,横截面上仍然只有一个值。一个不存在的维度,不会因为别的方向算得更细就长出来。

⛔ 由这一条派生出本课的一条硬纪律:本课不给任何「1D 比 3D 差几个百分点」式的精度对比。那种说法预设了两者在回答同一个问题,而这一节说的正是——它们回答的不是同一个问题。

教学算例 T1(全课共用的一组教学假设值) 本课后面每一个带数的结论,都出自同一组教学假设值。⛔ 它们不对应任何平台、任何车型、任何供应商产品,⛔ 禁止照抄取用。它们存在的唯一理由,是让本课的每一个数都能被你用一张纸复算出来、并且彼此严格自洽。本讲会用到其中这几项:主管当量总长 L = 6.0 m、控制体数 N = 20、散热器 UA = 350 W/K、空气侧质量流量 0.60 kg/s、热源功率 5.0 kW、进风温度 40 ℃;由泵曲线与系统阻力曲线联立出的稳态工作点为体积流量 18.80 L/min、管内流速 u = 0.638 m/s。完整输入清单与工作点怎么解出来,在第 4 讲。

工程量级。在 T1(⛔ 教学假设值,不对应任何平台)里,1D 给得出的是这三个平均量:散热器进口温度(也就是回路最高温)61.43 ℃、散热器出口温度 56.96 ℃、回路温升 4.47 K。而「芯体上哪一排扁管最热」「哪一片翅片先结霜」「水室里哪条流道少流」这三个问题,它一个也答不出来——⚠ 请注意措辞:不是答得偏了几度,是没有对应的变量可以拿来答

易错点。

  • 把「1D 算出来的温度偏低」当成精度问题去调 UA。真实情况往往是:你拿一个截面平均值去和一个局部峰值比。这是口径错不是参数错,调 UA 只会把一个本来对的参数调坏——这条线索直接连到第 6 讲的对不上原因分类。
  • 反过来,以为「1D 答不出局部」就等于「1D 不准」。于是在明明只需要平均量的系统级问题上也去做 3D,把开发周期从小时量级拖成周量级。
  • 把控制体数当成精度旋钮往上加,指望加够了就能看局部。它只会让沿流向更细、算得更慢,横截面上还是一个值(下一节把这笔账算给你看)。
一条自上而下的判据链:起点框写把要回答的问题写成一个响应量,向下进入判断菱形它是不是某个控制体上的平均量。判是的一支进入框 1D 是正解不是将就,其下再分两支——需要分辨空间不均则去准 2D 分层网络,不需要则用 0D 或 1D。判否的一支进入框加控制体、加时间步、压残差都无效,其下并列三个出口框:换方法去 3D 或共轭传热、换分工由 1D 出边界让 3D 看局部、换问题问平均量能不能替代。链的右侧另起一列典型响应量实例栏,左栏是平均量,列出回路流量分配、各节点平均温度、系统换热量、时间常数;右栏是非平均量,列出芯体最热扁管、翅片结霜位置、冷板内流量偏差、壁面膜温、分离与回流;两栏之间一条竖分隔线,标注右栏的问题不因为模型更细而变成左栏的问题。本图为判据流程示意,框的大小与位置不含任何数量含义,不得据图读取任何数值
图2 「1D 答不答得了这个问题」的判据闸。该读出三件事:① 判据只有一问——响应量是不是控制体平均量,⛔ 图上不出现「精度」「误差」这类词,因为本图判的不是精度;② 判「否」之后的三条出口没有一条是「把 1D 做细」,这是本图最要紧的一层;③ 右栏那五个问题在 1D 里永远无解,因为模型里没有承载它们的变量。「控制体平均量」在图上的定义与正文一字不差:某一个控制体上的截面平均压力、平均质量流量或平均温度。⚠ 判据与两栏实例均为本课归纳,⛔ 不来自任何标准条款。本图为判据流程示意,框的大小与位置不含任何数量含义,⛔ 不得据图读取任何数值。

1.3 加控制体只提高沿流向的分辨率——N 决定输运分辨率,不决定换热精度

是什么。控制体数 N 决定 Δx = L/N,而 Δx 决定两件事:温度锋面能被分辨到多细、管路的纯延迟能被表达成几个台阶。换热量准不准,则由降阶参数决定——由 UA 或 ε-NTU 关联式选得对不对、标得准不准决定,与 N 基本无关。⛔ 这两件事在 1D 里由完全不同的东西定,不得互相推导

为什么必须把它们分开。两个极端的例子就够说明问题:一个把整台散热器当成单个控制体的模型,只要 UA 标得准,稳态换热量就是准的;而一个把管路切成 200 段、UA 却直接抄了样本额定值的模型,稳态换热量照样错。前者的输运分辨率极差而换热准,后者的输运分辨率很好而换热不准——它们各自错在各自那一本账上,互相不能补救。

工程量级。在 T1(⛔ 教学假设值)里:主管当量总长 L = 6.0 m、管内流速 u = 0.638 m/s。取 N = 20 时,Δx = 0.300 m,单个控制体的输运时间 Δx/u = 0.470 s,而整条回路的输运纯延迟 L/u = 9.40 s。⚠ 这两个数的性质完全不同:Δx/u 随 N 变,L/u 不随 N 变——N 加倍则 Δx/u 减半,而 9.40 s 一动不动。⛔ 不要由此以为「加控制体让这个系统响应变快了」,物理上什么都没有变,变的只是你看它的分辨率。

把 N 从 20 加到 200,得到的是:输运分辨率提高 10 倍,温度锋面被数值抹开的尺度从 1.90 m 降到 0.600 m(同一条 CFL = 1 的口径下;这个抹平尺度是怎么算出来的、它的适用条件是什么,第 2 讲 2.6 给闭式与三条限定,第 4 讲 4.5~4.6 拿它算步长与控制体数的代价账)。而稳态的 61.43 ℃ 一点不动——因为稳态换热量由 UA 与两侧流量定,与 N 无关。⚠ 顺便记住 1.90 m 这个数的分量:它占整条回路 6.0 m 的 31.6%,也就是说 N = 20 的模型即使把步长压到 CFL = 1,锋面也已经被抹开了将近三分之一条回路。

易错点。

  • 用「加了网格结果没变」来论证「模型已经网格无关、所以准了」。这是把 3D 的网格无关性检验直接搬到 1D 的典型误植。在 1D 里,这个检验只证明了一件事——输运分辨率够了;它与降阶参数准不准毫无关系,而后者才是主要误差来源(下一节起会展开)。
  • 反过来,用「换热量对不上」去加控制体。加到多少都不会对上,因为那本账的钥匙在 UA 上不在 N 上。
  • 把 N 加大却不同步减小步长。Δx 变小而 Δt 不变,CFL 反而变大,输运分辨率不升反降——加了网格、算得更慢、结果更糊。这条等值线账在第 4 讲。
同一台液-气散热器的两种表达并排。左半是 1D 的表达:芯体被沿流向切成若干控制体,每个控制体上只画一个温度符号 T̄ᵢ 与一个流量符号 ṁ,下方画一条沿流向的一串离散平均值折线,底部标注变量是每个控制体一个值。右半是 3D 的表达:同一台芯体上画出扁管排布与一个二维温度分布的形态云带,并标出三处 1D 无法表达的对象——某一排扁管的局部高温、翅片区的结霜起始位置、水室内的流量不均,底部标注变量是一个场。两半之间一条竖线,标注同一台硬件两种自由度。右半的云带不带任何色标数值。本图为定性示意,色带无色标、折线无刻度,不得据图读取任何温度或倍数
图3 同一台散热器在 1D 与 3D 里的变量结构对照。该读出三件事:① 1D 的输出是一串离散的截面平均值,横截面方向上没有任何变量;② 右半标出的那三处对象,在左半的变量表里找不到对应项——不是值不对,是没有这一项;③ 因此「1D 算不准局部」这句话本身就说错了:它不是算得不准,是没有在算。左右两半用的是同一台假想散热器,⛔ 不对应任何实物型号。⚠ 本图不含任何计算结果:折线的起伏与云带的形态都是定性表达,⛔ 与教学算例 T1 的数值无关,也不得反过来用它去核对 T1。本图为定性示意图,色带无色标、折线无刻度,⛔ 不得据图读取任何温度或倍数。

1.4 三方分工:1D、3D 与试验之间只有两类东西在流动

是什么。1D、3D 与试验三者不是精度阶梯,是分工。把三者之间实际交换的东西列出来,只有两类

  • 降阶参数往 1D 里流。3D 给 1D 的是多孔系数、装车后的风扇曲线、分区换热系数;试验给 1D 的是 UA 与 ε、阻力系数、金属件热容、泵曲线。
  • 边界条件往 3D 与试验里流。1D 给 3D 的是各支路流量、进出口温度、要复现的工况时刻;1D 给试验的是试验矩阵与最危险工况、测点位置建议。

三方之间还有第三条链:3D 与试验互为验证。三者构成的是一个闭环,⛔ 不是一条「先 1D、再 3D、最后试验」的流水线。

为什么要把这两类分清。因为它们断掉时的症状完全不同,可以按症状反推是哪条链断了:

  • 参数链断(比如把样本额定点的 UA 直接拿来用),症状是「某一个工况凑得上、换一个工况就崩」——因为那个参数只在标定它的那个点附近成立。
  • 边界链断(比如 3D 用了一个拍出来的进口流量),症状是「算得很细但整体量级不对」。⚠ 这一条最贵:3D 那一侧看不出来,因为它内部完全自洽——网格漂亮、残差降到底、云图光滑,唯一的问题是它在回答一个不存在的工况。
  • 验证链断(3D 与试验各干各的、从不互相对),症状是「两边各自自洽,同一个错误无人发现」

工程量级。参数链的权重可以量出来。在 T1(⛔ 教学假设值)里,把四个输入各扰动 20%,回路最高温相对环境的那 21.43 K 会随之改变,改变量分别是:热源功率 +20% ⇒ +4.29 K;散热器 UA +20% ⇒ −2.41 K;空气侧风量 −20% ⇒ +1.09 K;系统阻力 K +20% ⇒ +0.18 K。⚠ 这一组数的口径是「解耦」——即扰动某一项时,不把 UA 与两侧流量之间的耦合算进去。⛔ 不得把它当成唯一的一组数:换成耦合口径,量值会明显不同,第 6 讲会把两套口径并列摆出来。这里只用它说明一件事:参数之间的权重差着一个数量级以上,因此「该标定谁」是有答案的

易错点。

  • 把 3D 当成「真值」去校 1D。在系统级瞬态这一层,3D 通常比 1D 不可信——它的边界是拍的,而且跑不起长瞬态(下一节与 1.7 各说一半)。
  • 把试验当成无误差的参照。试验侧本身有一整类不确定度:测点位置与安装、传感器精度与标定、工况控制误差、两边「稳态」的判据不同。本课只声明它存在且必须与仿真侧的取值口径对齐,具体的台架不确定度与测量方法前指 L7-01 测量技术:热电偶/红外/PIV/流量测量L2-03 冷却系统台架与流量分配试验
  • 三方各自迭代却不记录参数与边界的版本。三个月后没有人答得出「这条曲线用的是哪一版 UA、哪一套边界」。模型资产该留哪些要素才能重跑,前指 J7-07 仿真任务的输入、交付与模型资产管理
三个圆角框成三角布置,分别是 1D 系统模型、3D CFD、试验。三者之间六条有向箭头分成两族并用两种颜色区分:降阶参数族包含 3D 指向 1D 与试验指向 1D 两条,边界条件族包含 1D 指向 3D 与 1D 指向试验两条;3D 与试验之间另有一对箭头标注互为验证。每条箭头上写清运的是什么:3D 到 1D 写多孔系数、装车后风扇曲线、分区换热系数;试验到 1D 写 UA 与 ε、阻力系数、金属热容、泵曲线;1D 到 3D 写各支路流量、进出口温度、要复现的工况时刻;1D 到试验写试验矩阵与最危险工况、测点位置建议。三角形外侧各挂一个断链症状标签:参数链断是某工况凑得上换工况就崩,边界链断是算得很细但整体量级不对,验证链断是两边各自自洽无人发现同一个错。本图为关系示意,箭头粗细与框的大小不代表任何权重或数量,不得据图读取任何数值
图4 三方分工与两类流动物。该读出三件事:① 三方之间只有两类东西在流动,且方向是固定的——降阶参数只往 1D 流,边界条件只往 3D 与试验流;② 每条链断掉时症状不同,因此可以按症状反推是哪条链断了,这是本图最实用的一层;③ 3D 不是「真值」——它自己也需要边界,而边界来自 1D。图上「降阶参数」指代替一段物理的集总系数,「边界条件」指外部给定的状态量,与正文定义一致;⛔ 图上不出现任何参数的数值。⚠ 箭头内容与三条断链症状是本课归纳的,⛔ 不是标准或既有手册里现成的分类。本图为关系示意图,箭头粗细与框的大小不代表任何权重或数量,⛔ 不得据图读取任何数值。

1.5 第二根钉子:收敛的 1D 结果,差主要落在降阶参数上,不落在求解器上

是什么。上一节把「参数从哪来」摊开之后,本课的第二根钉子就是它的直接推论:一个已经收敛的 1D 结果,它的可信区间由降阶子模型的参数定,不由求解器定。

理由在 1D 的建模方式里:模型里的换热器不是被算出来的,是被 UA 或 ε-NTU 关联式代表的;管路不是被算出来的,是被阻力系数代表的;风扇不是被算出来的,是被一条曲线代表的;电芯生热不是被算出来的,是被一个功率或产热模型代表的。这些参数全部从别处借来——试验、3D、供应商数据表——而借来的东西不会比源头更准

由它推出本课的排查顺序:结果与试验对不上时,先查拓扑、再查参数、再查边界、再查口径,最后才查数值。⛔ 反过来做(先压残差、先减时间步、先加控制体)是本课最想拦住的那个动作——它把时间花在了通常最小的一本账上。这条顺序的完整展开与逐项做法在第 6 讲。

为什么数值误差通常更小——以及它的失效条件。1D 的方程规模小,隐式变步长积分器在给定容差内可以把局部误差压得很低。只要真的收敛,数值误差往往比参数误差小一个到几个数量级。⚠ 但这句话有一个明确的失效条件,必须一起记住:当步长或控制体数把某个快过程整个抹掉时,数值误差会一跃成为主导项。这不是理论上的边角情形,它恰恰是本课的反钉子所讲的那种形态,第 4 讲会把它当场算破。

工程量级。两本账各自的大小,在 T1(⛔ 教学假设值)里可以直接比:参数侧,四个输入各扰动 20% 造成的温度变化落在 0.18~4.29 K 之间(上一节那组解耦口径的数)。数值侧,同一个模型只要不踩进失效条件,误差远小于这个带;而一旦把步长取成 τ/10 这种看似合理的值(⚠ 这是一个错误动作,第 4 讲会把它当场点破),锋面被抹开的尺度就到了 7.08 m——比整条 6.0 m 的回路还长,瞬态形状整段失真。所以两本账谁大谁小,取决于你有没有踩进那个失效条件,⛔ 不能不加条件地说「数值误差不重要」。

易错点。

  • 结果不对就先加控制体、先压残差、先减时间步。三个动作都很容易做、都能立刻看到「我在干活」,而它们通常是在最小的那本账上花时间。
  • 反过来,认定「数值永远不是问题」。于是从不检查步长与离散事件,恰好落进上面那个失效条件里。
  • 把「收敛」与「正确」混为一谈。收敛只说明这组方程被解对了,它一个字都没说这组方程写对了——拓扑少一条支路、UA 抄错、边界拍成常数,残差照样降到底。

⚠ 两根钉子是配套的,但不得互相推导:第一根回答「这个问题该不该用 1D」,判据是响应量的类型;第二根回答「用了 1D,结果差在哪」,判据是参数的敏感度与来源。一个模型可以完全满足第一根(问的确实是平均量)而仍然错得离谱(参数全是拍的);也可以参数极准而答非所问。两者各由各的判据成立,谁也不是谁的推论。

1.6 第一个主场:单工况成本低到可以把「工况数」当成自由变量

前五节回答的是「该不该用 1D」与「用了差在哪」。剩下两节回到正面:判「是」之后,1D 强在哪里——而它强的原因,同样不是精度。

是什么。一次工程判断的总成本可以拆成两项:单工况成本 × 工况数 + 建模成本。3D 的单工况成本高到工况数只能取个位数,于是流程必然是先猜哪几个工况最危险、再去算;1D 的单工况成本低到工况数可以取几百上千,于是流程可以变成先扫出最危险的那个工况、再交给 3D 去看细节

为什么这是能力差异而不是效率差异。「先猜后算」的流程,其结论的上限就是猜的人的经验——清单之外的东西永远不会被算到。「先扫后看」的流程,可以让危险工况自己浮出来。真正咬人的往往是组合工况:高温怠速与快充叠加、模式切换的过渡段、长下坡回充与空调同开——这类工况通常不在任何人的清单上,因为它们不是任何单一部件的极限点,而是几件事同时发生时系统的极限点。⇒ 所以差别不在「同一件事 1D 做得快一点」,而在「有一类问题,只有工况数变成自由变量之后才问得出来」

工程量级。⛔ 1D 与 3D 的单工况机时都是平台相关量——它取决于模型规模、机器配置与许可数量,本课不给数,须在你自己团队的机器上实测一遍。⛔ 本课也不写「1D 比 3D 快 N 倍」这类倍数。本课只给结构结论:两者相差若干个数量级,而这个差距足以让「工况数」从一项约束变成一个自由变量——这才是可用的那个结论。

易错点。

  • 把低成本理解成「可以不挑工况、全扫一遍就行」。结果是产出几百条曲线而没有一条被读过。⚠ 扫描必须配一个判据:扫出来的产物是一份排序,不是一堆数据——排序才能被下一步用。
  • 把 3D 的高成本理解成「所以 3D 更准」。成本高低与误差来源没有关系,两者的误差来源根本不同(这正是上一节那根钉子)。
  • 堆工况却不记录每个工况的边界来源。最后无法追溯哪条曲线用的是哪套输入,整批结果作废(前指 J7-07 仿真任务的输入、交付与模型资产管理)。

1.7 第二个主场:系统级瞬态问的是「能量在时间上怎么在整套系统里搬」

是什么。pull-down、快充温升、模式切换、热浸后再启动——这一类问题的答案取决于三件事:哪些热容被充放、以什么速率、被谁限制。而热容与限制环节是分布在整个回路上的:冷却液、金属件、电芯、换热器芯体各有各的热容,泵、阀、风扇、压缩机各有各的限制。要回答它,就必须把全回路闭合起来,跑到分钟乃至小时量级。集总热容节点的能量方程 C·dT/dt = Q̇_in − Q̇_out 就是这类问题的基本积分单元——1D 做的事,就是把几十上百个这样的节点连成一张能量互相搬运的网。

为什么 3D 在这一层反而弱。3D 要闭合全回路,就得把所有部件都建进去;而它的时间步又被空间网格钉死,长瞬态的机时不可承受。于是 3D 的常规做法是取某一个时刻做稳态——而那恰恰把瞬态问题的全部内容丢掉了:瞬态问的就是「随时间怎么搬」,取一个时刻等于把时间维度删掉。⇒ 这是 1D 唯一在能力上强于 3D 的地方,而它强的原因是它要的是长时程与全回路闭合,不是局部分辨率

工程量级。在 T1(⛔ 教学假设值)里,回路的时间常数 τ = 126.10 s(⚠ 它的分母取的是 ε·C_min 而不是 UA,为什么必须这样取、以及取错会差多少,第 4 讲专门算一遍)。要看到稳态,得跑到 3~5τ,也就是 378~630 s。而同一个系统里最快的活跃模态——单个控制体的输运时间——只有 0.470 s。⇒ 若用显式定步长,光这一条回路就需要 8×10²~1.3×10³ 步;再把控制器采样周期算进来(10⁻²~10⁻¹ s,⚠ 这是行业典型带,⛔ 两个端点都不得当成可承诺值,具体取值须按项目的控制器规格确定,前指 K1-02 执行器驱动:电子阀/泵/风扇/压缩机 PWM),步数还要再涨一到两个数量级。同一个系统里同时活着相差两三个数量级的时间尺度——这就是第 4 讲要讲的「刚性」,也是隐式变步长积分器在 1D 里成为标配的原因。

易错点。

  • 把「瞬态」理解成「把稳态算例加一根时间轴」。加了时间轴却不去检查哪些热容被建进来了。⚠ 这是有代价的:稳态时 dT/dt = 0,热容完全不出现在方程里——所以一个漏了金属件热容的模型,稳态结果一分不差,瞬态却整段偏快。这条账第 2 讲会算给你看。
  • 用一段 3D 稳态结果去回答一个瞬态问题。典型形态是:拿最恶劣时刻的稳态温度当作峰值上报。而真实峰值受热容限制,在那个时刻根本到不了——这样报出来的是一个偏保守但也偏离事实的数,它会把设计推向不必要的余量。
  • 忘了瞬态的边界本身也在变。车速、热源功率、环温在整段瞬态里都是时间的函数,把它们拍成常数,等于用一个恒定边界去问一个变边界的问题。

1.8 第一根钉子推到验收那一层:取值断面必须在建模阶段就与测点对齐

是什么。第一根钉子还有一个直接后果,它落在验收上:仿真取的是控制体平均,试验测的是一个点。1D 里的「散热器进口温度」是某个控制体的截面平均值;试验里的「散热器进口温度」是某个位置、某一支热电偶的读数,它受插入深度、贴壁与否、传感器自身热惯量与辐射误差影响。两者在物理上不是同一个量。

为什么必须在建模阶段就对齐,而不是出报告时再解释。建完模再回头对齐,你多半会发现:测点落在一个弯头的下游、或者正好落在两个控制体的交界上、或者落在一段压根没建进模型的管路上。到那时改模型的代价,远大于建模时多切一个节点的代价。更麻烦的是,这时候人会倾向于走那条最省事的路——把仿真的取值位置挪到最接近实测的那个控制体上。⚠ 那是调参不是对齐,而且不可复现:换一个工况、换一个人,它就变了。

工程量级与做法。⛔ 传感器安装带来的偏差是平台相关量——它取决于插深、该处的流态与安装方式,本课不给数,须由实测标定,测量方法与不确定度前指 L7-01 测量技术:热电偶/红外/PIV/流量测量。本课只给做法:在拓扑图上把每一个测点画成一个节点,并注明它对应哪个控制体的哪一端。这件事零成本,而它决定了后面所有相关性数字有没有意义。相关性验证的完整流程(测点选择、误差指标定义、验收阈值)前指 J7-01 仿真—试验相关性验证(Correlation),⛔ 本课不给阈值。

易错点。

  • 报告里写「仿真 61.43 ℃、实测 59.00 ℃,偏差 2.43 K」,却不写两边取的是不是同一个断面。(⚠ 这里的「实测 59.00 ℃」是为讲清这件事而设的教学假设值,⛔ 不对应任何一次试验。)这个 2.43 K 在断面没对齐时,读者无法判断它是模型误差还是取值口径误差——⚠ 而这两者的处置方式完全相反:前者要改模型,后者一分钱不用花就能消掉。
  • 事后把仿真取值位置挪到最贴近实测的那个控制体上。如上,那是调参。
  • 只对齐温度不对齐流量。流量测点的位置同样决定了它测的是干路还是某一条支路——两边如果不是同一条路径上的流量,比出来的偏差没有含义。

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

会员专属

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

查看会员方案