J3-01 换热器芯体 CFD 与多孔介质模型

课程代码 J3-01 · 板块 J 仿真分析 CAE / J3 换热器/部件级 CFD 时长 约 3.5 小时（5 讲 + 1 次单元胞标定实操） 适合对象 CFD／仿真工程师；做前端模块、机舱与换热器性能仿真的人（必修） 前置 J2-01 CFD 基础、湍流模型与网格划分（Fluent/STAR-CCM+）；A3-01 换热器传热基础：ε-NTU 法与 LMTD 法；A2-05 空气侧压降、风量与风阻的工程处理 本讲义定位 讲师授课蓝本 / 学员自学讲义，是 J3-01 大纲的完整展开版

---

引言：残差照样掉到 1e-5、云图照样漂亮，而压降系统性地少了一半

一个纯电 SUV 的前端模块上层叠着低温散热器、冷凝器和中冷器，每个芯体里有成千上万片翅片，每片翅片上还有几十个百叶窗。要回答的问题其实很朴素——三个芯体各分到多少风、总压降多大、前一个的出风把后一个坑成什么样、格栅还能不能再关小一点。可只要把翅片一片片建进网格，这道题就不是算得慢，是根本排不上机器：本课教学假设值下的单个芯体就要 10⁷~10⁸ 量级的单元，三器层叠约 2.6×10⁸ 单元（第 1 讲 1.3 逐步代入），而这还没有算百叶窗锐边加密与棱柱层、没有算芯体之外的整个机舱，更没有算设计期要扫的十几版层叠顺序与格栅开度（这笔账第 1 讲会逐步拆开给你看）。于是整机 CFD 只剩一条可行路线：先在一个周期性单元胞上把这块芯体的空气侧特性标出来，再把它作为一组等效的体积阻力与体积热源植入整机，让求解器不去解翅片之间那些通道，只解它们合起来对流场做了什么。

⚠ 但等效并不是把问题变简单了，它换掉了出错的形态。逐翅片算不动的时候，你至少知道自己没算；等效之后，你会拿到一个残差掉到 1e-5、云图很漂亮、压降「看起来合理」的结果——它可以同时是错的，而且错得很有规律：偏在乐观侧。过芯风量偏高、层叠串扰偏轻、格栅可以关得更小，每一条结论都往好听的方向偏。本课整门要处理的就是这一类错：它不吵不闹，不以「不收敛」的形式暴露，只在整机与试验对不上时才露头，而那时最容易被归因到风扇模型或密封泄漏，于是去改一个本来没错的东西。学完这门课，你应当能做五件事：① 判断一个任务该用离散翅片，还是多孔介质／宏观换热器等效；② 从周期性单元胞 CFD 里提出空气侧的 Colburn j 因子与 Fanning f 因子；③ 把 f–Re 曲线反算成 Darcy-Forchheimer 的黏性与惯性阻力系数、按各向异性植入整机；④ 用双单元（dual-cell）模型给出芯体换热量与出口温度场；⑤ 把等效结果与焓差风洞试验对标，并判断什么时候允许修正、什么时候不允许。五讲依次落这五件事：第 1 讲定颗粒度并写清芯体黑箱要交付什么；第 2 讲在单元胞上标空气侧特性；第 3 讲把曲线反算成多孔系数并处理各向异性；第 4 讲用双单元把换热闭合起来并做守恒校核；第 5 讲对标焓差风洞并划清适用边界。其后的典型案例把五讲的次序合成一条判据链，动手做则用一段扁管加百叶窗翅片的芯体，把五条目标逐条挂到交付物上（⚠ 第五条只验到「七项口径核对表」这一层，真做对标要到 L1-01 换热器性能试验（风洞焓差法） 的台架，本课⛔ 不假装能验它）。

第一根钉子（全课下文简称「钉子①」）：芯体等效交付的不是「一个换热器」，而是「一组带口径的系数」。每一个系数背后都挂着五项口径（本课归纳）——① 速度定义（表观／迎面速度 u_face，还是最窄截面速度 u_m）② 参考面积（迎风面积 A_fr，还是最窄自由流通面积 A_c,min）③ 长度归一（压降除以哪一个厚度 L）④ 基准态（空气物性取在哪个温度、哪个压力）⑤ 有效窗口（标定到底覆盖了哪一段 Re 与迎面风速）。任何一项与求解器解释这组系数时所用的口径不一致、而中间又没有做换算，结果就是系统性偏差而不是随机误差——它不随网格加密收敛，也不随换一个更贵的湍流模型变好。⇒ 于是全课对每一个交付物只要求同一件事：数字旁边必须写着它的口径；口径写不出来的数字，本课认为它没有值，⛔ 不许交、也不许收。第 2 讲敲标定口径、第 3 讲敲换算口径、第 4 讲敲换热口径、第 5 讲敲对标口径，动手做要交的那张「标定卡」就是这根钉子的物化形式。

