[論文レビュー] Space-time least-squares finite elements for parabolic equations

Thomas Führer, Michael Karkulik|arXiv (Cornell University)|Nov 5, 2019

Advanced Numerical Methods in Computational Mathematics被引用数 2

ひとこと要約

本稿では、1階系の定式化におけるL²残差汎関数の最小化に基づき、熱方程式のための時空最小二乗有限要素法を提案する。この手法は一様安定性を保証し、対称正定値線形系を生成し、内蔵された後験的誤差推定器を備えた完全な時空適応メッシュを可能にし、単体メッシュ上での数値実験で最適収束率を達成する。

ABSTRACT

We present a space-time least squares finite element method for the heat equation. It is based on residual minimization in L2 norms in space-time of an equivalent first order system. This implies that (i) the resulting bilinear form is symmetric and coercive and hence any conforming discretization is uniformly stable, (ii) stiffness matrices are symmetric, positive definite, and sparse, (iii) we have a local a-posteriori error estimator for free. In particular, our approach features full space-time adaptivity. We also present a-priori error analysis on simplicial space-time meshes which are highly structured. Numerical results conclude this work.

研究の動機と目的

• 時間積分スキームの制限を回避する、放物型PDEのための安定的で適合的な時空有限要素法の開発。
• 最小二乗変分法を用いることで、任意の適合離散空間に対して一様安定性を保証。
• 最小二乗汎関数から導出された局所誤差指標を用いて、完全な時空適応を可能にする。
• 構造的単体時空メッシュ上での事前誤差推定を提供。
• 1次元および2次元問題における数値実験で、最適収束率と頑健な性能を実証。

提案手法

• 熱方程式を1階系系に定式化：∂tu − div σ = f, σ − ∇u = 0, 初期条件 u(0) = u₀。
• 最小二乗汎関数 j(u, σ) = ∫₀ᵀ ‖∂tu − div σ − f‖²_{L²(Ω)} dt + ∫₀ᵀ ‖σ − ∇u‖²_{L²(Ω)} dt + ‖u(0) − u₀‖²_{L²(Ω)} を定義。
• H¹(0,T;L²(Ω)) × L²(0,T;H¹(Ω)) の積空間上で、適切なトレースおよび正則性制約を満たすように j(u, σ) を最小化。
• 汎関数から対称的かつ強制的（coercive）な双線形形式を導出し、離散化において安定的かつスパースな剛性行列を保証。
• 汎関数を要素ごとの寄与に分解することで、局所後験的誤差推定器を構築。
• 誤差推定器に基づく適応的メッシュ refinement を実装し、局所的な時空メッシュの細分化を可能にする。

実験結果

リサーチクエスチョン

• RQ1 時空最小二乗有限要素法は、放物型問題において任意の適合離散空間に対して一様安定性を達成できるか？
• RQ2 時空領域全体でL²ノルムにおける残差最小化は、対称的かつ正定値でスパースな代数的線形系を生成するか？
• RQ3 最小二乗汎関数を用いて、適応的メッシュ refinement のための信頼性と効率性を兼ね備えた局所誤差推定器を導出できるか？
• RQ4 正則性が異なる問題に対して、単体時空メッシュ上でのこの手法が達成する収束率は何か？
• RQ5 均一メッシュ refinement と比較して、適応的メッシュ refinement は収束率および計算効率においてどのように異なるか？

主な発見

• 任意の適合離散化に対して一様安定性を達成し、任意の適合離散空間に対して頑健であることを保証。
• 得られる剛性行列は対称的かつ正定値であり、かつスパースであるため、プリコンディショニング付き共役勾配法（PCG）などの反復解法による効率的解法が可能。
• 最小二乗汎関数から導出された後験的誤差推定器は自然に局所的であり、効果的な時空適応を可能にする。
• 数値実験では、1次元で適応的メッシュ refinement を用いた場合の収束率が約0.45に達し、均一メッシュの0.25を上回る性能を示した。
• 特異性（例：再entrant角）を有する2次元問題では、均一メッシュの収束率約0.2から適応的メッシュの約0.24に向上し、正則性が低い場合の有効性が示された。
• 滑らかな解に対しては、L²ノルムにおける収束率が N⁻¹/² のオーダーに達し、最適な h-型メッシュ細分化の挙動と整合的である。