[논문 리뷰] Controlling the level of the ideal invariant fluxes for MHD turbulence using TURBO spectral solver

Bogdan Teaca, Cristian C. Lalescu외 2명|arXiv (Cornell University)|2011. 08. 12.

Solar and Space Plasma Dynamics참고 문헌 15인용 수 3

한 줄 요약

이 논문은 TURBO 스펙트럼 솔버에서 MHD 난류의 세 가지 이상의 보존량인 총 에너지, 크로스헬리시티, 자기헬리시티의 주입률을 독립적으로 제어할 수 있는 새로운 강제 메커니즘을 제시한다. 푸리에 공간에서의 나선형 분해를 통해 강제를 구현함으로써, 통계적 정상 상태 난류에서 플럭스 수준을 정밀하게 제어할 수 있으며, 이는 $5/3$ 스펙트럼 스케일링이 검증된 비균형성 및 나선형 MHD 상태의 체계적 연구를 가능하게 한다.

ABSTRACT

The ideal invariants present in the formalism of magnetohydrodynamics (MHD), i.e. global quantities that are conserved in the absence of sources and dissipative effects, play an important role in various theoretical and numerical studies of MHD turbulence. The fluxes of these ideal invariants represent separate channels that transfer the information across different scales in a turbulent system. Once a statistically stationary state of turbulence is reached, the amount of any ideal invariant quantity introduced in the system by a forcing mechanism equals the amount of the same quantity removed by the dissipative effects from the system. For highly developed turbulence, these two mechanisms act predominantly at different scales that are largely separated. Since the ideal invariant quantities cascade between scales, a constant flux is generated with great implication on the state of the system. Numerically, controlling the ideal invariant fluxes levels for a turbulent MHD system is important for the analysis of fundamental MHD turbulence properties. We propose a forcing mechanism that controls the three ideal invariants of MHD turbulence: the total energy, the cross-helicity and the magnetic helicity. This forcing is implemented in the freely available TURBO solver, that is also briefly presented.

연구 동기 및 목표

• MHD 난류의 세 가지 이상 보존량인 총 에너지, 크로스헬리시티, 자기헬리시티의 주입률을 독립적으로 제어할 수 있는 강제 메커니즘을 개발한다.
• 특히 비균형성 및 나선형 영역에서의 주입률 수준을 제어함으로써 MHD 난류의 수치적 연구를 가능하게 한다.
• 비압축성 MHD 난류의 고정밀 직접 수치 시뮬레이션(DNS)을 위한 자유로운 TURBO 스펙트럼 솔버에 강제를 구현한다.
• 통제된 주입 조건 하에서 의사에너지 스펙트럼에서 일관된 $5/3$ 스케일링을 보여줌으로써 방법을 검증한다.
• 크로스헬리시티 및 자기헬리시티의 통제된 주입을 통해 태양풍 난류 및 다이너모 물리학 연구를 지원한다.

제안 방법

• 강제는 엘사서어 변수를 사용해 푸리에 공간에서 수립되며, 속도 및 자기장이 $\mathbf{z}^{\pm} = \mathbf{u} \pm \mathbf{b}$로 분해되어 에너지 및 헬리시티 플럭스를 별도로 제어할 수 있도록 한다.
• 강제의 진폭 $\alpha^{u}_{\pm}(\mathbf{k})$, $\beta^{u}_{\pm}(\mathbf{k})$, $\alpha^{b}_{\pm}(\mathbf{k})$, $\beta^{b}_{\pm}(\mathbf{k})$는 목적한 주입률 $\varepsilon^{u}_{\pm}$, $\varepsilon^{b}_{\pm}$ 및 크로스헬리시티 제어 파라미터 $\sigma$에서 유도된다.
• 이 방법은 총 에너지, 크로스헬리시티, 자기헬리시티의 주입률을 $\varepsilon^{u}_{\pm}$, $\varepsilon^{b}_{\pm}$, $\sigma$를 통해 독립적으로 설정하며, 코시-슈와르츠 부등식을 위반하지 않도록 하기 위해 특수한 제약 조건을 적용한다.
• 강제는 통계적 정상성을 확보하고 이방성을 최소화하기 위해 많은 수의 $N_f$개의 푸리에 모드에 적용되며, 강제는 $\hat{\mathbf{f}}^{\pm}(\mathbf{k}) = \frac{\varepsilon_{\pm}}{N_f} \frac{\hat{\mathbf{u}}_{\pm}(\mathbf{k})}{E^{u}_{\pm}(\mathbf{k})}$로 표현된다.
• 이 방법은 비압축성 MHD 방정식을 주기적 경계 조건과 스펙트럼 공간 이산화를 사용해 푸는 TURBO 스펙트럼 코드에 구현된다.
• 강제는 통계적 정상 상태에서 스펙트럼과 플럭스를 계산하여 검증되었으며, 비균형 조건 하에서도 의사에너지 스펙트럼에서 $5/3$ 스케일링과 플럭스 분리가 확인되었다.

실험 결과

연구 질문

• RQ1 스펙트럼 시뮬레이션에서 MHD 난류의 세 가지 이상 보존량인 총 에너지, 크로스헬리시티, 자기헬리시티의 주입률을 어떻게 독립적으로 제어할 수 있는가?
• RQ2 비영인 크로스헬리시티 주입이 MHD 난류에서 의사에너지 플럭스와 스펙트럼의 분리에 어떤 영향을 미치는가?
• RQ3 제안된 강제 메커니즘이 비균형 난류 조건 하에서도 의사에너지 스펙트럼에서 $5/3$ 스케일링을 유지하는가?
• RQ4 이 방법은 비틀림이나 수치적 불안정성을 유발하지 않도록 스펙트럼 솔버에 구현될 수 있는가?
• RQ5 통제된 자기헬리시티 주입이 MHD 난류에서 다이너모 과정 연구에 얼마나 기여하는가?

주요 결과

• 제안된 강제 메커니즘은 MHD 난류 스펙트럼 시뮬레이션에서 총 에너지, 크로스헬리시티, 자기헬리시티의 주입률을 성공적으로 독립적으로 제어한다.
• 비균형 난류 조건에서 $\sigma = 0.4$일 경우, 두 의사에너지 플럭스($E^+$ 및 $E^-$)가 분리되어 양성 및 음성 헬리시티 모드에 대해 별개의 에너지 전달 역학이 존재함을 나타낸다.
• $E^+$ 및 $E^-$의 스펙트럼은 관성 영역에서 $5/3$ 스케일링을 보이며, 이는 콜모고로프 유사 스케일링과 일치한다. 비균형 강제 조건 하에서도 동일하게 유지된다.
• 코시-슈와르츠 부등식 위반을 방지하기 위해 제약 조건 $\sigma \varepsilon_{\pm} \leq \sqrt{4\varepsilon^{u}_{\pm}\varepsilon^{b}_{\pm}}$를 적용함으로써 수치적 안정성과 특이점 방지를 달성한다.
• 크로스헬리시티 주입은 두 의사에너지의 소산률에 영향을 미쳐 통계적 정상 상태에서 플럭스 수준에 유의미한 차이를 초래한다.
• TURBO 솔버에 구현된 이 방법은 보존량 플럭스를 통제한 고정밀 DNS를 가능하게 하여 태양풍 및 다이너모 현상 연구를 지원한다.