第二根钉子（简称「钉子②」，与钉子①配对，管的是「算出来的东西哪些能读」）：等效模型让「总量」对，⛔ 不让「场」对。多孔介质把翅片的真实阻力抹成了一个体积力源，双单元把翅片的真实换热抹成了一个体积热源——所以芯体内部的速度场、温度场、翅片壁温、局部换热系数全部是假的：它们只是一组「能给出正确体积平均量」的数值填充，不是物理解。能读的是芯体外的迎面速度分布、总风量与各芯体分风量、总压降、总换热量、出口截面在体积平均意义下的温度分布、层叠芯体之间的串扰（前一个的出风就是后一个的进风）；⛔ 不能读的是芯体内任一点的速度与温度、翅片表面温度、局部热流、翅片根部温差、芯体内的湍流量。⚠ 这里的「不能读」不是「读出来不准」，是读出来的那个东西不表示那个物理量——多孔区里每个单元都会有一个速度值和温度值，云图看起来完整无缺，再加密网格也长不出被抹掉的那层信息，这正是它骗人的地方。判据一句话：凡是要靠「解出翅片附近的流动」才能回答的问题，多孔介质都答不了，那类问题只能回到单元胞去问（第 2 讲）。⚠ 两根钉子各由各的判据成立：一根管「交出去的数有没有值」，一根管「算回来的场哪些能用」，谁也不是谁的推论。

围绕这两根钉子，本课建立三个判断力：① 按「你要交付的响应量」定颗粒度——颗粒度不是「离散／多孔／宏观」三选一，而是几何解析程度与换热闭合方式两个正交维度各取一档；算力只决定你能不能用某一档，不决定你该不该用（第 1 讲）。② 让标定与植入之间的每一次换算口径同源——单元胞里量到的量、拟合出来的系数、求解器面板收下的参数，三者必须说的是同一件事；四步换算链上每一步都有一个分叉（速度定义、除哪个厚度、窗口与残差、½ 归属与面板收 α 还是 1/α），任一步错都是确定倍数的系统偏差（第 2、3 讲）。③ 知道结果哪些能读、对不上时按什么顺序查、这条路线在哪里出界——先核口径与几何，再核边界与换热闭合，系数是最后才该动的东西（第 4、5 讲）。

第二条最容易被读反，本课要当场算破的就是这句话：「速度就是速度——单元胞里量到的压降–速度曲线，拟合出来的系数直接填进求解器就行。」它之所以必然被读反，诱因就藏在两条公式自己的形状里：Fanning 压降式 Δp＝f·(4L/D_h)·(½ρu_m²) 与多孔源项 −∂p/∂xᵢ＝(μ/αᵢ)·uᵢ＋C₂ᵢ·½ρ|u|uᵢ，两条都长着「½ρ 乘速度平方」的样子、都用一个叫 u 的字母，看起来是同一个量的两种写法。可它们不是同一个速度：紧凑式换热器一系的 f 与 j 按惯例以最窄自由流通截面速度 u_m（质量流速 G＝ρu_m）为参考，而多孔源项在主流求解器的常用表述里作用在表观速度上，也就是迎面速度 u_face＝体积流量÷迎风面积；二者相差 1/σ 倍，σ＝最窄自由流通截面比，其行业典型带为 0.4~0.7（标「典型」，承 J2-01 CFD 基础、湍流模型与网格划分（Fluent/STAR-CCM+） 同一条带，⛔ 本课不另立一套），对应的动压换算倍数 1/σ² 落在 2.04~6.25。读者会做错的那个具体动作是：把单元胞扫出来的 Δp 对 u_m 做二次拟合，再把拟合出来的两个系数直接填进求解器的多孔介质面板；这样得到的黏性阻力只有正确值的 σ 倍、惯性阻力只有 σ² 倍。在本课的教学算例里（σ＝0.60，教学假设值，⛔ 不对应任何芯体、⛔ 不得取用），同一芯体在迎面 2~8 m/s 各点的压降被低估 51%~59%；σ 落在典型带内时惯性项被压到 0.16~0.49 倍，最坏情形只剩六分之一。

