회전 Z 핀치에서 폭발적 펄스와 MHD 페디스털 주기

초록

전단 흐름이 자기장선과 평행하게 존재할 때, 급격한 압력 구배를 선형적으로 안정화시켜 Z‑핀치 플라즈마 내부에 “MHD 페디스털”이라 부르는 전이 장벽을 만든다. 이 장벽은 선형적으로는 안정하지만, 충분히 큰 교환 변위에 대해 메타스테이블해 급격히 붕괴한다. 저자들은 히팅과 전단 토크를 서서히 가함으로써 페디스털이 형성·붕괴·재구축을 반복하는 사이클을 수치 시뮬레이션으로 재현하고, 붕괴 시 방출되는 에너지와 최적 플럭스 튜브 교환을 조합 최적화 문제로 정량화한다. 이는 토카막의 ELM 현상과 직접적인 유사성을 가진다.

상세 분석

본 논문은 β≪1인 Z‑핀치 플라즈마에 전단 흐름이 존재할 경우, 총압 P=p+B²/2의 급격한 구배를 전단에 의해 생성된 원심력(ℓ²/χ)과 비교해 안정화할 수 있음을 수식(5)로 증명한다. 여기서 ℓ은 각운동량, χ는 비례 자기 플럭스이며, 전단 흐름이 초음속(Maϕ≈1) 수준에 이르면 ℓ²/χ의 급격한 상승이 s/χ(=p^{1/γ}/ρ) 구배를 역전시켜 “MHD 페디스털”을 형성한다. 선형 안정성 조건 L>0이 만족되면 작은 변위는 복원력을 갖지만, 비선형적으로는 식(4)의 s/χ 항이 ℓ²/χ 항을 압도해 결국 포텐셜 에너지 장벽을 넘어서는 교환 변위가 가능해진다. 이는 γ<2인 경우 메타스테이블성을 초래한다는 점에서 핵심이다.

시뮬레이션에서는 히팅(R_heat)과 azimuthal drag(τ_drag,ϕ)를 시간 의존적으로 적용해 두 가지 모드(L‑mode, H‑mode)를 구현한다. L‑mode에서는 전단 토크가 없으므로 교환 난류가 s/χ를 거의 일정하게 만들고 페디스털이 형성되지 않는다. 반면 H‑mode에서는 τ_drag,ϕ≈τ_out을 설정해 r<R_ped 구역에 강한 전단을 유지함으로써 ℓ²/χ 구배가 급격히 상승하고, 식(6)의 압력 제한 p<½ρu_ϕ²에 도달할 때까지 페디스털이 성장한다.

페디스털이 임계 높이에 도달하면, 플럭스 튜브들의 평형점이 두 개에서 하나로 합쳐지는 사다리꼭점(saddle‑node) 분기(그림 4a)가 발생한다. 이때 플럭스 튜브는 급격히 외부로 방출되며, 속도는 식(8)에 의해 지수적·지수적 증가 후 포화한다. 붕괴 과정은 레일리‑테일러형 불안정의 머시룸 모양 플룸을 형성하고, 외부 냉 플라즈마가 내부를 대체하면서 압력 프로파일이 급격히 낮아진다. 이후 교환 난류가 코어를 재혼합하고 s/χ가 다시 균일해지면서 페디스털이 재구축되고, 같은 사이클이 반복된다.

에너지 측면에서는 전체 자기·열·회전 에너지와 최적 플럭스 튜브 교환을 통한 최소 에너지 상태와의 차이를 “사용 가능한 포텐셜 에너지(APE)”라 정의하고, β→0 극한에서 NP‑hard인 비선형 할당 문제(NAP)를 선형 할당 문제로 근사한다. 이를 통해 메타스테이블 평형의 저장 에너지와 붕괴 시 방출되는 에너지의 상한을 첫 원리적으로 계산할 수 있음을 보였다.

이러한 메커니즘은 토카막의 엣지 펠릿(ELM) 사이클과 구조적으로 동일하며, Z‑핀치라는 단순한 기하학에서도 동일한 물리적 원리가 작동함을 입증한다. 따라서 페디스털 메타스테이블성, 사다리꼭점 붕괴, 그리고 APE 기반 에너지 예측은 고온 플라즈마 제어와 ELM 억제 전략에 새로운 이론적 기반을 제공한다.