热失控与热扩散传播仿真建模：化学源项、喷发边界与可信边界

课程大纲

• 能把自家电芯的 ARC/DSC/TGA 数据折算成多步放热反应源项的参数组（指前因子 A、活化能 E_a、反应级数、各步起始温度与放热量分配），并写明拟合区间、Φ 热惯量修正与「跨化学体系不可套用」的理由
• 能在「固体导热＋等效对流边界 / 加气相组分输运 / 再加颗粒相」三个喷发建模层级之间按关心量选层，并给出质量流率时程、气相温度、气相组分、颗粒相是否解析、持续时间五项边界条件的取值与来源
• 能判出算例的时间步卡点在化学时间尺度还是传热时间尺度，选出源项隐式/算子分裂/次循环中的处置方式，并用累计放热守恒校核证明推进方式没有吃掉或造出能量
• 能在包级 CHT、集总热网络、传播时序降阶三种口径之间按判据选型，并把隔热层、汇流排热短路与冷板横向热桥在三种口径下的等效方式互相换算、互相校核
• 能把传播时序的仿真—试验偏差按源项 / 接触热阻 / 喷发路径 / 触发实现四路归因，并用每一路的扰动量级判断该路是否足以解释实测偏差
• 能交付一份含源项来源与适用区间、建模层级与被略物理、改了物理的数值处置项、口径间校核残差、结论效力分级的可信边界声明，并指出哪些结论必须由整包试验定案

1. 多步放热反应源项的构造与参数来源
   • 源项形式是若干步 Arrhenius 反应的体积放热之和，决定包级传播时序的是「各步起始温度与放热量 ΔH_j 的分配」，而不是总放热量一个数：总量对、分配错，时序会整体偏快或偏慢，且换一套隔热组合后误差方向还可能反转。
   • 两级数据分工写死：DSC/TGA 定「有几步、各步落在哪个温度段、各步反应焓的相对占比」，ARC 定「整芯尺度上这几步合起来的速率与总量」。T1/T2/T3 的定义、判读与移动因素（体系、SOC、老化路径、电芯规格）见 D1-03，本课只取其数值做拟合输入，不重讲判读。
   • 进拟合之前必须先做 Φ 热惯量修正：未修正的 ARC 曲线被夹具热容摊薄、结果偏乐观（危险侧），直接拟合会把 A 与各步 ΔH 一起标低，传播算慢、隔热需求算低。Φ 的折算口径（真实自产热速率 ≈ Φ×实测 dT/dt、T3 按 T_start 起算折、TMR 按除以 Φ 折）沿用 D1-03，本课不重推。
   • 动力学拟合是欠定问题：A 与 E_a 在单条 dT/dt 曲线上高度相关（补偿效应），能得到多组等效参数，一外推到拟合区间之外的温度或升温速率就散开。工程做法是多升温速率/多 SOC 数据联合拟合，并把拟合区间外的使用显式声明为外推。
   • 跨化学体系一律不可套用：LFP、三元、钠离子、固态的反应步数、各步起始温度与产气行为都不同（差异见 D1-05），套别家参数正是 D5-01 已点名的「源项一错后面全错」；同体系不同规格也不能互抄，理由与 D1-03「T1/T2/T3 是实测量不是材料常数」同源。
   • 源项之外还有两类必备输入：电芯的 c_p 与各向异性 k（取值来源、测法与向电芯厂索取的验收清单见 D1-07，本课只用不测）；以及内短路电功率项——触发实现为过充/针刺时不可略，其量级判据（在 R_short ≈ R_int 处取极大，不是越小越危险）见 D1-03 的 P_short 式。
   • 参数不确定性怎么传到传播时序：把 A、E_a、ΔH 分配比与界面接触热阻列成输入分布表，交付「传播时间」的分位数而不是单值。本课只负责给出这几项的分布来源与相关性假设（同一次拟合导出的 A 与 E_a 强相关，独立抽样会把带宽算窄），传播方法、收敛判据与分位数裕度的框架本体见 J7-06，不重讲。
   • 落笔纪律：任何一组 A/E_a/ΔH 交出去必须同时附「由哪几条曲线、哪个温度区间、哪个 SOC 与老化档拟合、Φ 取值多少」，否则下游无法判断能不能用；本课不给任何典型值，与 D1-05「不给典型值、把量化输入的合法来源指向技术协议与实测报告」是同一条纪律。
