[논문 리뷰] A positivity-preserving, energy stable and convergent numerical scheme for the Poisson-Nernst-Planck system

Chun Liu, Cheng Wang외 3명|arXiv (Cornell University)|2020. 09. 17.

Gas Dynamics and Kinetic Theory참고 문헌 35인용 수 8

한 줄 요약

이 논문은 유일한 해법 보존, 양수성 보존, 에너지 안정성, 최적 수렴 속도를 동시에 확보하는 포isson-너른스트-플랑크(Poisson-Nernst-Planck, PNP) 시스템을 위한 유한차분 스킴을 제안한다. 에너지적 변분 접근법(Energetic Variational Approach, EnVarA)을 활용하여 이동도는 명시적으로, 로그/전자기적 항은 암시적으로 처리함으로써 볼록성과 특이 로그 잠재력의 특성을 활용하여 농도 폭주를 방지하며, 비정수 이동도 조건 하에서 최초로 이론적으로 수렴성과 안정성을 확보한다.

ABSTRACT

In this paper we propose and analyze a finite difference numerical scheme for the Poisson-Nernst-Planck equation (PNP) system. To understand the energy structure of the PNP model, we make use of the Energetic Variational Approach (EnVarA), so that the PNP system could be reformulated as a non-constant mobility $H^{-1}$ gradient flow, with singular logarithmic energy potentials involved. To ensure the unique solvability and energy stability, the mobility function is explicitly treated, while both the logarithmic and the electric potential diffusion terms are treated implicitly, due to the convex nature of these two energy functional parts. The positivity-preserving property for both concentrations, $n$ and $p$, is established at a theoretical level. This is based on the subtle fact that the singular nature of the logarithmic term around the value of $0$ prevents the numerical solution reaching the singular value, so that the numerical scheme is always well-defined. In addition, an optimal rate convergence analysis is provided in this work, in which many highly non-standard estimates have to be involved, due to the nonlinear parabolic coefficients. The higher order asymptotic expansion (up to third order temporal accuracy and fourth order spatial accuracy), the rough error estimate (to establish the $\ell^\infty$ bound for $n$ and $p$), and the refined error estimate have to be carried out to accomplish such a convergence result. In our knowledge, this work will be the first to combine the following three theoretical properties for a numerical scheme for the PNP system: (i) unique solvability and positivity, (ii) energy stability, and (iii) optimal rate convergence. A few numerical results are also presented in this article, which demonstrates the robustness of the proposed numerical scheme.

연구 동기 및 목표

• 이온 농도의 양수성과 에너지 감산을 유지하는 물리적 성질을 갖는 PNP 시스템을 위한 수치 스킴을 개발한다.
• 비정수 이동도와 특이 로그 잠재력 조건 하에서 유한차분 스킴의 엄밀한 수렴 분석을 수립한다.
• 유일한 해법/양수성, 에너지 안정성, 최적 수렴 속도라는 세 가지 핵심 이론적 성질을 하나의 수치 프레임워크에서 통합한다.
• 에너지 함수에서 비선형 포물형 계수와 특이 로그 항으로 인한 도전 과제를 극복한다.

제안 방법

• 에너지적 변분 접근법(EnVarA)을 활용해 PNP 시스템을 비정수 이동도를 갖는 H⁻¹ 경사하강 흐름으로 재구성한다.
• 이동도 함수를 명시적으로 처리하여 이산 시스템의 유일한 해법을 보장한다.
• 에너지 구조의 볼록성 덕분에 로그 항과 전기적 잠재력 확산 항에 대해 암시적 시간 이산화를 적용한다.
• 고차항 점근 전개(시간적으로는 3차까지, 공간적으로는 4차까지 정확도 확보)를 활용해 농도에 대한 거친 ℓ∞ 유계를 도출한다.
• 정교한 오차 추정을 통해 비선형 포물형 계수의 영향을 극복하고 최적 수렴 속도를 확보한다.
• 영점 근처에서 로그 잠재력의 특이성 특성을 활용하여 수치 해의 양수성을 이론적으로 보장한다.

실험 결과

연구 질문

• RQ1 PNP 시스템에 대해 유일한 해법과 농도의 양수성을 동시에 보장하는 유한차분 스킴을 구성할 수 있는가?
• RQ2 비정수 이동도와 특이 로그 잠재력 조건 하에서 에너지 안정성을 엄밀히 증명할 수 있는가?
• RQ3 이러한 복잡한 비선형성 조건 하에서 PNP 시스템의 수치 스킴의 최적 수렴 속도는 무엇인가?
• RQ4 양수성, 에너지 안정성, 최적 수렴 속도를 동시에 보장하는 통합된 이론적 프레임워크를 개발할 수 있는가?

주요 결과

• 특이 로그 잠재력의 특성 덕분에 수치 해의 폭주를 방지함으로써, 제안된 스킴은 이산 수준에서 이온 농도 n과 p의 유일한 해법과 양수성을 모두 보존한다.
• 에너지 함수의 볼록성과 로그 항 및 전기적 잠재력 항의 암시적 처리를 통해 에너지 안정성이 증명된다.
• 최적 수렴 속도가 확립되었으며, 시간적으로는 3차까지, 공간적으로는 4차까지 정확도를 확보한 고차항 점근 전개를 활용하였다.
• 비선형 포물형 계수를 다루기 위해 농도에 대한 거친 오차 추정을 통해 ℓ∞ 유계가 유도되었다.
• 정교한 오차 추정을 통해 최적 수렴 속도가 확인되었으며, 비정수 이동도와 특이 잠재력으로 인한 도전 과제를 극복하였다.
• 수치 결과는 이산 수준에서 질량 보존, 에너지 감산, 양수성 보존을 확인하였으며, 이는 이론 분석의 타당성을 검증한다.