燃料电池系统一维建模与水/热耦合仿真

课程大纲

• 能按「这个模型要回答什么」定电堆子模型的降阶粒度（整堆集总 / 沿叠层多节点 / 必须转三维），并逐档说清它给不了什么
• 能把极化曲线与产热源项 Q_stack=N_cell·(E_tn−V_cell)·I 正确落进 1D 模型，并在交付物上把 HHV/LHV 基准写死
• 能约定电堆与冷却回路的四个接口量，并把阴极尾气带走的显热＋潜热与 Q_stack 分账，避免把回路热负荷算成上界值
• 能按「是否进入输出量的因果链」「动态是否与热同尺度」两条判据，决定空压机 / 中冷 / 增湿器 / 氢路各降到第几档，以及何时可整体扔成边界条件
• 能在求解器里把水平衡闭合，并逐时间步做水量与能量两条守恒残差自检，按漏项清单定位残差不回零的原因
• 能表达自加热注入源项与冰相变源项，并说明相变温区宽度这个数值参数把冷启动结论推向哪个方向
• 能选出 FC 专属标定量与工况矩阵，把冷启动成功率算成模型侧通过率，并与环境舱样本按分布（而非单点）对标
• 能写出模型交付页的六项声明与可信域，使引用方知道这份模型在什么区间外必须降级使用

1. 电堆准稳态电化学—热子模型：降阶粒度怎么定、参数从哪来
   • 粒度不是越细越好，先问模型要回答什么：只算回路热负荷与出堆温度→整堆单热节点够；要评片间温差、最弱片电压→必须沿叠层方向分节点；要评片内脊-沟梯度、端板边缘冷点→一维根本给不了，必须转三维（判据同第 2 讲，固体侧本体见 O1-06，三维方法见 J2-01 / J2-04）
   • 产热源项一律按 Q_stack = N_cell·(E_tn − V_cell)·I 落在电堆热节点上，不许按电功率乘系数折算（式与两笔账的拆分承 O1-01）；⚠ 引 E_tn 前必须声明 HHV / LHV 基准，两套口径下同一台堆的产热差三成量级，换算的唯一归属是 O1-08
   • 极化曲线入模三种做法与各自的失效边界：查表插值最稳但外推区无效；分段半经验拟合（活化段 + 欧姆线性段 + 传质塌陷项）参数有物理含义、可随温度/压力/λ 修正，是全工况模型的默认选择；纯多项式拟合禁止用于外推。曲线一律是 BOL 口径，EOL 抬升另作一档（口径承 O1-07，本课不代讲）
   • 电化学取准稳态而非全动态是本课的默认假设：双电层与膜含水量的时间常数比电堆热时间常数小若干数量级，在热耦合尺度上可把 i-V 当代数关系用。⚠ 一旦输出量变成「电压跌落」本身（冷启动初期、秒级强变载），该假设失效，处置见第 5 讲
   • 热容必须分列而不是并成一个数：MEA + 双极板的片间热容与端板 / 集流板 / 拉杆的热容分开给——端板是升温最慢的一段（端部偏冷的机理承 O1-01、O1-08），并进整堆平均热容会把冷启动时间算短、成功率算高，误差方向偏危险
   • 参数来源与「拿不到数据时怎么办」：i-V 曲线、等效热容、堆侧 UA 取自电堆供应商规格与台架实测；⚠ 拿不到逐片数据时只能建整堆单节点，并在交付页写明「本模型不评片间温差与最低片」——不许用假设的分布凑出节点，那会让第 5 讲的成功率变成一个看不出错的假数
   • 本讲不重讲：三段极化的机理与曲线读法（前置 O1-08）、可逆熵热与不可逆极化热的拆分及膜内水传递三路径（前置 O1-01）、催化层—GDL—双极板的固体侧热阻链与冷却腔构型（O1-06）——本讲只把它们的结果降成 1D 模型里的节点、源项与参数
