[論文レビュー] Fully computable a posteriori error estimator using anisotropic flux equilibration on anisotropic meshes

Natalia Kopteva|arXiv (Cornell University)|Apr 14, 2017

Advanced Numerical Methods in Computational Mathematics参考文献 6被引用数 8

ひとこと要約

本稿では、非等方的メッシュ上での特異摂動付き半線形反応拡散方程式に対して、非等方的フラックス等衡化を用いた完全に計算可能な後験的誤差推定器を提示する。短い辺を共有する非等方的三角形のペアに対して発散なし補正を導入し、非等方的台形積分を採用している。誤差上界の定数は、要素の直径、アスペクト比、摂動パラメータ $ε$ に依存せず、ロバスト性と完全な計算可能性を保証する。

ABSTRACT

Fully computable a posteriori error estimates in the energy norm are given for singularly perturbed semilinear reaction-diffusion equations posed in polygonal domains. Linear finite elements are considered on anisotropic triangulations. To deal with the latter, we employ anisotropic quadrature and explicit anisotropic flux reconstruction. Prior to the flux equilibration, divergence-free corrections are introduced for pairs of anisotropic triangles sharing a short edge. We also give an upper bound for the resulting estimator, in which the error constants are independent of the diameters and the aspect ratios of mesh elements, and of the small perturbation parameter.

研究の動機と目的

• 非等方的メッシュ上での特異摂動付き半線形反応拡散方程式に対する完全に計算可能な後験的誤差推定器の開発。
• メッシュの非等方性および小さな摂動パラメータ $\varepsilon$ に対してロバストであることを保証すること。
• 誤差定数に未知の定数を含まず、明示的に計算可能であるようにすること。
• 非等方的要素における短辺ジャンプの残差の境界を精緻化することで、効率性を向上させること。
• 誤差定数が要素の直径、アスペクト比、および $\varepsilon$ に依存しない上界を提供すること。

提案手法

• フラックス等衡化の前段階として、短い辺を共有する非等方的三角形のペアに対して発散なし補正を導入する。
• 非等方的要素上で非等方的台形積分を採用し、精度を向上させるとともにメッシュ制約を軽減する。
• 局所的なメッシュ幾何に基づく非等方的重みをフラックス等衡化に用いることで、安定性と精度を向上させる。
• 補助局所問題の解法を避けるために、明示的な局所手順によるフラックス再構築を構築する。
• 修正されたフラックスを用いたパッチベースの等衡化戦略を適用し、局所的一致性とグローバル誤差制御を確保する。
• エネルギーノルムと重み付き残差を用いて上界を導出し、定数が $\varepsilon$、要素サイズ、アスペクト比に依存しないことを保証する。

実験結果

リサーチクエスチョン

• RQ1 非等方的メッシュ上での特異摂動問題に対して、未知定数を含まず完全に計算可能な後験的誤差推定器を構築可能か？
• RQ2 短い辺と高いアスペクト比を持つ非等方的要素に対応するため、フラックス等衡化をどのように適合化できるか？
• RQ3 $\varepsilon$ およびメッシュの非等方性に対してロバスト性を維持するため、フラックス再構築および台形積分にどのような修正が必要か？
• RQ4 短辺ジャンプ残差の取り扱いを精緻化することで、誤差推定器の効率性を向上させられるか？
• RQ5 摂動パラメータ $\varepsilon \to 0$ の極限においても、提案手法の誤差推定器はエネルギーノルムでロバスト性を保つのか？

主な発見

• 提案された誤差推定器は完全に計算可能であり、未知定数を含まず、キャリブレーションなしに直接数値的に評価可能である。
• 推定器の上界には、要素の直径、アスペクト比、摂動パラメータ $\varepsilon$ に依存しない定数 $C$ が含まれる。
• 局所的メッシュ幾何に基づく重みを用いた非等方的フラックス等衡化と、短辺に対する発散なし補正により、ロバスト性が達成される。
• 特定の非等方的メッシュ条件下で、等衡化フラックス $\boldsymbol{\tau}_z$ の精緻化版が、短辺ジャンプの残差境界を改善する。
• 数値的証拠から、標準的なバブル関数に基づく効率性境界は短辺残差を過大評価していることが示され、新たな残差処理の正当性が裏付けられる。
• 特に非等方的メッシュ条件 $\mathcal{A}1_{\rm ani}$ を満たすメッシュでは、$\varepsilon \to 0$ の極限においても推定器は効率的かつロバストである。