[논문 리뷰] Verification of Gyrokinetic codes: theoretical background and applications

N. Tronko, A. Bottino외 4명|arXiv (Cornell University)|2017. 03. 22.

Magnetic confinement fusion research참고 문헌 28인용 수 20

한 줄 요약

이 논문은 라그랑주 변분 형식과 교차 코드 기준 비교를 융합하여 기수역학 코드 검증을 위한 통합 이론적 및 수치적 프레임워크를 제시한다. ORB5 (PIC) 및 GENE (에일러리안) 코드가 전체 FLR 보정을 포함할 경우 동일한 선형 모델을 구현하고 있음을 입증하며, 전자기적 β-스캔 및 모드 구조 기준 비교에서 뛰어난 일치를 보이며, 이는 수치적 정밀도를 검증하고 미세 구조에 대한 정확한 포isson 해석기의 중요성을 부각시킨다.

ABSTRACT

In fusion plasmas the strong magnetic field allows the fast gyro-motion to be systematically removed from the description of the dynamics, resulting in a considerable model simplification and gain of computational time. Nowadays, the gyrokinetic (GK) codes play a major role in the understanding of the development and the saturation of turbulence and in the prediction of the subsequent transport. Naturally, these codes require thorough verification and validation. Here we present a new and generic theoretical framework and specific numerical applications to test the faithfulness of the implemented models to theory and to verify the domain of applicability of existing GK codes. For a sound verification process, the underlying theoretical GK model and the numerical scheme must be considered at the same time, which has rarely been done and therefore makes this approach pioneering. At the analytical level, the main novelty consists in using advanced mathematical tools such as variational formulation of dynamics for systematization of basic GK code's equations to access the limits of their applicability. The verification of numerical scheme is proposed via the benchmark effort. In this work, specific examples of code verification are presented for two GK codes: the multi-species electromagnetic ORB5 (PIC) and the radially global version of GENE (Eulerian). The proposed methodology can be applied to any existing GK code. We establish a hierarchy of reduced GK Vlasov-Maxwell equations implemented in the ORB5 and GENE codes using the Lagrangian variational formulation. At the computational level, detailed verifications of global electromagnetic test cases developed from the CYCLONE Base Case are considered, including a parametric $\beta$-scan covering the transition from ITG to KBM and the spectral properties at the nominal $\beta$ value.

연구 동기 및 목표

• 기하학적 이론적 모델 검증과 수치적 방법 검증을 통합하는 엄밀한 이중 검증 프레임워크를 구축하는 것.
• 다양한 코드에서 사용된 순서 및 근사치의 명확화를 통해 이론적 기수역학 모델과 코드 구현 간의 괴리를 해결하는 것.
• 전자기적 시험 케이스에 대한 상세한 교차 코드 기준 비교를 통해 ORB5 (PIC) 및 GENE (에일러리안) 코드의 수치적 정밀도를 검증하는 것.
• 특히 긴 파장 포isson 해석기로 인한 수치적 근사치가 전기적 잠재력의 미세 구조적 반경 방향 구조에 미치는 영향을 규명하는 것.
• 특히 유한 β 및 ρ∗ 효과 존재 시에 미세 불안정성 전이를 정확하게 시뮬레이션하기 위해 반경 방향 전역 모델이 필수적임을 보여주는 것.

제안 방법

• 동일한 라그랑주 변분 형식을 사용하여 기수역학 볼츠만-맥스웰 방정식을 유도함으로써 기초 이론적 모델과 그 적용 한계를系통적으로 규명한다.
• 동일한 변분 프레임워크를 사용하여 ORB5 (PIC) 및 GENE (에일러리안) 코드의 축소된 입자 역학 및 맥스웰-볼츠만 모델을 분석 수준에서 비교한다.
• 두 가지 시험 케이스에 대해 교차 코드 기준 비교를 수행한다: 고정된 토로이드 모드 수 n = 19에서의 전역 전자기적 β-스캔, 및 고정된 β에서의 반경 및 편심 모드 구조의 상세 비교.
• β-스캔에서는 선형 성장률과 실수 주파수를 관측량으로 사용하고, 모드 구조 비교에서는 전기적 잠재력 및 평행 자기장 잠재력의 전체 편심 및 반경 프로파일을 사용한다.
• 두 코드 간 일관성을 확보하기 위해 동일한 저-β 순서를 적용하고, GENE의 전체 FLR 해석기와 ORB5의 긴 파장 근사치를 대조한다.
• 전역 및 국소(유량 튜브) GENE 결과를 비교하여 유한한 크기(ρ∗) 효과를 정량화하며, 반경 방향 구조 효과로 인해 모드 전이 임계값이 다름을 보여준다.

실험 결과

연구 질문

• RQ1 ORB5 및 GENE 코드에 구현된 이론적 기수역학 모델이 동일한 변분 프레임워크에서 유도되었을 때, 어느 정도 일치하는가?
• RQ2 수치적 근사치—특히 포isson 해석기—는 전기적 잠재력의 미세 구조적 반경 방향 구조에 어떤 영향을 미치는가?
• RQ3 유한한 반경 크기(ρ∗)는 전역 전자기적 시뮬레이션에서 이온 온도 기울기(이온 온도 기울기, ITG) 모드와 운동학적 볼로우닝 모드(KBM) 간의 전이에 어떤 영향을 미치는가?
• RQ4 전역 전자기적 β-스캔에서 ORB5 및 GENE의 선형 성장률과 주파수는 어떤 범위의 β 값에서 얼마나 잘 일치하는가?
• RQ5 두 코드 간 전기적 잠재력 및 평행 자기장 잠재력의 모드 구조는 어느 정도 다를 수 있으며, 이러한 괴리의 원인은 무엇인가?

주요 결과

• 동일한 라그랑주 변분 형식을 통해 유도된 경우, ORB5 및 GENE의 이론적 모델은 선형 수준에서 동일하며, 기초 물리학의 일관성을 확인한다.
• n = 19에서의 전역 전자기적 β-스캔에서 ORB5 및 GENE 간 뛰어난 일치가 관찰되었으며, 두 코드 모두 β 범위 전반에서 유사한 선형 성장률과 실수 주파수를 보였다.
• 전역 모의와 국소(유량 튜브) 모의 간에 ITG에서 KBM 모드로의 전이가 다른 β 값에서 발생함을 확인하여, 유한한 크기(ρ∗ ≈ 1/182) 효과가 무시할 수 없고 코드에 따라 다름을 시사한다.
• 긴 파장 포isson 해석기를 사용하는 ORB5 코드는 GENE의 전체 FLR 해석기가 포착하는 전기적 잠재력의 미세 구조적 반경 방향 구조를 해석하지 못함을 확인하였으며, 특히 유리 표면 근처에서 두드러진다.
• ORB5에서 누락된 이러한 구조는 낮은 n에서 성장률과 주파수에 큰 영향을 주지 않지만, 높은 토로이드 모드 수(n ≳ 40)에서는 허용할 수 없는 편차로 나타나며, 전체 FLR 해석기가 필요함을 강조한다.
• 기준 비교 결과는 일관된 근사를 사용할 경우 기수역학 모델의 수치적 구현이 이론에 충실함을 확인하였으며, 핵 fusion 플라즈마 시뮬레이션에서의 신뢰성 있는 예측을 위해 교차 코드 검증이 필수적임을 입증한다.