[論文レビュー] Guaranteed optimal reachability control of reaction-diffusion equations using one-sided Lipschitz constants and model reduction

Adrien Le Coënt, Laurent Fribourg|arXiv (Cornell University)|Jul 25, 2019

Model Reduction and Neural Networks参考文献 39被引用数 4

ひとこと要約

本稿では、明示的オイラー時刻離散化と片側リプシッツ（OSL）定数を用いた1次元反応拡散方程式の保証された最適制御手法を提示する。OSL定数が負である場合、時間区間長に線形に依存する収束性と誤差境界を証明する。さらに、大規模なシステムにおける計算効率を向上させるためのモデル還元技術を導入する。

ABSTRACT

We show that, for any spatially discretized system of reaction-diffusion, the approximate solution given by the explicit Euler time-discretization scheme converges to the exact time-continuous solution, provided that diffusion coefficient be sufficiently large. By "sufficiently large", we mean that the diffusion coefficient value makes the one-sided Lipschitz constant of the reaction-diffusion system negative. We apply this result to solve a finite horizon control problem for a 1D reaction-diffusion example. We also explain how to perform model reduction in order to improve the efficiency of the method.

研究の動機と目的

• 非線形反応拡散PDEの有限区間最適制御における到達可能性を保証する課題に対処する。
• 負の片側リプシッツ定数を活用することで、時刻離散化された解の数値的安定性と収束性を保証する。
• 明示的オイラー積分法とモデル順序還元を組み合わせることで、計算効率の高い手法を開発する。
• 時間区間長に線形に依存する厳密な誤差境界を提供し、先行研究の指数的誤差境界を改善する。
• 射影に基づくモデル還元により、空間離散化が大きい場合の最適制御の適用を可能にする。

提案手法

• 反応拡散PDEを空間的半離散化により常微分方程式系に変換する。
• 明示的オイラー時刻積分を用いて数値解を求める。OSL定数が負である場合、収束が保証される。
• 片側リプシッツ（OSL）定数を用いて、時間区間長Tに線形に依存する事前誤差境界を導出する。
• 適切な直交分解（POD）を用いてモデル還元を行い、高次元系を低次元部分空間に射影する。
• 最適制御戦略を低次元系で合成し、誤差が有界であるように元の全系に適用する。
• 射影に基づく誤差推定を用いて、全系と還元系の軌道間の乖離を評価し、誤差は $ \frac{K_2 \sigma}{|\lambda_{h_1}|} $ に比例する。

実験結果

リサーチクエスチョン

• RQ1 特定の条件下で、明示的オイラー時刻離散化が反応拡散PDEの真の解に収束することが保証されるか？
• RQ2 片側リプシッツ定数は、PDEの有限区間最適制御における誤差境界にどのように影響するか？
• RQ3 反応拡散系にモデル還元を適用することで、計算コストを削減しつつ制御性能を維持できるか？
• RQ4 還元系と全系の軌道間の誤差を、保証された制御合成を可能にする形で定量的に境界づけることができるか？
• RQ5 提案手法は、先行手法とは対照的に、時間区間長に線形に依存する誤差スケーリングを達成するか？

主な発見

• 拡散係数が十分に大きく、片側リプシッツ定数が負である場合、明示的オイラースキームは反応拡散PDEの正確な解に収束する。
• 最適制御解の誤差境界は時間区間長Tに線形に依存するが、先行研究で得られた $ O(e^{L_f T}) $ の指数的依存とは対照的である。
• 射影に基づくモデル還元により、システムサイズを顕著に削減できる。例えば、$ M_2 = 10 $ から $ M_1 = 5 $ にまで縮小され、グリッドサイズは約 $ 15^5 \approx 7.6 \times 10^5 $ 倍にまで減少し、大規模制御が可能になる。
• 還元誤差は $ \frac{K_2 \sigma}{|\lambda_{h_1}|} $ で有界であり、$ \sigma $ はPODの切断誤差を表す。この境界により、射影された軌道が目標の近傍に常に存在することが保証される。
• 1次元二安定反応拡散系に対する実験結果では、還元モデルから導出された制御戦略が、全系を目標状態に的確に誘導することが確認され、シミュレーション結果も全系とよく一致する。
• 事後的シミュレーションでは、理論的誤差境界が保守的であることが示され、実際の軌道間距離は上界よりも著しく小さい。