[論文レビュー] MHD Equilibrium Equation in Symmetric Systems

M Y Kucinski, Iberê L. Caldas|arXiv (Cornell University)|Mar 25, 2011

Advanced Thermodynamics and Statistical Mechanics参考文献 1被引用数 3

ひとこと要約

本稿では、物理量が2つの変数にのみ依存する対称的プラズマ系に対して、曲線座標系における一般化されたMHD平衡方程式を導出する。無視可能な座標を仮定し、同一の対称性を持つ磁気表面関数を用いることで、非直交的な一般化座標に対しても適用可能な統一された平衡方程式を構築し、Grad-Shafranov方程式やねじれ系の式といった既知の式を特別な場合として導出可能にする。

ABSTRACT

In MHD symmetric systems the equilibrium physical quantities are dependent on two variables only. In this cases it is possible to find a magnetic surface function that has the same symmetry. Under the assumption that the metric determinant is also independent of a third, ignorable coordinate, a general MHD equilibrium equation in curvilinear coordinates is deduced. This equation is specially useful when non-orthogonal generalized coordinates are used.

研究の動機と目的

• 2つの独立変数を持つ対称的プラズマ系に適用可能な一般化されたMHD平衡方程式の開発。
• 非直交的な一般化曲線座標系における平衡方程式の定式化により、座標系選択の柔軟性を向上。
• 系の対称性を保つ磁気表面関数Ψと電流関数Iを用いた平衡方程式の導出。
• トロイダル系、ねじれ系、自然座標系における既知の平衡方程式を導出することで、一般方程式の適用可能性を示す。
• 計量テンソルの成分を用いてBおよびJ場をΨとIの関数として表現することにより、対称的システムにおけるMHD平衡の統一的取り扱いを実現。

提案手法

• u₃を無視可能な座標とし、u₂およびu₃で周期的であると仮定する曲線座標(u₁, u₂, u₃)を用いる。
• 磁気場成分B²から導かれる、u₃で境界づけられた表面を通過する横断的磁束Ψ(u₁, u₂)を定義する。
• B₃のu₃にわたる平均値を用いて、磁気場と類似した形でΨの関数として電流関数I(u₁, u₂)を導入する。
• ∇×B = μ₀J および ∇P = J×B を組み合わせ、計量成分を用いてP′をΨおよびIの関数として表現することで、一般化されたMHD平衡方程式を導出する。
• 一般方程式を用いて、トロイダル系（付録A）、ねじれ系（付録B）、自然座標系（付録C）における特定の形を導出する。
• ベクトル解析の恒等式および計量テンソルの関係式（例：g₃₃, √g）を用い、BおよびJを∇Ψおよび∇Iの関数として表現することで、対称性の保持を保証する。

実験結果

リサーチクエスチョン

• RQ1 非直交的な一般化座標を用いた対称的システムにおける統一的MHD平衡方程式をどのように導出できるか？
• RQ2 曲線座標系における磁気表面関数Ψおよび電流関数Iは、対称性の保持に果たす役割は何か？
• RQ3 無視可能な座標（u₃）の仮定が、平衡方程式の導出をどのように簡略化するか？
• RQ4 一般平衡方程式は、Grad-Shafranov方程式やねじれ系の式といった既知の結果を再現できるか？
• RQ5 任意の曲線座標系において、計量テンソルを用いて磁束および電流密度場をΨおよびIの関数としてどのように表現できるか？

主な発見

• 曲線座標系における一般化されたMHD平衡方程式は、∂²Ψ/∂u₁² + ∂²Ψ/∂u₂² = -μ₀P′ - μ₀²II′ として導出され、無視可能な座標u₃を持つ対称的システムに有効である。
• 磁束密度は B = (e₃/g₃₃)×∇Ψ + B₃e₃/g₃₃ として表現され、∇Ψ·B = 0 を満たし、Ψ = 定数 が磁気表面を定義することを確認する。
• 電流密度は J = (e₃/g₃₃)×∇I + J₃e₃/g₃₃ として与えられ、磁束密度の表現と類似しており、対称性を保つ。
• トロイダル座標系では、方程式は∂/∂ξおよび∂/∂ωを含む形に簡略化され、圧力および電流の微分に依存する源項を有する。
• 直線的な磁気軸を持つねじれ系では、(1+α²r²)/r 項を含む半径方向およびポリオイダル方向の微分形に変換される。
• 自然座標系では、方程式は P′ = -J₃χ′/g₃₃ - μ₀II′/g₃₃ に簡略化され、表面積分により P′V′ = -χ′J′ - Φ′I′ が得られ、磁束と電流の関係を結ぶ。