[論文レビュー] Variational Integration for Ideal Magnetohydrodynamics and Formation of Current Singularities

Yao Zhou|arXiv (Cornell University)|Aug 28, 2017

Advanced Numerical Methods in Computational Mathematics参考文献 116被引用数 3

ひとこと要約

本稿では、磁気抵抗性トポロジーを保存し、人工的な再接続を防ぐために、ラグランジュ座標と離散外微分計算を用いて、理想磁気流体力学（MHD）の変分積分法を開発した。3次元のライン・ティード・ハーム＝クルズルド＝テイラー問題に適用した結果、線形解は滑らかであるが、非線形解では解像度を上げるにつれて電流密度が発散する傾向を示した。これは、有限な系長さにおいて電流特異性が形成される可能性を示唆しているが、計算限界のため明確には確認されなかった。

ABSTRACT

Coronal heating has been a long-standing conundrum in solar physics. Parker's conjecture that spontaneous current singularities lead to nanoflares that heat the corona has been controversial. In ideal magnetohydrodynamics (MHD), can genuine current singularities emerge from a smooth 3D line-tied magnetic field? To numerically resolve this issue, the schemes employed must preserve magnetic topology exactly to avoid artificial reconnection in the presence of singular current densities. Structure-preserving numerical methods are favorable for mitigating numerical dissipation, and variational integration is a powerful machinery for deriving them. In this thesis, we develop variational integrators for ideal MHD in Lagrangian labeling by discretizing Newcomb's Lagrangian on a moving mesh using discretized exterior calculus. With the built-in frozen-in equation, the schemes are free of artificial reconnection, hence optimal for studying current singularity formation. Using this method, we first study a fundamental prototype problem in 2D, the Hahm-Kulsrud-Taylor (HKT) problem. It considers the effect of boundary perturbations on a 2D plasma magnetized by a sheared field, and its linear solution is singular. We find that with increasing resolution, the nonlinear solution converges to one with a current singularity. We then extend the HKT problem to 3D line-tied geometry, which models the solar corona by anchoring the field lines in the boundaries. The linear solution, which is singular in 2D, is found to be smooth. The nonlinear solution turns out smooth for short systems. Nonetheless, the scaling of peak current density vs. system length suggests that the nonlinear solution may become singular at a finite length. With the results in hand, we cannot confirm or rule out this possibility conclusively, since we cannot obtain solutions with system lengths near the extrapolated critical value.

研究の動機と目的

• 磁気トポロジーを保存することで人工的再接続を回避する構造保存型数値法を理想MHDに開発すること。
• 特に太陽コロナをモデル化するライン・ティード境界条件下で、真の電流特異性が3次元理想MHDにおいて形成可能かどうかを調査すること。
• 高解像度シミュレーションにより、ペイジャーの仮説（電流特異性が高エネルギーのナノフレアとコロナ加熱を引き起こす）を検証すること。
• 2次元モデルを超えて、3次元周期的およびライン・ティード幾何における電流特異性の兆候が残存するかどうかを特定すること。

提案手法

• 移動メッシュと離散外微分計算を用いてニューカンのラグランジアンを離散化し、変分積分法を構築する。
• ラグランジュ座標を用いて自然に「凍結された」条件を満たし、人工的再接続を排除する。
• 共鳴表面付近で解像度を細かくした四面体形状の非一様メッシュを実装し、特異的電流挙動を捉える。
• 非圧縮性をモデル化するため、剪断磁場と中程度の圧力を有するRMHD近似を適用する。
• 複数の波数モードを用いた境界駆動を導入し、共鳴を誘発し、電流集中を調査する。
• 解像度を段階的に向上させた際のピーク電流密度の収束性を検証し、特異性を示す非収束を検出する。

実験結果

リサーチクエスチョン

• RQ1 ペイジャーのナノフレア仮説が成立するためには、ライン・ティード境界条件下で3次元理想MHDにおいて真の電流特異性が形成可能か？
• RQ2 2次元のハーム＝クルズルド＝テイラー電流特異性は、3次元周期的およびライン・ティード幾何でも持続するか？
• RQ3 系長さが非線形3次元MHDにおける電流特異性形成にどのように影響するか？
• RQ4 正確なトポロジー保存を実現する変分積分法は、標準的手法が失敗する場面でも電流特異性を解明できるか？

主な発見

• 2次元では、解像度を上げるにつれて収束する解が得られ、電流特異性を示す。これは2次元予測の確認である。
• 3次元周期的幾何では、線形解は滑らかであるが、非線形解では解像度を上げるにつれてピーク電流密度が発散し、特異性の可能性を示唆する。
• 解像度 $N = 64, 96, 128, 160$ の間で、電流密度 $j_z$ の面内最大値は収束せず、平均値と標準偏差が増加するため、共鳴表面で非収束が示唆される。
• ピーク電流密度と系長さのスケーリングから、有限な臨界長さで特異性が発生する可能性があるが、計算限界のためその値に近い解は得られなかった。
• 本手法は磁気トポロジーを正確に保存し、人工的再接続を回避でき、特異的電流形成の信頼できる研究を可能にした。
• 共鳴表面 $x_0 = \pm k_z/k_y$ において、磁場線に沿った電流集中が観測され、3次元における磁場線に沿った整合性が確認された。