一维系统模型的数值内核与「跑飞」排查：DAE/代数环、步长与容差、一致初值与守恒自检

课程大纲

• 能判别一个 1D 回路模型的方程属 DAE 还是 ODE，按 ∂g/∂z 是否奇异判指标，并读懂非因果工具「非线性方程组求解失败／指标约简失败」与因果式「代数环」这两类报错各自指向什么
• 能用时间常数比估出模型的刚性比、据此选隐式变步长积分器，并给出分量化的 rtol/atol 取值与「收紧一个数量级复算」的自检结论
• 能构造一致初值（由代数约束反解代数变量），并在稳态求解器起步、瞬态 soak 起步、上一工况热启动三条路线间按判据选一条
• 能定位离散事件、判断何时必须事件重启与步长回缩，并把求解器层的事件抖振与控制策略层的 chattering 分开、分别处置
• 能对模型做质量/能量守恒残差自检与 1D 离散份数无关性检查，并按六层 verification 清单逐层定位跑飞源
• 能说清 verification 与 validation 的分工，判断一个异常结果该进本课清单还是进 J7-01 清单，并据此给出模型是否具备进 J7-01 相关性验证与 J7-06 UQ 的资格

1. 1D 系统方程的数学结构：DAE、指标，与「代数环」的两套表述
   • 1D 热管理回路的方程天然不是纯 ODE：部件热容给出微分方程 ẋ=f(x,z,t)，而不可压回路的节点压力—流量关系不含时间导数、是代数约束 g(x,z,t)=0，两者合成微分代数方程组（DAE）。判据只有一条——模型里有没有「无蓄能元件的节点」；只要有，DAE 已经存在，与用哪家工具无关
   • 指标（index）决定工具能不能直接推进：半显式 DAE 在 ∂g/∂z 非奇异时为指标 1，主流积分器可直接积分；∂g/∂z 奇异（阀全关后的孤立节点、两个理想压力源直接对接、一串无体积节点成环）则指标 ≥2，必须先做指标约简或加正则化。商用工具的通行做法是在节点引入小体积或人工可压缩性，把代数约束变成快速微分方程——这就是「模型明明是不可压液体，工具却要我填一个体积参数」的原因
   • 非因果物理网络工具（Amesim / GT / Dymola / Modelica）里，用户永远看不到「代数环」三个字：编译期做的是指标约简与方程撕裂（tearing），失败时的措辞形如「非线性方程组无法求解／初始化失败／指标约简失败」。按「代数环」去检索工具文档一无所获，正确的检索词是 tearing、index reduction 与 nonlinear solver failure
   • 因果式（Simulink 类信号流）建模里「代数环」才是原生术语：两个都带直通（direct feedthrough）的模块首尾相接即成环，典型现场是控制器—plant 联合仿真的直通反馈。识别靠工具的环诊断报告；打断办法（插单位延迟、改接口因果性、把一侧改成非直通）J1-04 §3 已写，本讲只补机理与代价——插单位延迟等于给该回路引入一个采样周期的相位滞后，强耦合回路上会把稳定性判断带偏
   • 两套表述必须并给，否则半数学员对不上号：同一个数学问题在两种范式下的报错措辞、诊断入口与修法完全不同。范式本身（因果 vs 非因果、四工具选型）见 J1-01 §2/§3，本讲不重讲，只接着讲它在求解层的后果
   • 雅可比与牛顿迭代是两套表述共同的底：非线性代数块每步都要解一次 Newton。求解器统计里「雅可比重算次数」「牛顿迭代不收敛次数」异常升高，是模型病态（状态量量级相差过大、除零、物性表外插）的早期信号，比等它报错早得多，可当日常体检指标读
