[論文レビュー] Buoyant magnetic flux ropes in a magnetized stellar envelope: Idealized numerical 2.5-D MHD simulations

Bertil F. Dorch|Sep 22, 2006

Solar and Space Plasma Dynamics参考文献 10被引用数 4

ひとこと要約

本研究では、2.5次元の理想化MHDシミュレーションを用いて、星の対流圏内におけるねじれたトロイダル磁気的フラックスロープが、多極磁場とどのように相互作用するかを調査している。その結果、周囲の多極磁場の極性および強度が、フラックスロープの浮上ダイナミクスと幾何学的形状に顕著な影響を及ぼすことが判明した。反平行にねじれたロープはより速く上昇し、平坦化するが、平行なねじれは、磁力線の再結合とシールド効果により、より安定で丸みを帯びた構造を形成する。

ABSTRACT

Context: The context of this paper is buoyant toroidal magnetic flux ropes, which is a part of flux tube dynamo theory and the framework of solar-like magnetic activity. Aims: The aim is to investigate how twisted magnetic flux ropes interact with a simple magnetized stellar model envelope--a magnetic "convection zone"--especially to examine how the twisted magnetic field component of a flux rope interacts with a poloidal magnetic field in the convection zone. Method: Both the flux ropes and the atmosphere are modelled as idealized 2.5-dimensional concepts using high resolution numerical magneto-hydrodynamic (MHD) simulations. Results: It is illustrated that twisted toroidal magnetic flux ropes can interact with a poloidal magnetic field in the atmosphere to cause a change in both the buoyant rise dynamics and the flux rope's geometrical shape. The details of these changes depend primarily on the polarity and strength of the atmospheric field relative to the field strength of the flux rope. It is suggested that the effects could be verified observationally.

研究の動機と目的

• 星の対流圏内における多極磁場が、ねじれたトロイダル磁気的フラックスロープの上昇および進化に与える影響を理解すること。
• 周囲の多極磁場の極性および強度が、フラックスロープの浮上ダイナミクスおよび幾何学的変形に影響を及ぼすかどうかを調査すること。
• 対流圏の磁場に対して反平行にねじれたフラックスロープが、磁気層をより効率的に貫通できるか、あるいはより容易に破壊されるかを特定すること。
• 磁気再結合が上昇中のフラックスロープのねじれおよびトポロジーをどのように変化させるかを評価すること。
• これらの相互作用が、太陽様の磁気活動周期およびフラックス出現パターンに観測的にどのように関連するかを探索すること。

提案手法

• 数値的シミュレーションにより、2048²の直交座標グリッド上で、6次精度の空間差分法および5次精度の補間を用いた高次精度有限差分法により、理想化された2.5次元圧縮性MHD方程式を解いている。
• ステaggeredメッシュ法により、∇·B = 0 を正確に保存し、磁束の整合性を保っている。
• 初期設定には、斜温対流圏モデルが含まれており、水平方向の磁場層と、浮上を促す初期条件を持つねじれたフラックスロープが埋め込まれている。
• フラックスロープは特定のねじれパラメータ χ で初期化され、磁化された層を通過する間にそのダイナミクスが追跡されている。
• 時間積分には3次精度の予測子・修正子スキームが用いられ、粘性およびジュール熱損失は Q_v および Q_J 項によってモデル化されている。
• 多項式状態方程式 P ∝ ρ^γ（γ = 5/3）を仮定し、フラックスロープ内の初期エントロピーは外部媒体と一致させることで浮上を誘発している。

実験結果

リサーチクエスチョン

• RQ1 周囲の多極磁場の極性が、ねじれたトロイダルフラックスロープの浮上速度および軌道にどのように影響を与えるか？
• RQ2 フラックスロープと上層の磁場層との間で発生する磁気再結合が、ロープのねじれおよびトポロジーをどのように変化させるか？
• RQ3 対流圏内の多極磁場の強度が、ねじれたフラックスロープの上昇を抑制するか、あるいは促進するか？
• RQ4 初期のねじれパラメータ χ が、磁気層を効率的に貫通するための臨界閾値にどのように影響を与えるか？
• RQ5 観測されたフラックスロープの幾何学的形状および出現ダイナミクスは、事前に存在する多極磁場との相互作用によって説明可能か？

主な発見

• 周囲の多極磁場に対して反平行にねじれたフラックスロープは、平行ねじれのものよりも速く上昇し、より高い高度に達することが、磁気的シールド効果の低減によるものである。
• 反平行ねじれのロープは、より平坦で傾斜がきつい軸方向磁場を示し、上部境界における磁場勾配が急になる。これは、活発領域で観測されるような迅速な出現を促進する可能性がある。
• 磁気層内での再結合により、反平行ねじれのロープのねじれが減少するが、これにより上昇が妨げられることはない。むしろ、層をより効率的に貫通する道筋を提供する。
• ねじれた磁場エネルギーと衝撃圧力のエネルギー収支に基づいて導かれた臨界ねじれパラメータ χ_c ≈ 5 が閾値として機能する。χ < χ_c のロープは強い浮上減衰を経験する。
• χ > χ_c のロープですら、対流圏内で停止する可能性があることから、臨界条件は χ のみに依存するのではなく、ねじれの進化および磁場配置にも依存することが示された。
• これらの相互作用により、電流シート内にプラズモイドが形成され、太陽および星の文脈で観測される再結合駆動ダイナミクスと整合的である。