[논문 리뷰] Electrostatic gyrokinetic simulation of global tokamak boundary plasma and the generation of nonlinear intermittent turbulence

Seung Hoe Ku, R.M. Churchill외 6명|arXiv (Cornell University)|2017. 01. 20.

Magnetic confinement fusion research인용 수 3

한 줄 요약

이 논문은 토카막 경계 플라즈마의 저배수 조건을 모의하기 위해 XGC1 코드를 사용하여 전자기적 기여도 기반의 기하학적 기반 시뮬레이션 프레임워크를 제시한다. 이는 척도 분리 없이 기여도 기반 이온과 드리프트 기반 전자를 포함한다. 비선형 간헐적 난류, 특히 블롭 유사한 구조가 자기 분리면 근처와 스케이프오프층에서 운동학적 미세난류에 의해 자연스럽게 발생함을 입증하여, 경계 플라즈마 물리학에 대한 모델의 정밀도를 검증한다.

ABSTRACT

Boundary plasma physics plays an important role in tokamak confinement, but is difficult to simulate in a gyrokinetic code due to the scale-inseparable nonlocal multi-physics in magnetic separatrix and open magnetic field geometry. Neutral particles are also an important part of the boundary plasma physics. In the present paper, noble electrostatic gyrokinetic techniques to simulate the flux-driven, low-beta electrostatic boundary plasma is reported. Gyrokinetic ions and drift-kinetic electrons are utilized without scale-separation between the neoclassical and turbulence dynamics. It is found that the nonlinear intermittent turbulence is a natural gyrokinetic phenomenon in the boundary plasma in the vicinity of the magnetic separatrix surface and in the scrape-off layer.

연구 동기 및 목표

• 스케일 분리가 불가능하고 국소적이지 않으며 다중 물리 현상이 복합된 난이도 있는 조건에서 토카막 경계 플라즈마를 위한 고정밀 기여도 기반 시뮬레이션 방법을 개발하기 위해.
• 기존의 핵심 영역 중심 기반 기여도 기반 코드의 한계를 극복하기 위해 자기 분리면을 넘어 스케이프오프층까지 실제 경계 조건을 반영하여 시뮬레이션을 가능하게 하기 위해.
• 비철도적 및 난류 동역학 간의 척도 분리 없이 기여도 기반 전자 및 이온을 사용하여 경계 플라즈마를 모델링함으로써 비맥스웰 분포 함수를 포착하기 위해.
• 경계 영역에서 비선형 간헐적 난류, 특히 블롭 구조를 포함한 구조의 자발적 생성을 조사하기 위해.
• 실험적으로 관측된 난류 특성, 예를 들어 블롭과 스트리머를 시뮬레이션 틀을 통해 자연스럽게 발생시켜 검증하기 위해.

제안 방법

• 기여도 기반 이온에 대해 5차원 전자기적 기여도 기반 보르츠만 방정식을 해결하기 위해 XGC1 코드를 사용하며, 계산 비용을 줄이기 위해 드리프트 기반 전자를 적용한다.
• 마커 입자와 위상공간 격자를 조합한 하이브리드 라그랑주 방식의 총-δf 기법을 사용하여 비맥스웰 분포 함수를 정확하게 해석한다.
• 라그랑주 운동 방정식을 사용하여 기여도 평균 전기장과 이온 기여도 중심의 운동을 계산하며, 자기 곡률 및 축방향 가속도 효과를 포함한다.
• 준중성 조건과 기여도 평균 포isson 방정식을 해결하여 전기적 위치 에너지를 결정하며, 이온 밀도는 기여도 중심 분포에서 유도한다.
• 시뮬레이션 영역은 X-점 기하학, 자기 분리면, 플라즈마 축 및 벽까지 포함되어 있어 글로벌 경계 조건을 가능하게 한다.
• 안정성과 정확도를 유지하면서도 빠른 전자 동역학을 해결하기 위해 하위주기 전자 기법을 사용한다.

실험 결과

연구 질문

• RQ1 비철도적 및 난류 동역학 간의 척도 분리 없이, 자기 분리면과 스케이프오프층을 포함한 전체 경계 플라즈마를 정확하게 모의할 수 있는 기여도 기반 시뮬레이션 프레임워크가 존재하는가?
• RQ2 저배수, 플럭스 주도 조건 하에서 기여도 기반 전자 및 이온은 어떻게 상호작용하여 경계 플라즈마에서 비선형 간헐적 난류를 생성하는가?
• RQ3 자기 분리면 근처에서 미세난류에 의해 블롭 유사한 구조가 시뮬레이션에서 자연스럽게 발생하는 정도는 어느 정도인가?
• RQ4 폴로이드 방향 ExB 유동이 경계 영역에서 블롭 운동과 물질 이동을 어떻게 유도하는가?
• RQ5 난류 강도와 구조는 경계 플라즈마 전반에서 어떻게 변화하는가? 특히 분리면 근처와 스케이프오프층에서 어떻게 변화하는가?

주요 결과

• 비선형 간헐적 난류, 특히 블롭 유사한 구조가 자기 분리면 근처와 스케이프오프층에서 기여도 기반 미세난류와 기여도 기반 전자에 의해 자연스럽게 발생한다.
• 난류 강도는 분리면에서 최고에 도달하며, 밀도 기울기가 강해지는 ψn = 0.96 부근에서 급격히 증가하기 시작한다.
• 블롭은 약 7cm의 폴로이드 크기와 1.5cm의 축방향 폭을 가지며, 밀도 변동은 평형 밀도의 10–20%에 이른다.
• 블롭은 약 15km/s의 속도로 폴로이드 방향으로 이동하며, 지역 폴로이드 ExB 유동과 매우 유사하여 이 구성에서 ExB 유동이 블롭 운동을 지배함을 시사한다.
• 블롭의 축방향 운동은 폴로이드 운동보다 훨씬 작아, 고ExB 유동 영역에서 물질 이동이 주로 폴로이드 유동에 의해 주도됨을 나타낸다.
• 시뮬레이션 시간 1.7ms 이후 경계 영역에서 안정 상태에 도달하였으며, 경계층에서 안정 상태의 열 이동이 달성되었다.