TMS BOOK · ACADEMY 讲义

一维/多维、稳态与瞬态导热的工程建模

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

A1-01 一维/多维、稳态与瞬态导热的工程建模

课程代码 A1-01 · 板块 A 基础理论与共性技术 / 传热与热力学基础 时长 约 3.0 小时(7 讲) 前置 体系外前置(需自备):高等数学(偏微分方程基础);体系外前置(需自备):大学物理传热概念 本讲义定位 讲师授课蓝本 / 学员自学讲义,是 A1-01 大纲的完整展开版


按任务导航

你手上的事 直接看 还要配
审一份别人交来的导热仿真报告 4.6 离散误差与模型误差是两本账,网格无关性只结清前一本 3.6 引用别人的 Bi 之前必须问出的两件事;7.6 交付要连着说清的三件事
定一条「这个件该不该上多维」的判据 2.1 判据是横向与纵向热阻之比,不是「件薄不薄」 4.1 必须上多维的三类触发,每条都指得出具体位置
选稳态还是瞬态求解 5.3 用 t_w 与 t_sys = max(τ, t_diff) 的比值分三区 5.1 τ 由谁主导;5.2 Fo 回答的是另一个问题
写建模输入规格 / 向供应商提模型要求 3.2–3.3 特征长度有两套口径,规格里必须写明用哪一套 7.4 接触热阻怎么进模型,并按上下限报区间,不报单值
把一块现成冷板缩放到新平台 6.2 两个时间尺度标度不同,缩放后「谁主导」会翻转 6.3 h 不随件等比缩放;6.5 相似准则不完备,不许外推

本课不覆盖:

【栏一 · 不讲哪些题目】h 从哪来(关联式、流态),见 A1-02 对流换热:自然对流、强制对流与换热系数估算;辐射那一支,见 A1-03 辐射换热与整车热害中的辐射屏蔽;k、ρ、c 的取值与温度相关性,见 A1-07 热管理常用工质的热物性与查用、H1-01 热管理用铝合金牌号、性能与选材;把判据接成热网络、标定 R 与 C、降阶与实时,见 A1-06 热阻网络法与集总参数建模入门(本课只交付判据,不搭网络);真实 CFD 的网格生成、湍流与共轭传热,见 J2-01 CFD 基础、湍流模型与网格划分(Fluent/STAR-CCM+)、J2-06 整车/机舱 CFD 前处理:从脏数模到可算网格(数模验收·删简·包面·计算域抽取);仿真与试验的相关性判定与不确定度,见 J7-01 仿真—试验相关性验证(Correlation)、J7-06 仿真不确定度量化(UQ)与稳健/可靠性设计:从确定性单点到带分布的裕度;TIM 选型与压力关系,见 D2-05 冷板与电芯的接触热阻与 TIM 应用、H5-01 导热界面材料(导热垫/凝胶/灌封)选型;热失控与快充的具体工况判据,见 D1-03 热失控触发机理与特征参数、D4-02 超充(4C/6C)热管理挑战与冷却强化。

【栏二 · 真正不给的数】① ASTM E1225 与 ASTM D5470 的条款内容、试样与限值——本项目未独立核对原文,本课只在 7.4 写这两本各测什么量、测出来的数落在模型的哪个输入上;引用前按现行版本核原文。导热系数与界面材料测试还有对应的国标,本项目对其中任何一本都没有做过适用范围判定(「这本管的是不是这件事」),所以本课不点名、不写内容;拿到它的路是:按你所在项目的引用标准清单确认,或按现行标准目录检索后核原文。② 本课的 1%(网格无关性)与「5 倍」(稳态分界)是本课自定的教学口径,不是验收阈值——真实项目的收敛判据与不确定度量化,见 J2-01 CFD 基础、湍流模型与网格划分(Fluent/STAR-CCM+)、J7-06 仿真不确定度量化(UQ)与稳健/可靠性设计:从确定性单点到带分布的裕度;本项目自己的验收阈值不在本课。③ 两处母排算例(7.5 与「典型案例详解」,它们用的是两组不同的输入)的截面、长度、电流、h 与环境温度都是教学假设值,不对应任何平台、不得取用;具体牌号的物性与电芯分方向导热系数,见供应商数据表、H1-01 热管理用铝合金牌号、性能与选材、D2-06 圆柱/方形/软包电池的冷却结构差异。

引言:先算三个比值,再决定建什么模型

教科书的导热章节回答的是「给定方程与边界条件,怎么把温度场解出来」。而工程上,在你打开前处理软件之前,还有一个更早、也更贵的问题没有答案:这个件,我该建一个多复杂的模型?