2. 刚性、积分器与容差：步长由什么定，rtol/atol 该怎么取
   • 刚性要按时间常数比读、不按「难算」读：刚性比 S ≈ τ_max/τ_min。1D 热管理回路里 τ_max 来自冷却液与金属大热容（10²~10³ s 量级），τ_min 来自小体积节点、两相段胞与快执行器（10⁻²~10⁻¹ s 量级），刚性比常在 10³~10⁴ 量级（均为量级估计，非实测值）。判据是：显式积分器的步长被 τ_min 钳住、与你关心的精度无关，命中即必须换隐式
   • 隐式变步长积分器（BDF/Gear 类）是系统仿真的通行选择：它按误差控制自动放大步长，代价是每步要解非线性方程组、要装配雅可比。这解释了两个反直觉现象——模型越大单步越贵；以及「步长明明放大了却更慢」（雅可比重算与牛顿失败的开销吃掉了收益），判别看统计量而不是看 CPU 曲线猜
   • rtol/atol 控的是每步的局部误差，不是最终结果的全局误差：全局误差随积分步数累积，收紧 rtol 只等价于「数值误差不再支配结果」，不等价于「结果更准」。建模误差与物性误差一分不减，那是 validation 的账（J7-01）
   • atol 必须按状态量分量给：温度以 K 计、压力以 Pa 计、质量以 kg 计，同一个标量 atol 对三者意义完全不同；对会过零或量级横跨几个数量级的状态量（切换瞬间的支路流量最典型），判据由 atol 支配而非 rtol——取大了等于放弃该量的精度，取小了等于把数值噪声当误差、步长被压死
   • 容差的工程取法有一条可执行的自检：把 rtol 收紧一个数量级重跑，比对真正要报的结果量（pull-down 时间、峰值温度、COP、能耗）；变化超过判读分辨力说明原设置下结果由数值误差支配，须再收紧后重跑；变化可忽略才允许把该设置写进仿真计划并入库（登记粒度见 J7-07）。这一步属 verification——它证明「结果不随数值设置漂」，不证明「结果对」
   • 事件定位精度是另一个独立容差，不受 rtol/atol 管：过零搜索容差决定阀切换时刻落在哪一步，取松了会把切换时刻挪动一整个宏观时间尺度、把压力冲击削平，机制见第 4 讲
   • 刚性比也是联合仿真宏时间步的判据来源：宏步须按快子系统（制冷剂侧）的 τ_min 取而不是按慢子系统（水路）取，取大了耦合失真。⚠ 宏步之外的耦合参数——耦合迭代次数与松弛因子——不在本课，其判据的全部前置（Jacobi/Gauss-Seidel 调度、显式/隐式耦合、FMI 状态回滚能力）都在 J1-04 §2 自己手里，归 J1-04
   • 「刚性」是同名不同物的重灾区：本讲指方程组的时间尺度分离，与 J5 板块的结构刚度、与 D5-07 三维 CHT 里化学源项的时间尺度分离都不是一回事，结论不可互相搬用（D5-07 已写死该划界）
3. 一致初值与起步：为什么「还没开始算就失败了」
   • 只给微分变量初值不构成完整初值：代数变量 z₀（节点压力、支路流量）必须由 g(x₀,z₀,t₀)=0 反解，各拍一个数就是不一致初值，工具在 t=0 直接报初始化失败或第一步就发散。这正是 J1-01 误区「稳态没收敛就切瞬态」下面那一层，也是 J1-04 误区「各子模型初值不一致」的机理
   • 一致初值与稳态初始化不是一回事：一致初值只要求满足代数约束 g=0，稳态初始化还要求 ẋ=0（即 f=0）。模型本就无定常解时（周期性输入、极限环、间歇热源）稳态求解器永远收敛不了——此时该给的是一致初值 + 一段瞬态 soak 起步，而不是继续加大稳态迭代上限
   • 指标 ≥2 时满足 g=0 仍不够：还须满足对 g 求导得到的隐含约束（工具做指标约简后自动生成的那几行）。现场表现是「初始化通过了，但第一步就把步长压到最小」，判据是看求解器日志里初始步长是否被反复拒绝——这个症状与容差取小完全同形，靠现象猜必然误判
   • 三条起步路线与取舍：① 稳态求解器起步——快、干净，但要求存在定常解且初猜足够近；② 瞬态 soak 起步——从明确的均温静止状态跑到准稳态，稳且可解释，费机时，且必须验证 soak 跑够了（判据是关心量的漂移速率进带，不是跑够某个约定时长）；③ 上一工况末状态热启动——最省，但会把上一工况未收敛的污染带进来，必须先过守恒残差自检（第 5 讲）才可用
   • 初猜要给物理量级而不是零：z 的初猜取泵在工作点附近的流量、回路名义压力，牛顿收敛域比全零初猜大得多。全零初猜在含 Δp ∝ Q² 型阻力的回路里会直接踩到雅可比奇异（Q=0 处该行导数为零），表现为「第一次迭代就报奇异矩阵」
   • 起步侧的边界：拓扑重构时刻（模式切换后）一致初值怎么在热泵换向域落地，属 J1-05 §5 的域专属项；燃料电池 -30℃ 冷启动的强瞬态起步只借本讲的刚性与一致初值两项，冷启动策略本体归 O1-04；soak 的物理侧初始条件（整车热浸是负初始条件、初温怎么给）归 B1-01，本讲只管数学侧的一致性