2. 喷发的建模层级判定与边界条件给法
   • 三个层级按「要回答什么」选、不按「算力允许什么」选：L1 固体导热＋等效对流边界，回答邻芯会不会被推过 T1、方案之间怎么排序；L2 加气相组分输运与温度场，回答烟气冲刷到哪些结构件、排烟通道温度；L3 再加颗粒相/多相，回答熔融物与高温颗粒的打穿与二次点火风险。判据是「关心量对喷发路径的敏感度」，敏感就必须升层。
   • L1 的两个关键取值必须成对标定：等效对流系数 h_eff 与被冲刷面积 A_wash 在结果里以乘积出现、会互相补偿，单标一个能把另一个的错误吸收掉，换个喷阀朝向立刻失效。A_wash 取实际被冲刷面积而不是总表面积（取全严重高估）这条来自 D5-01，本课补的是怎么定「实际」——按喷阀朝向与包内自由容积做射流覆盖判定。
   • L1 的失效边界要显式声明：等效对流边界把喷发压缩成一个时变面热流，同时丢掉组分、颗粒与流动路径三件事；一旦拿它去签排烟口温度或结构件热害结论，就已经越界——那类结论的通用方法体（长瞬态自然对流＋辐射、敏感件温度限、隔热罩整改）是 J2-05，本课只交源项接口。
   • 喷发源边界条件五项缺一不可（承 D5-02 的清单，本讲把清单变成建模输入）：质量流率时程 ṁ(t)、气相温度、气相组分、颗粒相是否解析、持续时间。ṁ(t) 的来源是 L1-11 的产气速率实测与 D5-02 的包内压升平衡反算，本课只讲怎么把它变成入口边界，不讲怎么测、也不讲泄压面积怎么定。
   • 时程形状比峰值单值更影响结果：把 ṁ(t) 简化成方波会抹平喷发段峰值功率，而峰值功率正是 D5-01 判「冷却不能当主防线」的那个量；至少要保留起始时刻、上升沿与主放热持续时间三个特征，起始时刻由 D1-03 交付的产气/喷发时刻定。
   • 颗粒相什么时候必须解析：结论落在「打穿 / 二次点火 / 熔融物沉积」这一类时必须解析；落在「邻芯温升」这一类时通常可并进气相的等效焓。多相方法族本体（VOF / 拉格朗日离散相 DPM / Eulerian wall film / SPH 各管哪一段液相、各自算力量级与失效边界）引 J3-06 不重讲，本课只做喷发源专属的层级判定。
   • 组分怎么给：H₂/CO/CH₄/CO₂ 与电解液蒸气的混合物，成分来源是 D1-03 的产气分析结果，采样与分析的落地方法见 L1-11。是否把包内燃烧算进去必须显式声明——这对应 D5-01 的两本能量账（Q_TR,ARC 与 Q_comb）；纪律是「算了就写明口径，没算就在可信边界里写明本次结论不含喷出可燃组分的燃烧放热」。
   • 本讲同时是 D5-02「喷发源怎么给、建到哪一层」两处前指的唯一落点：D5-02 保留设计侧判据（按被喷到设计、不按环境温度设计，以及排烟方向、IP 死结、乘员舱再侵入的设计结论），本课交建模层级与边界条件，两边不重叠、不互相搬结论。
