[논문 리뷰] Magnetic Fluctuations in Gyrokinetic Simulations of Tokamak Scrape-Off Layer Turbulence

Noah Mandell|arXiv (Cornell University)|2021. 03. 30.

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

한 줄 요약

이 학위논문은 토카막 스케일오프층(SOL)에서 개방된 자기장 선과 강한 자기장 변동이 존재하는 조건에서, 에너지 보존 및 앤퍼 법칙의 상쇄 문제를 방지하는 이산 갈레르킨(DG) 방법을 사용하여, 처음으로 전체-풀(fully-f) 전자기기구역 운동역학 시뮬레이션 능력을 개발한다. 이는 전자기적 불안정성과(blob 동역학 및 운반)의 영향을 안정적으로 시뮬레이션할 수 있으며, NSTX 유사 SOL 조건에서의 결과를 보여준다.

ABSTRACT

Understanding turbulent transport physics in the tokamak edge and scrape-off layer (SOL) is critical to developing a successful fusion reactor. The dynamics in these regions plays a key role in achieving high fusion performance by determining the edge pedestal that suppresses turbulence in H-mode. Additionally, the survivability of a reactor is set by the heat load to the vessel walls, making it important to understand turbulent spreading of heat as it flows along open magnetic field lines in the SOL. Large-amplitude fluctuations, magnetic X-point geometry, and plasma interactions with material walls make simulating turbulence in the edge/SOL more challenging than in the core region, necessitating specialized gyrokinetic codes. Further, the inclusion of electromagnetic effects in gyrokinetic simulations that can handle the unique challenges of the boundary plasma is critical to the understanding of phenomena such as the pedestal and ELMs, for which electromagnetic dynamics are expected to be important. In this thesis, we develop the first capability to simulate electromagnetic gyrokinetic turbulence on open magnetic field lines. This is an important step towards comprehensive electromagnetic gyrokinetic simulations of the coupled edge/SOL system. By using a continuum full-f approach via an energy-conserving discontinuous Galerkin (DG) discretization scheme that avoids the Ampere cancellation problem, we show that electromagnetic fluctuations can be handled in a robust, stable, and efficient manner in the gyrokinetic module of the Gkeyll code. We then present results which roughly model the scrape-off layer of the National Spherical Torus Experiment (NSTX), and show that electromagnetic effects can affect blob dynamics and transport. We also formulate the gyrokinetic system in field-aligned coordinates for modeling realistic edge and scrape-off layer geometries in experiments.

연구 동기 및 목표

• 개방된 자기장 선과 강한 자기장 변동이 존재하는 토카막 에지 및 스케일오프층(SOL)에서, 안정적이고 에너지 보존이 이루어지는 전자기기구역 운동역학 시뮬레이션 능력을 개발한다.
• 연속체 전체-풀(fully-f) 이산 갈레르킨(DG) 이산화 방법을 도입하여, 전자기기구역 운동역학 시뮬레이션에서 알려진 불안정성인 앤퍼 법칙의 상쇄 문제를 해결한다.
• 현실적인 SOL 기하구조에서 전자기적 효과가 난류 운반과 blob 동역학에 미치는 영향을 정확하게 모델링할 수 있도록 하며, 이는 페달 형성과 에지 국소화 모드(ELMs) 이해에 필수적이다.
• 현장에 따라 정렬된 좌표계를 사용하여 실험적 에지 플라즈마 구성의 더 나은 표현을 가능하게 한다.
• 분포 함수의 양수성 보존과 함께 질량, 운동량, 에너지 보존 법칙을 유지하는 새로운 DG 알고리즘 개발을 통해 수치적 안정성을 확보한다.

제안 방법

• 개방된 자기장 선에서 전자기적 변동을 다룰 수 있도록, 연속체 전체-풀(fully-f) 접근법과 이산 갈레르킨(DG) 공간 이산화를 사용한다.
• 앤퍼 법칙의 상쇄 문제를 피하는 에너지 보존 DG 방법을 구현하여, 전자기기구역 운동역학 코드에서 알려진 불안정성을 제거한다.
• 현실적인 에지 및 SOL 플라즈마 기하구조(예: X점 구성 및 플라즈마-벽 상호작용 포함)를 모델링하기 위해 현장에 따라 정렬된 좌표계를 사용한다.
• 전자기적 효과가 존재하는 상황에서도 질량, 운동량, 에너지를 보존하면서 분포 함수의 양수성을 유지하는 새로운 알고리즘을 개발한다.
• 단순화된 NSTX 유사 스케일오프층 구조를 대상으로 시뮬레이션을 수행하여 전자기적 효과가 blob 동역학과 운반에 미치는 영향을 연구한다.
• 충돌 연산자를 단순화된 Dougherty 모델로 통합하고, 향후 비선형 Fokker-Planck 모델 및 GPU 가속화로의 확장 가능성을 탐색한다.

실험 결과

연구 질문

• RQ1 토카막 에지 및 스케일오프층(SOL)에서 개방된 자기장 선과 큰 변동이 존재하는 조건에서, 전자기적 효과를 안정적이고 안정적으로 시뮬레이션할 수 있는 방법은 무엇인가?
• RQ2 실제 NSTX 유사 조건에서 전자기적 변동은 SOL 내에서 blob 동역학과 횡방향 운반에 어느 정도의 영향을 미치는가?
• RQ3 에너지 보존 및 양수성 보존을 갖춘 이산 갈레르킨(DG) 방법을 통해 개방된 자기장 선에서 안정적인 전체-풀(fully-f) 전자기기구역 운동역학 시뮬레이션이 가능할 수 있는가?
• RQ4 전자기적 효과는 페달 운반과 에지 국소화 모드(ELMs)에서 어떤 역할을 하는가? 그리고 페달과 SOL 간의 상호작용에 어떻게 영향을 미치는가?
• RQ5 수치적 프레임워크는 어떻게 확장되어 현실적인 충돌 연산자와 GPU 가속화를 포함할 수 있는가?

주요 결과

• 이산 갈레르킨(DG) 방법을 사용한 전체-풀(fully-f) 전자기기구역 운동역학 프레임워크는 앤퍼 법칙의 상쇄 문제를 효과적으로 방지하고, 개방된 자기장 선에서 안정적이고 에너지 보존이 이루어지는 시뮬레이션을 가능하게 했다.
• 시뮬레이션된 NSTX 유사 스케일오프층에서 전자기적 효과가 blob 동역학과 난류 운반에 상당한 영향을 미치는 것으로 나타났으며, 이는 에지 플라즈마 거동에서의 중요성을 시사한다.
• 양수성 보존 DG 알고리즘은 강한 전자기적 변동이 존재하는 상황에서도 질량, 운동량, 에너지를 보존하면서 수치적 안정성을 유지한다.
• 현장에 따라 정렬된 좌표계의 구현은 X점 구성과 플라즈마-벽 상호작용을 포함한 현실적인 에지 플라즈마 기하구조를 정확하게 모델링할 수 있게 하였다.
• 이 프레임워크는 고성능 계산을 위해 더 고급 충돌 모델과 GPU 가속화로의 확장이 가능하며, 초기 단계에서 NVIDIA GPU용 CUDA로 소자 커널을 이식하는 데 성공하였다.
• 이 능력은 고 confinement 모드(H-mode) 성능과 fusion 반응기의 전력 배출 문제를 이해하는 데 필수적인, 자발적 상호작용을 갖는 페달/SOL 시스템의 통합 연구를 가능하게 한다.