4. 离散事件：过零检测、事件重启与步长回缩，以及事件抖振
   • 阀切换、泵启停、饱和线穿越都是离散事件，变步长积分器不能「跨过」事件积分：必须给出事件函数（zero-crossing function），积分器检测其符号变化后在步内做二分/插值搜索定位到根，再在事件时刻截断当前步。没写事件函数的切换会被积分器当成一次剧烈非线性去硬啃，表现为步长塌陷而非报错
   • 事件重启是必需动作而非可选优化：事件后系统方程本身改变（拓扑重构、约束增减），积分器的历史阶数与步长必须重置、雅可比必须重装配，否则用旧历史外推会直接发散。J1-05 误区「切换时不重启／不缩步长」讲的就是这件事，本讲给的是机制与判据
   • 步长回缩的代价要量化再决定容忍不容忍：一次事件典型意味着丢掉当前步 + 从最低阶重新起步；事件密集时（带滞环的开关频繁翻转）积分器几乎退化成一阶小步长，墙钟时间可以塌一到两个数量级。判据看求解器统计里的事件计数与拒步数，不看 CPU 曲线猜
   • 事件抖振（数值层）是事件函数在阈值附近反复穿零导致的，模型侧处置是给事件函数加迟滞带或最小驻留时间。迟滞带宽的取法判据是双向的：宽到能盖住该信号在事件附近的数值噪声与容差带，窄到不改变关心的物理结论——切换时刻的偏移量必须小于判读分辨力，否则等于用数值手段改了物理答案
   • ⚠ 同形异义必须当堂划清：K3-01 的 chattering 是控制策略层的执行器抖动（成因是阈值单一、无迟滞），处置在控制器里加控制侧迟滞与延时；本讲的事件抖振是求解器过零检测层的现象，处置在模型侧。两者可同时存在，判别办法是把控制器换成开环固定指令重跑——还抖就是数值层的
   • 事件与物性的交互：跨饱和线是隐式事件（物性导数不连续），不当事件处理时的表现是步长被反复压到最小却不报错。两相/物性域的专属处置（关联式选择、两相段胞数、跨饱和线插值）归 J1-02 §5，本讲只给「把它当事件而不是当噪声」这一条判据
   • 事件与联合仿真宏步的关系：事件落在宏步内部时，若从属子系统不支持状态回滚（FMI 的 canGetAndSetFMUstate），事件只能被推到宏步边界、切换时刻带一个宏步的误差；协同机制本体见 J1-04 §2，本讲只给这一条后果
5. verification：守恒残差自检、1D 离散份数无关性与六层排查清单
   • 课首先划界：本讲面向 verification（方程解得对不对、结果是不是数值伪迹），J7-01 面向 validation（模型准不准、与试验对不对）。两份清单入口不同——算不出来／振荡／不守恒进本讲六层清单（单位→初值→边界→物性→离散→耦合），算得出来但与试验差进 J7-01 §5（试验数据有效性→边界条件→物性→数值→模型形式）。异常结果先判属哪一半再进清单。⚠ 措辞纪律：ASME V&V 20 正文重心在 validation 与不确定度量化，solution/code verification 主要作前置条件出现，故本讲清单是工程可操作化产物，不声称是该标准的对等半边
   • 质量守恒残差自检：对每个控制体与整个闭环算 R_m = |Σṁ_in − Σṁ_out − dM/dt| 并按特征流量归一化。判据不是「残差为零」（变步长积分器本就允许容差内误差），而是「归一化残差比关心量的分辨力小一个数量级以上」，且瞬态尖峰处短时抬高可接受、持续单调漂移不可接受。闭环充注量单调漂移而模型无泄漏路径，即为数值不守恒而非物理
   • 能量守恒残差自检：同法算 R_E 并按回路特征换热功率归一化。它抓「热量凭空产生或消失」，常见根因三条——换热器两侧换热量未严格互为反号、边界热流用了不同时刻的量、热容项与流项取了不同物性状态点。归一化基准落在零附近（零流量段、无热负荷段）时该判据失效，须改用绝对量判
   • 1D 的离散份数无关性：通用方法论（加密序列、Richardson 外推与 GCI、假收敛识别）引 A1-01 讲义 §3.4 与 J2-01 §5，本讲不重讲，只讲 1D 对象化的三件——管段离散份数（判据看温度锋面与输运时滞是否被抹平）、换热器胞数（判据看换热量与出口温度是否随胞数进带）、两相段胞数（最敏感，与 J1-02 误区「两相段只落一两个胞」是同一件事）。交付形式是一条「关心量—份数」曲线，不是一句「已加密」
   • 六层排查清单按代价从低到高逐层过：① 单位与量纲（最便宜、最常中，尤其跨工具接口的 bar/Pa、kW/W、℃/K）；② 初值一致性（第 3 讲）；③ 边界条件（给没给到位、是否互相矛盾、阀全关是否造出孤立节点）；④ 物性（调用是否外插、两相区导数是否连续）；⑤ 离散与容差（份数无关性、rtol/atol 收紧复算）；⑥ 耦合（宏步、接口方向、事件是否落在宏步内）。纪律与 J2-01 一致——一次只动一层，否则因果分不开
   • 排查前先做一次「最小化复现」：把控制器换成开环固定指令、工况砍成单点、子系统逐个换成理想边界，直到跑飞消失，最后被换掉的那块就是嫌疑对象。这比在整模型上调容差高效得多，也是把控制层 chattering 与数值层事件抖振分开的标准手段（第 4 讲）
   • 通过 verification 是下游三件事的准入条件而不是终点：进 J7-01 做相关性验证、进 J7-06 做 UQ（该课以「模型已能稳定求解」为前提，拿会跑飞的模型做 UQ 得到的只是数值噪声的分布）、进 J7-07 把积分器/步长/容差/初始化方式作为快照字段入库（本课只交付「凭什么这么取」，登记粒度归 J7-07）

