AXUV synthetic diagnostic for ASDEX Upgrade and its application for SPI simulations

기술 요약: ASDEX Upgrade 및 SPI 시뮬레이션을 위한 AXUV 합성 진단법

문제 정의
본 논문은 ASDEX Upgrade(AUG) 토카막의 파쇄 펠릿 주입(Shattered Pellet Injection, SPI) 실험 중 발생하는 AXUV(Absolute eXtended UltraViolet) 다이오드 측정값 해석의 어려움을 다룬다. AXUV 다이오드는 SPI와 같은 빠른 현상을 연구하기에 적합한 높은 시간 분해능($\sim\mu$s)을 제공하지만, 다음 두 가지 주요 요인으로 인해 절대 복사 출력(absolute radiated power)을 결정하는 데 한계가 있다:

1. 스펙트럼 반응성 저하: 극자외선(EUV) 및 진공 자외선(VUV) 노출로 인해 다이오드의 스펙트럼 민감도가 시간이 지남에 따라 저하되며, 이러한 저하는 스펙트럼 범위 전체에서 비균일하게 나타난다.
2. 시간에 따라 변하는 스펙트럼: 혼합 Ne/D$_2$ SPI 실험에서 플라즈마 방출 스펙트럼은 급격히 진화한다(예: 플라즈마 냉각에 따라 브렘스트랄룽(bremsstrahlung)에서 선 방사(line radiation)로 변화). 다이오드의 스펙트럼 반응성이 비균일하기 때문에, 유효 응답($\eta_{eff}$)은 방전 중에 동적으로 변화한다.

특히 네온(Neon)이 주입될 때(다이오드의 민감한 영역에서 강하게 복사함) 발생하는 이러한 불확실성은 절대 출력 재구성 및 토모그래피를 신뢰할 수 없게 만든다. 따라서 저자들은 절대적인 교정(calibration)에 의존하지 않고도 1차 원리 기반의 플라즈마 시뮬레이션과 실험적 AXUV 데이터 사이의 간극을 메울 수 있는 방법을 모색하였다.

방법론
저자들은 플라즈마 시뮬레이션 데이터를 기반으로 AXUV 카메라 신호를 시뮬레이션하기 위해 Cherab-Raysect 광학 모델링 프레임워크 내에서 3D 합성 진단 도구를 개발하였다.

• 프레임워크: 이 도구는 Raysect의 일반적인 가시광선 중심 범위를 넘어, AXUV 다이오드에 필요한 넓은 스펙트럼 범위(1 eV ~ 1.24 keV)를 처리할 수 있도록 조정된 Raysect 광선 추적 라이브러리를 활용한다.
• 기하 구조: 이 진단 도구는 두 개의 토로이달 섹터(SPI가 발생하는 Sector 16 및 약 $110^\circ$ 떨어진 Sector 5)에 위치한 4개의 AXUV 카메라(두 종류의 수직 및 수평 카메라, 총 192개 시선)를 모델링한다. 또한, 텅스텐 표면의 거칠기 매개변수 0.29를 가정하여 스펙트럼에 의존적인 반사를 고려하기 위해 플라즈마 대향 부품(PFC)의 전체 3D CAD 기하 구조를 포함한다.
• 시뮬레이션 입력: 합성 진단 프로세스는 비선형 확장 자기유체역학(MHD) 코드인 JOREK의 출력 데이터를 처리한다. JOREK으로부터 얻은 플라즈마 파라미터(전자 밀도, 온도, 불순물 전하 상태)는 2D 폴로이달 격자(2 cm $\times$ 2 cm 해상도) 상에 슈퍼샘플링(super-sampled)되며, 사전 계산된 다이오드 민감도 행렬과 컨볼루션(convolution)된다.
• 스펙트럼 처리: 계산 비용을 관리하면서 해상도를 유지하기 위해, 스펙트럼을 세 영역(1–10 eV, 10–100 eV, 100–1240 eV)으로 나누고 각 영역당 100개의 빈(bin)을 할당하여 총 300개의 빈을 구성하였다. 모델은 브렘스트랄룽과 네온 선 방사(OPEN-ADAS)를 포함하지만, 펠릿 방사 효과만을 분리하기 위해 시뮬레이션에서 배경 불순물은 제외하였다.
• 검증: 이 도구는 해석적으로 계산된 에텐듀(etendue) 값에 대해 검증되었으며(평균 상대 차이 2.3%), Cherab PowerPipeline Monte Carlo 방식과 교차 확인되었다(평균 상대 오차 0.017%).

주요 기여