建复杂了:几何一比一、网格上百万、跑一晚上还未必收敛,而评审日期不会因此推迟。建简单了:三分钟出一个温升数,可它恰恰在最该看的那几个位置——接头、拐角、贴合面——偏得最多,而你不知道它偏了,也说不出偏多少。

这门课交付的不是解法,是动手之前就能算出来的那条复杂度下限:这个件该按几维算、该稳态还是瞬态、哪些项不进模型就一定错。

第一条结论:该几维、该稳态还是瞬态,不是经验问题,是三个比值的问题

这三个比值都只要一次乘除就能估出来:

  1. 横向热阻与沿主热流方向总热阻之比——决定能不能把这个件降到一维;
  2. Bi = h·L_c/k——决定能不能把一整块当成一个温度;
  3. 工况特征时间 t_w 与系统时间尺度 t_sys = max(τ, t_diff) 之比——决定稳态结论够不够用。

由此落地成一条动作:建模之前先把这三个比值估出来,再决定建什么模型。不许先把模型建好、再回头解释它为什么够用——顺序一反,判据就退化成事后合理化:它永远不会否掉你已经做完的那件事。

三条判据线的性质要分开记,别一律当成「经验值」:Bi < 0.1 是工程惯例线,不是物理常数、也不是任何标准的条款;本课用来分稳态与瞬态的 t_w > 5·t_sys 里那个「5 倍」是本课自定的教学分界,不是标准也不是行业惯例,正文会说明越过它会错在哪一侧。

第二条结论:错得最狠的从来不是方程,是方程旁边那两个几何量

每一条判据都自带一个特征长度和一个面积,而这两个量本身没有唯一取法:

  • 长度取错:Bi 差 2~3 倍、Fo 差 4~9 倍(长度在 Fo 里是平方),判据当场翻转。本课的默认取法是 L_c = V/A_conv(体积 ÷ 参与对流的散热面积);而一项近似 / Heisler 图解那一套用的是另一个 L_H(平板取半厚、长圆柱与球取半径 R)。两套在长圆柱上差 2 倍、球上差 3 倍,在两面散热的平板上恰好相等——正因为「有时候相等」,混用才不会当场报错。
  • 面积取错:导热截面 A_c、参与对流的散热面 A_conv、接触界面 A_contact 混用一次,整段热阻的量级就错了。

由此立一条全课纪律:凡写下一个 Bi、Fo、τ 或一段热阻,必须同时写出它用的是哪一个长度、哪一个面积。别人递给你一个 Bi 时,先把这两件问清楚;问不出口径,就把它当成一个倍数未知的数,不许直接搬进自己的模型。

常见的读反:「α 大,所以导热好,所以稳态温升低」——按它选材,方向会整个做反

这句话把两件互不相干的事接成了一条因果链。真相是两条分工:稳态温度场只认 k——稳态导热方程 ∇·(k∇T) + q̇ = 0 里根本没有 ρc,因此也没有 α;α = k/(ρc) 只管瞬态,它决定温度扰动传播多快,而不决定最终温度多高。

一个当场可验算的反例。常温常压空气的热扩散率 α ≈ 2.2×10⁻⁵ m²/s(300 K、101.325 kPa,由物性库 CoolProp 算出);304 不锈钢的 α ≈ 4.1×10⁻⁶ m²/s(由 α = k/(ρc) 现算,输入 k ≈ 16.2 W/(m·K)、ρ ≈ 7900 kg/m³、c ≈ 500 J/(kg·K),来源等级=材料手册,不挂标准号)。空气的 α 约为不锈钢的 5 倍。而两者的 k 是 0.026 与 16.2 W/(m·K),差约 600 倍、方向相反。按 α 排序去挑「导热好的材料」,会把空气排在不锈钢前面——而没有人用空气去导热。

读反之后你会做的那个具体动作:拿两份材料数据表比 α,挑出 α 大的那一个去做散热路径的材料,或者据此断言「这块灌封胶散热不差」——而真正决定稳态温升的 k 可能低了两个数量级。

同型的第二面(τ 的读反):「这个件响应慢,是因为它热容大」。τ = ρcV/(h·A_conv) 里 h 与 A_conv 在分母,同一块件换一套冷却边界,τ 能差几倍。响应慢常常不是热容大、是散热差;照「减小热容」去改可能一点效果都没有,而改冷却边界立刻见效。因此判「τ 由谁主导」必须把分子分母一起看,不许只看热容那一半。

