[論文レビュー] Going with the Flow: a Lagrangian approach to self-similar dynamics and its consequences

D. G. Aronson, S. Betelú+1 more|ArXiv.org|Nov 23, 2001

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

ひとこと要約

本稿では、テンプレートベースの変換によって自己相似性を要因分解することで、非線形偏微分方程式から自己相似解を体系的に抽出するためのラグランジュ的で動的なピン留め手法——「MNダイナミクス」と呼ばれる——を導入する。この手法により、元の時間発展方程式が自己相似フレームにおける非局所的で時間に依存する方程式に変換され、第一種（バレンブラット型）および第二種（グラヴェロウ型）の解の両方の計算が可能になるとともに、自己相似指数も同時に得られる。

ABSTRACT

We present a systematic computational approach to the study of self-similar dynamics. The approach, through the use of what can be thought of as a ``dynamic pinning condition" factors out self-similarity, and yields a transformed, non-local evolution equation. The approach, which is capable of treating both first and second kind self-similar solutions, yields as a byproduct the self-similarity exponents of the original dynamics. We illustrate the approach through the porous medium equation, showing how both the Barenblatt (first kind) and the Graveleau (second kind) self-similar solutions arise in this framework. We also discuss certain implications of the dynamics of the transformed equation (which we will name "MN-dynamics"); in particular we discuss the discrete-time implementation of the approach, and connections with time-stepper based methods for the "coarse" integration/bifurcation analysis of microscopic simulators.

研究の動機と目的

• 非線形時間発展方程式における自己相似解を体系的に特定する計算手法の開発。
• 従来、移動波に用いられてきたテンプレートベースの動的ピン留め法を、自己相似ダイナミクスへと拡張すること。
• 連続時間における自己相似性の動的要因分解により、自己相似状態への到達を明らかにする変換された時間発展方程式を導出すること。
• 第一種（例：バレンブラット）および第二種（例：グラヴェロウ）の自己相似解を統一的な枠組み内で計算可能とする。
• 巨視的方程式の明示的導出を避けることで、微視的シミュレータの時間ステッパーに基づく粗視化された分岐解析を可能とする。

提案手法

• 自己相似性を要因分解するためのテンプレートベースの動的ピン留め条件を適用し、元のPDEを自己相似フレームにおける非局所的時間発展方程式に変換する。
• スケーリング変換を用いて、自己相似依存性を除去する新たな変数を定義し、これにより「MNダイナミクス」方程式が得られる。
• 変換された方程式を数値的に解くことで、自己相似解およびそれに関連する指数を抽出する。
• 連続時間および離散時間の両方の実装を可能とし、微視的ダイナミクスの時間ステッパーに基づくシミュレーションへの応用を可能にする。
• 変換されたダイナミクスを活用し、動的システムの手法（例：分岐理論、連続法）を自己相似解に適用する。
• 離散時間版を用いて、巨視的方程式を明示的に構築せずに、微視的シミュレータ（例：モンテカルロ法、格子ボルツマン法）における自己相似性を解析する。

実験結果

リサーチクエスチョン

• RQ1 移動波に用いられる動的ピン留めと類似する連続時間における自己相似性の動的要因分解はどのように可能か？
• RQ2 同じテンプレートベースのアプローチを、並進不変性から自己相似性へと一般化可能か？これにより非局所的時間発展方程式が得られるか？
• RQ3 第一種および第二種の自己相似解は、変換されたダイナミクスからどのように自然に出現するか？
• RQ4 自己相似解をMNダイナミクスの定常状態として扱うことで、数値的分岐解析にどのような影響を与えるか？
• RQ5 この手法は、巨視的方程式が得られない、またはあまりに複雑すぎて導出が困難な系の微視的シミュレータ（例：分子動力学、モンテカルロ法）に適用可能か？

主な発見

• MNダイナミクスフレームワークは、多孔質媒体方程式に対して、バレンブラット（第一種）およびグラヴェロウ（第二種）の自己相似解を効果的に計算できた。
• 解法プロセスの副産物として、自己相似指数が自動的に得られた。
• 変換された方程式により、連続法や分岐理論などの高度な動的システム技法を自己相似解に適用可能となった。
• 離散時間版の実装により、巨視的MN方程式を明示的に構築せずとも、微視的時間ステッパーの有効な分岐解析が可能となった。
• 本手法により、巨視的方程式が存在しない、あるいはあまりに複雑すぎて導出が困難な系における自己相似ダイナミクスのコンピュータ支援的解析への道筋が開かれた。
• 本手法はマルチスケール積分器（例：プロジェクティブ積分法）と互換性があり、微視的シミュレーションにおける物理的制約を組み込んだアンサンブルを定義するためにも利用可能である。