A $C^0$ weak Galerkin method with preconditioning for constrained optimal control problems with general tracking

この論文は、偏微分方程式（PDE）に支配された制約付き最適制御問題（OCP）を解くための新しい数値手法を提案しています。具体的には、一般的な追跡コスト関数（general tracking）と点ごとの状態制約（pointwise state constraints）を備えた問題に対して、**$C^0$弱ガラーキン法（$C^0$-WG）と加法的シュワルツ前処理（Additive Schwarz Preconditioner）**を組み合わせた手法を構築し、その理論的解析と数値的有効性を示しています。

以下に、論文の主要な内容を技術的に詳細に要約します。

1. 研究の背景と課題

• 問題設定: 物理法則（PDE）に従う空間・時間プロセスを制御し、所望の状態分布に近づける最適制御問題。特に、状態変数に対する点ごとの上下限制約（$\psi_- \le y \le \psi_+$）と、点、曲線、部分領域など多様な特徴に対する追跡コストを含む問題が対象です。
• 数学的難しさ:
  • 制御変数を PDE 制約を通じて消去すると、状態変数に関する4 階の変分不等式として定式化されます。
  • 一般的な追跡項や点ごとの制約により、最適解やラグランジュ乗数の正則性（滑らかさ）が低下し、測度値の乗数が現れる可能性があります。
  • 数値的には、4 階の問題を扱うために $C^1$ 連続な有限要素が必要になるか、あるいは内部ペナルティ法（C0-IP）などの弱定式化が必要となりますが、これらはパラメータの調整が困難だったり、要素ごとの行列組立が複雑であったりする課題があります。
  • 得られる線形システムは条件数が非常に大きく（ill-conditioned）、反復解法の収束が悪化する傾向があります。

2. 提案手法：$C^0$弱ガラーキン法（$C^0$-WG）

著者らは、4 階の問題に対して、大域的に連続な 2 次ラグランジュ要素を用いた新しい $C^0$-WG 法を提案しました。

• 弱ラプラシアン（Weak Laplacian）の導入:
  • 従来の C0-IP 法では、接線方向の微分の連続性をペナルティ項で弱く課すために、トレース項（境界積分）と大きなペナルティパラメータ $\rho$ が必要でした。
  • 提案手法では、要素内で定義された「弱ラプラシアン」$\Delta_w$ を用います。これは、要素内の関数値と隣接要素との境界での勾配の平均値（またはジャンプ）を用いて、要素ごとに局所的に計算される定数関数として定義されます。
• パラメータフリーの安定化:
  • 提案された双一次形式は、弱ラプラシアンの $L^2$ 内積と、勾配のジャンプに基づく安定化項（ペナルティ項に相当）で構成されます。
  • 最大の特徴: 従来の C0-IP 法と異なり、ペナルティパラメータの調整が不要であり、かつ要素ごとの行列組立が効率的に行えます（特定のトレース項を計算する必要がないため）。
• 定式化:
  • 離散化された最適制御問題は、$C^0$-WG 双一次形式 $a_w(\cdot, \cdot)$ を用いた変分不等式として記述されます。

3. 前処理手法：加法的シュワルツ法

4 階の問題から生じる巨大で条件数の悪い線形システムを効率的に解くため、**1 レベルの加法的シュワルツ前処理（Additive Schwarz Preconditioner, ASP）**を開発しました。

• ドメイン分解: 計算領域を重なり合う部分領域に分割し、各部分領域で局所問題を独立に解きます。
• 条件数評価: 理論解析により、前処理を施したシステムの条件数 $\kappa(B_h A_h)$ が、重なり幅 $\delta$ に対して $O(\delta^{-4})$ 以下に抑えられることを示しました（形状が規則的であればさらに改善可能）。
• 効果: 前処理なしでは $O(h^{-4})$ で発散する条件数が、前処理により $O(h^{-3})$ 以下に抑えられ、反復解法の収束性が大幅に向上します。

4. 収束解析

• 正則性: 最適解 $\bar{y}$ は、領域の幾何学的性質に依存して $H^{2+\alpha}(\Omega)$ ($0 < \alpha < 1$) の正則性を持つことが知られています。
• 誤差評価: 提案手法の誤差を、メッシュ依存ノルム $\|\cdot\|_w$ および最適制御問題特有のノルム $|||\cdot|||_W$ で評価しました。
• 結果: 理論的に、誤差が $O(h^\tau)$ のオーダーで収束すること（$\tau$ は正則性パラメータに依存）が証明されました。これは既存の C0-IP 法と同様の収束オーダーですが、パラメータフリーである点で優れています。

5. 数値実験結果

ビハルモニック問題と、制約付き最適制御問題（一般的な追跡項を含む）に対して数値実験を行いました。

• ビハルモニック問題:
  • 提案手法は、メッシュ依存ノルムにおいて理論通りの 1 次収束（あるいはそれ以上の超収束）を示しました。
• 制約付き最適制御問題:
  • 精度: 提案された $C^0$-WG 法は、既存の C0-IP 法と比較して、$L^\infty$、$L^2$、$H^1$ ノルムにおいてより高い精度と収束次数を示しました。特に、追跡項による特異性の影響をよりよく捉えていることが確認されました。
  • 前処理の効果: 前処理なしのシステム行列の条件数はメッシュ細分化に伴い急激に増加しましたが、加法的シュワルツ前処理を適用することで、条件数が 2 桁以上減少し、重なり幅を大きくするほどさらに改善されました。これにより、大規模な問題でも効率的に解けることが実証されました。

6. 結論と意義

• 主な貢献:
  1. 4 階の変分不等式問題に対して、パラメータフリーで効率的な行列組立が可能な新しい $C^0$-WG 法を提案し、その厳密な誤差解析を行ったこと。
  2. 4 階問題特有の ill-conditioning を克服するための、$C^0$-WG 法に特化した加法的シュワルツ前処理を設計し、その条件数評価を行ったこと。
• 意義:
  • この手法は、パラメータ調整の不要なロバストな実装を可能にし、並列計算にも適したスケーラブルなソルバーを提供します。
  • 点ごとの制約や一般的な追跡項を含む複雑な最適制御問題に対して、高精度かつ効率的な数値解法として確立されました。
  • 2 次元だけでなく、3 次元への拡張も容易であることが示唆されており、時間依存問題や非線形問題への応用への基盤となっています。

総じて、この論文は、制約付き最適制御問題の数値解法において、既存の C0-IP 法の課題（パラメータ依存性、計算コスト）を克服し、理論的保証と高い数値性能を両立した画期的なアプローチを提示しています。