第三面(网格无关性的读反):「网格无关性验证通过,结果就可信」。它只证明离散误差收敛了,并不证明收敛到的那个答案是对的。边界条件、接触热阻、各向异性物性取错时,网格越细,只会让模型更稳定地收敛到一个错的值。由此立一条全课纪律:离散误差与模型误差是两本账,网格无关性只结清前一本。

本课的适用范围:交付判据,不交付网络与取值

本课的目标是把量级判对、把失真位置事先指出来,不是把温升算到小数点后一位。判据用到的那些输入,本课当作已知的边界参数代入:h 的定义与估法见 A1-02 对流换热:自然对流、强制对流与换热系数估算;k、ρ、c 的取值与温度相关性见 A1-07 热管理常用工质的热物性与查用、H1-01 热管理用铝合金牌号、性能与选材;把判据接成热网络、标定参数、降阶与实时见 A1-06 热阻网络法与集总参数建模入门——本课第 2 讲只搭到「一维等效热阻的串并联」这一步,作为降维的落点。

也正因如此,本课能给你的东西是可以带走的:一张三个比值的估算单、两套特征长度口径的分辨方法、一份「这个模型答不了什么」的交付清单。它们不依赖你用哪个求解器。

七讲的路线

第 1 讲 场方程与接口:傅里叶定律的三个旋钮、导热微分方程四项各删掉一项等于做了哪条假设、找对称性这个降维动作、三类边界条件,以及那个不属于这三类的固体贴合面;末尾立起对这条常见读反的辨析。 第 2 讲 一维降维:热阻比判据、三条几何前提、肋片方程与肋效率、串并联搭法,以及「简化必然偏保守」为什么不成立。 第 3 讲 Bi 与它的两个用途:特征长度的两套口径、Bi 判不了的三件事、引用纪律。 第 4 讲 多维与离散:必须上多维的三类触发、扩散热阻、节点能量平衡、显式格式的稳定性上限、网格无关性怎么做,以及两本账。 第 5 讲 时间判据:τ、Fo、稳态/准稳态/瞬态三区、必须瞬态的三类场景、热穿透深度。 第 6 讲 尺度效应:散热面按 L² 而热容按 L³、两个时间尺度标度不同导致缩放后主导项翻转、h 本身也在变、相似准则不完备时不许外推。 第 7 讲 部件实例与隐藏项:母排焦耳热与导热的耦合、冷板肋片折算成等效对流热阻、叠片电芯各向异性、接触热阻怎么进模型,以及交付时必须连着说清的三件事。

第 1 讲 傅里叶定律、导热方程与三类边界

这一讲要建立的不是解题能力,而是一张能看懂自己模型的地图。整车热管理里绝大多数导热模型的致命伤都不在求解器里,而在三个地方:方程里被悄悄删掉的某一项、模型与外部世界接口处填错的边界类型、以及两块固体贴合的那个面被当成了"没有热阻"。这三处的错都有一个共同点——它们在结果里长得和"算对了"一模一样,不报错、不发散、云图很漂亮。

本讲末尾还会把全课最容易读反的一条一次讲透:热扩散率 α 大,不等于导热好。这一条做反,会直接改掉你的选材结论。

本讲开头先把几条全课口径说死(后面各讲一律照此,不再重复):

  • 三个"热流类"的量不许互写:q 是热流密度(W/m²)· Q 是热流功率(W)· q̇ 是体积产热率(W/m³)。同一个字母 q,带不带点、带不带面积,是三个不同量纲的量;本课任何一处写下它们,公式两侧单位必须自洽。
  • 温位一律写 ℃,温差一律写 K,两者同式出现时各自标注。唯一的例外是辐射相关表述(必须用绝对温标)——本课不做辐射定量,它与导热并联的处理方式与线性化前提归 A1-03 辐射换热与整车热害中的辐射屏蔽,这里只提一句、指明去处。
  • 面积只允许三个下标进公式:A_c(导热截面,热流穿过的那个面)· A_conv(参与对流的散热面)· A_contact(接触界面的名义面积)。禁止裸写 A。每写下一段热阻,你要能立刻说出"这里的 A 是哪一个"。
  • 特征长度一律写 L_c 或 L_H(第 3 讲定义两套口径),不裸写。注意,与它不是一回事的是导热热阻 R_cond = L/(k·A_c) 里的那个 L——那是这一段热路自身沿热流方向的长度,是几何尺寸,不是判据用的特征长度,两者不可互代。
  • 本课的 α 只表示热扩散率。部分文献用 α 表示对流换热系数(本课一律写 h)、A1-03 辐射换热与整车热害中的辐射屏蔽 用 α 表示太阳吸收比、母排那一讲的电阻温度系数写 α_R——四个量各自带自己的写法,本课不混用。

