[논문 리뷰] Splitting methods for Fourier spectral discretizations of the strongly magnetized Vlasov-Poisson and the Vlasov-Maxwell system

J. Ameres|arXiv (Cornell University)|2019. 07. 11.

Gas Dynamics and Kinetic Theory참고 문헌 40인용 수 5

한 줄 요약

이 논문은 강한 자기장에 의한 비틀림된 Vlasov-Poisson 및 Vlasov-Maxwell 시스템의 푸리에 스펙트럼 이산화에 대한 새로운 분할 방법을 제안하며, 4차원 및 3차원에서 고정밀 시뮬레이션을 가능하게 한다. 스트랑 분할과 해밀토니안 적분기를 사용하여 낮은 해상도에서도 에너지 오차가 10⁻¹² 이하로 유지되는 어려운 시험 케이스인 켈빈-헬름홀츠 및 웨이벨 스트리밍 불안정성에서 안정성과 기하학적 보존성을 입증한다.

ABSTRACT

Fourier spectral discretizations belong to the most straightforward methods for solving the unmagnetized Vlasov--Poisson system in low dimensions. In this article, this highly accurate approach is extended two the four-dimensional magnetized Vlasov--Poisson system with new splitting methods suited for strong magnetic fields. Consequently, a comparison to the asymptotic fluid model is provided at the example of a turbulent Kelvin--Helmholtz instability. For the three dimensional electromagnetic Vlasov--Maxwell system different novel charge conserving implementations of a Hamiltonian splitting are discussed and simulation results of the Weibel streaming instability are presented.

연구 동기 및 목표

• 강한 자기장 조건에서 고차원 Vlasov 시스템의 푸리에 스펙트럼 이산화에 대한 기하학적 분할 방법의 부족을 보완한다.
• 강한 자기장 하에서 4차원 자기장이 있는 Vlasov-Poisson 시스템에 대해 효율적이고 구조를 유지하는 시간 적분기를 개발한다.
• 해밀토니안 분할 기법을 3차원 Vlasov-Maxwell 시스템으로 확장하여 전하 및 에너지 보존을 보장한다.
• 켈빈-헬름홀츠 및 웨이벨 스트리밍과 같은 벤치마크 플라즈마 불안정성에 대해 저해상도 및 고해상도 조건에서 방법을 검증한다.
• 제안된 기법이 거친 격자와 강한 자기장 조건에서도 장기적인 안정성과 기하학적 정확성을 유지하는지 입증한다.

제안 방법

• 위치 공간과 속도 공간에서의 이동을 분리하기 위해 해밀토니안 분할을 적용하여 푸리에 변환을 통해 정확한 해를 구한다.
• 푸리에 스펙트럼 방법을 공간 및 속도 공간 이산화에 사용하며, 효율적인 계산을 위해 FFTW를 활용한다.
• 일정 계수 이동을 푸리에 공간에서 정확하게 통합하기 위해 지수 분할 기법(예: 스코벨의 방법 및 지수 보리스)을 구현한다.
• Vlasov-Maxwell 시스템에서 전하 보존을 확보하기 위해 기하학적 적분기를 사용하여 기저 해밀토니안 구조를 유지한다.
• 2차 정밀도를 가지며 에너지 및 운동량 보존 성질을 유지하는 2차 정밀도를 가지는 스트랑 분할을 시간 적분에 적용한다.
• 베셀 함수와 가우스 구적법을 사용하여 앙페르 방정식의 기하 평균화를 처리하여 스코벨의 방법에서 정확한 시간 적분을 가능하게 한다.

실험 결과

연구 질문

• RQ1 새로운 분할 기법을 사용하여 푸리에 스펙트럼 방법을 강한 자기장 조건에서 4차원 자기장이 있는 Vlasov-Poisson 시스템으로 효과적으로 확장할 수 있는가?
• RQ2 제안된 분할 기법은 강한 자기장 조건에서 켈빈-헬름홀츠 불안정성과 같은 난류 불안정성을 시뮬레이션할 때 어떻게 성능을 발휘하는가?
• RQ3 해밀토니안 분할 기법을 3차원 Vlasov-Maxwell 시스템에 적응시켜 전하 및 에너지 보존을 보장할 수 있는가?
• RQ4 장기 시뮬레이션에서 웨이벨 유형의 불안정성에 대해 해상도와 시간 간격이 에너지 및 운동량 오차에 어떤 영향을 미치는가?
• RQ5 기존의 의사 스펙트럼 방법과 비교해 새로운 기법이 기하학적 구조와 안정성을 얼마나 잘 유지하는가?

주요 결과

• 제안된 스코벨의 분할 방법은 점근적 극한이 없더라도 광범위한 자기장 강도에서 고정밀도와 안정성을 유지한다.
• 켈빈-헬름홀츠 불안정성의 경우, 강한 자기장 조건에서의 운동론적 모델은 유체 모델과 밀도적으로 일치하며 점근적 영역을 검증한다.
• 해밀토니안 분할을 사용한 Vlasov-Maxwell 솔버는 고해상도(128×128×128 격자점)에서 에너지 오차가 1.71×10⁻¹² 이하로 유지되며, 고차원 GEMPIC 방법과 유사한 성능을 보인다.
• 낮은 해상도(32×32×32)에서도 장기간의 시뮬레이션 기간 동안 안정성을 유지하며 에너지 오차가 10⁻¹² 이하로 유지되어 강건성을 입증한다.
• 웨이벨 스트리밍 불안정성의 시뮬레이션은 운동론적 에너지와 자기장의 일관된 증가를 보이며, 대칭 및 비대칭 케이스에서 각각 다른 필라멘트 패턴을 해석할 수 있었다.
• 베셀 함수와 가우스 구적법의 사용은 스코벨의 방법에서 앙페르 방정식의 정확한 시간 적분을 가능하게 하여 장기적 보존성의 핵심 요소가 되었다.