⚠ 它比一般的错更危险，因为这个偏差在标定点上同样存在——「我做过单元胞标定了」这句话保护不了你，错的不是数据，是数据与系数之间那一次解释；而它也不会被残差、被网格无关性研究、被压降收敛曲线抓到。⇒ 正确的判据不是「哪个速度才对」，而是：拟合时用的速度定义，必须与求解器解释这组系数时用的速度定义是同一个。合法做法有两条，选哪条都行、⛔ 不许混：(a) 单元胞后处理直接按迎面速度整理 Δp–v 曲线再拟合，什么都不用换算；(b) 按 u_m 整理（这样才能与 j／f 图册对上），拟合后再按 σ、σ² 换算成表观口径。⚠ 反向同样要说清：如果单元胞本来就是按迎面速度扫的，那就不需要任何换算，硬乘一个 σ 反而把对的改错——这条反钉子点破的是「不看口径就换算」与「不看口径就不换算」两个方向的同一个病。第 3 讲会把两个方向逐点算破。

本课的定义域是面向前端模块与机舱、来流大体均匀的空气侧芯体等效：颗粒度选型的判据、单元胞的建域与后处理口径、j 与 f 的提取、Darcy-Forchheimer 系数的反算与各向异性赋值、双单元换热闭合与守恒校核、等效结果与焓差试验的口径对齐。⚠ 下面这些事不在本课射程内，本课只写去向、⛔ 不给限值、不给等级号、不给条款内容（本课归纳的射程声明）：翅片强化机理与通用关联式 → A3-02 翅片管/扁管/微通道换热形式与传热强化；ε-NTU／LMTD 与 UA 热阻串联 → A3-01 换热器传热基础：ε-NTU 法与 LMTD 法；紧凑度与换热-压降权衡 → A3-03 换热器性能评价：换热量、压降与紧凑度权衡；湍流模型选型、壁面处理与 y⁺、网格无关性与网格收敛指标 → J2-01 CFD 基础、湍流模型与网格划分（Fluent/STAR-CCM+）；整机计算域、风扇建模与热风回流 → J2-02 机舱热管理与前端模块 CFD；空气侧阻力链与整车进风量 → A2-05 空气侧压降、风量与风阻的工程处理；焓差风洞台架、喷嘴标定与两侧热平衡判据 → L1-01 换热器性能试验（风洞焓差法）；空调箱内的斜射流与极短入口段（本课的标定路线在那里会失真）→ J4-04 HVAC 空调箱总成级 CFD：箱内建模、开度扫描与内漏标定；两相与相变芯体 → J3-04 两相流、沸腾与冷凝仿真；液冷流道内的流量均匀性 → J3-02 电池包液冷流道 CFD 与流量均匀性；湿空气与结露侧的物性修正 → A2-04 湿空气物理与焓湿图：空调负荷计算基石；辅助侧温升由一维系统模型给出的那条链 → J1-01 1D 热管理仿真导论（KULI/AMESim/GT/Flowmaster）。

⚠ 由此有三种读法会落空。一是把本课当标准条文课读——正文出现的 ASME V&V 20、ASHRAE 33、AHRI 410、AIAA G-077 四份文件，本课均未取到原文，因此只写标准号与「它管什么、要用时去哪查」，⛔ 一个条款、一个限值、一个判据数都不写；版次与实施年份以现行目录为准，引用前须核。二是来找可以直接抄进输入面板的数——本课凡带代入值的算例，输入与结果一律是教学假设值（含芯体几何、σ，以及为让全课算例可逐点复算而自设的教学关联式），不对应任何芯体、任何平台，⛔ 禁止取用、⛔ 禁止外插；孔隙率、面积密度、三个方向的阻力比、可接受的热平衡失衡幅度与对标偏差这类平台相关量，本课不给取用值，只给「须由本项目量取、标定或实测确定」这条纪律与相应的自检动作。三是把「等效」当成「精度低一档的离散」来读——它不是精度问题，是能回答的问题不同，这正是钉子②要钉住的东西。本课能给你的不是数，是推出那个数的那条链，以及判断这条链走没走对的判据。

