[논문 리뷰] Runko: Modern multi-physics toolbox for simulating plasma

Joonas Nättilä|arXiv (Cornell University)|2019. 06. 14.

Magnetic confinement fusion research참고 문헌 1인용 수 7

한 줄 요약

Runko는 다중 물리 및 다중 스케일 모델링을 지원하는 오픈소스이자 하이브리드 C++/파이썬 프레임워크로, 천체물리학적 플라즈마를 시뮬레이션하는 데 사용됩니다. 최신 C++14를 활용해 성능을 확보하고 파이썬을 통해 고수준 제어를 가능하게 하여 빠르고 확장 가능한 시뮬레이션을 가능하게 하며, 상대론적 PIC 모듈을 통해 자기장 지배 쌍입자 플라즈마에서의 붕괴하는 상대론적 난류가 열적 난류 캐스케이드와 별개의 비열적 입자 집단을 생성하고, 이들이 에너지 수납체로 작용하는 것으로 밝혀졌습니다.

ABSTRACT

Runko is a new open-source plasma simulation framework implemented in C++ and Python. It is designed to function as an easy-to-extend general toolbox for simulating astrophysical plasmas with different theoretical and numerical models. Computationally intensive low-level kernels are written in modern C++14 taking advantage of polymorphic classes, multiple inheritance, and template metaprogramming. High-level functionality is operated with Python3 scripts. This hybrid program design ensures fast code and ease of use. The framework has a modular object-oriented design that allow the user to easily add new numerical algorithms to the system. The code can be run on various computing platforms ranging from laptops (shared-memory systems) to massively parallel supercomputer architectures (distributed-memory systems). The framework also supports heterogeneous multi-physics simulations in which different physical solvers can be combined and run simultaneously. Here we report on the first results from the framework's relativistic particle-in-cell (PIC) module. Using the PIC module, we simulate decaying relativistic kinetic turbulence in suddenly stirred magnetically-dominated pair plasma. We show that the resulting particle distribution can be separated into a thermal part that forms the turbulent cascade and into a separate decoupled non-thermal particle population that acts as an energy sink for the system.

연구 동기 및 목표

• 복잡한 천체물리학적 플라즈마 시스템을 모델링하기 위한 유연하고 확장 가능하며 고성능의 시뮬레이션 프레임워크를 개발하기 위해.
• 다양한 물리적 해법이 하나의 계산 환경 내에서 공존하고 결합될 수 있도록 하여 다중 물리 플라즈마 현상을 시뮬레이션하는 데 도전하는 것.
• 공유 메모리 및 분산 메모리 아키텍처를 모두 지원하여 랩톱부터 슈퍼컴퓨터까지의 이식성과 확장성을 확보하는 것.
• 특히 비열적 입자 집단의 발생을 포함한 자기장 지배 쌍입자 플라즈마에서의 상대론적 운동학적 난류를 연구하는 것.

제안 방법

• 프레임워크는 하이브리드 아키텍처를 사용하며, 저수준 계산 커널은 최신 C++14를 활용해 다형성 클래스, 다중 상속, 템플릿 메타프로그래밍을 통해 성능과 확장성을 확보합니다.
• 고수준 제어 및 워크플로우 관리는 파이썬3 스크립트를 통해 처리되어 사용성과 빠른 프로토타이핑을 용이하게 합니다.
• 모듈러하고 객체지향적인 설계 덕분에 사용자는 쉽게 새로운 수치 알고리즘과 물리 모델을 시뮬레이션 파이프라인에 통합할 수 있습니다.
• 다양한 해법(예: 상대론적 PIC 모듈 포함)을 병렬로 통합하고 실행할 수 있도록 하여 이질적 다중 물리 시뮬레이션을 지원합니다.
• 상대론적 PIC 모듈은 입자-격자 방법을 사용하여 자기장이 가해진 플라즈마 내에서 상대론적 전자 및 양성전자 동역학을 시뮬레이션합니다.
• 랩톱부터 슈퍼컴퓨터까지 다양한 컴퓨팅 플랫폼에 배포되며, 공유 메모리 및 분산 메모리 병렬 처리를 모두 활용합니다.

실험 결과

연구 질문

• RQ1 갑작스럽게 교반된 자기장 지배 쌍입자 플라즈마에서 상대론적 운동학적 난류는 어떻게 진화하는가?
• RQ2 그러한 난류 시스템에서 나타나는 별개의 입자 집단 구성 성분은 무엇인가?
• RQ3 난류 캐스케이드 동안 열적 및 비열적 입자 집단 간의 에너지 분할은 어떻게 이루어지는가?
• RQ4 비열적 입자 집단이 상대론적 플라즈마 시스템에서 안정적인 에너지 수납체로 작용할 수 있는가?
• RQ5 Runko 프레임워크는 복잡한 플라즈마 현상의 확장 가능한 다중 물리 시뮬레이션을 어느 정도 지원할 수 있는가?

주요 결과

• 시뮬레이션은 열적 입자 집단이 난류 에너지 캐스케이드에 참여하는 데 반해, 별개의 비열적 집단이 존재하는 명확한 분리를 보여줍니다.
• 비열적 입자 집단은 독립적으로 형성되며, 난류 캐스케이드에서 에너지를 빼가는 주요 에너지 수납체로 작용합니다.
• 열적 성분은 콜모고로프 유사 에너지 스펙트럼을 보이며, 발달된 난류 캐스케이드를 나타냅니다.
• 비열적 집단은 공간적·운동학적으로 열적 성분과 분리되어 있어 별개의 에너지 증가 메커니즘이 존재할 가능성을 시사합니다.
• 프레임워크는 랩톱부터 슈퍼컴퓨터까지 다양한 하드웨어 플랫폼에서 성능과 사용성을 유지하며 성공적으로 확장 가능합니다.