[论文解读] A fast solver for the pseudo-two-dimensional model of lithium-ion batteries

Rachel Han, Colin B. Macdonald+1 more|arXiv (Cornell University)|Nov 17, 2021

Advanced Battery Technologies Research参考文献 20被引用 7

一句话总结

本文提出了一种基于有限差分法、后向欧拉时间积分和牛顿法的快速、高精度求解器，用于锂离子电池的伪二维（P2D）模型，采用Python实现。通过重构线性代数结构并利用JAX的自动微分技术，该方法在不损失精度的前提下显著降低了计算成本，实现了对复杂电化学行为的高效模拟。

ABSTRACT

The pseudo-two-dimensional (P2D) model is a complex mathematical model that can capture the electrochemical processes in Li-ion batteries. However, the model also brings a heavy computational burden. Many simplifications to the model have been introduced in the literature to reduce the complexity. We present a method for fast computation of the P2D model which can be used when simplifications are not accurate enough. By rearranging the calculations, we reduce the complexity of the linear algebra problem. We also employ automatic differentiation, using an open source package JAX for robustness, while also allowing easy implementation of changes to coefficient expressions. The method alleviates the computational bottleneck in P2D models without compromising accuracy.

研究动机与目标

• 解决锂离子电池全伪二维（P2D）电化学模型计算成本过高的问题。
• 克服单颗粒模型等简化模型的局限性，后者在某些应用场景中精度不足。
• 开发一种数值高效的求解器，在保持全P2D模型保真度的同时降低计算负担。
• 实现在电池管理系统（BMS）中对P2D模型的实际部署，支持实时控制与状态估计。
• 通过自动微分实现模型系数的便捷修改，提升模型的灵活性与鲁棒性。

提出的方法

• 在x方向（电极厚度）和径向r方向（粒子扩散）的一维空间域上，采用有限差分法（FDM）对P2D模型进行空间离散化。
• 采用后向欧拉时间积分方法，确保微分-代数系统时间推进的无条件稳定性。
• 在每个时间步使用牛顿法求解得到的非线性代数系统，确保快速收敛。
• 利用JAX实现自动微分，精确高效地计算雅可比矩阵和敏感性，支持快速原型设计与模型更新。
• 重构计算结构，降低线性代数求解的复杂度，提升性能而不引入近似。
• 将热效应整合至P2D模型中，包括欧姆、反应和可逆过程产生的热量，并考虑温度依赖的材料属性。

实验结果

研究问题

• RQ1 能否开发一种快速且高精度的求解器，用于全P2D模型，避免单颗粒模型等简化模型带来的不准确性？
• RQ2 如何在不损失物理保真度的前提下，降低求解P2D模型微分-代数方程组的计算负担？
• RQ3 通过JAX实现的自动微分在多大程度上能提升P2D模型实现的鲁棒性与可维护性？
• RQ4 重构后的线性代数运算在多大程度上能显著减少求解P2D系统所需的时间，同时保持精度？
• RQ5 所提出的求解器在模拟关键电池行为（如浓度分布、电压响应和热动力学）方面表现如何？

主要发现

• 通过重构线性代数问题并优化牛顿求解过程，该求解器在速度上显著优于标准黑箱DAE求解器（如DASSL、SUNDIALS）。
• 利用JAX进行自动微分可实现雅可比矩阵的精确且高效计算，减少实现错误，并简化模型更新。
• 该方法在所有状态变量（包括电解质浓度、固相浓度、电势和温度）上均保持高精度，即使在不同电流负载下亦然。
• 求解器能够准确处理温度依赖的材料属性，包括扩散系数、电导率和反应动力学，且采用阿伦尼乌斯型温度依赖关系。
• 由于JAX的函数式编程接口，该计算框架支持对模型表达式（如扩散系数、反应速率）的便捷修改，显著提升模型可扩展性。
• 该方法能够稳定且准确地模拟包含热效应的全P2D模型，涵盖电极和隔膜区域中欧姆、反应和可逆过程产生的热量。