1.1 傅里叶定律:整个热管理设计只有三个旋钮

傅里叶定律把导热写成一组比例关系:

q = −k·∇T

三个量各自的角色很清楚:∇T 是驱动力(温度梯度,指向温度升高的方向);k 是材料的导热能力;q 是单位面积上的热流。负号不是记号约定,它是方向——热自发地由高温流向低温,与梯度方向相反。把负号漏掉或写反,整条热路的方向就反了,而这种错在一维手算里往往要到对不上实测才被发现。

把它在一段等截面直路上积分,就得到工程上真正天天用的那个形式:

Q = k·A_c·ΔT_∥ / L ⇔ R_cond = L/(k·A_c)

(ΔT_∥ = 这一段两端、沿热流方向的温差,K;本课温差一律带下标写明是哪两点之间,不裸写 ΔT。L = 这一段热路自身沿热流方向的长度,m。)

所以这个式子说明:降低一段导热热阻,你手上只有三个旋钮——换材料(提高 k)· 缩短路径(减小 L)· 加大截面(增大 A_c)。没有第四个。所有听起来花样繁多的降温手段,最后都落回其中之一:铝换铜是第一个,把热源挪到离冷面更近是第二个,加肋片、加扩热板是第三个。当有人给你一个"新方案"而你说不出它动的是哪个旋钮时,多半是它其实动的是对流那一侧(h 或 A_conv),那是另一本账,归 A1-02 对流换热:自然对流、强制对流与换热系数估算。

k 的跨度有几个数量级,这正是热管理里"选材"能起决定作用的原因(下列均为常温值,来源等级 = 材料手册 / 物性库,未挂任何标准号;空气与水由物性库 CoolProp 在 300 K、101.325 kPa 算出):

介质 k(W/(m·K)) 备注
纯铜 ≈ 398 不同牌号与状态(纯铜/铜合金、退火/冷加工)差异可观,按你的牌号复核 → H1-02 铜、钢及复合材料的应用场景
铝合金 6061-T6 ≈ 167 注意,1 系纯铝可达 200 以上,压铸铝合金明显更低,因此不得把 6061 的值当作"铝"的通用值 → H1-01 热管理用铝合金牌号、性能与选材、H1-04 铸造铝合金(压铸为主)的材料特性与热管理壳体选材
304 不锈钢 ≈ 16.2 1.6 反例对照的一端(见 1.6)
工程塑料 典型 0.2~0.5 只给区间,牌号与填充体系跨度大;本课不用它做任何算例,选材归 H1-03 轻量化材料与薄壁化工艺趋势
静止空气 ≈ 0.026 300 K、101.325 kPa

今天就能用的一步:把你手上那条热路按段拆开,每一段算出它的单位面积热阻——导热段算 L/k、对流段算 1/h——然后排序。排在最前面的那一段,才是值得花钱改的那一段;顺手去把一块本来就不是瓶颈的铝板加厚,是热管理里最常见的无效投入。

拿两个数当场做一遍,结论会让很多人意外:

  • 0.1 mm 的静止空气层:L/k = 1×10⁻⁴ m ÷ 0.026 W/(m·K) ≈ 3.85×10⁻³ K·m²/W = 38.5 K·cm²/W
  • 5 mm 厚的 6061 铝板:L/k = 5×10⁻³ m ÷ 167 W/(m·K) ≈ 3.0×10⁻⁵ K·m²/W = 0.30 K·cm²/W
  • 两者之比 ≈ 130 倍

(换算口径:1 K·cm²/W = 1×10⁻⁴ K·m²/W;工程上报出来的界面热阻多数用 K·cm²/W,而公式里要代 K·m²/W,这一步换算是本课最容易丢一个 10⁴ 的地方。)

这就意味着,一条头发丝粗细的空气隙,热阻可以比几毫米厚的铝还大两个数量级。这不是一个修辞,它就是接触热阻的物理根源——1.5 会回到这里。

本节三个易错点:① 把 q 和 Q 混写(差一个面积,量纲都不同);② 以为负号只是记号约定;③ 在各向异性材料上把 k 当标量用——叠片电芯上这么做会把厚度方向的传热高估数倍,第 7 讲展开。


1.2 导热微分方程:你删掉哪一项,就等于做了哪一条假设

完整的导热微分方程有四项:

ρc·∂T/∂t = ∇·(k∇T) + q̇

