[논문 리뷰] Gyrokinetic Continuum Simulation of Turbulence in Open-Field-Line Plasmas

E. L. Shi|arXiv (Cornell University)|2017. 08. 24.

Magnetic confinement fusion research참고 문헌 184인용 수 12

한 줄 요약

이 논문은 토카막 스케이프오프 레이어에서 개방된 자기장 선에서 난류를 모델링할 수 있는 최초의 기하학적 연속성 시뮬레이션 코드를 제시한다. 에너지 보존 불연속 갈레르킨 방법과 새로운 도전성 쉬트 경계 조건을 사용한다. 주요 기여는 안정적이고 모멘트 보존 알고리즘을 통해 경계 영역에서의 플라즈마 난류, 전자기 효과 및 쉬트 물리 현상을 정확하게 시뮬레이션할 수 있도록 한 것이다. 직선 및 나선형 자기장 선 기하구조에서의 검증을 통해 성능을 입증하였다.

ABSTRACT

The properties of the boundary plasma in a tokamak are now recognized to play a key role in determining the achievable fusion power and the lifetimes of plasma-facing components. Accurate quantitative modeling and improved qualitative understanding of the boundary plasma ultimately require five-dimensional gyrokinetic turbulence simulations, which have been successful in predicting turbulence and transport in the core. The additional challenges of boundary-plasma simulation necessitate the development of new gyrokinetic codes or major modifications to existing core gyrokinetic codes. In this thesis, we develop the first gyrokinetic continuum code capable of simulating plasma turbulence on open magnetic field lines, which is a key feature of a tokamak scrape-off layer. In contrast to prior attempts at this problem, we use an energy-conserving discontinuous Galerkin discretization in space. To model the interaction between the plasma and the wall, we design conducting-sheath boundary conditions that permit local currents into and out of the wall. We start by designing spatially one-dimensional kinetic models of parallel SOL dynamics and solve these systems using novel continuum algorithms. By generalizing these algorithms to higher dimensions and adding a model for collisions, we present results from the first gyrokinetic continuum simulations of turbulence on two types of open-field-line systems. The first simulation features uniform and straight field lines, such as found in some linear plasma devices. The second simulation is of a hypothetical model we developed of the NSTX scrape-off layer featuring helical field lines. These developments comprise a major step towards a gyrokinetic continuum code for quantitative predictions of turbulence and transport in the boundary plasma of magnetic fusion devices.

연구 동기 및 목표

• 토카막 스케이프오프 레이어와 같은 개방된 자기장 선 플라즈마에서 난류를 시뮬레이션할 수 있는 기하학적 연속성 코드 개발.
• 경계 플라즈마 시뮬레이션에서의 큰 진폭의 변동, 전자기 효과, 복잡한 자기 기하학적 구조에 대응하는 데 목표.
• 정확한 장기 시뮬레이션을 위해 에너지 보존 및 모멘트 보존 수치적 방법 구현.
• 벽과의 플라즈마 상호작용을 도전성 쉬트 경계 조건을 통해 모델링하여 지역적 전류 교환을 허용.
• 직선 및 나선형 자기장 선 기하구조에서의 코드 검증 — 선형 장치 및 NSTX 유사 기하구조를 대표함.

제안 방법

• 공간 이산화에 불연속 갈레르킨 유한요소 방법을 사용하여 에너지 보존 및 위상공간에서 고차 정확도 확보.
• 벽과의 상호작용을 모델링하기 위해 지역 전류를 벽으로 주입 및 추출할 수 있는 도전성 쉬트 경계 조건 개발.
• 기존 다항식 기저에 비해 해상도 효율을 향상시키기 위해 속도 공간에서 지수 가중 다항식 도입.
• 분포 함수의 모든 모멘트를 보존하는 수정된 불연속 갈레르킨 알고리즘 제안 — 기존 방법의 결함 보완.
• 일차원 병렬 SOL 동역학의 키네틱 모델을 고차원으로 일반화하고 현실적인 난류 시뮬레이션을 위해 충돌 모델 통합.
• 선형화된 방정식을 사용한 기하학적 연속성 프레임워크를 적용하여 이온 및 전자의 동역학 기술 — 쉬트 전류 효과는 수정된 병행 전류 모델을 통해 통합.

실험 결과

연구 질문

• RQ1 기하학적 연속성 시뮬레이션 코드는 토카막 스케이프오프 레이어의 개방된 자기장 선에서 난류를 정확하게 시뮬레이션할 수 있는가?
• RQ2 도전성 쉬트 경계 조건은 개방된 자기장 선 시스템에서 난류의 안정성과 역학에 어떤 영향을 미치는가?
• RQ3 기존 다항식에 비해 속도 공간 이산화에서 지수 가중 다항식이 얼마나 높은 효율성을 제공하는가?
• RQ4 수정된 불연속 갈레르킨 알고리즘은 분포 함수의 모든 모멘트를 보존하는가? 그리고 이는 장기 시뮬레이션에서 왜 중요한가?
• RQ5 쉬트 전류는 개방된 자기장 선 기하구조에서 난류 모드의 성장률에 어떤 영향을 미치는가?

주요 결과

• 에너지 보존 성질을 지닌 불연속 갈레르킨 방법을 사용하여 개방된 자기장 선에서의 기하학적 연속성 난류 시뮬레이션을 최초로 성공적으로 수행.
• 도전성 쉬트 경계 조건을 통해 실제 벽과의 전류 교환을 모델링할 수 있었으며, 이는 스케이프오프 레이어 물리학에 있어 핵심적이다.
• 수정된 불연속 갈레르킨 알고리즘은 표준 구현 방식과 달리 분포 함수의 모든 모멘트를 성공적으로 보존.
• 속도 공간에서의 지수 가중 다항식은 특히 강한 기울기가 있는 영역에서 분포 함수의 해상도 효율을 뛰어나게 했다.
• 강한 쉬트 전류 조건(εₜ ≫ 1)에서 난류 모드의 성장률은 1/√εₜ의 비율로 감소하여 상당한 안정화 효과를 보였다.
• 직선 및 나선형 자기장 선 기하구조에서의 시뮬레이션을 통해 코드의 강건성과 실제 토카막 경계 구조 적용 가능성 확인.