[論文レビュー] Weak Galerkin Finite Element Methods for Optimal Control Problems Governed by Second Order Elliptic Partial Differential Equations

Chunmei Wang, Junping Wang+1 more|arXiv (Cornell University)|Oct 2, 2023

Advanced Numerical Methods in Computational MathematicsEngineering被引用数 3

ひとこと要約

本稿は、2階楕円型偏微分方程式に従う最適制御問題に対する弱ガレルキン（WG）有限要素法を開発した。不連続有限要素を用いることで、柔軟性と精度を向上させた。エネルギーノルムおよび$L^2$ノルムにおける最適な誤差推定を確立し、正多角形メッシュを用いた収束解析を通じて数値的に検証した。

ABSTRACT

This paper is concerned with the development of weak Galerkin (WG) finite element method for optimal control problems governed by second order elliptic partial differential equations (PDEs). It is advantageous to use discontinuous finite elements over the traditional $C^1$ finite elements here. Optimal order error estimates are established and confirmed by some numerical tests.

研究の動機と目的

• 2階楕円型PDEに従う最適制御問題に対する弱ガレルキン有限要素法の開発。
• 従来の$C^1$有限要素法の制限を克服するため、不連続近似を用いて柔軟性を向上。
• 数値解におけるエネルギーノルムおよび$L^2$ノルムの最適収束率の確立。
• 構造的および非構造的正多角形メッシュ上での広範な数値実験を通じた理論的結果の検証。
• 従来、関連する結果が得られていなかった最適制御設定に、弱ガレルキン手法を拡張。

提案手法

• 状態方程式とノイマン境界制御を含むコスト関数 functional を含む制約付き最小化問題として最適制御問題を定式化。
• 多角形または多面体要素上で定義された弱ラプラシアン作用素を用いた弱ガレルキン変分形式を導入。
• 状態変数およびその弱勾配に対して不連続有限要素空間を採用し、一般のメッシュタイプと高次精度を実現。
• 離散WGスキームを導出し、インフラ・サブ安定性解析を用いて離散解の存在および一意性を証明。
• 射影作用素を用いて正確解と離散解を比較することで誤差方程式を導出。
• WG射影の安定性および近似性質を用いて、エネルギーノルムおよび$L^2$ノルムにおける事前誤差推定を適用。

実験結果

リサーチクエスチョン

• RQ1 2階楕円型PDEに従う最適制御問題に対して、弱ガレルキン有限要素法を効果的に拡張できるか？
• RQ2 このような制御問題において、WG法のエネルギーノルムおよび$L^2$ノルムにおける収束挙動はいかなるものか？
• RQ3 正則化パラメータ$\alpha$は、数値解の条件数および精度にどのように影響するか？
• RQ4 一般の正多角形または多面体メッシュ上でも、WG法が最適の収束順序を達成できるか？
• RQ5 WG有限要素空間における多項式次数の違いは、精度および収束率にどのように影響するか？

主な発見

• $P_k$-$P_{k-2}$ WG有限要素ペアでは、エネルギーノルムで$O(h^{k})$、$L^2$ノルムで$O(h^{k+1})$の最適収束率が達成された。
• $P_4$-$P_2$ WG法では、エネルギーノルム誤差が$2.8$～$3.0$の収束順、$L^2$誤差が$2.7$～$4.9$の収束順となり、メッシュの細分化および$\alpha$の値に依存して変動した。
• $P_5$-$P_3$ WG法では、エネルギーノルム誤差が$5.9$～$6.0$の収束順、$L^2$誤差が$6.0$～$6.2$の収束順となり、最適収束が確認された。
• 数値結果により、$\alpha \to 0$ または $\alpha \to \infty$ のとき、解が解析的極限に近づき、メッシュの細分化に伴い誤差が一貫して減少することが確認された。
• 極端な値（$10^{-8}$ から $10^9$）を含む広い$\alpha$の範囲においても、メソッドは安定的かつ高精度を維持し、ロバスト性を示した。
• 数値解の可視化により、$\alpha$の増加に伴い、解が平坦または構造的プロファイルに滑らかに収束することが確認され、理論的予想と整合的であった。