Distributed optimal control problems governed by poroelasticity equations

技術要約：多孔質弾性方程式によって支配される分散型最適制御問題

問題設定
本論文は、ビオの圧密モデル（Biot's consolidation model）に基づく、多孔質弾性における分散型最適制御問題の数値近似について扱う。この物理系は、間隙流体の流れと多孔質媒体の変形との相互作用を記述するものであり、地球物理学、生体組織力学、および炭素貯留などの応用において極めて重要なモデルである。特定の制御問題は、固体変位（$u$）および流体圧（$p$）を目標値（$u_C, p_C$）から偏差させる度合いを測定し、同時に制御変数の大きさを罰則化することを目的とする。制御変数は、領域内に分布する流体のソース／シンク項（$m$）として特定される。本研究は、加速度効果が無視できる準静的な圧密に焦点を当てており、特に、比貯留係数（$s_0$）がゼロに近い場合に、流体と弾性媒体がほぼ非圧縮となる「ポロエラスティック・ロッキング（poroelastic locking）」に関連する数学的および数値的な課題に対処している。

手法
著者らは、主要な未知数として固体変位と流体圧を用いる、新しい二変数対称定式化を提案し、分析を行っている。その手法は以下の段階を経て進められる：

1. 連続定式化およびウェル・ポーズドネス（適当性）：

   • 著者らは、運動量バランスの時微分を取ることで、系を一連の微分方程式へと書き換えることにより、支配方程式を再定式化する。
   • 適切な関数空間（変位のための$V$、圧力のための$Q$）を定義し、鞍点理論およびLadyzhenskaya–Babuška–Brezzi（LBB）条件を用いて、弱解のウェル・ポーズドネスを確立する。
   • コスト関数の凸性と状態演算子の線形性を利用して、最適制御の存在と一意性を、変分法の直接法を用いて証明する。
   • 随伴状態アプローチを通じて、一次の必要かつ十分な最適条件を導出し、制御に関する変分不等式を得る。

2. 離散スキーム：

   • 空間離散化： 領域は、inf-sup条件を満たす有限要素（$V_h \times Q_h$）（例：MINI要素またはTaylor-Hood要素）を用いて離散化される。
   • 時間離散化： 後退オイラー法（具体的には時間に関するdG(0)法）が採用され、これにより完全離散化された系が得られる。
   • 制御の離散化： 論文では、制御変数に対して変分離散化アプローチ（Hinzeによって提案されたもの）を利用している。これは、制御空間が明示的に離散化されるのではなく、離散制御が離散的な随伴状態を許容制御集合へと投影することによって得られることを意味する。これにより、制御空間を離散化することによって通常導入される近似誤差を回避できる。

3. 誤差解析：

   • 著者らは、完全離散スキームに対するア・プリオリ誤差評価を導出する。
   • 解析では、総誤差を、離散化誤差（固定された制御に対するもの）と、連続的な最適制御と離散的な最適制御の差による摂動誤差に分離する。
   • 前方状態系および後方随伴系の両方に対して安定性評価を確立し、これらの誤差を抑え込む。

主な貢献
本論文は、文献に対して主に3つの貢献を行っている：

• 新しい対称定式化： ビオモデル内の最適制御問題に特化した、二変数対称定式化（変位と圧力）を導入している。これは、従来の四変数または混合定式化とは異なるものである。
• 存在と一意性： 新しい定式化のウェル・ポーズドネスを厳密に証明し、最適制御の存在と一意性を確立するとともに、一次の最適条件を導出した。
• ア・プリオリ誤差評価： 完全離散スキームに対する完全なア・プリオリ誤差解析を提供した。特筆すべきは、$L^\infty(0, T; H^1(\Omega))$ および $L^2(0, T; L^2(\Omega))$ ノルムにおける変位と圧力の誤差、および $L^2(0, T; L^2(\Omega))$ ノルムにおける制御変数の誤差を導出し、標準的な正則性の仮定の下で最適な収束率を達成している点である。

結果

• 理論的側面： 本論文は、数値スキームが $O(h^s + \Delta t)$ のオーダーで収束することを証明している（ここで $h$ は空間メッシュサイズ、$\Delta t$ は時間ステップ、$s$ は解の正則性に関連する）。具体的には、解が十分に正則である場合（$s=1$）、空間と時間の両方において線形な収束（$O(h + \Delta t)$）が得られる。
• 数値的側面： 単位正方形領域における数値実験により、理論的知見を検証している。オープンソースライブラリであるFEniCSを用い、前方向・後方向アルゴリズム（forward-backward algorithm）を用いて最適系を解いた。結果は、提案されたスキームが、最低次のMINI要素およびTaylor-Hood要素を使用した場合に、変位、圧力、および制御変数の予測された最適収束率を達成することを実証している。

意義
著者らは、本研究を、多孔質弾性によって支配される最適制御問題の数値解析における基礎的な一歩として位置づけている。彼らは、流体の流れ（貯留係数が消失するという仮定の下でのBociuらによる研究など）の最適制御は研究されてきたものの、一般的なビオモデルにおけるこのような問題の数値近似については、これまで調査されてこなかったと指摘している。新しい対称定式化のウェル・ポーズドネスを確立し、完全離散スキームに対する厳密な誤差評価を提供することで、本論文は文献の空白を埋め、生体組織モデリングや地球物理学的シミュレーションなどの応用における、多孔質弾性プロセスの制御のための堅牢な理論的枠組みを提供している。