[論文レビュー] A systematic construction of finite element commuting exact sequences

Bernardo Cockburn, Guosheng Fu|arXiv (Cornell University)|Apr 30, 2016

Advanced Numerical Methods in Computational Mathematics参考文献 13被引用数 3

ひとこと要約

本稿では、1次元、2次元、3次元の一般の多面体要素に対して、安定な混合有限要素法を保証する、有限要素可換完全系列を構築する体系的な手法を提示する。多項式形状関数と段階的次元別構築を活用することで、標準要素の既知の系列を再発見し、多角形、角柱、錐台に対して新たな族を特定。特に、先行の有理関数に基づく手法と比較して、錐台では顕著な次元削減を達成した。

ABSTRACT

We present a systematic construction of finite element exact sequences with a commuting diagram for the de Rham complex in one-, two- and three-space dimensions. We apply the construction in two-space dimensions to rediscover two families of exact sequences for triangles and three for squares, and to uncover one new family of exact sequence for squares and two new families of exact sequences for general polygonal elements. We apply the construction in three-space dimensions to rediscover two families of exact sequences for tetrahedra, three for cubes, and one for prisms; and to uncover four new families of exact sequences for pyramids, three for prisms, and one for cubes.

研究の動機と目的

• 1次元、2次元、3次元の任意の多面体要素に対して、有限要素空間の可換完全系列を構築する一般的で体系的なフレームワークの開発。
• 四面体、立方体、角柱、錐台などの標準要素に対する既存の構築を、単一の整合的な手法から導出することで統一・拡張すること。
• 一般の多角形および多面体要素に対して、これまでに未知であった族を同定し、完全系列を構築すること。特に正方形、角柱、錐台に対して未発見の族を特定。
• 特に錐台要素に対して、先行の有理関数に基づく構築と比較して、高次系列において顕著な次元削減を達成すること。
• 偏微分方程式の解法に依存せず、明示的かつ計算可能な基底関数を提供する、仮想有限要素法や有限要素系に対する実用的で構成的な代替手法を提供すること。

提案手法

• 1次元から出発し、2次元多角形要素へ、さらに3次元多面体要素へと段階的に拡張する再帰的構築により、de Rham複体全体にわたる連続性と可換性を保証する完全系列の構築。
• テンソル積および階層的構成を用いた多項式形状関数により、基底要素上の有限要素空間を定義。自由度と多項式次数の精密な制御を実施。
• 次元数え上げと核の解析（例：ker(∇×) ⊂ ∇H の証明）を用い、特に多変数多項式空間の分解を通じて、系列の完全性を検証。
• 標準基底要素（四面体、立方体、角柱、錐台）にこの手法を適用し、既知の系列を再現するとともに、一般要素に適用して新たな系列を発見。
• テンソル積および混合多項式空間（例：P_k|k, P_k|k+1, RT_k）を用い、de Rham複合体内の異なる微分形式にわたる整合性のある空間を構築。
• 明示的な多項式解析と係数一致を用いて、補間作用素が勾配、回転、発散作用素と可換であることを証明することで、可換図式性を検証。

実験結果

リサーチクエスチョン

• RQ1 任意の多面体要素に対して、多次元における有限要素可換完全系列を構築する統一的で体系的な手法を開発可能か？
• RQ2 四面体、立方体、角柱などの標準要素に対して既知の有限要素系列族（例：四面体、立方体、角柱）は、この体系的アプローチで再現可能か？
• RQ3 一般の多角形および多面体要素に対して、これまでに発見されていなかった新たな完全系列族は存在するか？
• RQ4 特に困難な要素（例：錐台）に対して、得られた有限要素空間の次元は、既存の構築と比較してどのように異なるか？
• RQ5 この手法は、偏微分方程式に基づく基底関数に依存せず、実用的で構成的な方法で最小互換有限要素系（mcFES）を生成可能か？

主な発見

• 2次元において、三角形に対して2つの系列族、正方形に対して3つの系列族を再発見。さらに、正方形に対して1つの新族、一般多角形要素に対して2つの新族を同定。
• 3次元において、四面体に対して2つの系列族、立方体に対して3つの系列族、角柱に対して1つの系列族を再発見。さらに、錐台に対して4つの新族、角柱に対して3つの新族、立方体に対して1つの新族を発見。
• 錐台要素に対して、本手法は[22,23]の有理関数に基づく系列と比較して顕著な次元削減を達成し、より効率的で実用的な有限要素空間を実現。
• 多項式空間の分解と係数解析を用いて、ベクトル値空間における回転作用素の核が、スカラー空間上の勾配作用素の像に丁度一致することを証明することで、系列の完全性を厳密に示した。
• 系列に含まれる空間の次元の公式を明示的に導出し、検証した。具体的には、2次元系列 S_{6,k}^{2d} に対して dim H₃ = 3(k+1)(k+2)(k+3)/2 という正確な数え上げが得られた。
• 偏微分方程式の解法や複雑な調和形式サブシステムに依存せず、明示的かつ構成的な方法で最小互換有限要素系（mcFES）を構築可能であることを示した。これにより、既存のFESおよびVEMフレームワークに対する構成的代替手法を提供した。