本课的三处符号约定（全课统一，⛔ 不许中途换写法） · α 一律指渗透率（Darcy-Forchheimer 黏性阻力项里的那个 α，单位 m²）——⛔ 本课不用 α 表示热扩散率。 · ε 一律指换热效能（承 A3-01 换热器传热基础：ε-NTU 法与 LMTD 法 的定义）；孔隙率一律写 γ——文献里孔隙率常写 ε，本课把这个字母让给效能，⛔ 正文不出现表示孔隙率的裸 ε。 · u 必须带下标：u_face＝表观（迎面）速度＝体积流量÷迎风面积，u_m＝最窄自由流通截面速度＝质量流量÷(密度×最窄自由流通面积)。⛔ 正文禁止裸写不带下标的 u——这门课一半的坑就长在这个下标上。同理，全课的 f 一律指 Fanning f_F，与管路水力学的 Darcy f_D 之间是定义式层面的恒等关系 f_D＝4·f_F，⛔ 二者不可互算、不可比大小。

第 1 讲　颗粒度不是一条轴：先写下要交什么，再定两维各取哪一档

这一讲不讲多孔介质怎么填参数，那是第 3 讲的事。它要先把两件更靠前的事定下来：这个芯体在你这次计算里被建成什么，以及它作为一个黑箱，要向下游交出哪些东西。这两件事定错了，后面四讲每一步都是在一个错的底座上做精细化——单元胞标定得再仔细、多孔系数拟合得再漂亮，也补不回来。

这一讲要建立三个判断：① 建模颗粒度是两个正交维度（几何解析程度 × 换热闭合方式）的组合谱，选型的问题从来不是「用哪一档」，而是「这两维各取哪一档」；② 定档的依据是你要交付的响应量，算力只决定你能不能用某一档，不决定你该不该用；③ 芯体黑箱要交的是五组特性，而且每一组都必须附着口径——全课的两根钉子都在这一讲第一次立起来。

1.1 建模颗粒度是两个正交维度的组合谱，不是「离散／多孔／宏观」三选一

是什么。把一个芯体放进 CFD 模型里，你实际上要做两个互相独立的选择，而不是一个。

维度 A：几何解析程度——芯体的几何被解到什么程度，自细到粗六档：

• A1 全离散翅片：每一片翅片、每一个百叶窗都建出来；
• A2 周期性单元胞：只建一个代表性周期，其余方向用周期边界代替（第 2 讲整讲讲它）；
• A3 分区多孔：把芯体切成若干块，每块给自己的一组系数，用来表达芯体内部的不均；
• A4 单块均质多孔：整个芯体是一块均质多孔体，一组系数管到底；
• A5 零厚度阻力面（porous jump／压降跃变面）：厚度为零，只在一个面上加一个 Δp(v) 关系；
• A6 完全不建几何：芯体在一维系统模型里只是一个部件（→ J1-01 1D 热管理仿真导论（KULI/AMESim/GT/Flowmaster））。

维度 B：换热闭合方式——热是怎么交换的，五档：

• B1 共轭传热真解固体翅片；
• B2 双单元（dual-cell）：流体胞与辅助胞占据同一块体积，逐单元按当地温差闭合（第 4 讲）；
• B3 宏观换热器（macro HX）：按 macro 分区给 UA 或 ε，逐 macro 对闭合；
• B4 给定壁温或给定热流边界；
• B5 绝热：完全不换热，只算流阻。

实际工程里用的是这两维的组合：A1×B1 是单元胞标定用的（第 2 讲）；A4×B2 与 A4×B3 是整机的两个主力（第 3、4 讲）；A5×B5 是只关心风量分配时的快扫；A3×B2 是迎面严重不均时的分区做法（A3 这一档的系数怎么给在第 3 讲，迎面不均本身的代价在第 5 讲 5.3）。

为什么必须写成两维。把它压扁成一维的「离散／多孔／宏观」三档，会同时丢掉三样东西：

• 会让人以为「选了多孔介质就得自己给换热系数」——而 A4×B3（动量用均质多孔、换热交给宏观分区）恰恰是最常见的组合之一；
• 会让人看不见「只算风量分配时可以整层不换热」这条最省的路（B5），也看不见 A5 零厚度阻力面——而这两个正是机舱 CFD 做风量分配快扫时最常用的；
• 会让「同一次计算里动量取一档、换热取另一档」这种完全正常的做法，显得像是在混用两套模型。

★ 这里值得停一下，因为它是本课后面反复出现的一个形态：沿着一条最熟悉的线索往下枚举，漏掉的不是一两条，是整类。只沿「几何解析程度」这条线数，B5 与 A5 会成建制地不出现，而你数出来的清单看起来是完整的——它不会以「少了一条」的形式提醒你。自检只有两问：我是沿着什么线索数出来的？有没有不在这条线索上、但同样属于这个集合的东西？答不出第二问，说明全集还没建起来。