2. 冷却回路与电堆的水力—热耦合接口约定
   • FC 回路有两条量级特征直接改写建模取舍：散热量与电功率同量级、可用温差只有内燃机的一半（承 O1-01、O3-01），于是回路走的是「大流量、小温差」。小温差意味着出堆温度对流量的相对误差极敏感——流量测不准的误差几乎全额落在进出口温差上，这是 FC 模型标定与 BEV 最大的手感差异
   • 接口量必须一次约定死并随模型交付，共四个：进堆冷却液温度、冷却液流量 ↔ 出堆冷却液温度、进入冷却液的放热量 Q_cool。⚠ Q_cool ≠ Q_stack——HHV 口径的 Q_stack 隐含「生成水全部在堆内冷凝」，实际大量水以蒸汽随阴极尾气带走显热＋潜热，Q_stack 只是进入冷却液热量的上界（该上界性质承 O1-01）；模型里必须显式建一条尾气带热支路，否则回路热负荷与散热器面积被系统性算大，并逼着标定去改 UA 补偿
   • 回路水力本体不重讲：管段沿程与局部阻力的 Δp-Q、泵曲线与系统阻力曲线求工作点、并联支路流量分配、多通阀切换的离散事件、空气侧阻力链与气—液闭环迭代，全部承 J1-03；本讲只补 FC 专属的三件——DI 旁通支路、电堆低阻大流量特性、板换 / 散热器与电堆的边界交换量定义
   • DI（去离子）支路在 1D 里的接法与简化判据：它是并联的小流量旁通，水力上只占百分之几量级、热上可忽略，做热与水力仿真时可简化成一个固定分流比的旁通阻力；但它决定的电导率是不可标量的硬边界（限值来源与在线管理承 O3-03、B2-04）——一旦要评电导率随运行时间的爬升或树脂寿命，就必须加一条电导率状态方程，那属 O3-03 的口径，本课只留接口不建模型
   • 电堆在水力网络里怎么等效：几百片并联流道整体等效为一条 Δp-Q 特性 + 一个（或多个）热节点。⚠ 该特性只能由供应商数据或整堆实测给——用单流道 CFD 结果乘片数外推会漏掉进出口集管的分配不均，而分配不均恰恰是片间温差的主因（本体见 O1-03）
   • 「何时必须上三维」的判据（本课只给判据、不讲三维方法）：输出量里出现片间流量分配、脊-沟局部温度、端板边缘冷点、局部液态水堵塞位置这四类中的任何一类，一维即失效；网格与湍流方法本体见 J2-01，固体—流体共轭传热见 J2-04，电堆内部热阻链见 O1-06
3. 空气路与氢路子模型的粒度取舍：什么时候可以整体扔成边界条件
   • 统一取舍判据（本讲的核心，三档）：某部件若不进入本次仿真输出量的因果链，降为边界条件；进入因果链但其动态比电堆热时间常数快一个数量级以上，取准稳态 map；既进因果链、动态又与热同尺度的（中冷、增湿器、堆入口温湿度），必须建动态。判据要按「这一次的输出量」重新过一遍，不是按部件重要性一次定死
   • 空压机取准稳态 map 查表（流量—压比—效率—转速），出口温度由压比与等熵效率折算，它同时是中冷段的进口边界。⚠ map 的实测口径、喘振线与修正折算见 L1-12（本课不代讲测法）；部件冷却本身归 O2-01
   • 中冷段必须建、不能省：它直接定堆入口空气温度，而堆入口温度是水平衡的起点（饱和蒸气压随温度陡升，是水热耦合的物理根源，承 O1-01）。1D 里降成 ε-NTU 单元 + 空气侧压降即可；⚠ 空气侧风量与进风温度不许当外部常数灌进来，接法承 J1-03 与 A2-05 的阻力链—风扇工作点求交
   • 增湿器降阶分三档，按输出量选：① 固定有效度 ε_m——只要堆入口 RH 的量级；② ε_m 随工况修正的查表——要评全工况水平衡；③ 双通道分布式传质——要评沿程结露与低温冻堵位置。ε_m 的定义、选型与实测口径承 O2-02、L1-12，本讲只定取哪一档、每档失效在哪；⚠ 三档共同的失效边界是把常温常压标定的 ε_m 直接外推到高原与低温（该误区承 O2-02）
   • 氢路的简化最激进，也最容易简化过头：稳态热仿真里氢路常可整体降成「按法拉第定律算出的氢耗 + 一个入堆温度边界」，因为氢的热容流率比空气与冷却液小一个量级。⚠ 三种情况下这条简化立即失效——要算引射器能否卷吸够（λ_H2 的系统级口径 ≈1 与堆入口口径 1.2~2 不可混用，两个口径的分野承 O1-08，换算与选型归 O2-03）、要算 purge 引起的压力与流量脉冲、要算低温下阳极再循环环路的积水与冻堵（承 O2-02、O1-04）
   • 排放 / purge 阀在 1D 里是离散事件，它同时改写阳极压力、氮与水的积累和水平衡账本：建模上要么按周期性边界注入、要么按积累量触发，两种做法对水量结果影响不同，取哪一种必须写进交付页。⚠ 离散事件的过零检测与步长回缩属通用数值机制，前指 J7-08（本课不代讲）
   • 本讲不重讲部件热设计本体：空压机冷却与进气中冷承 O2-01、膜增湿器承 O2-02、氢循环泵 / 引射器与供氢温控承 O2-03——本讲只回答「在系统级 1D 里它降到第几档、什么时候能整体扔成边界条件」
