[논문 리뷰] A fast solver for the pseudo-two-dimensional model of lithium-ion batteries

Rachel Han, Colin B. Macdonald외 1명|arXiv (Cornell University)|2021. 11. 17.

Advanced Battery Technologies Research참고 문헌 20인용 수 7

한 줄 요약

이 논문은 유한차분법, 역오일러 시간 적분, 뉴턴 방법을 사용하여 리튬이온 배터리의 가상이중차원(P2D) 모델을 빠르고 정확하게 푸는 해법을 제시한다. 파이썬으로 작성된 이 방법은 선형 대수의 재구성과 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 모델의 실용적 구현을 가능하게 하기 위해.
• 자동미분를 통해 모델 계수의 수정을 쉽게 할 수 있도록 하여 모델의 유연성과 견고성을 향상시키기 위해.

제안 방법

• 1D 도메인에서의 x방향(전극 두께)과 반경 r방향(입자 확산)에 대해 유한차분법(FDM)을 사용하여 P2D 모델을 구현한다.
• 미분-代수 시스템의 시간 스텝에서 무조건적 안정성을 확보하기 위해 역오일러 시간 적분을 적용한다.
• 각 시간 단계에서 유도된 비선형 대수 시스템을 뉴턴 방법으로 풀어 빠른 수렴을 보장한다.
• JAX를 사용하여 자동미분를 통해 야코비안과 민감도를 정확하고 효율적으로 계산함으로써 빠른 프로토타이핑과 모델 업데이트를 가능하게 한다.
• 선형 대수 해법의 복잡도를 줄이기 위해 계산 구조를 재구성하여 정확도를 유지하면서 성능을 향상시킨다.
• 열역학 효과를 P2D 모델에 통합하여 오믹, 반응, 가역 과정에서 발생하는 열 생성을 포함하고, 온도에 따라 변화하는 물성값을 반영한다.

실험 결과

연구 질문

• RQ1 단일입자모델과 같은 단순화된 모델의 부정확성을 피하면서도 전체 P2D 모델에 대해 빠르고 정확한 해법을 개발할 수 있는가?
• RQ2 P2D 모델의 미분-대수 방정식을 해결하는 데 드는 계산 부담을 물리적 정밀도를 훼손하지 않고 어떻게 줄일 수 있는가?
• RQ3 JAX를 통한 자동미분는 P2D 모델 구현의 견고성과 유지보수성에 어느 정도 기여하는가?
• RQ4 다시 구조화된 선형 대수 연산은 정확도를 유지하면서도 P2D 시스템을 푸는 데 소요되는 시간을 상당히 줄일 수 있는가?
• RQ5 제안된 해법은 농도 프로파일, 전압 반응, 열역학적 거동과 같은 핵심 배터리 거동을 얼마나 잘 시뮬레이션하는가?

주요 결과

• 선형 대수 문제의 재구성과 뉴턴 해법 최적화를 통해 표준 블랙박스 DAE 해법(DASSL, SUNDIALS 등)에 비해 상당한 속도 향상을 달성하였다.
• JAX를 통한 자동미분를 활용하여 야코비안의 정확하고 효율적인 계산이 가능해져 구현 오류가 감소하고 모델 업데이트가 간편해졌다.
• 다양한 전류 하중 조건에서도 전해질 농도, 고체상 농도, 전위, 온도 등 모든 상태 변수에서 높은 정확도를 유지하였다.
• 배터리 물성값(확산도, 전도도, 반응kinetics 등)의 온도 의존성(아레니우스 형식)을 정확히 반영하여 온도 변화에 따른 물성 변화를 효과적으로 처리하였다.
• JAX의 기능형 프로그래밍 인터페이스 덕분에 모델 식(예: 확산도, 반응속도)의 수정이 용이하여 모델의 확장성이 향상되었다.
• 열역학 효과를 포함한 전체 P2D 모델의 안정적이고 정확한 시뮬레이션을 전극 및 분리막 영역 전반에 걸쳐 성공적으로 수행하였다. 이는 오믹, 반응, 가역 과정에서 발생하는 열 생성을 포함한다.