[논문 리뷰] High Lundquist Number Resistive MHD Simulations of Magnetic Reconnection: Searching for Secondary Island Formation

C. S. Ng, Srivathsan Ragunathan|arXiv (Cornell University)|2011. 06. 02.

Ionosphere and magnetosphere dynamics인용 수 3

한 줄 요약

이 연구는 룬드슈트 수가 $2 \times 10^5$에 이르는 고해상도 저항성 MHD 시뮬레이션을 통해 자기 재결합 시 2차 도시(플라스모이드) 형성을 조사한다. 외부 노이즈가 없는 조건에서는 재결합이 수렴-파커 척도($S^{-1/2}$)를 따르며 플라스모이드는 형성되지 않지만, 소규모 노이즈를 추가하면 플라스모이드가 나타나고 재결합 속도는 $S$에 의존하지 않게 된다. 이는 노이즈 유도 흐트러진 turbulent 재결합 메커니즘이 저항성 MHD에서 빠른 재결합을 가능하게 할 수 있음을 시사한다.

ABSTRACT

Recently, secondary island formation due to the tearing instability of the Sweet-Parker current sheet was identified as a possible mechanism that can lead to fast reconnection (less sensitive dependence on Lundquist number $S$) both in numerical simulations using Particle-in-Cell (PIC) method [Daughton et al. 2009], as well as using resistive magnetohydrodynamics (MHD) [Lapenta 2008; Bhattacharjee et al. 2009]. This instability is thought to appear when $S$ is greater than a certain threshold. These recent results prompt us to perform more resistive MHD simulations of a basic reconnection configuration based on the island coalescence instability, using much higher resolutions and larger $S$. Our simulations are based on a fairly standard pseudo spectral code, which has been tested for accuracy, convergence, and compared well with codes using other methods [Ng et al. 2008]. In our simulations, formation of plasmoids were not found, except when insufficient resolution was used, or when a small amount of noise was added externally. The reconnection rate is found to follow the Sweet-Parker scaling when no noise is added, but increases to a level independent of $S$ with noise, when plasmoids form. Latest results with $S$ up to $2 imes 10^5$ will be presented.

연구 동기 및 목표

• 고루드슈트 수가 큰 저항성 MHD 시뮬레이션에서 2차 도시(플라스모이드) 형성이 천연적으로 발생하는지 조사하기 위해.
• 이전 연구에서 관측된 빠른 재결합 속도가 수치적 오류인지 물리적 메커니즘인지 판단하기 위해.
• 수치적 해상도와 외부 노이즈가 플라스모이드 형성과 재결합 속도 향상에 미치는 역할을 평가하기 위해.
• 노이즈 없는 시뮬레이션과 노이즈가 포함된 시뮬레이션 간의 재결합 역학을 비교하여 빠른 재결합이 나타나는 조건을 분리하기 위해.
• 통제된 조건에서 고해상도 저항성 MHD에서 수렴-파커 척도의 타당성을 평가하기 위해.

제안 방법

• 시뮬레이션은 2차원 비압축성 저항성 MHD 방정식을 정규화된 형태로 푸는 가상 스펙트럼 코드를 사용하며, $\partial_t \Omega + [\phi, \Omega] = [A, J] + \nu \nabla_\perp^2 \Omega$ 및 $\partial_t A + [\phi, A] = \eta \nabla_\perp^2 A$ 를 포함한다.
• 루드슈트 수 $S$는 점성도 $\eta$를 통해 제어되며, $\eta \sim S^{-1}$을 사용하여 $S$ 값이 $2 \times 10^5$에 이르게 한다.
• 수렴성 확보와 부적절한 플라스모이드 형성 방지를 위해 최대 $8192^2$ 격점의 수치적 해상도를 확보한다.
• 스펙트럼 영역에서 제어된 진폭(에너지 비율 ~$10^{-6}$)과 상관 시간(~0.25)을 가진 외부 무작위 노이즈를 도입하여 난류 편향을 모의한다.
• 재결합 속도는 $dA_0/dt$를 통해 측정되며, $A_0$는 최대 플럭스 함수 크기이고 재결합 단계 동안 평균화된다.
• 다양한 $S$와 해상도에서 노이즈 유무에 따라 시뮬레이션을 비교하여 노이즈가 플라스모이드 형성과 재결합 속도에 미치는 영향을 분리한다.

실험 결과

연구 질문

• RQ1 고해상도 저항성 MHD 시뮬레이션에서 고루드슈트 수($S \sim 10^5$) 조건에서 수렴-파커 전류판이 갈라지는 진동 모드에 불안정해지는가?
• RQ2 외부 교란 없이 고해상도 저항성 MHD 시뮬레이션에서 2차 도시(플라스모이드) 형성이 자연스럽게 발생할 수 있는가?
• RQ3 수치적 노이즈 또는 외부 교란이 플라스모이드 형성과 빠른 재결합을 유도하는 데 어떤 역할을 하는가?
• RQ4 노이즈가 추가된 상태와 추가되지 않은 상태에서 재결합 속도는 $S$에 따라 어떻게 스케일링되는가?
• RQ5 시뮬레이션 결과가 해상도와 소규모 노이즈의 포함 여부에 얼마나 의존하는가?

주요 결과

• 외부 노이즈가 없는 고해상도 시뮬레이션에서는 재결합 속도가 수렴-파커 척도 $S^{-1/2}$를 따르며, $S = 2 \times 10^5$ 조건에서도 플라스모이드 형성이 관찰되지 않는다.
• 플라스모이드 형성과 배출은 소규모 외부 노이즈가 추가된 경우에만 관찰되며, 이는 불안정성에 대한 임계값 유사한 촉발 메커니즘을 시사한다.
• 노이즈가 추가된 경우 재결합 속도는 $S$에 거의 의존하지 않으며, 큰 $S$ 조건에서 수렴-파커 예측치보다 훨씬 높은 수준에 도달한다.
• 노이즈 추가 시 진동하는 그러나 향상된 재결합 속도를 보이며, 평균 $dA_0/dt$는 $S$ 값이 $1.25 \times 10^4$에서 $2.5 \times 10^4$에 이르는 범위에서 일정하게 유지된다.
• 결과는 플라스모이드 형성이 고 $S$ 조건에서의 수렴-파커 전류판의 내재적 특성이 아니라, 외부 교란 또는 수치적 노이즈의 결과임을 시사한다.
• 이 연구는 해상도와 소규모 교란의 존재가 시뮬레이션에서 빠른 재결합 메커니즘이 나타나는지 여부를 결정하는 데 있어 극히 중요한 역할을 한다고 강조한다.