4. 水平衡在一维求解器里的闭合与守恒校核
   • 守恒的状态量一律用湿度比 ω 或水的质量流，不用 RH：RH 是温度与压力的函数，压力/温度一变就失真，拿它做质量平衡必错（该误区承 O2-02）；RH 只出现在判据侧与输出侧（干膜 / 水淹的判断口径承 O1-01）
   • 膜水通量三项在 1D 里的表达与最小节点要求：电渗拖曳 J_drag = n_d·(I/F) 与电流同步、可作代数源项；反扩散需要膜两侧含水量之差，因此至少要给阳极、阴极两个水节点，只建一个水节点等于把模型钉死在「生成水全在阴极」这个错误结论上；生成水由法拉第定律直接从电流给（三项本体承 O1-01、O2-02）。⚠ n_d 随膜含水量与电流密度变，全工况取常数是明确误区（承 O1-01）
   • 气液两相在系统级的标准处理是平衡假设：按局部温度与压力算饱和含湿量，超饱和部分即刻析出为液水。它的失效边界是过饱和与液滴携带——低温段、高电流密度段、快速降温段会明显偏离，此时析水量被算多、堵塞风险被算早（方向偏保守但会误导控制策略）；要更细就必须上三维（判据同第 2 讲）
   • 守恒自检是本课的硬交付、不是可选项，每步核三本账：① 水的质量——生成 + 进气带入 + 增湿转移 = 排出 + 液态积存变化 + 膜储水变化；② 能量——Q_stack = Q_cool + 尾气显热 + 尾气潜热 + 蓄热变化 + 对环境散热；③ 元素 / 摩尔——H 与 O 的进出与电流一致。判据是残差不随时间发散、且量级远小于最小的真实项；残差持续爬升就是模型有漏项，不能当「数值噪声」放过
   • 漏项就那几处，按这个顺序查：进气侧起点湿度被当成干空气（环境湿度这一路在高湿工况可与生成水同量级，承 O2-02）→ 增湿器只记了传向干侧的水、没记湿侧被带走的水（必须双侧同时守恒）→ purge / 排放阀那一路没记账 → 膜储水这个状态被省掉，导致动态工况下水量凭空出现或消失
   • 通用数值内核不重讲，全部前指 J7-08：DAE 与代数环、rtol/atol 的工程取法、一致初值、步长与刚性积分器选择、离散事件定位。⚠ 本课不得代讲，只保留占位前指。本讲只保留 FC 水域专属的数值项——低温下含湿量趋零导致的分母病态、饱和判据在相变切换点的不连续与迟滞处理、以及 λ 与含水量互相反馈形成的代数环该在哪一侧打断
   • 与部件课的分界线：水通量表达式、ε_m 的定义与选型、水平衡的物理判据全部承 O1-01 与 O2-02（O2-02 已声明「本课给出的水通量表达式即该模型的输入」）；本讲只负责把它们接成一个在求解器里能闭合、能逐步自检的账本
5. FC 专属强瞬态与冷启动成功率：源项表达、粒度失效点与环境舱分布对标
   • 强瞬态的三个 FC 专属机制：① 自加热注入——物理与策略承 O1-04，模型侧只表达为「V_cell 被限压 / 限流曲线钳住后 (E_tn−V_cell)·I 抬升」的源项；② 冰相变——残留水在 0℃ 的潜热平台，用焓法或等效热容法表达；③ 低温段可解性——含湿量趋零与膜电导率极低带来的病态，处置见第 4 讲
   • 变载瞬态里「谁跟得上谁」决定粒度：电流可在毫秒级跳变，空压机与增湿器是百毫秒到秒级，冷却回路热是几十秒到分钟级。所以强变载下电压跌落的主因是气路响应滞后（λ_air 瞬时跌破目标、传质段提前塌陷，承 O1-08 第 3 讲），不是热。判据：输出量是温度就可把电化学取准稳态；输出量一旦是电压跌落，第 1 讲的准稳态假设必须放弃
   • 冷启动成功率在模型上怎么算：它是统计量、不是单次结果——把 O1-04 的失败判据（超时、单片电压跌落触发保护）原样搬进模型作判据，对残留水量、初始温度分布、吹扫质量三个不确定输入抽样，通过率即模型侧成功率。⚠ 样本量、置信与启动耗时分位数的统计口径本课不另立，承 O1-05 与 O3-07 第 5 讲；输入散差怎么给分布、怎么向输出传播属 UQ 方法本体，见 J7-06
   • 与环境舱数据的对标，对的是分布不是单点：模型给出的成功率与耗时分布，要与环境舱同一「环温 × 冷存时长 × 吹扫参数」格子里的实测样本比。⚠ 台架与环境舱不可互换——台架给的是可复现单回路边界下的热标定量，只有环境舱才有 −30℃ 的真实初始状态分布（归层判据承 O3-07 第 5 讲）；试验规程、DV/PV 编排与统计报告口径归 L3-06（本课不代讲）
