ASDEX 업그레이드에서 관측된 국소화된 ELM 시뮬레이션

초록

JOREK 비선형 감소 MHD 코드를 이용해 ASDEX Upgrade 토카막의 타입‑I ELM을 재현하였다. 시뮬레이션에서 얻은 토로이달 모드 번호, 폴로이달 필라멘트 크기, 그리고 SOL로의 방사형 전파 속도가 실험과 일치했으며, 최근 보고된 “단일 자기 교란”(solitary magnetic perturbations)과 유사한 국소화된 불안정성을 보여준다. 이러한 현상은 높은 토로이달 해상도를 갖는 계산에서만 포착된다.

상세 분석

본 연구는 비선형 감소 MHD(Reduced‑MHD) 모델을 기반으로 하는 JOREK 코드를 활용해 ASDEX Upgrade 토카막의 가장 흔한 유형인 타입‑I ELM(Edge Localized Mode)의 동역학을 정밀하게 재현하였다. 감소 MHD는 전류와 플라즈마 압력의 변화를 중심으로 전자기와 유체 운동을 결합한 방정식 체계로, 고전압 플라즈마 경계층에서 발생하는 급격한 비선형 현상을 효율적으로 계산한다. JOREK는 3‑D 토로이달 구조를 구현할 수 있는 유한요소 기반 코드이며, 특히 토로이달 방향의 고해상도 모드 전개가 가능하도록 설계되었다.

시뮬레이션 설정은 실제 ASDEX Upgrade 장치의 기하학적 파라미터(주축 반경 R₀≈1.65 m, 플라즈마 전류 Iₚ≈0.8 MA 등)를 그대로 반영했으며, 초기 평형은 EFIT 재구성 데이터를 사용하였다. 토로이달 모드 전개는 n=520 범위의 정수 모드들을 포함했으며, 특히 n≈812 영역에서 가장 큰 성장률을 보였다. 이는 실험에서 보고된 ELM의 토로이달 구조와 일치한다.

폴로이달 방향에서는 필라멘트 폭이 약 2~4 cm, 길이는 전체 폴로이달 각도(≈180°)에 걸쳐 분포하는 형태를 나타냈다. 이러한 필라멘트는 ELM 발발 직후 급격히 팽창하면서 SOL(스크래프‑오프‑레이어)로 진입한다. 전파 속도는 평균 1.5 km/s 정도로 측정되었으며, 이는 실험적 레이더·베이컨 측정값과 거의 동일하다. 특히, 필라멘트가 SOL에 진입하면서 발생하는 전류 단락 현상은 “solitary magnetic perturbations”(SMP)이라 불리는 국소화된 자기 교란 신호와 형태·시점이 일치한다.

핵심적인 발견은 이러한 SMP와 유사한 국소화된 불안정성이 높은 토로이달 해상도(Δφ≈2° 이하)에서만 재현된다는 점이다. 낮은 해상도(Δφ≈10°)에서는 전반적인 ELM 전개는 잡히지만, 국소적인 자기 교란이 흐릿해져 실험과의 정량적 차이가 발생한다. 이는 토로이달 방향의 비선형 상호작용이 ELM 초기 단계에서 중요한 역할을 함을 시사한다.

또한, 시뮬레이션은 ELM 전개 과정에서 압력 피크가 급격히 감소하고, 전류 층이 재구성되면서 플라즈마 경계가 불안정하게 변하는 메커니즘을 상세히 보여준다. 이러한 과정은 기존의 선형 안정성 분석으로는 설명하기 어려운 비선형 급변 현상을 포착한다는 점에서 의미가 크다.

한계점으로는 감소 MHD 모델이 전자 온도와 이온 온도 차이를 무시하고, 전자 전도성 및 비등방성 효과를 단순화한다는 점이다. 따라서 고에너지 입자 손실이나 중성 입자 재결합 효과는 별도 모델링이 필요하다. 또한, 시뮬레이션에 사용된 초기 잡음 수준이 결과에 민감하게 작용할 수 있어, 실제 플라즈마의 잡음 스펙트럼을 보다 정밀히 반영할 필요가 있다.

요약하면, 본 연구는 고해상도 토로이달 전개를 통해 ASDEX Upgrade에서 관측된 타입‑I ELM의 핵심 특성(모드 번호, 필라멘트 크기, 전파 속도)과 최신 실험에서 보고된 SMP 현상을 성공적으로 재현했으며, 향후 실시간 ELM 제어 및 예측 모델 개발에 중요한 물리적 인사이트를 제공한다.