逐项读一遍,每一项都对应一个你迟早要做的建模决定:

  • ρc·∂T/∂t:蓄热项。材料把热存起来的能力乘以温度变化速率。删掉它 = 假设稳态——你声明"在我关心的时间尺度里,温度已经不再变了"。
  • ∇·(k∇T):扩散项。热量在材料内部的输运。把 k 提到散度外写成 k∇²T = 假设常物性且各向同性——你声明"这块材料各个方向导热一样,且 k 在我这个温区里不随温度变"。
  • q̇:内热源项。删掉它 = 假设体内不产热。

初学者常把"稳态方程""常物性方程""拉普拉斯方程"当成三个彼此独立的公式去背。它们不是——它们是同一个式子逐项删减的结果:

你删了什么 得到的式子 你因此声明了什么
删蓄热项 ∇·(k∇T) + q̇ = 0 稳态:我关心的时间尺度内温度不再变
删蓄热项 + 删 q̇ ∇·(k∇T) = 0 稳态且体内不产热
k 提到散度外 + 常物性 ∂T/∂t = α·∇²T + q̇/(ρc),α = k/(ρc) 各向同性、k 不随温度变
全删只剩扩散 ∇²T = 0 稳态、无源、常物性、各向同性(拉普拉斯方程)

为什么这一步值得当成一个动作来做:反过来用这张表,只要你说得出自己删了哪几项,就说得出自己的模型在什么条件下失效。这句话是可交付的——它正是第 7 讲要求的"这个模型答不了什么"的第一块内容。

整车热管理里最危险的一删是 q̇。母排的焦耳热、电芯的产热、电机的铜损铁损,全都是体内产热,它们进的是 q̇ 这一项,而不是某个面上的热流边界。把一个体内均匀产热的件写成"一个面上给定热流",温度分布的形状从一开始就错了(真实的最高温在体内,而不在加热面上)。

本节三个易错点:① 件是叠片电芯(各向异性)却把 k 提到散度外——此时 k 必须留在散度里按张量赋值,方向写作 k_∥(面内)与 k_⊥(厚度方向),取值与倍数在第 7 讲;② 物性随温度变化显著时仍按常物性算——金属在热管理温区内变化有限,塑料与界面材料不然;③ 写稳态方程时顺手把 q̇ 也删了,这是两条独立的假设,别打包。

导热方程逐项删减的决策树:完整式 ρc·∂T/∂t = ∇·(k∇T) + q̇ 的三项分别标出蓄热项、扩散项、内热源项,四层判断分别问时间尺度内温度是否还在变、各方向导热是否一样、物性随温度变化是否可忽略、体内是否产热,每条分支的终端方框下都写出这一步所做的假设;稳态分支旁特别标出该式里没有 ρ、没有 c,因此也没有 α。本图为决策示意,不含任何数值,不得据图取值。
图2 该读出两件事:那些"稳态方程""常物性方程""拉普拉斯方程"不是彼此独立的公式,而是同一式逐项删减的结果;每删一项就绑定一条假设,说得出删了哪几项,就说得出模型在什么条件下失效。特别注意稳态分支那一格里没有 α——它是 1.6 那条被读反的说法的出处。本图为决策示意,框的大小与位置不含任何定量含义,不得据图取值。

1.3 选坐标系不是"按件形套公式",是找对称性

这是全课第一次降维机会,也是"能手算"与"必须上有限元"之间的那条分界线。

三套坐标系各自的主场:

  • 直角坐标——板、壁、方块:冷板基板、电池壳壁、母排的矩形截面。
  • 柱坐标——管、圆棒、圆柱:冷却液管路、圆柱电芯、圆形连接器。
  • 球坐标——球形件,工程中少见。

真正起作用的不是"件长得像什么",而是坐标系能不能把某个方向上的变化整个消掉。一段圆管,热量沿径向由内壁流向外壁,周向处处相同,因此在柱坐标下直接写成一维径向导热,有解析解;同一个问题硬套直角坐标,就是一个二维甚至三维问题。一次成功的坐标选择,常常把求解规模从上万节点降到一个解析式——这是本课性价比最高的一步,而它发生在打开任何软件之前。

但对称性必须几何与边界同时成立,这是本节唯一真正会咬人的地方。圆管几何上是轴对称的,可一旦它一侧贴着冷板、另一侧对着空气,周向就不再对称了——此时"径向一维"当场失效,而几何看上去毫无变化。工程件上这种情形是常态而不是例外:贴合、遮挡、单侧迎风、单侧受热,每一个都会打破边界的对称性。

今天就能用的两问清单(顺序不能反):

  1. 几何在哪个方向上不变?(绕轴旋转?沿长度平移?镜像?)
  2. 这个方向上,边界条件也处处相同吗?——把每一个面逐个点名:它贴着什么、朝向哪里、有没有被挡住。

