[논문 리뷰] Linear MHD stability studies with the STARWALL code

P. Merkel, E. Strumberger|arXiv (Cornell University)|2015. 08. 20.

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

한 줄 요약

이 논문은 3차원 이상의 이상적이고 저항성 등가 구조를 갖는 복잡한 다중 연결 구조를 가진 복잡한 도전성 물체를 고려한 선형 MHD 안정성 분석을 위한 수치 도구인 STARWALL 코드를 제시한다. 이는 OPTIM과의 연계를 통해 저항성 벽 모드(RWMs)의 견고한 피드백 안정화를 가능하게 하며, 축대칭 모드를 정확하게 모델링하기 위해 Lüst-Martensen 항을 포함하고 있다. ASDEX Upgrade 및 ITER 구성에 대한 시뮬레이션을 통해 고유모드 및 전류 반응 분석을 통해 그 능력을 입증하고 있다.

ABSTRACT

The STARWALL/CAS3D/OPTIM code package is a powerful tool to study the linear MHD stability of 3D, ideal equilibria in the presence of multiply-connected ideal and/or resistive conducting structures, and their feedback stabilization by external currents. Robust feedback stabilization of resistive wall modes can be modelled with the OPTIM code. Resistive MHD studies are possible combining STARWALL with the linear, resistive 2D CASTOR code as well as nonlinear MHD simulations combining STARWALL with the JOREK code. In the present paper, a detailed description of the STARWALL code is given and some of its applications are presented to demonstrate the methods used. Conducting structures are treated in the thin wall approximation and depending on their complexity they are discretized by a spectral method or by triangular finite elements. As an example, a configuration is considered consisting of an ideal plasma surrounded by a vacuum domain containing a resistive wall and bounded by an external wall. Ideal linear MHD modes and resistive wall modes in the presence of multiply-connected walls are studied. In order to treat the vertical mode self-consistently the STARWALL code has been completed by adding the so-called Luest-Martensen term generated by a constant normal displacement of the plasma. The appendix contains the computation of the 2D Fourier transform of singular inductance integrals, and the derivation of an asymptotic expansion for large Fourier harmonics.

연구 동기 및 목표

• 3차원 토카막 구성에서 복잡한 도전성 구조를 고려한 선형 MHD 안정성 연구를 위한 융통성 있고 고정밀 수치 프레임워크 개발
• OPTIM 코드와의 연계를 통해 STARWALL을 활용한 저항성 벽 모드(RWMs)의 견고한 피드백 안정화 기능 제공
• 축대칭 불안정성에 대한 수치 모델링을 향상시키기 위해 Lüst-Martensen 항을 도입하여 축대칭 모드의 안정성 처리 확장
• JOREK 및 CASTOR3D와 같은 비선형 및 저항성 MHD 코드와의 연계를 통해 종합적인 안정성 분석 지원
• 다음 세대 fusion 장치인 ITER와 같은 장치에서 RWMs 및 수직 이탈 사건(VDEs)을 연구하기 위한 계산 기반 제공

제안 방법

• STARWALL은 박막 벽 근사법을 사용하며, 도전성 구조를 스펙트럼 방법 또는 삼각형 유한요소로 이산화한다.
• 이론적 MHD 에너지 원리를 변분 프레임워크 내에서 해결하며, 경계 적분 방법을 통해 진공 및 벽 기여를 통합한다.
• 이 코드는 3차원 등가 구조의 고유모드와 성장률을 계산하며, 선형 저항성 2차원 CASTOR 코드와의 연계를 통해 저항성 벽 효과를 포함한다.
• 축대칭(n=0) 모드에 대해 일정한 법선 플라즈마 이격을 갖는 진공 반응을 모델링하기 위해 Lüst-Martensen 항을 구현한다.
• 피드백 안정화는 OPTIM 코드를 사용하여 모델링되며, 이는 불안정 모드와 안정 모드를 모두 고려하여 코ils의 전압을 최적화한다.
• STARWALL은 CAS3D(이론적 MHD), NEMEC(3차원 등가 구조), JOREK(비선형 MHD)를 포함한 더 큰 코드 세트에 통합되어 다중 스케일 시뮬레이션을 가능하게 한다.

실험 결과

연구 질문

• RQ1 3차원 등가 구조에 다중 연결 저항성 및 이상적 도전성 벽이 존재할 경우 선형 MHD 안정성을 어떻게 정확하게 계산할 수 있는가?
• RQ2 벽 기하학적 복잡성이 저항성 벽 모드(RWMs)의 성장률과 구조에 어떤 영향을 미치는가?
• RQ3 Lüst-Martensen 항의 포함이 축대칭 구성에서 수직 모드 안정성 모델링에 어떤 개선을 이끌어내는가?
• RQ4 실제 토카막 벽 구성에서 피드백 제어가 RWMs를 얼마나 효과적으로 안정화할 수 있으며, 최적 제어는 어떻게 달성되는가?
• RQ5 STARWALL이 JOREK, CASTOR3D와 같은 비선형 및 저항성 MHD 코드와 연계될 경우, 붕괴 및 VDE 예측 능력은 어떻게 향상되는가?

주요 결과

• STARWALL은 저항성 및 이상적 벽을 갖는 3차원 등가 구조의 고유모드와 성장률을 성공적으로 계산하였으며, 분석 기준과 양호한 일치를 보였다.
• Lüst-Martensen 항의 포함으로 축대칭 수직 모드를 정확하게 모델링할 수 있었으며, ASDEX Upgrade 구성에서 n=0일 경우 성장률 26 s⁻¹, n=1일 경우 67 s⁻¹, n=2일 경우 314 s⁻¹를 기록하였다.
• 다중 조화 모드(n=0,1,2)에 대해 비대칭 벽 기하학으로 인해 표면 전류가 명확히 결합되고 분리되는 경향을 보였으며, 고유값은 단일 n 해와 유일하게 관련되어 있었다.
• CAS3D에서 비균일 간격의 축방향 격자 구조는 합리적인 점 수(931점)로도 유비의 표면(q=3/2 등) 근처에서 고해상도를 확보하였으며, 이는 균일 격자 구조가 약 40,000점이 필요로 하는 것에 비해 뛰어난 성능을 보였다.
• OPTIM을 통한 RWM 피드백 안정화 모델링은 안정 모드의 의도하지 않은 자극을 방지하면서도 견고하게 수행되었으며, 코ils 전압 최적화를 가능하게 하였다.
• JOREK와의 연계를 통해 VDE, 붕괴 및 QH-모드의 비선형 연구가 가능해졌으며, 초기 결과는 저항성 벽 및 홀로 전류 효과에서 일관된 행동을 보였다.