1. 3D 합성 진단 도구 개발: 복잡한 기하 구조, 스펙트럼 반응성 저하 및 PFC 반사를 고려하여 3D MHD 시뮬레이션으로부터 합성 AXUV 신호를 생성할 수 있는 견고한 도구를 개발하였다.
2. 방법론적 전환: 저자들은 혼합 Ne/D$_2$ SPI의 경우 스펙트럼 불확실성 때문에 절대 복사 출력 재구성이나 토모그래피를 시도하는 것에 반대한다. 대신, 합성 신호와 실험 신호 사이의 질적 특징 및 시간 진화 패턴을 비교할 것을 제안한다.
3. JOREK 시뮬레이션에 적용: 이 진단 도구는 2022년 SPI 캠페인의 두 가지 특정 AUG 방전 데이터에 적용되었다:
   • #40673: 높은 네온 함량 (10 mol% Ne).
   • #41007: 낮은 네온 함량 (0.17 mol% Ne).

결과
합성 신호는 2022년 SPI 캠페인의 실험적 AXUV 데이터와 비교되었다:

• 높은 네온 함량 (#40673):

  • 일치도: 방사 프런트(radiation front)의 시간 진화 측면에서 특히 좋은 일치도를 보였다. 시뮬레이션에서의 방사 구름의 내측 방향 반경 침투는 $t \approx 2.3285$ s까지 속도와 진폭 모두에서 실험과 일치하였다.
  • 특징: 합성 및 실험 데이터 모두에서 방사 구름으로부터 고자기장 측(high-field side)으로 이동하고 외측 플럭스 표면을 따라 퍼지는 미세한 방사 구조를 보여주었다.
  • 불일치: 수평 카메라 데이터에서 수직 확산의 타이밍과 강도에서 미세한 차이가 관찰되었다. 합성 데이터는 단일한 급격한 확산 이벤트를 보여준 반면, 실험에서는 여러 번의 작은 이벤트가 나타났다. 저자들은 이를 부분적으로 펠릿 파쇄의 확률적 특성(시뮬레이션은 펠릿 크기를 과대평가할 수 있는 Parks 모델을 사용함)과 시뮬레이션에서의 배경 불순물 부재 때문으로 돌린다.
  • 토로이달 확산: Sector 5($\sim110^\circ$ 떨어진 곳)에서 시뮬레이션은 실험에 비해 방사 발생이 지연되었는데, 이는 시뮬레이션이 평형 토로이달 흐름(equilibrium toroidal flow)을 반영하지 못했기 때문으로 보인다.

• 낮은 네온 함량 (#41007):

  • 불일치: 일치도가 현저히 낮았다. 실험에서는 펠릿 도착 직후 미세한 수직 방사 구조가 아래로 드리프트하는 현상이 나타났으나, 시뮬레이션에는 이러한 현상이 없었다.
  • 원인: 저자들은 이를 낮은 네온 펠릿에서 두드러지게 나타나는 **로켓 효과(rocket effect)**와 플라즈모이드 드리프트(plasmoid drift) 때문이라고 분석하였으며, 이는 JOREK 시뮬레이션에 포함되지 않았다.
  • 배경 복사: 실험에서는 배경 불순물 방사(예: 텅스텐)가 펠릿 방사와 맞먹는 수준이었으나, 시뮬레이션(고유 불순물 없음)은 반사에 의한 배경만을 보여주었다.
  • 결론: 드리프트 불일치에도 불구하고, 시뮬레이션과 실험 모두 열에너지가 10%로 떨어질 때까지($t_{0.1W}$) 방사가 폴로이달 방향으로 확산되는 모습을 보였다.

의의 및 주장
본 논문은 개발된 합성 진단 도구가 절대적인 출력 교정에 의존하는 대신, 결여된 물리 모델을 식별함으로써 SPI 시뮬레이션을 실험적 AXUV 데이터와 검증하고 해석할 수 있는 유능한 도구임을 주장한다.

• 물리 검증: 이 도구는 특정 물리 모델의 부재를 성공적으로 탐지하였다. 저네온 사례의 경우, 방사 드리프트 패턴의 불일치는 향후 시뮬레이션에 로켓 효과와 플라즈모이드 드리프트를 포함시켜야 한다는 필요성을 직접적으로 나타냈다.
• 강건성: 본 연구는 채널 민감도와 저하에 있어 1차 규모(order-of-magnitude)의 불확실성이 존재하더라도, 신호 특징(공간 구조 및 시간 진화)의 질적 비교가 여전히 유효하고 유익하다는 것을 입증한다.
• 향후 권고 사항: 분석을 바탕으로 저자들은 향후 SPI 시뮬레이션에 다음 사항을 포함할 것을 권고한다:
  1. 로켓 효과 및 플라즈모이드 드리프트 모델 (특히 미량 네온 시나리오의 경우).
  2. 배경 불순물 방사.
  3. 평형 토로이달 흐름.

본 논문은 이러한 접근 방식이 연구자들이 어떤 물리 모델이 실험적 관측을 재현하는 데 필요한지 평가할 수 있게 함으로써, 결과적으로 ITER를 위한 파괴 완화 전략 개발을 지원할 수 있다고 결론짓는다.