3. 强放热源的数值实现：刚性、时间步与守恒校核
   • 刚性来自 Arrhenius 的指数项——温度小幅上移让反应速率跳一个台阶（机理式见 D1-03），化学特征时间可比传热特征时间小几个数量级。这与 J2-04 讲的「固体慢、流体快」是两回事：那里的处置是把慢的加速（稳态 CHT、冻结流场、固体时间步缩放），这里必须把快的解出来，把固体时间步缩放搬过来会直接把反应跨过去。
   • 开算前先判卡点：分别估化学时间尺度与传热时间尺度并比出刚性比，说明时间步被哪一侧卡死。卡在化学侧时加密网格不解放时间步、只是更慢——这一条与 J3-06「卡在转速侧时加密气隙网格不解放时间步」同型，都属于「先说卡在哪一侧再动手」的纪律。
   • 工程可行的三条路都不改物理：源项隐式/半隐式处理（把 ∂q̇/∂T 进对角）、算子分裂（化学步与传热步分开、各自用合适步长）、对源项做次循环。改了物理的做法——人为限幅源项、把放热摊到更长时间、强行加大亚松弛——必须逐项写进可信边界声明。
   • 亚松弛与限幅是本讲最隐蔽的失效模式：为压住振荡松弛过头等于人为把放热速率降低，表现为传播时序整体偏慢，而残差曲线还很好看；「跑完了没发散」不构成稳定性证据。
   • 守恒校核是唯一兜底：把全域累计化学放热对时间积分，与由源项参数解析算出的总放热量比对，残差超过百分之几量级就说明推进方式在吃能量或造能量。带喷发出流的算例必须把出流带走的焓单列，否则残差会被误判成源项问题。这一步要逐算例做，不能靠「上次做过」。
   • 稳定性判据写成可检查的形式：单步内单元最大温升有上界、化学子步内温度不越过下一步反应的起始温度、能量残差随时间的漂移量有上界。三条都要给数值并写进算例记录，判据的具体取值由本项目的网格与源项参数确定，本课不给通用数。
   • 与 J7-08 划清边界：J7-08 管一维系统模型的 DAE 指标与代数环、步长与容差、一致初值与积分器选型，对象是 1D 求解器；本讲管三维 CHT 里化学源项与传热的时间尺度分离，对象是有限体积推进。两门的「刚性」同名不同物，结论不可互相搬用。
   • 前提声明：本讲默认网格与步长收敛已按 J7-01 引的口径做过。拿一个未做收敛研究的模型讨论源项数值实现，分不清误差来自源项还是来自离散，后面第 5 讲的四路归因也就无从做起。
4. 包级模型架构选型与分工：CHT、集总热网络、传播时序降阶
   • 三种口径各回答不同问题：包级 CHT 给空间分辨的温度场与冲刷路径，最贵，用来定隔热层与冷板热桥这类局部设计；集总热网络给分钟级的传播趋势与措施排序，用来在方案期扫大量组合；传播时序降阶模型给「第 n 颗在什么时刻被触发」，用来与试验对标、与整车报警时序对齐。D5-01 只并列了前两种且不给选型判据，本讲补的正是这条判据。
   • 选型判据三条，任一条不满足就不许用低阶口径：关心量是空间分布还是时序标量；隔热层与邻芯的温度梯度是否需要分辨（用毕渥数判能否集总，Bi 定义见 J2-04）；要扫的方案数与单次算例成本之比。
   • 三个口径必须互相校核而不是只跑一个：低阶口径的参数（等效热阻、等效热容、触发判据阈值）由高阶口径标定，标定后用低阶口径扫方案，选中的方案再回高阶口径复算一遍。只跑集总网络就签隔热厚度，等于把标定这一步整个跳了。
   • 隔热层的三种等效必须能互相换算：CHT 里按实体解并给高温段的 λ(T)（常温导热系数不算数，高温段收缩/烧穿/失去支撑的选型判据见 D5-02）；热网络里折成面热阻 δ/λ；降阶模型里只保留「延迟时间」一个参数，用傅里叶数标定。换算对不上，说明某一层的假设有问题。
   • 汇流排热短路的等效：CHT 里必须建实体（铝铜把热直接送到邻芯极柱，是 D5-01 点名的第三条传播路径）；热网络里是一条低热阻支路，最容易因为「截面看起来很小」被漏建——判据是把该支路的热导与芯间隔热层的热导比一比，同量级就不能省。汇流排未熔断时并联支路持续灌流是第四条路径，建模上是一个额外的电功率源项，不是热阻问题。
   • 冷板热桥的等效：冷板既是散热件也是横向热桥（D5-01 的案例正是冷板把热沿铝板铺给同排邻芯）。CHT 里靠实体几何自然体现；热网络里必须显式建「邻芯—冷板—邻芯」的横向支路，只建「电芯—冷板—冷却液」的纵向支路会把热桥整个漏掉，得出「再加厚隔热就够」的错结论。
   • 传播期间冷板到底搬不搬热要分两问：一问冷却液还在不在（管路是否被烧穿、失控瞬间泵停还是转，该接口判据属 D5-01 的设计侧，本课只把两种情形都建成算例）；二问它的搬热能力与喷发段峰值功率差几个量级（D5-01 已给出差一到两个量级的结论，本课不重算）。建模上的常见错是默认冷却液全程存在且流量不变。
   • 与既有建模课划清边界：D4-01 的模型选型是快充常规工况口径、J7-02 是控制导向降阶、D1-06 的分区热网络是量产包在线估计，三门都不是故障传播口径；本讲的降阶对象是传播时序，标定数据来自一次性破坏试验而不是常规工况，模型结构与标定纪律都不能互相搬用。