工程量级。两维必须分开选，压力来自算力这一侧：本课教学假设值下的一块迎风 0.6 m × 0.4 m、厚 0.020 m 的芯体，A1 全离散约 8.6×10⁷ 单元，A2 周期性单元胞约 4×10³ 单元，两者相差四到五个数量级（逐步代入与算式见 1.3）。⚠ 这是教学假设值下的量级估算，⛔ 不对应任何平台的芯体、不得取用。

易错点。

• 把颗粒度当成精度等级，认为「越离散越准」。A1 若热边界条件用错，会比一个口径对的 A4 错得更远。颗粒度决定的是能回答什么问题，不是答得多准。
• 在同一次计算里两维取了互不兼容的档。典型形态：几何取 A5 零厚度面，却要求它给出出口温度分布——零厚度面里没有体积，也就没有地方放体积热源（1.7 把这条讲透）。
• 忘了还有 A3 分区多孔这一档。于是迎面严重不均时只能在「整块均质」与「全离散」之间二选一，而这两条路一条答不了、一条算不动。

1.2 定档的判据是「你要交的响应量需要解到哪一层」，算力只决定你能不能用

是什么。动手选档之前先做一件事：把这次计算要交出去的响应量逐条写下来。常见的有——总风量与各芯体分风量 / 芯体总压降 / 总换热量与出风温度 / 出口断面的温度分布 / 芯体内的翅片壁温 / 局部结霜风险 / 气动噪声。写完之后回到图1 的矩阵，在右栏「这一档能交什么」里找到能交出这些量的最低档，两维各定一档，最后检查这对组合有没有被禁止符号标住。

为什么是这个顺序。这与 J2-01 CFD 基础、湍流模型与网格划分（Fluent/STAR-CCM+） 贯穿全课的那条主线——先判这个响应量由什么支配，再决定动哪个旋钮——是同一条思路落在颗粒度这一层上；本课不重讲那条判据，只把它接到颗粒度选型上（本课归纳）。反过来做（先定颗粒度、再想要交什么）在工程上极其常见，代价是算完才发现最关心的那个量这套模型根本给不出——而那时网格、边界、算例矩阵都已经做完了，重来一遍的不是一次计算，是一整轮。

工程量级。总风量、总压降、总换热量这三个量在 A4×B2 或 A4×B3 上就能交；出口断面的温度分布要 B2 或 B3 才有（B4、B5 给不了）；翅片壁温、局部热流、结霜起始位置必须回到 A1×B1——多孔介质在这三个量上给出的数没有物理含义（1.6 讲这一条）。⛔ 本课不给各档之间的算力比：它取决于网格策略与所用机器，属平台相关量，须按本项目自己估。

易错点。

• 先定颗粒度，再想要交什么。这条错误的代价不在计算本身，在日程——发现得越晚，能改的越少。
• 把「云图上有温度」当成「温度可读」。多孔区里每一个单元都有一个温度值，云图完整无缺，但那是体积平均意义下的填充值，不是翅片附近的真实温度。
• 为了「保险」一律往细里选。前端模块在设计期要扫十几版层叠顺序与格栅开度，每版都重划一次全离散网格的结果是方案扫不完——扫不完的方案本身就是一种不保险。

1.3 「整机算不动」不是判据，把网格量级的账拆完才是

是什么。「逐翅片建模算不动」这句话人人会说，但它当不了判据——判据要能换一组参数重算。把账拆开，用本课的一组教学假设值逐步代入（⚠ 以下几何全部是教学假设值，⛔ 不对应任何芯体、不得取用）：

• 迎风面 0.6 m × 0.4 m，芯体厚 0.020 m；翅距 F_p ≈ 1.4 mm，翅高 F_h ≈ 8 mm；
• 通道数 ＝ (0.6/0.0014) × (0.4/0.008) ＝ 428 × 50 ＝ 21 400；
• 每通道单元数：断面布 10×10、沿厚度布 40 层 ⇒ 10×10×40 ＝ 4×10³；
• 单芯体 ＝ 21 400 × 4 000 ≈ 8.6×10⁷ 单元；
• 三芯体层叠（前端模块常态）≈ 2.6×10⁸ 单元。

而这还没有算百叶窗锐边的局部加密、壁面棱柱层，以及芯体之外的整个机舱。对照之下，一个周期性单元胞约 4×10³ 单元——与单芯体全离散相差约 2×10⁴ 倍，与三芯体层叠相差约 6×10⁴ 倍。这就是「先标定、再等效」这条路线成立的全部理由。

