[論文レビュー] Control-constrained parabolic optimal control problems on evolving surfaces - theory and variational discretization

Morten Vierling|arXiv (Cornell University)|Jun 3, 2011

Advanced Numerical Methods in Computational Mathematics参考文献 33被引用数 3

ひとこと要約

本稿は、R^{n+1} 内の変化する超曲面上での制御制約付き線形二次最適制御問題に対して、ベクトル値分布および変分的離散化手法を用いて弱解の存在および一意性を確立する。最小限の正則性仮定のもとで、完全離散近似の最適収束率を証明し、数値例を通じて L² 誤差では収束順序が約 2、L∞ 誤差では 2.0 に近い実験的収束順序が観察された。

ABSTRACT

We consider control-constrained linear-quadratic optimal control problems on evolving surfaces. In order to formulate well-posed problems, we prove existence and uniqueness of weak solutions for the state equation, in the sense of vector-valued distributions. We then carry out and prove convergence of the variational discretization of a distributed optimal control problem. In the process, we investigate the convergence of a fully discrete approximation of the state equation, and obtain optimal orders of convergence under weak regularity assumptions. We conclude with a numerical example.

研究の動機と目的

• 変化する超曲面上における制御制約付き放物型最適制御問題の厳密な関数解析的枠組みを確立すること。
• ベクトル値分布および W(0,T)-型の解空間を用いて、状態方程式の弱解の存在および一意性を証明すること。
• 弱正則性仮定のもとで最適収束率を達成する完全離散変分的離散化スキームを構築・解析すること。
• 理論を制御制約付き最適制御問題に拡張し、変分的離散化を用いて実装可能な数値アルゴリズムを導出すること。
• 理論的結果を数値実験により検証し、最適収束順序を示すこと。

提案手法

• 変化する $ C^2 $-滑らか超曲面 $ \Gamma(t) \subset \mathbb{R}^{n+1} $ 上での物質微分および表面勾配を用いて、状態方程式を弱形式で定式化する。
• 固定領域へのプルバックを適用し、[LM68] の意味でのベクトル値分布を用いて弱解を定義し、やや弱い正則性のもとで存在および一意性を保証する。
• 制御制約付き最適制御問題に [Hin05] の意味での変分的離散化を適用し、完全離散的かつ実装可能な最適化スキームを可能にする。
• 弱正則性のもとで、状態方程式の誤差推定を $ L^2(0,T;L^2(\Gamma(t))) $-類似ノルムにおいて行い、最適収束率を達成する。
• 時間方向に不連続ガレルキン法、空間方向に有限要素法を用い、エネルギー技法および補間推定を用いて誤差境界を導出する。
• 制御のボックス制約を扱うために、離散最適制御問題にための半スムーズニュートン法を実装する。

実験結果

リサーチクエスチョン

• RQ1 制御制約付きの変化する曲面上における放物型状態方程式に対して、弱解が存在する条件は何か？
• RQ2 最小限の正則性仮定のもとで、状態方程式の完全離散近似に対して最適収束率を達成できるか？
• RQ3 変分的離散化は、変化する曲面上における制御制約付き最適制御問題の実装可能なアルゴリズムの構築をどのように可能にするか？
• RQ4 高正則性および低正則性データの下で、離散最適制御問題の実験的収束順序（EOC）はどのように変化するか？
• RQ5 離散設定において、$ L^2 $ と $ L^\infty $ ノルムの間で誤差推定および収束挙動はどのように異なるか？

主な発見

• 本稿は、分布的枠組みを用いて、変化する曲面上における放物型状態方程式の弱解の存在および一意性を証明し、[DE07] の先行研究を拡張した。
• 弱正則性仮定のもとで、状態近似における $ L^2(0,T;L^2(\Gamma(t))) $-ノルム誤差に対して、$ \mathcal{O}(h^2) $ の最適収束率が達成された。
• 数値実験により、$ L^2 $-誤差では実験的収束順序（EOC）が約 2.0、$ L^\infty $-誤差では 2.01～2.03 に近い値が得られ、最適収束挙動が確認された。
• 低正則性データ（例：$ y_T = 1/(x+y)^{0.45} $）に対しては、$ L^2 $-誤差の EOC は約 1.0 に低下し、最小可積分性のもとで予想される $ \mathcal{O}(h) $ 収束と整合的であった。
• 変分的離散化アプローチにより、完全に実装可能な最適化アルゴリズムが得られ、半スムーズニュートン法が離散制御問題を効果的に解いた。
• 誤差解析は、表面の時間的・空間的変化を完全に反映しており、三角形分割された表面におけるリフトおよび射影技法を用いて一貫性のある誤差境界が導出された。