[論文レビュー] Model dynamics on a multigrid across multiple length and time scales

A. J. Roberts|ArXiv.org|Feb 8, 2008

Advanced Mathematical Modeling in Engineering参考文献 29被引用数 5

ひとこと要約

本稿では、有限なスケール分離を持つマルチグリッド階層において、多スケールの対流拡散ダイナミクスを厳密にモデル化するための動的システムフレームワークを提示する。中心多様体理論を用いて、各レベルにおける一貫性のある粗スケールモデルを導出する。この手法により、安定性と精度を保ちながら、2倍の要因による反復的粗化が可能となり、対流に起因する拡散の増強を含む解析的に導出された有効係数が得られる。

ABSTRACT

Most methods for modelling dynamics posit just two time scales: a fast and a slow scale. But many applications, including many in continuum mechanics, possess a wide variety of space-time scales; often they possess a continuum of space-time scales. I discuss an approach to modelling the discretised dynamics of advection and diffusion with rigorous support for changing the resolved spatial grid scale by just a factor of two. The mapping of dynamics from a finer grid to a coarser grid is then iterated to generate a hierarchy of models across a wide range of space-time scales, all with rigorous support across the whole hierarchy. This approach empowers us with great flexibility in modelling complex dynamics over multiple scales.

研究の動機と目的

• 極端なスケール分離（$\epsilon \to 0$）を必要とせず、有限なスペクトルギャップを有する複数の長さスケールおよび時間スケールにおけるダイナミクスを数学的に厳密にモデル化する手法を開発すること。
• 従来の二スケールモデリングを、各レベルが2倍に粗化される階層的マルチグリッドフレームワークに拡張すること。
• 特に中心多様体理論を用いた厳密な動的システム理論により、粗スケールモデルの安定性と精度を保証すること。
• 制限演算子と延長演算子を用いて、スケール間で初期条件を体系的に射影し、ダイナミクスを保存する方法を提供すること。
• 線形および弱非線形系（バーガース方程式を含む）への適用可能性を示すこと。

提案手法

• 中心多様体理論を用いて、細粒度グリッドのダイナミクスを粗いグリッドに簡約することの正当性を示し、ゆっくりな多様体への指数的吸引を保証する。
• 幾何的アプローチを用いて、粗いグリッド上での有効な対流および拡散係数を導出し、明示的な式 $\bar{c} = \frac{1}{2}c$ および $\bar{d} = \frac{1}{4}d + \frac{c^2}{16d}$ を得る。
• コンピュータ代数を用いてゆっくりな多様体を近似し、マルチスケール展開における高次補正を計算する。
• 複数のレベルにわたり、粗化変換を反復的に適用し、細かいスケールから粗いスケールへのモデルの階層を生成する。
• 細粒度グリッド値から粗粒度グリッド値への写像を $\bar{u}_j = u_{2j}$ で定義する制限演算子を用い、細分化には延長を用いる。
• 摂動解析を用いて非線形ダイナミクスへ拡張し、各レベルで制御された誤差項を有する有効方程式を導出する。

実験結果

リサーチクエスチョン

• RQ1 有限なスペクトルギャップを有する複数スケールにおけるダイナミクスを、$\epsilon \to 0$ を要求しない厳密な数学的フレームワークとして開発できるか？
• RQ2 離散的対流拡散方程式のダイナミクスを、安定性と精度を保持したまま2倍に一貫して粗化できるか？
• RQ3 細粒度グリッドのダイナミクスから導かれる粗スケールモデルを安定化するための拡散の増強の役割は何か？
• RQ4 制限演算子を用いて、スケール間で初期条件を一貫して射影する方法は何か？
• RQ5 この手法を多数のレベルにわたり反復適用し、制御された誤差を有する完全なマルチグリッド階層のモデルを生成できるか？

主な発見

• グリッド間隔 $\bar{h} = 2h$ の対流拡散の粗スケールモデルは、$\frac{d\bar{u}_j}{dt} = -\bar{c}\frac{1}{2}(\bar{u}_{j+1} - \bar{u}_{j-1}) + \bar{d}(\bar{u}_{j+1} - 2\bar{u}_j + \bar{u}_{j-1})$ で与えられ、$\bar{c} = \frac{1}{2}c$ および $\bar{d} = \frac{1}{4}d + \frac{c^2}{16d}$ を満たし、安定性と精度が保証される。
• 拡散項 $\frac{c^2}{16d}$ は微細スケールのダイナミクスを解明することに起因し、循環削減法と一致するため、安定化の物理的根拠を提供する。
• この手法は反復的に適用可能である：レベル0からレベル1へ、さらにレベル2へと変換することで、有効な散逸が増加し、非線形項が $\alpha^2 u^2$ のスケーリングを示す一貫性のあるモデルが得られる。
• レベル2では、非線形項 $\frac{5}{64}\alpha^2 u^2$ が線形散逸項を上回り、より大きなスケールでのダイナミクスに本質的な変化が生じることを示している。
• 反復処理の間で誤差が蓄積せず、中心多様体理論により、すべてのレベルのモデルが指数的吸引性を示す。
• 制限を用いた初期条件の射影が可能であり、内部要素構造はウェーブレットに類似しており、ダイナミクスのマルチスケール表現を可能にする。