[論文レビュー] Curl-curl conforming elements on tetrahedra

Qian Zhang, Zhimin Zhang|arXiv (Cornell University)|Jul 20, 2020

Advanced Numerical Methods in Computational Mathematics参考文献 26被引用数 4

ひとこと要約

本稿では、3次元クワッドカール問題を解くための四面体メッシュ上でのH(curl²)-適合有限要素を構築する。これは電磁気学および磁気流体力学に現れる4階微分カール系である。提案された要素は最適な補間誤差推定を達成し、一般のリプシッツ領域において適合有限要素離散化を可能にし、数値的検証によりその精度と収束性が確認されている。

ABSTRACT

In [24], we proposed H(curl^2)-conforming elements on both a triangle and a rectangle. This family of elements provides a brand new method to solve the quad-curl problem in 2 dimensions. In this paper, we turn our focus to 3 dimensions and construct an H(curl^2)-conforming tetrahedral finite element. The newly proposed element has been proved to have the optimal interpolation error estimate. Having tetrahedral elements, we can solve the quad-curl problem in any Lipschitz domain by conforming finite element method. We also provide several numerical examples of using our element to solve the quad-curl problem. The results of the numerical experiments show the effectiveness and correctness of our element.

研究の動機と目的

• 3次元クワッドカール問題を解くための四面体メッシュ上でのH(curl²)-適合有限要素の開発。
• H¹(curl)要素よりも弱い正則性要件で、3次元におけるカール-カール適合要素を構築する課題の克服。
• 最低次要素で315の自由度を有するが、それでも最適な補間誤差推定を保証すること。
• H²-適合法に内在する射影誤差を回避し、正確な解への収束をH(curl²)空間で達成する適合有限要素法の提供。
• 任意のリプシッツ領域におけるクワッドカール問題の高精度な数値解法を、四面体メッシュを用いて実現すること。

提案手法

• ベクトル場の点値、そのカール、および辺と面上での法線微分を含む自由度（DOFs）を定義することで、四面体上でのH(curl²)-適合有限要素を構築する。
• 参照要素と変換行列Dを用い、DOFsを参照要素から物理的四面体に写像することで、一意性と適合性を保証する。
• 最適な補間誤差推定の証明を容易にするために、参照要素に関連する中間的要素を導入する。
• 従来のラグランジュ型基底関数の計算が非現実的であるため、[8]の手法を用いて高自由度要素（315 DOFs）の基底関数を計算する。
• ヤコビ行列および辺・面上の法線ベクトルを含む変換行列を用いて、物理的要素上のすべてのDOFsを参照要素上の関数値の線形結合として表現する。
• DOFsを用いて補間作用素を表現し、一貫性と安定性を検証することで、H(curl²)ノルムにおける補間誤差が最適に収束することを証明する。

実験結果

リサーチクエスチョン

• RQ1 3次元四面体メッシュ上にH(curl²)-適合有限要素を構築し、最適な収束率でクワッドカール問題を解くことは可能か？
• RQ2 自由度に辺および面上での法線微分が含まれる場合、これらが参照要素と物理的要素の間で自然に写像されない状況において、カール-カール適合要素の一致性と適合性をどのように保証できるか？
• RQ3 四面体上での提案されたH(curl²)-適合有限要素空間における最適な補間誤差推定は何か？
• RQ4 提案された手法は、H²-適合法に内在する射影誤差を回避し、H(curl²)空間において正確な解への収束を達成できるか？
• RQ5 自由度が315に達する高次H(curl²)-適合要素の基底関数を、どのようにして効率的に計算できるか？

主な発見

• 提案された四面体上でのH(curl²)-適合有限要素は、H(curl²)ノルムにおいて最適な補間誤差推定を達成しており、理論的信頼性が裏付けられている。
• 要素の構築は数学的に厳密であり、自由度とその変換に関する詳細な解析を通じて、一意性と適合性が確認されている。
• 最低次要素は315の自由度を有しており、標準的なネデレック要素よりも顕著に高い。これは、カールおよび2階微分を捉える必要があるためである。
• 数値実験により、手法の正しさと収束性が確認され、3次元クワッドカール問題における最適な収束率が得られている。
• H²-適合法や非適合法の制限を克服し、任意のリプシッツ領域においてクワッドカール問題の適合有限要素離散化を可能にしている。
• https://github.com/QianZhangMath/3D-curl-curl-conforming-FE にコードを公開した計算フレームワークにより、自由度の高さにもかかわらず基底関数の効率的計算が可能となっている。