为什么要自己算一遍。这串数不是拿来记的，是拿来换参数重算的：翅距减半，通道数翻倍、单元数跟着翻倍；换一块 400 mm 宽的芯体，代进去就有答案。记住「8.6×10⁷」这个数没有用，记住「(迎风尺寸/翅距) × (迎风尺寸/翅高) × 每通道单元数」这条算式才有用——以后遇到别的芯体，你要能自己判它落在哪个量级。

工程量级。数量级是 10⁷~10⁸（单芯体）与 10⁸（三芯体层叠）。⛔ 具体值随几何与网格策略变，本课给的是教学假设值下的算例，不对应任何平台；真实项目的网格量须按本项目的几何与网格策略自行估算。

易错点。

• 把这个数当成「某某芯体的网格量」去引用。它是一组教学假设值算出来的，换一组几何就是另一个数。
• 只算通道数、不算每通道的单元数。少乘那一步，量级差三个数量级，得出的结论会是「才两万个单元，能算」。
• 以为「买更大的机器就行」。真正的墙不是内存：前端模块在设计期要扫十几版层叠顺序与格栅开度，每改一版几何就要重划一次这种网格——挡住你的是人和日历，不是核数。

1.4 芯体黑箱要交的是五组特性，不是「阻力＋换热」两组

是什么。把芯体做成黑箱，要交给下游的东西一共五组：

• ① 动量特性：Δp–v 曲线或 f–Re 曲线，三个主方向各自一组（第 3 讲）。
• ② 换热特性：j–Re 曲线，或 Q–ṁ_air–ṁ_aux 的性能面，或 UA–ṁ 关系（第 2、4 讲）。
• ③ 几何与拓扑特性：孔隙率 γ、面积密度 a_v、水力直径 D_h、最窄自由流通截面比 σ、芯体厚度 L 与迎风面尺寸、主流向在整机坐标系下的方向余弦。
• ④ 热物性与热容：固体材料、芯体金属质量与比热，以及管内存液的那一份。
• ⑤ 有效性边界：标定覆盖的 Re 窗口与迎面风速窗口、物性基准态、干态还是湿态（表冷器结不结露）、是否含结霜、来流均匀性假设。

⚠ 符号纪律（全课统一，从这里开始）：孔隙率一律写 γ，⛔ 正文中不裸写 ε 表示孔隙率——本课的 ε 只表示换热效能；α 一律指渗透率，⛔ 不指热扩散率。另外，γ（体积比）与 σ（最窄断面的面积比）是两个不同的量，⛔ 不得互相顶替，第 3 讲会把它们分开讲透。

为什么是五组。写成「阻力＋换热」两组是最常见的交接遗漏，而漏掉的三组各自会在一个不同的下游步骤上出事：

• 缺 ③：双单元换热闭合要的 a_v 没人给；速度口径换算要的 σ 没人给；整机模型里那块多孔体的厚度 L 与迎风面积无从与标定同源（第 3 讲会把「同源」这件事量化）。
• 缺 ④：热浸、冷启动这类瞬态算出来的芯体响应快得离谱，因为求解器眼里那块体积就是一团空气（1.8 讲这一条）。
• 缺 ⑤：系数被外插到标定窗口之外，而外插不报错——没有任何一条残差曲线会告诉你「这个风速点不在标定范围里」。

★ 又是同一个形态：沿「阻力＋换热」这条最熟的线索数下去，漏掉的三组不会以「少一条」的形式出现，它们是整类不出现，清单看起来是齐的。

工程量级。五组里通常只有前两组会出现在供应商样本上，后三组多半要自己量、自己算。⇒ 交接清单必须显式列出五组，而不是等对方给什么就用什么（本课归纳）。

易错点。

• 拿到 Δp–v 与 Q–ṁ 两张表就开工。等到要算热浸才发现没有芯体金属质量，而那时供应商已经进入下一轮报价。
• 把 A3-03 换热器性能评价：换热量、压降与紧凑度权衡 的气侧面积密度量级线当成某芯体的 a_v 填进求解器。「气侧面积密度达 10³ m²/m³ 量级即属紧凑式」（典型）是判「紧不紧凑」的量级线，⛔ 不是取用值；双单元里要填的 a_v 须按本芯体的换热面积与体积算，属平台相关量。
• ⑤ 只写了 Re 窗口，没写干态还是湿态。于是把干态标定的 j 用到会结露的表冷器上，少算的是潜热那一大块（去向：J3-04 两相流、沸腾与冷凝仿真、L1-01 换热器性能试验（风洞焓差法））。