两问都答"是",才能降维;只有第 1 问答"是"就降维,是本课最常见的一次错。

本节三个易错点:① 对称性只看几何不看边界(上面那条);② 多热源件上误以为"形状规则就能降维"——热源的分布也是边界的一部分;③ 把"用柱坐标"当成"一定是一维",这是两个动作,选坐标系只是给降维创造了条件。

三行对照:直角坐标对应板/壁/方块并标出主热流方向,柱坐标对应管/圆棒/圆柱并用大括号标出周向对称因而径向一维,球坐标对应球形件;每行右侧的警示框写明对称性要同时看边界,第二行下方另画反例——同一段圆管下半周贴冷板、上半周对空气,周向温度分布不均匀,打禁止符号标注此时不可按径向一维。本图为几何示意,尺寸比例不代表任何真实件,不得据图取值。
图3 该读出两件事:选坐标系的动作是找对称性,选对了问题直接降维;而对称性必须几何与边界同时成立——几何对称、边界不对称时降维就是错的,这正是工程件上最常见的情形。本图为几何示意,尺寸比例不代表任何真实件,不得据图取值。

1.4 三类边界条件:模型与外部世界的唯一接口

方程管的是件的内部,边界条件管的是它与世界的接口。内部算得再准,接口错了结果就是错的——而三类边界的区别,说到底就是一句话:你手上已知哪个量。

类别 你已知的量 写法 典型场合
第一类(定温 / Dirichlet) 边界温度 T = T₀ 贴着一个温度被牢牢控住的大热沉;恒温槽试验
第二类(定热流 / Neumann) 边界热流密度 q 给定;绝热面是其特例(q = 0) 已标定的加热片;对称面;判为可忽略漏热的面
第三类(对流 / Robin) h 与流体温度,壁温由求解得出 q = h·(T_w − T_∞) 冷板对冷却液、壳体对空气、换热器对气流

整车热管理里绝大多数面是第三类。因此这类边界的可信度,直接由 h 的取值精度决定;h 怎么估、它随流态与几何怎么变,是 A1-02 对流换热:自然对流、强制对流与换热系数估算 的内容,本课只把它当作已知的边界参数代进判据。

本节要带走的一条结论:把一个实际上在对流的面写成定温边界,等于给模型喂进一条现实中不存在的无限散热能力。定温边界的含义是"无论多少热量流到这个面上,它的温度都不变"——那只有在 h 无穷大时才成立。后果有两层,第二层更麻烦:温升被系统性低估,而且低估了多少你估不出来(低估幅度取决于那个被你抹掉的、未知的真实 h)。这不是一个"偏保守一点"的近似,它是一个方向确定、幅度不可知的乐观误差。

今天就能用的一步——定温边界的反解自检:如果你已经用了定温边界(有时确实合理),就把求解结果里流过该面的总热流 Q 取出来,再反算"要维持这个壁温,需要多大的 h":

h = Q / (A_conv·(T_w − T_∞))

把这个 h 拿去和 A1-02 对流换热:自然对流、强制对流与换热系数估算 给出的量级对照:如果它落在自然对流的量级里,定温边界是荒唐的;如果它高到只有相变或强液冷才能达到,也是荒唐的。判出来之后的下一个动作:把该面改成第三类边界,用一个真实可达的 h 重算一遍,两次结果的差,就是这条假设值多少钱。

本节三个易错点:① 图省事把散热面写成定温;② 把"环境温度"直接当成壁温——第三类边界里 T_w 是解出来的,不是给定的,这两者差的正好就是对流那一段的温差;③ 忘了"绝热"也是一个假设——真实的绝热面往往在漏热,只是漏得少;当那个面恰好是你的最热点所在时,这条假设就该被单独核一遍。


1.5 两块固体贴合的那个面,不属于这三类

这是本课各种"隐藏项"里最常被漏掉的一项,而且漏掉的方向是确定的乐观侧。

母排接头的两段导体、电芯与冷板之间的界面材料层——它们的两侧都是固体,没有流体温度可对,因此既不是定温、也不是定热流、更不是对流。它由界面条件表达:

q = h_c·ΔT_c = ΔT_c / R''_c , h_c = 1/R''_c

其中 h_c 是接触热导(W/(m²·K))、R''_c 是面积接触热阻(K·m²/W)、ΔT_c 是界面两侧的温度跳变(K)。注意这里说的是跳变:在理想化的界面模型里,温度在穿过这个面时不连续地掉一个台阶,台阶高度就是 ΔT_c。这个台阶正是接触热阻在温度场里唯一可见的形式。

