[論文レビュー] Theoretical Considerations on the Computation of Generalized Time-Periodic Waves

Dirk Pauly, Tuomo Rossi|arXiv (Cornell University)|May 20, 2011

Numerical methods in inverse problems参考文献 42被引用数 7

ひとこと要約

本稿では、一般化された波動方程式の時間周期的解を計算するための新しい最小二乗制御法を提案する。この手法は、微分形式で定式化された一階系を用いる。$L^2$-型エネルギーノルムを、一階のマクスウェル型系を介して活用することで、前処理を要しない共役勾配法の最小化が可能となり、$O(n)$の計算複雑性と、2次誤差低減を伴う数値的収束性を達成した。

ABSTRACT

We present both, theory and an algorithm for solving time-harmonic wave problems in a general setting. The time-harmonic solutions will be achieved by computing time-periodic solutions of the original wave equations. Thus, an exact controllability technique is proposed to solve the time-dependent wave equations. We discuss a first order Maxwell type system, which will be formulated in the framework of alternating differential forms. This enables us to investigate different kinds of classical wave problems in one fell swoop, such as acoustic, electro-magnetic or elastic wave problems. After a sufficient theory is established, we formulate our exact controllability problem and suggest a least-squares optimization procedure for its solution, which itself is solved in a natural way by a conjugate gradient algorithm operating in the canonical Hilbert space. Therefore, it might be one of the biggest advances of this approach that the proposed conjugate gradient algorithm does not need any preconditioning.

研究の動機と目的

• 不均質媒質における一般化された波動方程式の時間周期的解を効率的かつ安定に計算するための手法を開発すること。
• 音響波、電磁波、弾性波を統一的に扱える、微分形式を用いた波動方程式の制御ベースのアプローチの一般化。
• 問題を$ L^2$-型ヒルベルト空間に再定式化することで、共役勾配法の最小化における前処理の必要性を排除すること。
• 良好に条件付けられた最小二乗法の定式化により、現代のマルチコアおよび分散アーキテクチャにおける効率的な並列実装を可能にすること。
• 正確な解が得られる3次元散乱問題を用いた数値的検証を通し、収束性と精度に焦点を当てる。

提案手法

• 音響波、電磁波、弾性波問題を統一的に取り扱えるように、交互に配置された微分形式を用いた一階系として波動方程式を定式化する。
• 正確な可制御性理論を適用し、1周期後の初期状態と最終状態の差を測る二次的汎関数の最小化問題に時間周期的解を再定式化する。
• $ L^2$-型ヒルベルト空間における最小二乗法の定式化により、良好に条件付けられた最小化問題を保証し、前処理の必要性を回避する。
• 自然に前処理を要しない$ L^2$空間で直接作用する共役勾配法アルゴリズムを実装し、並列計算に適した構造を有する。
• 適切に中心をもつ単体メッシュと、外接円中心に基づく双対メッシュを用いた離散外微分計算（DEC）により、系を離散化し、対角行列を有する離散ホッジ作用素を保証する。
• ドメインの切断のため、吸収境界条件（ABC）、一般化されたディリクレ・ツー・ノイマン（DtN）作用素、および完全一致層（PML）を微分形式の枠組みに組み込む。

実験結果

リサーチクエスチョン

• RQ1 微分形式で定式化された波動方程式の一階系は、前処理を要しない良好に条件付けられた最小二乗問題を導くことができるか？
• RQ2 提案された共役勾配法は、メッシュサイズに依存せず、実際には$ O(n)$の計算複雑性を達成するか？
• RQ3 正確な解が得られる3次元時間周期的波動散乱問題において、本手法の収束性と精度はいかがなっているか？
• RQ4 非構造的単体メッシュを用いて、可変係数媒質や複雑な幾何形状を扱えるように、本フレームワークを拡張可能か？
• RQ5 異なるメッシュの細分化において、共役勾配法の残差と誤差の収束特性はどのように変化するか？

主な発見

• 共役勾配法は140反復後に残差を$10^{-5}$の要因で低減し、メッシュ要素のサイズに依存しない収束性を示した。
• シミュレートされた電場における相対誤差は、メッシュの細分化に伴い減少し、平均辺長に対する2次収束を示唆した。
• $ L^2$-型エネルギーノルムを用いた定式化により、最小化問題が良好に条件付けられ、共役勾配ソルバーにおける前処理の必要性が排除された。
• 正確な解が得られる3次元電磁波散乱問題に対する数値結果は、6回のメッシュ細分化において一貫した収束性と精度を示した。
• 外接円中心に基づく双対メッシュを用いた離散外微分計算（DEC）の実装は、非構造的3次元メッシュと非同期時間刻みを効果的に処理できた。
• 本フレームワークは、音響波、電磁波、弾性波といったさまざまな波動問題を、統一的な微分形式の定式化により自然に扱える。