[논문 리뷰] On Runge-Kutta discontinuous Galerkin schemes for Vlasov-Poisson systems

Yingda Cheng, Irene M. Gamba외 1명|arXiv (Cornell University)|2012. 09. 28.

Gas Dynamics and Kinetic Theory참고 문헌 49인용 수 3

한 줄 요약

이 논문은 충돌 없는 플라즈마를 모델링하는 일차원 Vlasov-Poisson 시스템을 해결하기 위한 고차수 Runge-Kutta 불연속 Galerkin (RKDG) 스킴을 제시한다. 이 방법은 보존 성질을 보장하며, 리미터를 통해 분포 함수의 양성도 유지하고, 푸리에 분석과 다항식 근사법을 통해 Landau 감쇠 및 두 개의 스트림 불안정성과 같은 벤치마크 문제에 대해 정확한 해를 도출한다.

ABSTRACT

In this paper we consider Runge-Kutta discontinuous Galerkin (RKDG) schemes for Vlasov-Poisson systems that model collisionless plasmas. One-dimensional systems are emphasized. The RKDG method, originally devised to solve conservation laws, is seen to have excellent conservation properties, be readily designed for arbitrary order of accuracy, and capable of being used with a positivity-preserving limiter that guarantees positivity of the distribution functions. The RKDG solver for the Vlasov equation is the main focus, while the electric field is obtained through the classical representation by Green's function for the Poisson equation. A rigorous study of recurrence of the DG methods is presented by Fourier analysis, and the impact of different polynomial spaces and the positivity-preserving limiters on the quality of the solutions is ascertained. Several benchmark test problems, such as Landau damping, two-stream instability and the KEEN (Kinetic Electrostatic Electron Nonlinear) wave, are given.

연구 동기 및 목표

• 일차원 충돌 없는 플라즈마 역학에서 Vlasov-Poisson 시스템을 해결하기 위한 고차수, 보존적인 수치 스킴을 개발하는 것.
• 양성 보존 리미터를 사용하여 운동학적 시뮬레이션에서 분포 함수의 양성을 확보하는 것.
• 푸리에 분석을 통한 DG 방법의 재귀성 및 안정성 성질을 분석하는 것.
• Landau 감쇠 및 두 개의 스트림 불안정성과 같은 표준 벤치마크 문제에 대한 스킴의 검증.
• 다항식 차수와 리미터 선택이 해 품질과 정확도에 미치는 영향을 조사하는 것.

제안 방법

• Runge-Kutta 불연속 Galerkin (RKDG) 방법을 Vlasov 방정식에 적용하여 고차수 정밀도와 국소 보존성을 확보한다.
• 정확한 잠재력 재구성에 위해 포아송 방정식의 그린 함수 표현을 전기장 계산에 사용한다.
• 부정적인 분포 함수를 방지하기 위해 RKDG 프레임워크에 양성 보존 리미터를 통합한다.
• DG 이산화의 재귀성 및 스펙트럼 성질을 연구하기 위해 푸리에 분석을 활용한다.
• 해 운동 방정식의 수렴성과 정확도 평가를 위해 다양한 차수의 다항식 공간을 사용한다.
• 안정성, 정확도, 보존 성질을 평가하기 위해 표준 벤치마크 문제에 대한 수치 실험을 수행한다.

실험 결과

연구 질문

• RQ1 다양한 다항식 공간은 Vlasov-Poisson 시스템에 대한 RKDG 스킴의 정확도와 안정성에 어떤 영향을 미치는가?
• RQ2 양성 보존 리미터는 장기 시뮬레이션에서 물리적 일관성을 어느 정도 유지하는가?
• RQ3 DG 이산화에서 재귀성의 역할은 무엇이며, 해 품질에 어떤 영향을 미치는가?
• RQ4 RKDG 스킴은 Landau 감쇠 및 두 개의 스트림 불안정성과 같은 주요 플라즈마 현상을 얼마나 잘 포착하는가?
• RQ5 고차수 공간 및 시간 이산화가 해 정확도와 보존 성질에 어떤 영향을 미치는가?

주요 결과

• RKDG 스킴은 공간 및 시간에서 고차수 정밀도를 보이며, 매끄러운 해에 대해 이론적 기대와 일치하는 수렴 속도를 나타낸다.
• 양성 보존 리미터는 모든 테스트된 벤치마크에서 음이 아닌 분포 함수를 성공적으로 유지한다.
• 푸리에 분석 결과, RKDG 방법은 유리한 재귀 성질을 보이며, 부정확한 진동을 최소화한다.
• 전기장의 기대되는 지수 감쇠를 유지하면서, 전기장에 대한 수치적 소산이 최소화된 상태에서 Landau 감쇠를 정확히 포착한다.
• 두 개의 스트림 불안정성의 경우, 스킴은 불안정성의 성장과 비선형 포화를 고해상도로 해석한다.
• KEEN 파동 해는 정확하게 시뮬레이션되어, 이 스킴이 복잡한 비선형 운동학 현상을 다룰 수 있음을 확인한다.