[論文レビュー] An Adaptive Finite Element DtN Method for Maxwell's Equations in Biperiodic Structures

Xue Jiang, Peijun Li+4 more|arXiv (Cornell University)|Nov 29, 2018

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

ひとこと要約

本稿では、3次元二周期的グレーティング構造におけるマクスウェル方程式を解くために、ディリクレからノイマン（DtN）透明境界条件を用いた適応型有限要素法を提案する。有限要素離散化誤差とDtN作用素の切断誤差の両方を考慮した後験的誤差推定を導出し、メッシュと切断パラメータを同時に適応的に調整することで、準最適収束率を達成し、従来のPMLに基づく手法に比べて計算効率と精度で優れている。

ABSTRACT

Consider the diffraction of an electromagnetic plane wave by a biperiodic structure where the wave propagation is governed by the three-dimensional Maxwell equations. Based on transparent boundary condition, the grating problem is formulated into a boundary value problem in a bounded domain. Using a duality argument technique, we derive an a posteriori error estimate for the finite element method with the truncation of the nonlocal Dirichlet-to-Neumann (DtN) boundary operator. The a posteriori error consists of both the finite element approximation error and the truncation error of boundary operator which decays exponentially with respect to the truncation parameter. An adaptive finite element algorithm is developed with error controlled by the a posterior error estimate, which determines the truncation parameter through the truncation error and adjusts the mesh through the finite element approximation error. Numerical experiments are presented to demonstrate the competitive behavior of the proposed adaptive method.

研究の動機と目的

• 3次元二周期的グレーティング構造における時間調和マクスウェル方程式を解くための適応型有限要素法を開発すること。
• 非有界計算領域における非局所的DtN境界作用素の切断という課題に対処すること。
• 有限要素誤差とDtN切断誤差の両方を同時に制御する適応アルゴリズムを設計すること。
• 計算領域サイズを小さくし、効率を向上させたPML手法の代替手段を提供すること。
• 数値実験を通じて、準最適収束性と競争力のある性能を実証すること。

提案手法

• 透明境界条件をDtN作用素から導出し、有界領域における境界値問題として問題を定式化する。
• 非局所的DtN作用素をパラメータ$N$を用いて有限和に切断し、切断誤差は$N$とともに指数的に減少する。
• 双対性アプローチを用いて、有限要素誤差と切断誤差を統合した信頼性と効率性に優れた後験的誤差推定を導出する。
• 後験的誤差推定が適応的メッシュ・リファインメントを導く：有限要素誤差に基づきメッシュを細分化し、切断誤差に基づきパラメータ$N$を調整する。
• 適応アルゴリズムにより両誤差成分を動的にバランスさせ、最適収束率を達成する。
• ベクトル場のカール適合離散化にNédélecエッジ有限要素を用いて数値解を計算する。

実験結果

リサーチクエスチョン

• RQ1 3次元二周期的マクスウェル問題における、有限要素誤差とDtN切断誤差の組み合わせに対する後験的誤差推定が導けるか？
• RQ2 適応的メッシュ・リファインメントと適応的DtN切断を同時に制御することで、最適収束を達成できるか？
• RQ3 計算コストと精度の観点から、適応型DtN法は適応型PML法に比べてどのように性能を発揮するか？
• RQ4 準最適収束率$\|\mathbf{E} - \mathbf{E}_h^N\|_{H(\mathrm{curl},\Omega)} = O(N_k^{-1/3})$を達成できるか？
• RQ5 PMLが適用できない場合を除き、非有界領域における波動伝播問題に対して、DtN法はPMLの有効な代替手段となり得るか？

主な発見

• 提案手法は、数値実験により確認されたように、自由度$N_k$に関して準最適収束率$\|\mathbf{E} - \mathbf{E}_h^N\|_{H(\mathrm{curl},\Omega)} = O(N_k^{-1/3})$を達成する。
• 後験的誤差推定$\eta_h$は真の誤差を正確に捉え、メッシュと切断パラメータの適応的調整を効果的に導く。
• 平板例では、適応型DtN法は、自由度を少なくし、計算領域もより小さくすることで、適応型PML法と同等の精度を達成する。
• チェッカーボードグレーティング例では、複雑な散乱パターンを高精度に捉え、適応的リファインメントの過程でも誤差推定が信頼性を保つ。
• PML層を人工的に設ける必要がなくなるため、計算領域を散乱体に近づけることができ、計算コストを削減する。
• 適応型PML法と同等の性能を発揮し、特にPMLが適用できない場合に顕著な代替手段として有効である。