6. FC 专属标定量、工况矩阵与模型交付页：可信域怎么写、下游怎么接
   • FC 专属标定量只有四类，其余全走通用相关性方法学（承 J7-01：验证计划与覆盖矩阵、标定集与独立验证集分开、误差口径按标度类型选、E=S−D 与 ±u_val 的读法）：① 单片电压——判据必须落在最低片与平均片的偏差上，总电压会掩盖单片失效（承 O1-08 第 4 讲）；② 堆进出口温差——ΔT 本就小，用相对误差必须声明基准是 ΔT 而非绝对温度；③ 堆入口 / 出口 RH 或露点；④ 空气路与冷却路压降
   • 工况矩阵至少铺四维，缺一维就会出现「标好的模型换个工况就崩」：负载点 × 环温 × λ_air × 增湿档，BOL / EOL 另作一档。⚠ 负载点必须含低负荷高电压点——可逆熵热在低负荷段占产热的比例可达四成以上（承 O1-01），只标额定点会掩盖产热模型的系统性偏差；EOL 口径承 O1-07
   • 模型交付页必须带六项声明，缺一项引用方就会误用：① 所用产热 / 效率基准（HHV 还是 LHV）；② 电堆节点数与「能不能评片间温差与最低片」的明示；③ 气路各部件所取的降阶档；④ 水平衡是否含液态积存与膜储水；⑤ 标定所用工况矩阵与独立验证点清单；⑥ 可信域与已知失效工况（集管分配不均、过饱和携液、purge 密集段、超出标定 λ_air 区间）。超区间使用一律降级为趋势判断
   • 本课的出口：这份子模型接进整车联合仿真的接口与 FC 支路接法归 J1-04（该课已声明「子模型本身的建模粒度与标定不在本课」，两侧口径就此对齐）；实时化与控制导向降阶见 J7-02；相关性方法学本体见 J7-01；系统级模式管理与功率—热仲裁见 O3-07；前端散热能力核算见 O3-01

• 以为把电功率乘一个系数就能当电堆热源，其实必须按 N_cell·(E_tn−V_cell)·I 算，且 E_tn 要声明 HHV / LHV——基准混用一处就差三成量级（承 O1-01、O1-08）
• 以为按 HHV 算出的 Q_stack 就是冷却回路负荷，其实它只是上界：大量热随阴极尾气以显热＋潜热排走；不建尾气支路会把散热器选大，还会逼着标定去改 UA 做补偿
• 以为整堆一个热节点足够用，其实只要输出量里出现「片间温差」或「最低片电压」，单节点就给不出来——此时算出的温差结论与冷启动成功率是假数，须沿叠层分节点或在交付页明示不评
• 以为水平衡用 RH 记账更直观，其实 RH 随温度与压力变，做质量守恒必失真；守恒一律用湿度比 ω 或水的质量流，RH 只出现在判据侧（承 O2-02）
• 以为电渗拖曳系数 n_d 取个常数够用，其实它随膜含水量与电流密度变化，全工况常数会让阳极侧水量系统性算错（承 O1-01）
• 以为只建一个水节点就能表达水平衡，其实反扩散需要膜两侧含水量之差，单水节点等于把模型钉死在「生成水全在阴极」上，低电流与厚膜工况必偏
• 以为水量残差小就说明模型对，其实残差只证明账本闭合，不证明 n_d、ε_m 与析水假设取对——那要靠独立验证点（承 J7-01）
• 以为冰相变温区宽度 ΔT_mush 是可以调着去对实测的参数，其实它是数值参数：调宽了冷启动时间算短、成功率算高，方向偏危险；正确动作是敏感性扫描并写进交付页
• 以为模型能跑通一次冷启动就算能预测成功率，其实成功率是统计量，必须对残留水量 / 初温分布 / 吹扫质量抽样后按通过率给，并与环境舱同格子的样本比分布（统计口径承 O1-05、O3-07，规程归 L3-06）
• 以为氢路永远可以简化成一个边界条件，其实一旦要算引射器卷吸比、purge 脉冲或低温阳极积水，这条简化立即失效；且 λ_H2 的系统级口径 ≈1 与堆入口口径 1.2~2 不可混用（承 O1-08、O2-03）
• 以为空气侧风量与进风温度是可以拍常数的外部输入，其实它由阻力链与风扇工作点求交解出、且受层叠串扰影响（接法承 J1-03、A2-05），拍常数会让中冷与堆入口温度整条链失真
• 以为在台架上标好的门限可以直接用于冷启动结论，其实台架给的是可复现单回路边界下的热标定量，−30℃ 的真实初始状态分布只有环境舱有（归层判据承 O3-07 第 5 讲）

需求区 · 想听众筹

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