[論文レビュー] Analysis and Computation of a Discrete KdV-Burgers Type Equation with Fast Dispersion and Slow Diffusion

Zvi Artstein, C. W. Gear+3 more|ArXiv.org|Aug 19, 2009

Numerical methods for differential equations参考文献 11被引用数 4

ひとこと要約

本稿では、角変数が明示的に与えられていないため古典的平均化法が失敗する、高速分散と低速拡散を有する離散KdV-Burgers方程式の長時間挙動を解析するための新しい計算フレームワークを構築する。不変測度とLax不変量の平均化を用いて、明示的な角変数を必要とせずに有効な低速ダイナミクスを導出する。この手法により、高速系の不変量とヤング測度平均化を活用したプロジェクト型積分が可能となり、直接シミュレーションと比較して40:1を超える高速化を達成する。

ABSTRACT

The long time behavior of the dynamics of a fast-slow system of ordinary differential equations is examined. The system is derived from a spatial discretization of a Korteweg-de Vries-Burgers type equation, with fast dispersion and slow diffusion. The discretization is based on a model developed by Goodman and Lax, that is composed of a fast system drifted by a slow forcing term. A natural split to fast and slow state variables is, however, not available. Our approach views the limit behavior as an invariant measure of the fast motion drifted by the slow component, where the known constants of motion of the fast system are employed as slowly evolving observables; averaging equations for the latter lead to computation of characteristic features of the motion. Such computations are presented in the paper.

研究の動機と目的

• 古典的平均化法が角変数の明示的表記が不足するため失敗する、高速分散と低速拡散を有する離散KdV-Burgers型方程式の長時間挙動を解析すること。
• 高速・低速成分の明示的分離や高速ダイナミクスの周期軌道の知識を必要としない計算フレームワークを構築すること。
• 高速系のLax不変量と不変測度を低速観測量として用い、低速拡散過程を支配する有効方程式を導出すること。
• これらの低速観測量に基づくプロジェクト型積分が、直接数値積分と比較して顕著な計算高速化（40:1以上）を実現できることを示すこと。
• ヤング測度平均化と方程式フリー手法の両方を用いて手法を検証し、直接シミュレーションと区別できない結果が得られることを示すこと。

提案手法

• 系はKdV-Burgers PDEの半離散近似としてモデル化され、非線形対流項による高速分散と粘性項による低速拡散が特徴である。
• 高速ダイナミクスはLaxペア表現を用いて解析され、第一積分（Lax不変量）の豊富な族が得られ、これが低速観測量として機能する。
• 平均化は状態変数そのものではなく、これらの不変量の時間平均値に対して実行され、高速フローの不変測度を用いて有効なドリフトを計算する。
• ヤング測度理論を用いて平均化プロセスを厳密に正当化し、極限測度が高速振動の統計的挙動を捉える。
• プロジェクト型積分は低速観測量に適用される：各不変量の時間微分は1つの高速周期内での数値差分で近似され、その後大ステップで時間発展が行われる。
• 方程式フリー手法は、1つの高速周期にわたって全系をシミュレートし、低速観測量の傾きを計算して、低速方程式を明示的に知らずに時間前方に投影する。

実験結果

リサーチクエスチョン

• RQ1 角変数が明示的に与えられていない場合、高速分散と低速拡散を有する離散KdV-Burgers系の有効ダイナミクスをどのように導出できるか？
• RQ2 高速系のLax不変量は、長時間挙動を支配する有効方程式を構築するための信頼できる低速観測量として機能するか？
• RQ3 不変測度と観測量の数値差分に基づくプロジェクト型積分は、直接数値積分と比較してどの程度計算速度で優れるか？
• RQ4 ヤング測度平均化と方程式フリー手法は、系の真の低速ダイナミクスをどの程度正確に再現するか？
• RQ5 直接シミュレーションと比較して、不変量に基づくプロジェクト型積分によって達成可能な定量的高速化率はどの程度か？

主な発見

• 直接積分と比較して、慢性的拡散係数νを10倍小さくした場合、計算高速化率が40:1を超える。
• 1周期分の時間微分を用いた方程式フリー手法によるプロジェクト型積分は、直接シミュレーションおよびヤング測度平均化とプロット精度で区別できない結果を生成した。
• ν = 0.001の場合、低速観測量の減衰を10⁻⁵のスケールで追跡するには、直接シミュレーションで約23,000ステップを要したが、プロジェクト型積分では約5,000ステップで十分であった。
• 低速観測量の時間微分は1つの高速周期内での数値差分により近似され、t = 1,492.08で10³スケールの誤差がv₃に対して5.9未満であった（12回のプロジェクトステップ使用時）。
• 高速系の不変測度（式3.1）はトーラス上に台を持つことが判明し、その低速拡散下での減衰は図7に可視化され、系の統計的構造の緩やかな進化が確認された。
• Lax不変量v₁からv_{N/2+1}が低速観測量として効果的に使用され、平均化法と方程式フリー手法の両方でその低速拡散下での時間発展が正確に捉えられた。