[論文レビュー] Data Assimilation In Large Prandtl Rayleigh-benard Convection From Thermal Measurements

Aseel Farhat, Nathan Glatt-Holtz+3 more|arXiv (Cornell University)|Mar 4, 2019

Fluid Dynamics and Turbulent Flows参考文献 59被引用数 38

ひとこと要約

本稿では、温度測定値のみを用いて、連続的データ同化スキームを大スケールのプラントル数を持つレイノルズ＝ベナール対流に適用する手法を提案している。ノッチングアルゴリズムを用いてモデルと観測値を同期化するが、理論的条件下では同期化が保証されるものの、数値結果では理論的上限が過剰に保守的であることが判明しており、緩和パラメータや観測分解能の閾値よりもはるかに低い値でも同期化が達成されることがある。

ABSTRACT

This work applies a continuous data assimilation scheme-a framework for reconciling sparse and potentially noisy observations to a mathematical model-to Rayleigh-Benard convection at infinite or large Prandtl numbers using only the temperature field as observables. These Prandtl numbers are applicable to the earth's mantle and to gases under high pressure. We rigorously identify conditions that guarantee synchronization between the observed system and the model, then confirm the applicability of these results via numerical simulations. Our numerical experiments show that the analytically derived conditions for synchronization are far from sharp; that is, synchronization often occurs even when sufficient conditions of our theorems are not met. We also develop estimates on the convergence of an infinite Prandtl model to a large (but finite) Prandtl number generated set of observations. Numerical simulations in this hybrid setting indicate that the mathematically rigorous results are accurate, but of practical interest only for extremely large Prandtl numbers.

研究の動機と目的

• 無限大または大きなプラントル数におけるレイノルズ＝ベナール対流に対して、温度観測値のみを用いた数学的に厳密なデータ同化フレームワークの構築を目的とする。
• モデルと観測系との間の同期化を保証するための、ノッチングパラメータおよび観測モード数の最小条件を同定すること。
• 理論的同期化境界が実用的適用性を有するかどうかを、無限大プラントル数および有限プラントル数の両領域における数値シミュレーションと比較することで評価すること。
• 観測が有限プラントル数系から得られるが、同化モデルが無限大プラントル数を仮定するというモデル不一致が、同期化精度に与える影響を調査すること。
• 無限大プラントル数モデルが有限プラントル数の現実にどの程度収束するかを評価し、モデルと観測の不一致が生じる状況での誤差を定量化すること。

提案手法

• 観測データを緩和項を通じて組み込むことで、モデル方程式を修正するノッチングアルゴリズムを用いた連続的データ同化スキームを定式化する。
• 温度場測定値のみを観測変数として用い、無限大プラントル数を仮定した2次元レイノルズ＝ベナール系にノッチングスキームを適用する。
• 緩和パラメータµと射影観測モード数Nに依存する理論的同期化条件を導出し、時間の経過とともにモデル状態が真の状態に収束することを保証する。
• 理論的境界を検証するための数値的シミュレーションを実施し、中程度の乱流状態を高精度ソルバおよび適応的時間刻み法で解明する。
• 観測が有限プラントル数シミュレーションから生成され、同化モデルが無限大プラントル数を仮定するハイブリッド設定を構築し、モデルと観測の不一致を検証する。
• 誤差の減衰率を定量化し、理論的予測と数値的結果を比較する。特に、ハイブリッド設定におけるO(Pr⁻¹)誤差境界の分析を行う。

実験結果

リサーチクエスチョン

• RQ1 無限大プラントル数極限において、緩和パラメータµおよび観測モード数Nがどのような条件下で、モデルと真の系との間で同期化が保証されるか。
• RQ2 理論的に導出された同期化条件はどの程度鋭く、数値実験はそれらを確認するか、あるいは超えるか。
• RQ3 有限プラントル数領域において、温度測定値のみを用いても同期化は達成可能か。また、無限大プラントル数の場合と比較して収束速度はいかがなっているか。
• RQ4 同化モデルが無限大プラントル数を仮定するが、観測が有限プラントル数系から得られるというモデルと観測の不一致が、同期化精度および誤差ダイナミクスに与える影響は何か。
• RQ5 無限大プラントル数モデルと有限プラントル数の現実との間の誤差は、プラントル数とともにどのようにスケーリングされるか。また、この誤差は動的同期化誤差を上回る主要因となるか。

主な発見

• 無限大プラントル数の場合の同期化の理論的条件は、相対的に大きな緩和パラメータµと十分な観測モード数Nを要するが、数値実験ではこれらの条件を満たさない場合でも同期化が達成されることが示された。
• 理論的同期化に要するµの値は、数値シミュレーションで成功した同期化が得られる値と比較して、少なくとも1桁以上大きいことが判明しており、理論的境界が過剰に保守的であることを示している。
• 有限プラントル数領域においても、温度測定値のみを用いた同期化は依然として達成可能であるが、収束は遅く、無限大プラントル数の場合と比べてパラメータ関係が明確でない。
• ハイブリッド設定（観測が有限プラントル数シミュレーションから、モデルが無限大プラントル数を仮定）では、同期化誤差がプラントル数の不一致によって支配され、プラントル数が増加するに従いO(Pr⁻¹)で減少する。
• 数値結果はハイブリッド設定における理論的O(Pr⁻¹)誤差境界を確認しており、無限大プラントル数を仮定することによるモデル誤差が、同期化失敗とは別に主要因であることが示された。
• 理論的厳密性は保たれているが、真のモデルが正確に分かっていない現実的状況、特にプラントル数や境界条件の不一致が生じる状況では、このデータ同化スキームの実用的適用性は制限される。