把它接进一维热阻链时,还要再除一次面积:

R_c = R''_c / A_contact

R''_c 与 R_c 不许互代——两者差一个 A_contact,而工程上报出来的数多数是 R''_c。引用任何一个"接触热阻"数值之前,先看它的单位里有没有面积:有面积(K·cm²/W、K·m²/W)就是 R''_c,没有(K/W)才是某一个具体界面的 R_c。

为什么这一段最容易被删掉:因为删它不需要任何动作。求解器里两个零件贴在一起,默认的装配关系是"绑定 / 理想接触",而理想接触的含义就是 R''_c = 0——界面温差被算成零,热路上凭空少了一段。没有警告、没有报错,模型跑得比加了界面条件还顺。

它有多大?面积接触热阻 R''_c 常见落在 0.2~0.6 K·cm²/W,多数量产平台在 0.35 K·cm²/W 附近。典型值,不是验收依据,你的平台须实测确认;本项目未取得可核出处——它来自行业工程实践的常见范围,不是任何标准的条款数据。→ 拿到它的路:按界面热阻/界面材料的相关测试方法实测(测试方法的去处在 7.4 一并给出),取值与压力的关系、界面材料选型归 D2-05 冷板与电芯的接触热阻与 TIM 应用、H5-01 导热界面材料(导热垫/凝胶/灌封)选型。

把它和 1.1 那笔账放在一起看,能读出一件对建模很有用的事:连续 0.1 mm 空气层的 38.5 K·cm²/W,比典型真实界面的 0.35 K·cm²/W 大约 100 倍。由此可见,真实贴合面并不是一层均匀的空气膜——它靠稀疏的真实接触点和极薄的间隙导热。这直接解释了本节最后那个易错点。

今天就能用的三步(把界面合法地放进模型):

  1. 场模型里:在两个零件的接触面上显式设置 h_c 或 R''_c,而不是留着默认绑定;
  2. 一维网络里:在热路上串进一段 R_c = R''_c/A_contact,A_contact 用名义接触面积;
  3. 无论哪一种:按区间的上下限各算一遍,报区间、不报单值。因此,判出来之后的下一个动作:如果上下限给出的最热点温度落在同一个结论里(都合格或都不合格),这个参数就不必再花钱去测;如果结论跨了线,那它就是你下一个必须实测的量。

本节三个易错点:① 把接头面当成第三类边界,凑一个虚构的流体温度——那等于给模型编了一条现实中不存在的散热路径,比不建界面还坏;② 直接用求解器默认绑定(把这一段设成零);③ 按界面材料的厚度和 k 建实体网格,就以为算到了接触热阻——接触热阻主要来自两侧微观不平整造成的真实接触点稀疏,不是填隙材料本体的厚度热阻,两者是两笔账,要分开计。

四格并排剖面示意:第一格定温边界,壁面温度被钉死在 T₀;第二格定热流边界,给定 q 并标出绝热面是其特例;第三格对流边界,壁面外侧是流体 h 与 T_∞,温度由壁温平缓过渡到远场,式为 q = h·(T_w − T_∞);第四格是界面条件,两块固体贴合、贴合处画锯齿表示微观不平整、缝隙里标注空气,界面两侧温度出现台阶跳变 ΔT_c,式为 q = h_c·ΔT_c = ΔT_c/R''_c,并标明两侧都是固体、没有流体温度可对。本图为机理示意,壁厚、缝隙与温度线斜率均不代表任何真实尺度,不得据图取值。
图1 该读出四件事:前三格的区别只在于你手上已知哪个量(温度、热流、还是 h 与流体温度);第三格的壁温是解出来的而不是给定的;第四格的温度是跳变的,这个跳变就是接触热阻的可见形式;第四格不在前三类里,把它当成对流边界去凑一个流体温度,等于编了一条不存在的散热路径。本图为机理示意,不得据图取值。

1.6 热扩散率 α 只管瞬态,稳态温度场里根本没有它

这是全课最容易被读反的一条。先把那句错话原样摆出来,再拆掉它:

误区:"热扩散率 α 大,导热就好,稳态温升就低。"

这句话把两件互不相干的事接成了一条因果链,按它去选材、去判断谁更能散热,方向可以整个做反。

真相是两条泾渭分明的分工:

  • 稳态温度场只认 k。把 1.2 那张表往回翻一格:稳态导热方程是 ∇·(k∇T) + q̇ = 0——这个式子里没有 ρ、没有 c,因此也没有 α。最终的温度分布由 k、几何和边界三样东西决定,与材料存不存得住热完全无关。
  • α = k/(ρc) 只管瞬态。它出现在瞬态项里,决定温度扰动传播多快,不决定最终温度有多高。

