[논문 리뷰] Rubrics for Charge Conserving Current Mapping in Finite Element Particle in Cell Methods

Zane D. Crawford, Scott O’Connor외 2명|arXiv (Cornell University)|2021. 01. 28.

Plasma Diagnostics and Applications인용 수 4

한 줄 요약

이 논문은 de Rham 복합체 사상과 이산 하우지 연산자를 활용하여 유한요소 입자-장(Fe-PIC) 방법에서 전하 보존을 보장하기 위한 일반적인 수학적 기준을 수립한다. 전류 맵핑이 모서리 기반과 면 기반 유한요소 기저 간의 쌍대성과 일치할 때 전하 보존이 달성되며, 수치 결과를 통해 검증된 새로운 본질적으로 전하 보존 성질을 갖는 FE-PIC 방법을 제시한다.

ABSTRACT

Modeling of kinetic plasmas using electromagnetic particle in cell methods (EM-PIC) is a problem that is well worn, in that methods developed have been used extensively both understanding physics and exploiting them for device design. EM-PIC tools have largely relied on finite difference methods coupled with particle representations of the distribution function. Refinements to ensure consistency and charge conservation have largely been an ad-hoc efforts specific to finite difference methods. Meanwhile, solution methods for field solver have grown by leaps and bounds with significant performance metrics compared to finite difference methods. Developing new EM-PIC computational schemes that leverage modern field solver technology means re-examining analysis framework necessary for self-consistent EM-PIC solution. In this paper, we prescribe general rubrics for charge conservation, demonstrate how these are satisfied in conventional finite difference PIC as well as finite element PIC, and prescribe a novel charge conserving finite element PIC. Our effort leverages proper mappings on to de-Rham sequences and lays a groundwork for understanding conditions that must be satisfied for consistency. Several numerical results demonstrate the applicability of these rubrics.

연구 동기 및 목표

• 전자기 입자-장(EM-PIC) 시뮬레이션에서 전하 보존을 보장하기 위한 일반적이고 해법에 종속되지 않는 기준을 개발하기 위해.
• 유한요소 방법에서 전류 맵핑이 전하 보존을 유지하는 데 필요한 수학적 조건을 규명하기 위해.
• 기존의 유한차분 PIC에서 현대적인 유한요소 장 해법으로 전하 보존 원리를 확장하기 위해.
• FE-PIC 체계에서 일致하고 발산이 없는 전류 침착을 체계적으로 구성하기 위한 프레임워크를 제공하기 위해.
• 수치적 검증을 통해 제안된 기준의 적용 가능성과 정확성을 입증하기 위해.

제안 방법

• de Rham 복합체 내 이산 미분형식 간의 쌍대성에 기반하여 전하 보존을 요구 조건으로 설정함으로써, 전류 맵핑이 하우지 스타 연산자와 코일 관계를 존중하도록 보장한다.
• 전기장과 자기장에 대해 모서리 기반(1형식) 및 면 기반(2형식) 유한요소 기저 함수를 사용하여 이산 맥스웰 방정식을 유도한다.
• 입자 궤적의 선적분을 통해 전류 자유도를 정의함으로써 암페어의 법칙과 패러데이의 법칙과의 일致성을 확보한다.
• 이산 하우지 스타 연산자를 적용하여 원시 메쉬와 쌍대 메쉬 간의 맵핑을 수행함으로써 정확한 장 및 전류 보간을 가능하게 한다.
• de Rham 수열을 사용하여 이산 코일 및 발산 연산자가 연속 경우와 동일한 대수적 항등식을 만족하도록 보장한다.
• 시간에 따른 전하 보존과 맥스웰 방정식과의 일致성을 확인하는 수치 실험을 통해 체계를 검증한다.

실험 결과

연구 질문

• RQ1 유한요소 PIC 방법에서 전류 맵핑이 전하 보존을 보장하기 위해 충족해야 할 수학적 조건은 무엇인가?
• RQ2 임의의 유한요소 장 해법에 대해 전하 보존 성질을 갖는 전류 침착을 체계적으로 유도할 수 있는 방법은 무엇인가?
• RQ3 de Rham 복합체 사상과 이산 하우지 연산자는 FE-PIC에서 일致적인 전하 보존을 어떻게 가능하게 하는가?
• RQ4 기존의 FDTD-PIC와 비교할 때 제안된 FE-PIC 체계는 전하 보존 및 일치성 측면에서 어떤가?
• RQ5 제안된 기준은 유한요소 방법 외의 다른 장 해법 유형으로 일반화될 수 있는가?

주요 결과

• 논문은 de Rham 복합체를 통해 모서리 기반과 면 기반 유한요소 기저 간의 쌍대성이 존중될 경우 FE-PIC에서 전하 보존이 보장된다고 규명한다.
• 제안된 방법은 모든 시간 단계에서 이산 연속 방정식이 만족됨을 보장하여 부적절한 전하 생성을 제거한다.
• 수치 결과는 장기간의 시뮬레이션 동안 기계 정밀도 수준에서 전하 보존이 유지됨을 확인한다.
• 구조적 격자에 적용했을 때 기존의 FDTD-PIC 결과를 재현함으로써 기존 체계와의 일致성을 검증한다.
• de Rham 수열과 이산 하우지 연산자의 사용은 코일 및 발산 연산자가 연속 경우와 동일한 대수적 항등식을 만족하도록 보장한다.
• 이 프레임워크는 후처리 또는 발산 정리가 필요 없이 새로운 본질적으로 전하 보존 성질을 갖는 FE-PIC 체계를 구성할 수 있도록 한다.