[논문 리뷰] Electromagnetic gyrokinetic instabilities in STEP

Daniel Kennedy, M. Giacomin외 5명|arXiv (Cornell University)|2023. 07. 04.

Magnetic confinement fusion researchPhysics and Astronomy인용 수 3

한 줄 요약

이 연구는 STEP 구형 토카막 설계에서 전자기적 마이크로불안정성에 대한 선형 기하학적 기반 분석을 제시하며, 이온 기울기와 갇힌 전자기인으로 인해 유도되는 하이브리드 기체 복합파동 모드(KBM)가 이온 이편방향 척도에서 지배적임을 규명한다. 충돌성 미세갈라짐 모드(MTMs)는 보조적이나 유의미하며, 고유한 균형 척도를 가진다. 결과는 다중 코드 기준 비교 및 해상도 수렴성 검증을 통해 검증되어, 후속 비선형 연구를 위한 신뢰를 확보하였다. 동반 논문에서 후속 비선형 연구를 위한 기반을 마련한다.

ABSTRACT

We present herein the results of a linear gyrokinetic analysis of electromagnetic microinstabilites in the conceptual high-$β$, reactor-scale, tight-aspect-ratio tokamak STEP (Spherical Tokamak for Energy Production, https://step.ukaea.uk). We examine a range of flux surfaces between the deep core and the pedestal top for two candidate flat-top operating points of the prototype device. Local linear gyrokinetic analysis is performed to determine the type of microinstabilities that arise under these reactor-relevant conditions. We find that the equilibria are dominated at ion binormal scales by a hybrid version of the Kinetic Ballooning Mode (KBM) instability that has significant linear drive contributions from the ion temperature gradient and from trapped electrons, while collisional Microtearing Modes (MTMs) are sub-dominantly also unstable at similar binormal scales. The hybrid-KBM and MTM exhibit very different radial scales. We study the sensitivity of these instabilities to physics parameters, and discuss potential mechanisms for mitigating them. The results of this investigation are compared to a small set of similar conceptual reactor designs in the literature. A detailed benchmark of the linear results is performed using three gyrokinetic codes; alongside extensive resolution testing and sensitivity to numerical parameters providing confidence in the results of our calculations, and paving the way for detailed nonlinear studies in a companion article.

연구 동기 및 목표

• 고-β, 반응로 규모, 높은 비율의 구형 토카막(STEP) 평형 상태에서 전자기적 마이크로불안정성 특성을 평가하기 위해.
• 핵심에서 페달까지의 플럭스 표면에 걸쳐 두 가지 후보 STEP 운영 조건에서 지배적 및 보조적 불안정성을 규명하기 위해.
• 미래 비선형 시뮬레이션을 위한 다수의 기하학적 기반 코드를 사용한 수치 수렴성 및 해상도 요구 조건 평가하기 위해.
• 세 종류의 기하학적 기반 코드(GS2, CGYRO, GENE) 간 결과 기준 비교 및 선형 안정성 예측 검증을 위해.
• 후속 논문(Paper II)에서 비선형 난류 시뮬레이션을 위한 기반을 마련하기 위해 해상도 및 매개변수 민감도를 확립하기 위해.

제안 방법

• 두 후보 STEP 평형 상태인 STEP-EC-HD 및 STEP-EB-HD에 대해 GS2 코드를 사용하여 局부 선형 기하학적 기반 분석을 수행한다.
• 핵심에서 페달 영역까지를 포함하는 다수의 플럭스 표면(Ψₙ = 0.49에서 0.95까지)에 걸쳐 시뮬레이션을 수행한다.
• 선형 성장률 및 모드 주파수의 교차 검증 및 기준 비교를 위해 세 종류의 기하학적 기반 코드(GS2, CGYRO, GENE)를 활용한다.
• 수치적 안정성을 확보하기 위해, 균형, 병렬(볼로노이 공간), 속도 공간 그리드에 대한 엄격한 해상도 수렴성 테스트를 수행한다.
• 충돌 모델(Lorentz, Sugama, 및 완전 선형화 Fokker-Planck)과 입자 종류 수(2–5)에 대한 민감도를 분석하여 성장률 안정성과의 관계를 조사한다.
• 비선형 시뮬레이션을 위한 최소 해상도 요구 조건을 도출하기 위해 세부 매개변수 스캔(nθ, nperiod, npass, nε)을 수행한다.

실험 결과

연구 질문

• RQ1 고-β, 반응로 관련 STEP 평형 상태에서 이온 이편방향 척도에서 지배적인 전자기적 마이크로불안정성은 무엇인가?
• RQ2 플럭스 표면 계층에서 하이브리드 기체 복합파동 모드(KBM)와 충돌성 미세갈라짐 모드(MTM)의 성장률과 균형 척도를 비교하면 어떻게 되는가?
• RQ3 비선형 시뮬레이션 준비를 위한 정확한 선형 안정성 예측을 위해 필요한 수치적 해상도는 어느 정도인가?
• RQ4 선형 성장률과 모드 주파수는 충돌 연산자 및 플라즈마 입자 종류 수의 선택에 얼마나 민감한가?
• RQ5 GS2, CGYRO, GENE의 다중 코드 기준 비교는 불안정성 예측의 강건성을 어느 정도 확인하는가?

주요 결과

• 평형 상태는 이온 이편방향 척도(k_yρ_i ≲ 1)에서 이온 기울기와 갇힌 전자기인으로 인한 선형 구동력이 큰 하이브리드 기체 복합파동 모드(KBM)에 의해 지배된다.
• 충돌성 미세갈라짐 모드(MTMs)는 보조적이지만 유사한 이편방향 척도에서 불안정하며, 하이브리드-KBM보다 훨씬 넓은 균형 척도를 가진다.
• 낮은 k_y 모드의 경우 nθ ≥ 32에서 해상도 수렴이 달성되나, k_yρ_s > 0.4인 경우 nθ ≥ 64가 필요하여 고파동수 모드에 더 높은 해상도가 요구됨을 시사한다.
• 속도 공간 해상도는 k_y 스펙트럼 전반에서 성장률 및 주파수에 대해 약한 의존성을 보이며, npass 및 nε의 테스트 범위에서 수렴이 관찰된다.
• 선형 스펙트럼은 충돌 모델 선택에 대해 약한 민감도를 보이며, CGYRO에서 Sugama 연산자를 사용할 경우 성장률이 Lorentz 또는 완전 Fokker-Planck 모델 대비 약 20% 낮게 나타난다.
• GS2, CGYRO, GENE 간의 다중 코드 기준 비교를 통해 성장률 및 모드 주파수에 대해 강력한 일치가 확인되어, 비선형 시뮬레이션에 활용하기 위한 결과의 검증이 완료되었다.