[論文レビュー] Continuous data assimilation for the magnetohydrodynamic equations in 2D using one component of the velocity and magnetic fields

Animikh Biswas, Joshua Hudson+2 more|arXiv (Cornell University)|Apr 7, 2017

Advanced Mathematical Modeling in Engineering参考文献 47被引用数 3

ひとこと要約

本稿では、2次元磁気流体力学（MHD）方程式に対する連続的データ同調法を提案し、速度または磁場の1成分という部分観測に基づくフィードバック制御を用いる。適切な制御パラメータおよび観測解像度の条件下で、誤差が $L^2$ および $H^1$ ノルムにおいて指数的収束することを証明している。これは、フィードバック制御が1成分または水平成分にのみ適用されても成立する。

ABSTRACT

We propose several continuous data assimilation (downscaling) algorithms based on feedback control for the 2D magnetohydrodynamic (MHD) equations. We show that for sufficiently large choices of the control parameter and resolution and assuming that the observed data is error-free, the solution of the controlled system converges exponentially (in $L^2$ and $H^1$ norms) to the reference solution independently of the initial data chosen for the controlled system. Furthermore, we show that a similar result holds when controls are placed only on the horizontal (or vertical) variables, or on a single Elsässer variable, under more restrictive conditions on the control parameter and resolution. Finally, using the data assimilation system, we show the existence of abridged determining modes, nodes and volume elements.

研究の動機と目的

• 2次元MHD方程式に対する部分観測に基づく連続的データ同調法の開発。
• 同調された解が $L^2$ および $H^1$ ノルムにおいて真の解へ指数的収束することの確立。
• 例えば水平速度または1つのエールサッター変数に制限された1成分へのフィードバック制御でも収束が保証されることの示唆（ただし、より制限の厳しい条件を要する）。
• データ同調とMHD系における決定的モード、ノード、体積要素との関連の確立。
• ナビエ＝ストークス方程式におけるフィードバック制御フレームワークをMHD方程式へ拡張し、厳密な収束解析を実施。

提案手法

• 偏微分方程式（PDE）レベルでフィードバック制御を用い、部分観測に基づきMHD系の解を補正する。
• アズアニ、オルソン、ティティが提唱した連続的データ同調フレームワークを2次元MHD方程式へ応用する。
• 解析を簡素化し、系のダイナミクスを分離するためにエールサッター変数を用いる。
• 速度または磁場の1成分、あるいは1つのエールサッター変数にのみフィードバック制御を適用する。
• エネルギー推定と対数的ソボレフ不等式を用い、誤差ダイナミクスにおける非線形項を評価する。
• 指数的収束率を確立するために、$L^2$ および $H^1$ ノルムにおける事前推定を導出する。

実験結果

リサーチクエスチョン

• RQ1 速度または磁場の1成分へのフィードバック制御を用いた連続的データ同調が、2次元MHD方程式の真の解へ指数的収束を保証できるか。
• RQ2 部分観測の下で指数的収束を達成するための最小限の観測要件（解像度および制御パラメータ）は何か。
• RQ3 MHD用のフィードバック制御フレームワークは、既存のカルマンフィルターや変分法と比較して、収束保証においてどのように異なるか。
• RQ4 データ同調システムを用いて、MHD系における決定的モード、ノード、体積要素の存在を厳密に示すことができるか。
• RQ5 部分観測の下で $H^1$ ノルム収束を保証するための制御パラメータおよびデータ間隔の条件は何か。

主な発見

• すべての変数が観測される場合、制御パラメータおよび観測解像度に適切な条件が満たされれば、フィードバック制御系において $L^2$ ノルムでの指数的収束が確立される。
• フィードバックが速度および磁場の水平成分にのみ適用される場合、制御パラメータおよび解像度により制限の厳しい条件が要求されるが、依然として $L^2$ ノルムでの指数的収束が成立する。
• 標準的な仮定の下で、全系について $H^1$ ノルムでの指数的収束が証明され、水平成分のみの観測の場合でも、補間関数に関する弱い仮定の下で同様に成立する。
• データ同調フレームワークを用いて、短縮版の決定的モード、ノード、体積要素の存在を確立した。
• 収束速度は指数的であり、誤差は $e^{- u t}$（$ u > 0$）の形で減少する。この係数は制御パラメータおよび系のパラメータに依存する。
• 初期データの不確実性に対しても結果は頑健であり、制御系の初期条件に依存しない。