[論文レビュー] Multisymplectic Geometry Method for Maxwell's Equations and Multisymplectic Scheme

Hongling Su, Mengzhao Qin|ArXiv.org|Feb 24, 2003

Numerical methods for differential equations参考文献 5被引用数 3

ひとこと要約

本稿は、可変なεおよびμを有する不均一で等方的かつ損失のない媒質におけるマクスウェル方程式の多重シンプレクティック幾何的枠組みを構築する。ラグランジアンのルジャンドル変換により多重シンプレクティックハミルトニアン形式を導出し、その多重シンプレクティック保存則を保つ九点プレッサマンスキームを提案する。このスキームは、電磁波伝搬の数値シミュレーションにおいて長期にわたり安定性と精度を示す。

ABSTRACT

In this paper we discussed the self-adjointness of the Maxwell's equations with variable coefficients $ε$ and $μ$. Three different Lagrangian are attained. By the Legendre transformation, a multisymplectic Bridge's (Hamilton) form is obtained. Based on the multisymplectic structure, the multisymplectic conservation law of the system is derived and a nine-point Preissman multisymplectic scheme which preserve the multisymplectic conservation law is given for the Maxwell's equations in an inhomogeneous, isotropic and lossless medium. At last a numerical example is illustrated.

研究の動機と目的

• 可変なεおよびμ係数を有するマクスウェル方程式の自己随伴な変分形式を確立すること。
• ラグランジアン密度からルジャンドル変換を用いて多重シンプレクティックハミルトニアン形式を導出すること。
• 多重シンプレクティック保存則を保存する構造保存型数値スキームを構築すること。
• 正確な解を用いた波動伝搬の数値シミュレーションを通じてスキームを検証すること。

提案手法

• 逆変分問題技法を用いて、可変係数媒質におけるマクスウェル方程式の3つの異なるラグランジアン密度を導出する。
• ルジャンドル変換を適用し、ラグランジュ系を多重シンプレクティックハミルトニアン形式に変換することで、系の背後にある幾何的構造を明らかにする。
• ハミルトニアン系の幾何的構造から多重シンプレクティック保存則を導出する。
• ハミルトニアン形式を直接離散化することで、多重シンプレクティック保存則を保つ九点プレッサマン多重シンプレクティックスキームを構築する。
• 中間段階で補助変数（U, V, P, Q）を用いて最終的な九点スキームを導出し、数値実装のため行列形式に表現する。
• 空間および時間方向にステaggered leapfrog型離散化を採用し、対称性および保存性を維持する。

実験結果

リサーチクエスチョン

• RQ1 可変係数媒質におけるマクスウェル方程式が、変分形式を可能にするために必要な自己随伴性の条件は何か？
• RQ2 マクスウェル方程式のラグランジアンからどのように多重シンプレクティックハミルトニアン構造を導出できるか？
• RQ3 空間的に変化するεおよびμを有するマクスウェル方程式の文脈において、多重シンプレクティック保存則の形は何か？
• RQ4 多重シンプレクティック構造を保存しつつ、電磁波伝搬を正確にシミュレートできる数値スキームを構築できるか？
• RQ5 九点多重シンプレクティックスキームの性能は、長時間シミュレーションにおいて正確な解と比較してどうなるか？

主な発見

• 本稿は、可変なεおよびμを有する1次形式のマクスウェル方程式が自己随伴であるための必要十分条件を確立し、変分形式の可能性を示した。
• 場の表現の選択に応じて、マクスウェル方程式に対して3つの異なるラグランジアンが導出された。
• ルジャンドル変換により、複素多重シンプレクティックハミルトニアン形式が得られ、系の本質的幾何的構造が明らかになった。
• 九点プレッサマンスキームが多重シンプレクティック保存則を正確に保存することが示され、長期的安定性が保証された。
• ε = μ = 1およびJ = 0の条件における数値結果は、正確な解と良好に一致しており、1000ステップの間、波動的挙動が維持され、誤差の増大も最小限に抑えられた。
• t = 1000における絶対誤差は、正確な解の形状と密接に一致する規則的で周期的なパターンを示し、高い精度と構造的忠実性を示した。