[論文レビュー] Computation of Maxwell's equations on Manifold using DEC

Zheng Xie, Yujie Ma|ArXiv.org|Aug 31, 2009

Electromagnetic Simulation and Numerical Methods参考文献 17被引用数 4

ひとこと要約

本稿では、時空ピラミッド格子を用いて曲がった多様体上でのマクスウェル方程式を解くための新しい離散外微分計算（DEC）フレームワークを提示する。この手法により、2次元および3次元幾何構造における電磁波伝搬の安定的かつ高精度なシミュレーションが可能になる。本手法は重要な幾何的構造を保持し、非構造格子をサポートしており、時間的および空間的精度は1次であり、構造格子では2次精度を達成する。

ABSTRACT

In this paper, the method of discrete exterior calculus for numerically solving Maxwell's equations in space manifold and the time is discussed, which is a kind of lattice gauge theory. The analysis of its stable condition and error is also accomplished. This algorithm has been implemented on C++ plateform for simulating TE/M waves in vacuum.

研究の動機と目的

• 一般の2次元および3次元多様体上でマクスウェル方程式を解くための幾何的数値法を開発すること。
• Yeeスキームを時空ピラミッド格子に拡張し、非構造格子および可変時間ステップを可能にする。
• 切り捨て誤差および固有値条件の厳密な解析を通じて、スキームの安定性と収束性を保証すること。
• C++を用いた実装と検証を行い、複雑な表面および3次元多様体上での真空中のTE/TM波モードのシミュレーションを行う。

提案手法

• 時空を1次元時間セグメントと2次元／3次元空間の三角形分割または単体分割の積として離散化し、ピラミッド格子を形成する。
• 空間のリーマン計量（辺の長さ、面積、体積）と負の時間計量を組み合わせることで、離散ローレンツ計量を定義する。
• 外接円中心の双対セルを用いて、原始セルと双対セルの体積をスケーリングする離散Hodgeスター作用素を定義する。
• k次単体と(k+1)次単体の間の接続行列の転置として、離散外微分作用素$d$を実装する。
• 離散1形式上の変分原理を用いて離散マクスウェル方程式を導出し、源項および連続の式を得る。
• 局所的な幾何的データ（辺の長さ、面積など）を用いて、TEおよびTMモードのための明示的反復スキームにシステムを分割する。

実験結果

リサーチクエスチョン

• RQ1 非構造格子を用いた空間多様体上でのマクスウェル方程式の整合的離散化は、どのように実現可能であり、その幾何的・位相的構造をどのように保持できるか？
• RQ2 時空ピラミッド格子上での可変時間ステップを有するDECベースの時間領域スキームの安定性条件は何か？
• RQ3 提案されたDECスキームの電磁波伝搬における収束性および精度の順序は何か？
• RQ4 干渉や曲がった表面および3次元多様体上での伝搬といった複雑な波動現象をシミュレートできるか？
• RQ5 非長方形格子上での本スキームの性能は如何であり、構造格子と非構造格子における精度はどのように異なるか？

主な発見

• DECスキームは、自動的に離散ビアンキ恒等式$dF = 0$を満たし、電磁場曲率の局所的保存を保証する。
• 一般の非構造格子では時間的および空間的精度が1次であり、時間ステップと空間ステップが等価な長方形格子では2次精度を達成する。
• 時間的および空間的固有値解析から導かれたCFL型条件の下で、安定性が保証される：$\Delta t \leq \min_{P_{123}} \left( \frac{1}{c} \sqrt{ \frac{2P_{123}}{ \left( \frac{l_{23}}{l_{A0}} + \frac{l_{34}}{l_{B0}} + \frac{l_{45}}{l_{C0}} \right) } } \right)$。
• 連続マクスウェル方程式と整合的であり、切り捨て誤差が$O((\Delta t)^2 + \Delta t |*e|)$であるため、ラクスの等価定理により収束性が保証される。
• C++を用いた数値実装により、スタンフォード・バニーおよび球面上でのガウスパルスの伝搬と干渉を効果的にシミュレートできた。
• 本手法は、鋭角三角形、長方形、正多角形を含む任意の三角形分割をサポートしており、複雑な幾何構造上での柔軟なシミュレーションを可能にする。