[논문 리뷰] Unconditional positivity-preserving and energy stable schemes for a reduced Poisson-Nernst-Planck system

Hailiang Liu, Wumaier Maimaitiyiming|arXiv (Cornell University)|2019. 09. 28.

Quantum, superfluid, helium dynamics참고 문헌 30인용 수 6

한 줄 요약

이 논문은 이온 채널을 모델링하는 축소된 1차원 포아송-너른스트-플랭크 시스템을 위한 이차순, 조건부로 양성 보존 및 에너지 안정적인 유한체적 스킴을 제시한다. 비로그형 라운드 포맷의 반익명적 시간 이산화를 사용함으로써, 방법은 어떠한 시간 간격 크기나 포아송 해법 선택과도 무관하게 이온 농도의 양성 보존을 보장하며, 영류 경계 조건 하에서 질량 보존과 자유 에너지 감소를 유지한다.

ABSTRACT

The Poisson-Nernst-Planck (PNP) system is a widely accepted model for simulation of ionic channels. In this paper, we design, analyze, and numerically validate a second order unconditional positivity-preserving scheme for solving a reduced PNP system, which can well approximate the three dimensional ion channel problem. Positivity of numerical solutions is proven to hold true independent of the size of time steps and the choice of the Poisson solver. The scheme is easy to implement without resorting to any iteration method. Several numerical examples further confirm the positivity-preserving property, and demonstrate the accuracy, efficiency, and robustness of the proposed scheme, as well as the fast approach to steady states.

연구 동기 및 목표

• 축소된 포아송-너른스트-플랭크 시스템에서 이온 농도의 물리적 양성 보존을 보장하는 수치 스킴을 개발한다.
• 시간 간격 크기나 포아송 해법 선택과 무관하게 무조건적인 양성 보존을 보장한다.
• 영류 경계 조건 하에서 질량 보존과 자유 에너지 감소를 유지한다.
• 이온 이동 시뮬레이션을 위한 단순하고 반복이 필요 없으며 쉽게 구현 가능한 수치 방법을 제공한다.
• 다양한 초기 조건과 경계 조건 하에서 스킴의 정확도, 강건성 및 안정 상태 수렴 속도를 검증한다.

제안 방법

• 양성 보존 이산화를 가능하게 하기 위해 축소된 PNP 시스템을 비로그형 라운드 포맷으로 변환한다.
• 이차 정확도를 갖는 반익명적 유한체적 방법을 시간에 대해 적용한다.
• 비균일 메쉬 위에서 세포 중심 기반의 유한체적 공간 이산화를 사용한다.
• 이온 농도와 전기적 위치 에너지 방정식을 결합하기 위해 비스태그너드 격자 배열을 활용한다.
• 반복 보정을 피하기 위해 유량 재구성 기반의 양성 보존 업데이트 전략을 구현한다.
• 디리클레 및 로빈 유형 경계 조건(영류 및 전위 반사 조건 포함)을 모두 통합한다.

실험 결과

연구 질문

• RQ1 축소된 PNP 시스템에서 어떠한 시간 간격 크기나 포아송 해법 선택과도 무관하게 이온 농도의 양성을 무조건 보존하는 이차순 수치 스킴을 구성할 수 있는가?
• RQ2 영류 경계 조건 하에서 스킴이 질량 보존과 자유 에너지 감소를 유지하는가?
• RQ3 스킴은 초기 자료, 시간 간격 크기, 포아송 해법 선택에 대해 얼마나 강건한가?
• RQ4 다양한 이온 채널 구조에서 안정 상태 해수에 대한 수렴 거동은 어떠한가?
• RQ5 고정 전하 밀도와 초기 농도 프로파일은 안정 상태 이온 분포에 어떤 영향을 미치는가?

주요 결과

• 제안된 스킴은 어떠한 시간 간격 크기나 포아송 해법 선택과도 무관하게 이온 농도의 양성을 무조건 보존하며, 정리 3.1에서 엄밀히 증명되었다.
• 스킴은 공간과 시간에 대해 이차 정확도를 달성하며, 수치 수렴률이 이론적 기대와 일치한다.
• 예제 4.4의 초기 자료 (4.8)를 사용한 경우, τ=10⁻³로 설정했을 때 시스템은 2.7410 시간 단위(2741 반복) 내에 안정 상태에 도달하며 잔차 오차 ||ψⁿ−ψⁿ⁻¹||∞ ≈ 10⁻⁷이다.
• 임의의 초기 자료를 사용한 경우, 시스템은 2.0440 시간 단위(2044 반복) 내에 안정 상태에 도달하여 초기 조건과 무관한 강건성을 확인한다.
• 예제 4.4에서 두 자료 세트 모두 동일한 프로파일을 보이며, 안정 상태 이온 분포가 초기 자료에 의존하지 않음을 확인한다.
• 영류 및 로빈 경계 조건이 적용된 예제 4.5에서 스킴은 총 질량을 보존하고 시간에 따라 단조롭게 자유 에너지가 감소함을 보이며, 정리 3.2를 검증한다.