5. 与整包热扩散试验的时序对标与四路归因
   • 对标量是时序而不是温度峰值：整包试验能稳定测到的是「第 n 颗被触发的时刻」与「传播是否停止」，喷发段的峰值温度受测点位置与存活状况支配（测点在喷发环境下怎么活到试验结束、线缆烧断时数据怎么保全见 L1-11），拿峰值温度当对标量会把模型往测点误差上标。
   • 传播时序误差按四路归因，且查的顺序不能乱：先查触发实现（加热片功率与贴装、触发时刻定义，落地见 L1-11）与 Φ 修正做没做——这两项会整体平移时序；再查源项的 ΔH 分配与界面接触热阻（TIM/间隙折成界面热阻层的做法见 J2-04）——这两项改变时序的形状；喷发路径（h_eff 与 A_wash、喷阀朝向）与形状和幅值都有关，放在最后。
   • 归因必须写成可检验的形式：每一路单独（或按相关性联合）扰动一次，读第 2 颗触发时刻被推向哪边、推多少；若某一路的扰动量级不足以解释实测偏差，就不许把偏差归给它。这是本讲区别于「调参调到对上」的关键——参数调到对上但归因说不出来，模型换一个方案就不成立。
   • 单次破坏性试验的统计性质要说清：一次试验一个样本，重复性/再现性框架在此失效（该失效边界的完整论述见 L1-11）。因此对标产出不是一个「误差百分比」，而是「在这一次试验的条件下，模型时序落在哪个区间、哪些条件一变结论就不成立」。
6. 结论效力分级、法规接口与可信边界声明的交付
   • 仿真结论按效力分三级并写进交付：可用——方案之间的排序（哪个隔热/冷板组合更晚传播、传播能否被挡住的定性趋势）；有条件可用——传播延迟时间的量级（须同时声明源项拟合区间与喷发建模层级）；不可用——直接判「是否满足法规判据」。最后一条是 D5-01 已写死的「判据由整包试验定案」，本课不松口。
   • 法规接口按版本落笔、不单值化：GB 38031 的热扩散要求已由旧版「热事件报警后留出逃生时间」演进到不起火不爆炸方向；现行版对新申请型式批准车型与已获批车型分档实施，各档生效窗口按现行版本逐档确认，禁止写成一个日期。UN GTR 20 与 SAE J2464 同样按现行版本确认。合规判定与型式认证证据链本体见 M4-03，本课只提供仿真侧的输入与效力声明。
   • 可信边界声明的最小要素六条：源项的拟合数据来源与适用温度/SOC 区间、喷发建模层级与本次被略掉的物理、数值处置中改了物理的项、口径间校核的残差、对标的试验条件、本次结论的有效范围与效力分级。缺任何一项，下游都无法判断这份结果能用到什么程度——这与 J7-06 的评审级 UQ 结果表是同一条纪律，本课只补热扩散专属的这几行。
   • 交付纪律：「传播被挡住」这个结论必须同时给出它在哪些输入分布下成立（分位数口径承 J7-06），单点算例得出的「挡住了」不得进设计评审结论；同样，本课的算例结论不能反过来当作 D5-01/D5-02 的设计判据来源，那两门的判据各有其试验与设计侧出处。

