전류 흐르는 제트의 자기유체역학 불안정성

초록

본 논문은 축방향과 방사형(azimuthal) 자기장이 동시에 존재하는 천체 제트의 선형 MHD 안정성을 분석한다. 서브열적(β≫1) 자기장 환경에서 비축대칭(m≠0) 섭동에 대해 언제나 불안정함을 보이며, 축방향 필드가 약할 때는 모든 m에 대해 성장률이 포화하고, 축방향 필드가 강할 때는 높은 m값에서만 불안정이 나타난다.

상세 분석

본 연구는 무한히 긴 원통형 제트 흐름을 가정하고, 기본 상태에서 축방향 자기장 B_z는 일정하고 방사형 자기장 B_φ은 반경 s에 따라 변하는 형태(B_φ=B_φ0 ψ(s))를 채택하였다. 플라즈마 β≫1인 고β 환경을 전제로 비압축성 근사를 적용했으며, 이는 알벤 파동 주파수보다 낮은 성장률을 갖는 느린 모드에 적합하다. 선형화된 MHD 방정식(1‑4)을 전개하고, 섭동을 exp(σt‑ik_z z‑im φ) 형태로 가정함으로써 방사형 속도 섭동 v_{1s}에 대한 2계 미분 방정식(10)을 도출하였다. 이 방정식은 σ(성장률)와 파라미터(α=∂lnB_φ/∂lns, λ=1+m²/(s²k_z²) 등) 사이의 비선형 고유값 문제이며, 경계조건으로는 축축면에서 유한성, 외부에서는 v_{1s}=0(강체 전도벽) 를 적용하였다.

수치해석에서는 ψ(s)=x^n e^{-p(x‑1)} 형태를 이용해 다양한 B_φ 프로파일을 모사하였다. n=1은 축 중심에 전류가 존재하는 경우, n=2는 전류가 외부 경계에 집중되는 경우를 의미한다. 파라미터 ε=B_z/B_φ0, q=k_z s_1, Γ=σ/ω_{B0} 등을 도입해 무차원화된 방정식(13‑14)을 풀었다. 결과는 다음과 같다.

1. 축방향 필드가 약한 경우(ε≪1): 모든 azimuthal 모드 m에 대해 불안정이 발생한다. 특히 m이 작을수록 성장률이 급격히 증가하지만, m≈5~10 정도에서 포화값에 도달한다. 이는 전류 구동 불안정(current‑driven instability)이 주도함을 의미한다.

2. 축방향 필드가 강한 경우(ε≫1): 불안정은 주로 높은 m(≥10)에서만 나타난다. 이때 성장률은 m이 증가함에 따라 다시 감소하는 경향을 보이며, 최적 m는 B_φ 프로파일에 따라 달라진다. 즉, 강한 축방향 필드가 낮 m 모드(특히 m=1 kink)를 억제하지만, 고‑m 압력 구동 불안정(pressure‑driven) 혹은 혼합 모드가 여전히 존재한다.

3. B_φ 프로파일 의존성: n과 p 값에 따라 B_φ의 최대 위치가 경계에 있거나 내부에 있을 때 성장률 차이가 나타난다. B_φ가 외부 경계에 집중될수록(큰 n, 작은 p) 불안정이 더 강해지며, 이는 전류가 외부에 몰려 있을 때 라모어 모드가 쉽게 성장함을 시사한다.

4. 경계조건 민감도: 외부 경계조건을 v_{1s}=0 외에 v’_{1s}=0 등으로 바꾸어도 정성적 결과는 크게 변하지 않는다. 이는 제트와 주변 매질 사이의 얇은 전단층이 실제 경계 역할을 하며, 전단층 두께가 성장률에 미치는 영향이 미미함을 확인한다.

전체적으로 이 연구는 전류 구동 불안정과 압력 구동 불안정이 동시에 작용하는 복합적인 불안정 메커니즘을 제시한다. 특히 β≫1인 서브열적 제트에서도 비축대칭 섭동이 언제든지 성장할 수 있음을 보여주며, 관측된 제트의 구조적 복잡성(와인딩, 꼬임, 국부적인 밝기 변동 등)이 이러한 MHD 불안정에 기인할 가능성을 뒷받침한다.