제한된 지각판 모델을 통한 코로나 가열 메커니즘

초록

본 논문은 초기 균일 자기장을 선결(line‑tied)된 두 면 사이에 두고, 한 방향(x)으로만 변하는 광구 발판 움직임을 가정한 제한된 지각판(tectonics) 모델을 제시한다. 저항성·점성에 의한 에너지 소산을 장기적으로 추적한 결과, 발판 움직임의 상관시간이 저항성 확산시간보다 훨씬 짧을 때 가열률이 플라즈마 저항성에 무관함을 보였다.

상세 분석

이 연구는 태양 코로나 가열 문제를 다루기 위해, 기존의 복잡한 3차원 MHD 시뮬레이션 대신, 강력한 기하학적 제한을 가한 2차원 RMHD 모델을 구축하였다. 초기 상태는 B = B₀ ẑ이며, z = 0, L 두 면에 선결(boundary‑tied) 조건을 부여한다. 발판(footpoint) 움직임을 x‑좌표에만 의존하도록 가정함으로써, 전류와 속도는 y‑방향으로만 존재하게 되고, RMHD 방정식(1‑2)은 선형화된 형태(4‑5)로 축소된다. 이때 전류 밀도 J와 속도 포텐셜 φ는 각각 Fourier 모드 Aₙ(z,t), φₙ(z,t) 로 전개되며, 각 모드는 z‑방향 2차 미분 방정식을 만족한다.

첫 번째 분석에서는 발판 속도가 시간에 대해 상수인 경우를 고려한다. 정적 해를 구하면 φₙ(z)∝sinh(η kₙ² z)·φₙᴸ/ sinh(η kₙ² L)와 Aₙ(z)∝cosh(η kₙ² z)·φₙᴸ·ν/η·1/ sinh(η kₙ² L) 가 얻어진다. 저항성이 매우 작을 때 η ν kₙ² L≪1 이면 Aₙ(z)≈φₙᴸ/(η kₙ² L) 로 z에 거의 의존하지 않으며, 이는 B⊥∼v_L τ_r/L 형태의 스케일링을 의미한다. 여기서 τ_r = w²/η 는 저항성 확산시간이며, 전류와 전자기 에너지의 축적이 저항성에 역비례함을 보여준다.

두 번째 분석에서는 실제 태양 광구의 무작위 발판 움직임을 모사하기 위해, 코히런스 시간 τ_coh 를 갖는 랜덤 워크 형태의 φ₀(t) 를 도입한다. φ_L(x,t)는 다중 사인 모드의 합으로 표현되며, 각 모드의 위상은 τ_coh 간격마다 무작위로 변한다. 수치 시뮬레이션은 스펙트럼 분해와 암시적 시간 적분법을 사용해 장시간(수천 τ_r) 동안 통계적 정상 상태에 도달하도록 설계되었다. 결과적으로, τ_coh ≪ τ_r 인 경우 전자기 에너지와 전류는 랜덤하게 변동하지만 평균값은 일정하게 유지된다. 가장 중요한 발견은, 이 평균 가열률 W_d ≈ B₀² v_L l_r/L 가 저항성 η 에 독립적이라는 점이다. 즉, 발판 움직임이 충분히 빠르게 무작위화될 때, 전류 시트가 형성되더라도 저항성에 의한 소산이 포화되며, 전체 가열은 플라즈마의 전도성에 좌우되지 않는다.

이러한 스케일링은 실제 코로나에서 η가 매우 작아 Lundquist 수가 10¹² ~ 10¹³ 수준일 때도 적용될 수 있음을 시사한다. 따라서, 미세한 나노플레어나 전류 시트가 연속적으로 발생하더라도, 전체 에너지 공급은 발판 운동의 통계적 특성에 의해 결정되고, 저항성은 2차적인 역할만을 수행한다는 결론에 도달한다.