[論文レビュー] 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

ひとこと要約

この論文では、有限差分法、後退オイラー時刻積分、およびニュートン法を用いて、リチウムイオン電池の擬似2次元（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.

研究の動機と目的

• リチウムイオン電池の完全な擬似2次元（P2D）電気化学的モデルの高い計算コストに対処すること。
• 特定の用途において十分な精度を欠く単一粒子モデルなどの簡略化モデルの限界を克服すること。
• 完全なP2Dモデルの忠実性を保ちつつ、計算負荷を軽減する数値的に効率的なソルバーを開発すること。
• バッテリマネジメントシステム（BMS）におけるP2Dモデルの実用的導入を可能にし、リアルタイム制御および状態推定を実現すること。
• 自動微分によるモデル係数の容易な変更を可能にし、モデルの柔軟性と頑健性を高めること。

提案手法

• 空間方向の離散化に有限差分法（FDM）を用い、x方向（電極厚さ）および径方向r（粒子拡散）の1次元領域にP2Dモデルを実装する。
• 微分代数方程式系の時刻積分において無条件安定性を確保するため、後退オイラー法を適用する。
• 各時刻ステップで得られる非線形代数方程式系をニュートン法で解き、高速収束を実現する。
• JAXを用いてヤコビ行列と感度を正確かつ効率的に計算し、迅速なプロトタイピングとモデル更新を可能にする。
• 線形代数のソルブの複雑さを低減するための計算構造の再編成を実施し、近似を加えずに性能向上を図る。
• 熱的効果をP2Dモデルに統合し、オーム損失、反応、および逆転可能プロセスに起因する発熱を含め、温度依存性材料特性を考慮する。

実験結果

リサーチクエスチョン

• RQ1 単一粒子モデルなどの簡略化モデルの誤差を回避しつつ、完全なP2Dモデルに対して高速かつ高精度なソルバーを開発できるか？
• RQ2 P2Dモデルの微分代数方程式の解法にかかる計算負荷を、物理的忠実性を損なわずどのように低減できるか？
• RQ3 JAXを用いた自動微分が、P2Dモデル実装の頑健性と保守性をどの程度向上できるか？
• RQ4 再構築された線形代数演算により、P2Dシステムの解法に要する時間が著しく短縮され、精度を保持できるか？
• RQ5 提案されたソルバーは、濃度プロファイル、電圧応答、熱的ダイナミクスといった重要なバッテリ挙動をどの程度正確にシミュレートできるか？

主な発見

• 線形代数問題の再構築とニュートン法の最適化により、標準的なブラックボックスDAEソルバー（例：DASSL、SUNDIALS）と比較して顕著な高速化を達成した。
• JAXを用いた自動微分により、ヤコビ行列の正確かつ効率的な計算が可能になり、実装エラーの低減とモデル更新の簡素化が実現された。
• 電解質濃度、固体相濃度、電位、温度といったすべての状態変数において、高い精度を維持した。これは、変動する電流負荷下でも同様に成立する。
• 温度依存性材料特性（拡散率、電導率、反応キネティクスなど）を、アレニウス型の温度依存性を含めて適切に扱うことができた。
• JAXの関数型プログラミングインタフェースのおかげで、モデル式（例：拡散率、反応速度）の容易な変更が可能となり、モデルの拡張性が向上した。
• オーム損失、反応、および逆転可能プロセスに起因する発熱を含め、全電極およびセパレータ領域において、熱的効果を含む完全なP2Dモデルの安定的かつ高精度なシミュレーションが可能となった。