1.5 ★钉子①：交付的不是「一个换热器」，是一组带口径的系数

是什么。五组特性里的每一组数字，背后都挂着五项口径：

• ① 速度定义：这条曲线里的 u 是表观／迎面速度 u_face，还是最窄截面速度 u_m？
• ② 参考面积：迎风面积 A_fr，还是最窄自由流通面积 A_c,min？
• ③ 长度归一：压降除的是哪个厚度 L——芯部的，还是含边框的？
• ④ 基准态：空气物性取在哪个温度、哪个压力？（芯体前后有温升，取进口态还是平均态，本身就是口径的一部分。）
• ⑤ 有效窗口：标定覆盖了哪一段 Re、哪一段迎面风速？

把这五项做成一张固定格式的标定卡，每一组系数必须附卡交付——这就是本课「动手做」最后要交的那个东西。⛔ 口径写不出来的数字，本课认为它没有值：不许交，也不许收。

为什么这是全课第一根钉子。因为口径不一致造成的不是随机误差，是系统性偏差；而它不会以「不收敛」的形式暴露——残差照样掉到 1e-5，云图照样漂亮，压降照样是个「看起来合理」的数。本课后面每一次出错，都能回到这五项里的某一项：第 2 讲的标定口径（①②④⑤）、第 3 讲 α 与 C₂ 的换算（①③）、第 4 讲把空气侧当成了总的（参考对象错）、第 5 讲与试验对标（②④⑤）。⇒ 全课每一个交付物都要求同一件事：数字旁边必须写着它的口径。

工程量级。口径错的量级由口径本身决定，而不是由网格或湍流模型决定。本课教学算例里，只把速度口径这一项弄错（拿按 u_m 整理的数据拟合出系数，直接填进按表观速度解释系数的求解器），迎面 2/4/6/8 m/s 处的压降就分别被低估 51%/55%/57%/59%（教学假设值 σ ＝ 0.60 下的算例，逐点复算见第 3 讲；分母一律取正确口径下的压降值）。⚠ 这些百分比不是任何芯体的性质，⛔ 不得当作可乘的修正系数。它的量级远大于换一个湍流模型带来的差别——⇒ 先把口径核对了，再谈模型精细化。

易错点。

• 把口径写在邮件正文里，而不是写在数据文件里。转手一次就丢，而系数还在继续往下传。
• 只写「按 SI 单位」。SI 不解决「u 是哪个 u、A 是哪个 A」——两个速度都是 m/s，两个面积都是 m²，单位一模一样。
• 收到别人的系数时不问口径就直接填，出了问题再去追责。等到整机与试验对不上，几何已经改了三版，追不回来了。

1.6 ★钉子②：等效模型让「总量」对，⛔ 不让「场」对

是什么。多孔介质把翅片的真实阻力抹成了一个体积力源，双单元把翅片的真实换热抹成了一个体积热源。省下来的算力正是从这里省的，而被抹掉的，正是芯体内部那一层信息。所以两张清单必须在交付时就写死：

可读的：芯体外的迎面速度分布 · 总风量与各芯体分风量 · 芯体总压降 · 总换热量 · 出口截面在体积平均意义下的温度分布 · 层叠芯体之间的串扰（前一个的出风就是后一个的进风）。

⛔ 不可读的：芯体内部任一点的速度与温度 · 翅片表面温度 · 局部热流 · 翅片根部温差 · 芯体内的湍流量。

判据可以压成一句话：凡是要靠「解出翅片附近的流动」才能回答的问题，多孔介质都答不了；那类问题只能回到单元胞去问（第 2 讲）。

为什么必须写死，不能让读者自己悟。因为多孔区里每一个单元都有一个速度值和一个温度值，云图看起来完整无缺、颜色连续、没有任何异常。⇒ 「不可读」不是「读出来不准」，是读出来的那个数不表示那个物理量——它只是一组「能给出正确体积平均量」的数值填充。这两件事的处置完全不同：前者可以靠加密网格、换模型去改善；后者加密到天荒地老也长不出来，因为那一层信息在模型里根本不存在。

工程量级。⛔ 本课不给「多孔区内部温度偏离真实翅片温度多少」这类数——它不是一个偏差，是一个没有定义的比较。这一条给的是判据本身，不是一个数。

易错点。

