플라스미드 이동 메커니즘 3D MHD 모델링

초록

본 연구는 JOREK 3D 비선형 MHD 코드를 이용해 JET H‑mode 플라즈마에 저자장면(LFS) 중간면에 순간적인 중성 가스원을 주입함으로써 대량 물질 주입(MMI) 후 플라스미드의 이동 메커니즘을 조사한다. 시뮬레이션은 전단 알레븐 파동(SAW) 패킷 전파에 의한 “SAW 브레이킹”과 외부 저항 전류에 의한 “페구리에 브레이킹”이 전하 분리를 억제하고 E×B 드리프트를 제한한다는 기존 이론을 재현한다. 작은 소스 진폭에서는 마이크로초 수준의 SAW 브레이킹이 속도를 제한하고, 큰 토로이달 폭을 가진 경우에는 더 긴 시간 스케일의 페구리에 브레이킹이 조기에 작용한다. 또한 E×B 흐름 영역의 크기가 플라스미드 속도에 결정적인 영향을 미치며, SAW 브레이킹이 지배할 때는 유효 압력을 고려한 이론값과 좋은 일치를 보인다.

상세 분석

이 논문은 플라스미드 드리프트를 이해하기 위해 두 가지 주요 브레이킹 메커니즘을 정량적으로 검증한다. 첫 번째는 플라스미드 양단에서 발생하는 전단 알레븐 파동(SAW) 패킷이 알레븐 속도 CA 로 전파되면서 전하를 플라스미드 외부로 전달하는 과정이다. SAW 브레이킹은 전류 jA = –∇⊥²Φ/(μ0 CA) 형태로 표현되며, 이는 전하 분리를 빠르게 상쇄시켜 초기 가속을 억제한다. 이 메커니즘에 의해 도출된 가속 방정식(식 1)에서 두 번째 항이 SAW 브레이킹을 나타내며, 압력 차이(p0–p∞)에 비례하는 첫 번째 항과 경쟁한다. 작은 중성 가스 주입(소스 진폭 ≤ 1.4×10²⁷ s⁻¹)에서는 플라스미드 압력이 충분히 낮아 SAW 브레이킹이 지배적이며, 시뮬레이션 결과 V_D,lim 이 이론식 V_D,lim ≈ μ0 CA p0 Z0 R / Bφ² 와 거의 일치한다. 여기서 Z0 은 플라스미드 길이, R 은 대형 반경, Bφ 는 토로이달 자기장이다.

두 번째 메커니즘인 페구리에 브레이킹은 플라스미드 양단을 연결하는 외부 저항 전류 j∥ 가 형성되는 현상이다. SAW 패킷이 플라스미드 양쪽에서 반사·재결합하면서 전류 루프가 완성되고, 이는 전하 분리를 장기적으로 소멸시킨다. 페구리에 브레이킹은 연결 길이 L_con 에 역비례하는 시간 스케일 τ_P≈L_con/(2 CA) 로 나타나며, 토로이달 소스 폭 Δφ 가 커질수록 L_con 가 짧아져 τ_P 가 감소한다. 논문에서는 Δφ 를 0.5 rad에서 2 rad 로 확대했을 때 V_D 가 급격히 감소하는 현상을 관찰했으며, 이는 페구리에 브레이킹이 조기에 작동함을 명확히 보여준다.

또한, E×B 흐름 영역의 크기가 플라스미드 속도에 미치는 영향도 강조한다. 플라스미드 내부에서 발생하는 전위 U 는 E×B 흐름을 유도하고, 이 흐름이 차지하는 부피가 클수록 유효 압력 p_eff 가 증가해 SAW 브레이킹에 의한 포화 속도가 상승한다. 시뮬레이션에서는 플라스미드 중심에서 측정한 n_e 와 p_e 가 소스 진폭에 비선형적으로 변함에도, V_D,lim 은 p_eff 를 이용한 이론식과 일치함을 확인했다.

마지막으로, 시뮬레이션에서 실제 SAW 존재를 확인하기 위해 전자기 파동 구조와 전류 분포를 분석하였다. 전류 패턴이 플라스미드 양단에서 대칭적으로 전파되고, 전위 Φ 가 알레븐 속도로 전파되는 모습이 관측되었으며, 이는 기존 이론에서 가정한 SAW 패킷 전파와 일치한다. 전반적으로 JOREK 코드가 3D 비선형 MHD 현상을 충분히 포착하고, 실험적 플라스미드 드리프트와 이론적 예측을 연결하는 중요한 교량 역할을 함을 입증한다.