一个当场可以验算的反例(α 由 α = k/(ρc) 现算,三个输入一并列出,你可以自己按一遍;常温 20 ℃ 固体值来源等级 = 材料手册,空气与水由物性库 CoolProp 在 300 K、101.325 kPa 算出,均未挂任何标准号):

介质 k(W/(m·K)) ρ(kg/m³) c(J/(kg·K)) α = k/(ρc)(m²/s)
纯铜 398 8930 385 ≈ 1.16×10⁻⁴
铝合金 6061 167 2700 896 ≈ 6.9×10⁻⁵
304 不锈钢 16.2 7900 500 ≈ 4.1×10⁻⁶
常温常压空气 0.026 — — ≈ 2.2×10⁻⁵
水(300 K) 0.61 — — ≈ 1.5×10⁻⁷

看这两行:空气的 α(≈ 2.2×10⁻⁵)约为 304 不锈钢(≈ 4.1×10⁻⁶)的 5 倍;而两者的 k 是 0.026 与 16.2 W/(m·K),相差约 600 倍,方向完全相反。

因此,按 α 排序去选"导热好的材料",会把空气排在不锈钢前面。而没有人用空气去导热。

第二组对照把这件事钉得更死:水的 k 约为空气的 23 倍,而水的 α 反而比空气小约两个数量级。同一对介质,两个排序不但不同,而且相反。(这两个数在本课只用于此处的对照,不参与任何换热计算——那归 A1-02 对流换热:自然对流、强制对流与换热系数估算、A1-07 热管理常用工质的热物性与查用。)

今天就能用的两句判据:

  • 有人问"这个件最后会到多少度",就看 k(以及几何和边界)。换一种 ρc 完全不同、k 相同的材料,稳态温度一点都不会变。
  • 有人问"这个件多久才热起来/多久才降下去",就看 α,以及第 5 讲的 τ。
  • 因此,听到"导热好"这三个字,先反问一句:你说的是最终温升低,还是响应快?这两件事可以指向不同的材料,而"导热好"这个词把它们混成了一个。

本节三个易错点:① 用 α 判稳态温升(本节主线);② 用 k 判响应快慢——铜的 k 高、α 也高,容易让人误以为两者总是同向,而空气与水的对照立刻推翻它;③ 把 α 与对流换热系数、与太阳吸收比共用一个符号(本讲开头的口径已把四个量分列)。

顺带说破一个同型的错,第 5 讲会正式处理:"这个件响应慢,是因为它热容大"——时间常数 τ = ρcV/(h·A_conv) 里,h 与 A_conv 在分母。同一块件换一套冷却边界,τ 能差几倍;响应慢常常不是热容大,是散热差。判"τ 由谁主导"必须把分子分母一起看。

左右两张并排的对数刻度水平条形图,同一组五种介质(纯铜、铝合金 6061、304 不锈钢、水、空气)在两张图上各自按自己的量从大到小重排:左图按 k 排序(单位 W/(m·K)),右图按 α 排序(单位 m²/s),每根条端写出数值标签,两图之间用细引线连接同一种介质使排序交叉可见;空气那两根用醒目色标出,右图旁标注空气的 α 比不锈钢还大约 5 倍。固体为材料手册常温值,空气与水由物性库在 300 K、101.325 kPa 算出,未挂任何标准号。
图4 该读出三件事:同一组材料在两张榜上的顺序不同,尤其空气与不锈钢的相对位置发生了对调;因此"导热好不好"这句话必须先问"你在说稳态还是瞬态";稳态温升看 k、响应快慢看 α,两者不能互相替代。条长按对数刻度绘制,不得据条长直接读比值(要读数值标签);材料牌号与状态不同时数值会变,不得作为选材依据。

本讲小结:傅里叶定律给你三个旋钮(k · 路径 · 截面),没有第四个;导热方程的四项里,你删掉哪一项就等于做了哪一条假设,说得出删了什么才说得出模型何时失效;选坐标系是在找对称性,而对称性必须几何与边界同时成立;三类边界的区别只在"你手上已知哪个量",把对流面写成定温=喂进一条无限散热能力;两块固体贴合的面不属于这三类,它是界面条件,而求解器的默认理想接触会把这一段热阻整个删掉;最后,稳态只认 k,α 只管瞬态——按 α 选"导热好的材料",会把空气排在不锈钢前面。

下一讲把"能不能降到一维"从"件薄不薄"这种几何描述,换成一条可以真正判真假的判据。

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

会员专属

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

查看会员方案