• 以为「1D 没有网格所以没有离散误差」——其实管段份数、换热器胞数、两相段胞数就是 1D 的网格，份数不足会把温度锋面与输运时滞抹平，方法论与 J2-01 的网格无关性同源，只是对象换了
• 以为 rtol 收得越紧结果越准——其实 rtol/atol 只控每步局部误差，收紧只能让数值误差不再支配结果；建模误差与物性误差一分不减，那是 validation 的账（J7-01）
• 以为工具没报「代数环」就没有代数环问题——其实非因果工具（Amesim/GT/Dymola）在编译期做指标约简与撕裂，失败时的措辞是「非线性方程组无法求解／指标约简失败」，按「代数环」去检索工具文档一无所获
• 以为给了各部件温度初值就算给了初值——其实节点压力与支路流量是代数变量，必须由 g=0 反解成一致初值，各拍一个数就是起步即不可解
• 以为稳态求解器不收敛就说明模型建错了——其实稳态求解要求 ẋ=0，模型本就无定常解时（周期输入、极限环、间歇热源）它永远收敛不了；此时该改走一致初值 + 瞬态 soak 起步，而不是加大迭代上限
• 以为切换处抖动一定是控制策略的锅——其实求解器过零检测层的事件抖振与控制层 chattering 现象几乎同形，判别办法是把控制器换成开环固定指令重跑，还抖就是数值层的（控制侧归 K3-01）
• 以为守恒残差不为零就是模型错——其实变步长积分器本就允许容差内的残差，判据是归一化残差相对关心量分辨力的量级、以及是否持续单调漂移，不是「等于零」
• 以为跑通了就可以交——其实跑通只过了 verification 的一半：还须过守恒自检与份数无关性，才有资格进 J7-01 的相关性验证与 J7-06 的 UQ；J7-06 已明写以「模型已能稳定求解」为前提，拿会跑飞的模型做 UQ 得到的只是数值噪声的分布
• 以为本课的「刚性」可以搬到别处用——其实它指方程组的时间尺度分离，与 J5 板块的结构刚度、与 D5-07 三维 CHT 里化学源项的时间尺度分离同名不同物，D5-07 已写死该划界
• 以为跑飞了就该先调容差——其实容差是六层清单里的第 ⑤ 层，前面还有单位量纲、初值一致性、边界条件、物性四层更便宜也更常中；从最贵的一层开始调是最常见的时间浪费
• 以为一个标量 atol 可以全模型通用——其实 atol 必须按分量给：温度以 K 计、压力以 Pa 计、质量以 kg 计，同一数值对三者意义完全不同，且对会过零的流量类状态量是 atol 而非 rtol 在支配判据
• 以为「步长塌到最小」只有一种根因——其实指标 ≥2、事件抖振、容差取小、物性表外插四者的现场表现完全同形，只能靠六层清单逐层过、一次只动一层来分开，靠现象猜根因必然反复返工

需求区 · 想听众筹

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