• 拿多孔区内部的温度云图去判断「哪片翅片会结霜」「哪一排管子先沸腾」。
• 看到多孔区里速度场很平滑，就断定「芯体内流动很均匀」。那是体积源抹出来的平滑，不是物理。
• 把多孔区网格加密，指望「看清芯体内部」。加密改变不了体积源的闭合关系——第 3 讲末尾（3.6）会把这条的机理讲透（⛔ 那里同样不给「几层就够」的层数，属平台相关量）。

1.7 零厚度阻力面与有厚度多孔块不是精粗之分，是能回答的问题不同

是什么。A5 零厚度阻力面（porous jump／压降跃变面）只在一个面上施加一个 Δp(v) 关系，在网格上不占厚度；A4 有厚度多孔块占据真实芯体的那块体积，能承载体积动量源与体积热源。

为什么不是精粗之分。两者的差别是有没有体积，而体积决定三样东西能不能放进去：体积热源（换热）、固体热容（瞬态）、以及沿厚度方向的速度变化。所以——只关心风量分配与总压降时，零厚度面省掉整块芯体的网格、还不引入多孔源项带来的数值刚性，是最经济的做法；一旦要算换热、要出风温度、要算瞬态，就必须有体积，没有第二条路。

工程量级。⛔ 本课不给两者的算力比——它取决于网格策略，属平台相关量。判据是一句话：你要交的量里，有没有需要「体积」才能承载的量；有，就必须用多孔块。

易错点。

• 用零厚度面算完风量，再想在同一套网格上加换热。发现无处可加，只好重划网格重算——省下的那点时间全赔进去。
• 在零厚度面上按芯体厚度做归一化。它没有厚度，Δp 是直接给的，⛔ 不除以 L。
• 把两种做法的系数混着填。多孔块的系数是 Pa/m 量纲的压降梯度式，零厚度面收的是 Pa 量纲的压降式；量纲不同，混填的后果是量级错，不是精度差。

1.8 稳态可以不给芯体热容，但省了就不能拿这套模型去算瞬态

是什么。芯体的金属质量 × 比热（再加上管内存液的那一份）是它的热惯量。稳态解与热惯量无关——稳态方程里根本没有这一项；瞬态则与它强相关：热浸、冷启动、快充后停车、风扇启停，全在它的射程里。

为什么这条线要画清楚。因为「省」与「错」在这里只隔一步。多孔区若只定义了流体、没有固体域，求解器眼里那块体积就是一团空气；空气的热容比真实芯体（铝 ＋ 管内液体）小得多，于是算出来的温度响应会快得离谱——升温快、降温也快，热浸峰值被低估。⚠ 方向必须写死：偏差落在乐观侧，它看起来像好消息，因此不会有人去查它。

工程量级。⛔ 本课不给热容比的具体倍数——它取决于芯体材料、壁厚、管内存液量，属平台相关量，须按实物称重与几何统计确定。给的是一条判据：凡是时间尺度进入分钟量级的工况，动手之前先问一句「芯体热容进模型没有」（本课归纳）。

易错点。

• 拿稳态标定好的多孔区直接跑瞬态，从不检查固体域是否存在。稳态那一版跑得很对，这正是它能一路传下去的原因。
• 给了固体热容，却忘了管内液体的那一份。
• 把热浸算出的「响应很快」当成好消息，拿去支撑一个降配决定。它是模型缺件，不是散热好。

本讲小结：这一讲交出五件可以直接写进建模方案的东西。①颗粒度是两个正交维度（几何解析程度六档 × 换热闭合五档）的组合谱，工程上真正在用的是组合——A4×B3、A5×B5 这些都在被压扁成一维之后整类消失；②定档的判据是你要交的响应量，先写清单再查矩阵，算力只决定能不能用、不决定该不该用；③「整机算不动」要拆成网格量级的账才当得了判据，教学假设值下单芯体全离散约 8.6×10⁷ 单元、周期性单元胞约 4×10³ 单元，相差约 2×10⁴ 倍——记算式，不要记数；④芯体黑箱要交五组特性（动量／换热／几何拓扑／热容／有效性边界），后三组通常没人给你，缺哪一组会在哪一步以什么形式出事，要一一对得上；⑤两根钉子在这里立起来——交付的是一组带五项口径的系数（口径写不出来的数字没有值），而等效模型只让总量对、⛔ 不让场对（要靠翅片附近流动才能回答的问题，只能回到单元胞去问）。第 2 讲就去建那个单元胞，把 j 与 f 两条曲线连同它们的口径一起提出来。