https://jaywillow0.github.io/tags/%E7%87%83%E6%96%99%E7%94%B5%E6%B1%A0/Recent content in 燃料电池 on LiuYang LabHugozh-CN© 2026 LiuYangFri, 09 Oct 2026 00:00:00 +0800https://jaywillow0.github.io/posts/pemfc-cold-start-part2-symbols-units/Wed, 07 Oct 2026 22:30:00 +0800https://jaywillow0.github.io/posts/pemfc-cold-start-part2-symbols-units/<p>上一篇结尾把 Part2 的清单按优先级列了八条，前两条就是建符号转换表和全模型换 SI 摩尔基准。真到写实施规格的时候，这两条占掉的篇幅比预想多得多，干脆从原计划的第二篇里分出来单独成篇。原计划第二篇直接写脚本框架和离散化，动笔才发现那部分内容每一步都踩在符号和单位上，先把地基打实，网格留到下一篇。</p> <p>分出来的理由不体面也不复杂。逐式核对时留下的那些坑，量纲冲突、单位混用、正负号不明，没有一样能靠求解器兜住。$EW$ 印作 1100 $\mathrm{kg/kmol}$，配上 $F=96487\ \mathrm{C/mol}$ 直接代进去，误差一千倍，程序不报错，曲线照画。符号反向更隐蔽，曲线形状大体还在，物理早就反了。</p> <ol class="solitude-tag series-items tag-series"><li><a href="https://jaywillow0.github.io/posts/pemfc-cold-start-part1-mechanisms/" title="燃料电池冷启动复现札记（一）· 解构 Huo 2019 的水热机理与公式">燃料电池冷启动复现札记（一）· 解构 Huo 2019 的水热机理与公式</a></li><li><a href="https://jaywillow0.github.io/posts/pemfc-cold-start-part2-symbols-units/" title="燃料电池冷启动复现札记（二）· 符号转换表与 SI 摩尔基准">燃料电池冷启动复现札记（二）· 符号转换表与 SI 摩尔基准</a></li></ol> <h2 id="一三套正负号一张转换表">一、三套正负号，一张转换表</h2> <p>上一篇第二节记过这笔账。式(1)把 Nernst 电压、活化、欧姆三项当带符号的加数，式(9)又明确 $V_{\rm ohm}\le0$，到了热源那边，式(35)要求耗散非负。三个地方三套隐含约定，单看每处都能自圆，合在一起写程序就要出事。式(7)、(8)把两极过电位写成相反符号的相电势差，还没说明是否扣掉了平衡电位参考。这种账，光靠小心防不住，程序里换一套从第一步就自洽的约定，再建表把原文逐条接进来。</p> <p>程序内部只认一套。定义体积法拉第电流 $j_F$，阳极氧化为正，阴极还原为负，电子与质子电流沿全局 $+x$ 方向为正</p> $$ \partial_x i_s=-j_F,\qquad \partial_x i_H=j_F, $$$$ i_s=-\sigma_e^{\rm eff}\partial_x\phi_s,\qquad i_H=-\sigma_H^{\rm eff}\partial_x\phi_m. $$<p>每个 CL 控制体的两个电流散度使用同一个 $j_F$，两式相加直接得到 $i_s+i_H=I$ 逐点成立。恒流加载求解完，回头检查阳极端电流是否等于给定值，这是不花钱的正确性检查。原文式(7)、(8)里那些容易混淆的正负号，从此不再进入程序。</p> <p>过电位只认常规定义</p> $$ \eta_e=\phi_s-\phi_m-U_e. $$<p>配合阳极 $j_F\ge0$、阴极 $j_F\le0$，乘积 $j_F\eta_e$ 天然非负，活化热 $Q_{\rm act}=j_F\eta_e$ 直接可写。端电压取两端电子电势之差，端面值由半单元通量外推，CL 活化损失和各层欧姆损失都已经含在电势解里，不能再减一次 $IR$。各层损失照样可以输出，但只做诊断分解，不进电压公式。</p> <p>反应若出现明显反向，程序报告超出冷启动单向反应协议，不取绝对值硬跑。指数项用双曲正弦形式计算，大过电位下的指数溢出触发代数失败重试，不被当成物理停机。</p> <p>转换表本身不长，贴在这里。</p> <table> <thead> <tr> <th>原文表达</th> <th>程序内约定</th> </tr> </thead> <tbody> <tr> <td>式(1) 三项带符号相加</td> <td>端电压取两端电子电势差，损失只做诊断</td> </tr> <tr> <td>式(3)、(4) $J_i$ 分别作源</td> <td>同一个 $j_F$ 驱动两相散度，逐点闭合</td> </tr> <tr> <td>式(7)、(8) 相电势差过电位</td> <td>$\eta=\phi_s-\phi_m-U_e$，$U_e$ 有明确参考</td> </tr> <tr> <td>式(9) $V_{\rm ohm}\le0$</td> <td>耗散 $\sigma\|\nabla\phi\|^2$ 天然非负，不进电压</td> </tr> </tbody> </table> <p>这张表往后比正文更常用。任何时候和原文对不上，先查表，别凭记忆翻译符号。</p>https://jaywillow0.github.io/posts/pemfc-cold-start-part1-mechanisms/Mon, 28 Sep 2026 15:30:00 +0800https://jaywillow0.github.io/posts/pemfc-cold-start-part1-mechanisms/<p>读博那几年，我的研究方向是质子交换膜燃料电池冷启动的跨尺度建模，从催化层微结构一直做到单电池的温度场。毕业进了储能行业，这套手艺日常用在磷酸铁锂电池上，燃料电池的线偶尔也会想起，读博时的一些想法也没机会想完。现在纯粹出于兴趣，想把这些补上，做法是挑一篇当年给过我较多启发的文献，用 Python 从零把它的模型实现一遍，当年怎么算的，这回一行一行重新走。</p> <p>就是 Huo、Jiao 和 Park 2019 年发表在 Applied Energy 上的那篇冷启动模型论文（<a href="https://doi.org/10.1016/j.apenergy.2018.10.068">原文</a>）。读博搭自己模型的时候，这篇文章在机理组织和参数选择上给过我不少启发，这次换了身份重读，感受完全不同。当读者和当复现者要的东西不一样，读者要读懂主线，复现者要让每个符号都能落进代码，量纲、单位、参考态、分支条件，一处含糊后面都是坑。</p> <p>先把总体判断放在这里。这篇文章的模型结构足以支撑一次严肃的复现，五种水状态、六类相变通道、电与热的双向耦合都摆在明面上；参数表和符号体系却不能直接抄进代码。我逐页核对了全部 36 个编号公式和 5 张参数表，量纲冲突、单位混用、数值排版问题都能找到，还有些闭合关系原文没有给。这篇札记就是这份对账单，公式编号与原文一致，方便对照。</p> <p>按计划写六篇。本篇拆文献的水热机理和公式；第二篇先立符号与单位的规矩，再往后是离散骨架、物性与加载的坑、脚本框架与对照表；最后一篇记录调参和结果对标。</p> <ol class="solitude-tag series-items tag-series"><li><a href="https://jaywillow0.github.io/posts/pemfc-cold-start-part1-mechanisms/" title="燃料电池冷启动复现札记（一）· 解构 Huo 2019 的水热机理与公式">燃料电池冷启动复现札记（一）· 解构 Huo 2019 的水热机理与公式</a></li><li><a href="https://jaywillow0.github.io/posts/pemfc-cold-start-part2-symbols-units/" title="燃料电池冷启动复现札记（二）· 符号转换表与 SI 摩尔基准">燃料电池冷启动复现札记（二）· 符号转换表与 SI 摩尔基准</a></li></ol> <h2 id="一一块电池和五种水">一、一块电池和五种水</h2> <p>质子交换膜燃料电池低温启动的核心矛盾很清楚。零下环境里电池发生氧化还原反应产水，水在离聚物、孔隙和冰之间反复分配，冰占住孔隙，氧气就到不了反应位点，电压跟着垮。启动成败取决于产水结冰和发热融冰谁跑得快。Huo 这篇的价值在于把过冷液水也纳入模型，用五种水状态把这条因果链完整摆开。</p> <p>模型沿膜电极组件厚度方向建一维，完整层序是</p> $$ \mathrm{BP_a\,|\,Channel_a\,|\,GDL_a\,|\,MPL_a\,|\,CL_a\,|\,Membrane\,|\,CL_c\,|\,MPL_c\,|\,GDL_c\,|\,Channel_c\,|\,BP_c} $$<p>下标 a、c 表示阳极和阴极，BP 是双极板，GDL 是气体扩散层，MPL 是微孔层，CL 是催化层。层序明确不等于边界闭合，原文没有给出每层网格数，也没有讲清流道在一维传热传质里如何等效，复现时要自己决定流道是规定值边界还是独立控制体。</p> <table> <thead> <tr> <th>水状态</th> <th>符号</th> <th>主要位置</th> <th>迁移方式</th> </tr> </thead> <tbody> <tr> <td>未冻膜水</td> <td>$\lambda_{\rm nfm}$</td> <td>膜与两侧 CL 的离聚物</td> <td>扩散与电渗拖曳</td> </tr> <tr> <td>冻结膜水</td> <td>$\lambda_{\rm fm}$</td> <td>膜内</td> <td>就地冻结与融化</td> </tr> <tr> <td>水蒸气</td> <td>$C_{\rm vap}$</td> <td>CL、MPL、GDL 孔隙</td> <td>气相扩散</td> </tr> <tr> <td>过冷液水</td> <td>$s_{\rm lq}$</td> <td>CL、MPL、GDL 孔隙</td> <td>压力驱动流动</td> </tr> <tr> <td>孔隙冰</td> <td>$s_{\rm ice}$</td> <td>CL、MPL、GDL 孔隙</td> <td>就地累积</td> </tr> </tbody> </table> <p>膜水这个名字容易误导。催化层里的离聚物同样装着膜水，名字说的是介质，不能据此只在膜层设变量。冻结膜水和孔隙冰是两个库存，前者影响还能导质子的未冻水量，后者直接占用气体通道，混在一起会把两条物理反馈搅乱。</p>