• 以为源项只要总放热量对了就行，其实决定传播时序的是各步起始温度与 ΔH 的分配——总量对、分配错，时序会整体偏快或偏慢，且换一套隔热组合后误差方向还可能反转（判据是闭合校核只校总量、校不出分配，分配错只能靠第 5 讲的时序对标抓）。
• 以为 ARC 曲线可以直接拿去拟合动力学参数，其实未做 Φ 热惯量修正的曲线被夹具热容摊薄、偏乐观（危险侧不是保守侧），拟合出的 A 与 ΔH 系统性偏低，传播算慢、隔热需求算低；Φ 的折算口径见 D1-03。
• 以为动力学参数拟合出来就是材料常数、可以跨体系甚至跨规格套用，其实 A 与 E_a 在单条 dT/dt 曲线上存在补偿效应、本身只是一组等效解，换体系、换规格、换 SOC、换老化档都得重标（体系差异见 D1-05，「实测量不是材料常数」的理由见 D1-03）。
• 以为喷发建模层级越高越好、算力不够就退到最低层，其实层级要按关心量选：拿只做等效对流边界的结果去签排烟口温度或结构件热害结论，等于把被略掉的组分、颗粒与流动路径当成不存在（那类结论的方法体是 J2-05）。
• 以为等效对流系数 h_eff 与被冲刷面积 A_wash 可以分别标定，其实两者在结果里以乘积出现、会互相补偿——单标一个能把另一个的错误整个吸收掉，换个喷阀朝向立刻失效；判据是两者必须成对标定、成对交付。
• 以为算例「跑完了没发散、残差曲线很平」就算数值处置到位，其实亚松弛或源项限幅过头会把放热速率悄悄压低、残差反而更好看；唯一能抓住它的是累计放热守恒校核残差。
• 以为 J2-04 那套「固体慢、流体快」的瞬态加速能处理化学源项，其实那里的处置是把慢的加速，而化学项是必须解出来的快过程；把固体时间步缩放搬过来会直接跨过反应，得到一个不起火的假结论。
• 以为 J7-08 讲的刚性可以直接搬来用，其实那门课的对象是一维系统模型的 DAE 指标与积分器选型，本课第 3 讲的对象是三维有限体积推进里化学源项与传热的时间尺度分离，两个「刚性」同名不同物。
• 以为集总热网络里冷板就是散热件，其实它同时是横向热桥（D5-01 的案例正是冷板把热沿铝板铺给同排邻芯）；只建纵向支路会得出「再加厚隔热就够」的错结论，判据是把横向支路热导与芯间隔热层热导比一比。
• 以为汇流排截面小、可以在热网络里省掉，其实铝铜热短路是 D5-01 点名的第三条传播路径，且汇流排未熔断时并联支路持续灌流是第四条——后者在建模上是一个额外的电功率源项，不是热阻问题，两者都省掉会同时低估传入邻芯的热量。
• 以为仿真时序调到与试验对上就算标定完成，其实参数调到对上但四路归因说不出来，模型换一个方案就不成立；判据是每一路误差源的扰动量级能否解释实测偏差，解释不了就不许把偏差归给它。
• 以为对标要盯着喷发段的峰值温度，其实那个量受测点位置与存活状况支配（测点怎么活到试验结束见 L1-11），对标量应取「第 n 颗被触发的时刻」与「传播是否停止」这两个时序量。
• 以为一次整包试验对上了就能报「误差 X%」，其实一次试验只有一个样本、重复性/再现性框架在此失效（论述见 L1-11）；能报的只有「在这一次试验的条件下模型时序落在哪个区间、哪些条件一变就不成立」。
• 以为传播仿真做得够细就能判「是否满足法规判据」，其实喷发相态与颗粒轨迹标不准（D5-01 已写死），仿真只交方案排序与趋势，法规判据必须由整包试验定案。
• 以为把包级 CHT 跑通了就不必再跑低阶口径（或反过来只跑集总网络就签隔热厚度），其实三口径的分工是「高阶标参数、低阶扫方案、选中方案回高阶复算」，跳掉标定这一步等于把低阶模型的等效热阻当成了物理常数。

需求区 · 想听众筹

想听/点名均为需求登记，不涉及任何付款（意向金仅登记不收款）。满 5 人点名 → 自动向平台专家发邀约。