ELM 스펙트럼 비정상 분석으로 본 플라즈마 변동의 새로운 통찰

본 논문은 토카막 플라즈마에서 발생하는 엣지 로컬라이즈드 모드(ELM)의 비정상 스펙트럼 특성을 정밀하게 분석하기 위해 최신 신호 처리 기법들을 적용한 연구이다. 실험은 TFTR(토카막 플라즈마 테스트 리액터)에서 H‑mode 방전 #49035를 대상으로 진행되었으며, 마이크로파 산란 진단을 이용해 플라즈마 가장자리(반경 r≈1.0 cm)에서의 전자 밀도 변동을 측정하였다. 두 개의 검출기(θ‑파수 kθ=3.3 cm⁻¹와 8.5 cm⁻¹)를 사용해 전자와 이온 드리프트 방향으로 전파되는 변동을 구분하고, 5 MHz 샘플링, 524 288 포인트(≈0.1 s)의 데이터를 확보하였다. 1. **정상 스펙트럼 추정** 다중 테이퍼(Multitaper) 방법을 적용해 20개의 직교 테이퍼와 20 kHz 대역폭 커널 스무딩을 수행, 전체(정상+비정상) 스펙트럼을 추정하였다. 전자 드리프트 쪽(주파수 >1 MHz)에서는 300–700 kHz 대역의 넓은 밴드가 지배적이며, 이는 플라즈마 가장자리 근처에 음전기장이 존재함을 의미한다. 반면 이온 드리프트 쪽(주파수 <1 MHz)에서는 급격히 감소하는 지수형 스펙트럼이 관측되었다. 두 검출기의 스펙트럼을 비교하면, 전자 스펙트럼은 검출기 #2(큰 kθ)에서 약 3배 강하게 나타나며, 이는 파수에 따라 전자 난류가 강해짐을 시사한다. 2. **비정상 성분 강조 – 자동회귀(AR) 필터** 비정상 신호(ELM 및 프리커서)를 강조하기 위해 AR(p) 필터(p≈10)를 적용하였다. AR 모델은 기존 시계열을 과거 샘플들의 선형 결합으로 표현하고, 잔차(z_t)를 혁신 시퀀스로 추출한다. 이 과정에서 1 MHz IF 피크와 같은 정적 배경이 크게 감소하고, 폭넓은 스펙트럼을 가진 ELM과 고주파 프리커서가 잔차에 두드러지게 나타난다. 필터링 후 잔차는 ELM 발생 시 약 60 포인트(≈11 µs) 동안 급격히 상승하고, 1 ms 전후에 10 µs 정도 지속되는 250–650 kHz 버스트가 두드러진다. 3. **시간‑주파수 스펙트로그램 및 특이값 분해(SVD)** 멀티테이퍼 스펙트로그램을 계산하고, 이를 시간‑주파수 행렬로 전개한 뒤 특이값 분해를 수행하였다. 첫 번째 특이벡터는 ELM 동안의 스펙트럼을 포착하며, 정상 상태 스펙트럼에 비해 더 넓고 대칭적인 형태를 보인다. 두 번째 특이벡터는 배경 잡음 및 잔여 변동을 나타낸다. 이 분석을 통해 ELM이 발생할 때 스펙트럼이 급격히 넓어지고, 전자와 이온 드리프트 방향 모두에서 에너지 분포가 변함을 정량화하였다. 4. **ELM 사이 스펙트럼 진화** ELM 종료 직후 전자 드리프트 스펙트럼은 0.5 ms 동안 50–80 % 상승한 뒤 포화한다. 반면 이온 드리프트 스펙트럼은 5 ms에 걸쳐 지수적으로 8배 증가한다. 이 지수 성장률은 ELM 발생 시점과 강하게 연관되어, 이온 스펙트럼이 일정 임계값에 도달하면 ELM이 촉발되는 메커니즘을 제시한다. 5. **고주파 프리커서 탐지** 고해상도 사전백색화(prewhitened) 멀티테이퍼 진화 스펙트럼을 이용해 프리커서를 탐지하였다. 500 kHz 부근에서 약 10 µs 지속되는 짧은 버스트가 ELM 1 ms 전부터 반복적으로 나타났으며, 이는 이온 드리프트 방향으로 전파된다. 두 번째 검출기에서는 이러한 프리커서가 관측되지 않아, 프리커서 발생이 파수와 측정 위치에 민감함을 보여준다. 프리커서가 나타난 순간을 기준으로 ELM 발생을 90 % 이상 정확도로 예측할 수 있음을 보고하였다. 6. **결론 및 의의** 자동회귀 필터와 멀티테이퍼 기반 SVD 시간‑주파수 분석은 플라즈마 비정상 변동을 고해상도로 분리·정량화하는 강력한 도구임을 입증한다. 특히, 짧은 고주파 프리커서는 ELM 예측에 실용적인 신호로 활용될 가능성을 제시한다. 향후 연구에서는 프리커서 메커니즘을 이론적으로 규명하고, 실시간 제어 시스템에 적용하여 ELM에 의한 플라즈마 손상을 최소화하는 방안을 모색할 수 있다.