[논문 리뷰] A numerical analysis focused comparison of several Finite Volume schemes for an Unipolar Degenerated Drift-Diffusion Model

Clément Cancès, Claire Chainais-Hillairet외 2명|arXiv (Cornell University)|2019. 01. 01.

Advanced Mathematical Modeling in Engineering참고 문헌 35인용 수 1

한 줄 요약

이 논문은 화학적 위치 에너지 h(c) = log(c/(1−c))를 갖는 단극성 비가역 드리프트-디퓨전 모델에 대해 네 가지 유한체적 스킴의 비교 수치 분석을 제시한다. 스킴들은 중심형, Sedan, Bessemoulin-Chatard, 활성도 기반의 서로 다른 유량 형식에 기반하며, 이동도 η(c) = c를 사용한다. 연구는 두 스킴에 대해 안정성, 존재성, 수렴성을 입증하였으며, 수치 실험을 통해 Sedan 및 Bessemoulin-Chatard 스킴이 양성과 엔트로피 구조를 보존하는 데 뛰어난 성능을 보였다.

ABSTRACT

In this paper, we consider an unipolar degenerated drift-diffusion system where the relation between the concentration of the charged species $c$ and the chemical potential $h$ is $h(c)=\log \frac{c}{1-c}$. We design four different finite volume schemes based on four different formulations of the fluxes. We provide a stability analysis and existence results for the four schemes. The convergence proof with respect to the discretization parameters is established for two of them. Numerical experiments illustrate the behaviour of the different schemes.

연구 동기 및 목표

• h(c) = log(c/(1−c))를 갖는 단극성 비가역 드리프트-디퓨전 모델에 대해 네 가지 유한체적 스킴을 분석하고 비교하기.
• 모든 네 스킴에 대해 안정성 및 존재성 결과를 확립하기.
• 두 스킴에 대해 이산화 매개변수에 대한 수렴성을 증명하기.
• 양성과 엔트로피 구조 같은 물리적 성질을 유지하는 데 있어 스킴의 수치 성능 평가하기.
• 다종성 이온 액체 및 고체 산화물 전해질 모델로의 프레임워크 확장 기초 제공하기.

제안 방법

• 중심형, Sedan, Bessemoulin-Chatard, 활성도 기반의 서로 다른 유량 형식에 기반한 네 가지 유한체적 스킴 제안.
• h(c) = log(c/(1−c)) 및 η(c) = c를 사용하여 유량 항식 J = −c∇(h(c) + Φ) 적용.
• 엔트로피 구조 분석을 통해 모든 스킴에 대한 안정성 및 존재성 결과 유도.
• 이산 엔트로피 추정과 컴팩턴스 추론을 통해 중심형 및 Sedan 스킴의 수렴성 확립.
• 면 농도 기능성 C를 도입하고, 일관성 평가를 위해 기준 기능성 eC와 비교.
• 수치 실험을 수행하여 농도 보존, 수렴 속도, 엔트로피 감쇠 평가.

실험 결과

연구 질문

• RQ1 어느 유한체적 스킴이 농도 c가 1−로 수렴하는 비가역 영역에서 농도의 양성을 가장 잘 유지하는가?
• RQ2 다양한 유량 형식은 수치 스킴의 수렴성과 안정성에 어떤 영향을 미치는가?
• RQ3 Sedan 및 Bessemoulin-Chatard 스킴은 비가역 한계에서 엔트로피 구조를 유지하고 비물리적 진동을 피할 수 있는가?
• RQ4 면 농도 기능성과 스킴의 일관성 간의 관계는 무엇인가?
• RQ5 수렴 속도와 엔트로피 감쇠 측면에서 스킴의 수치 성능은 어떠한가?

주요 결과

• Sedan 및 Bessemoulin-Chatard 스킴은 중심형 및 활성도 기반 스킴보다 농도의 양성 보존과 엔트로피 감쇠에서 뛰어난 성능을 보였다.
• 이산 엔트로피 추정과 컴팩턴스 추론을 통해 중심형 및 Sedan 스킴에 대해 수렴성이 엄밀히 증명되었다.
• Sedan 스킴의 면 농도 기능성 C는 기준 기능성 eC에 대해 균일하게 유계이므로 일관성이 보장된다.
• 활동도 기반 스킴은 cK/cL → ∞ 근처에서 비한계적인 비율으로 인해 조건 (C.2)를 충족하지 못한다.
• 수치 실험을 통해 Sedan 및 Bessemoulin-Chatard 스킴이 엔트로피 구조를 유지하고 비물리적 진동을 피하는 것으로 확인되었다.
• 중심형 스킴은 단순하지만, c = 1 근처에서 나쁜 조건으로 인해 비가역 영역에서 최적의 성능을 내지 못한다.