[論文レビュー] Computing Coherent Sets using the Fokker-Planck Equation

Andreas Denner, Oliver Junge+1 more|arXiv (Cornell University)|Dec 11, 2015

Statistical Mechanics and Entropy参考文献 18被引用数 4

ひとこと要約

本稿では、確率的摂動下での転送作用素を直接近似するためにFokker-Planck方程式を解くことで、動的系における一貫集合を計算する新しい手法を提案する。スペクトルコロケーションと指数時刻差分法を用いることで、高精度かつ高効率な計算が実現される。本手法は軌道積分を回避し、Fokker-Planckに基づく転送作用素の特異ベクトル解析により、一貫集合の検出を安定的に行う。

ABSTRACT

We perform a numerical approximation of coherent sets in finite-dimensional smooth dynamical systems by computing singular vectors of the transfer operator for a stochastically perturbed flow. This operator is obtained by solution of a discretized Fokker-Planck equation. For numerical implementation, we employ spectral collocation methods and an exponential time differentiation scheme. We experimentally compare our approach to the more classical method by Ulam that is based on integration of the exact transfer operator.

研究の動機と目的

• 有限次元の動的系における一貫集合を特定する数値的手法を開発すること。この手法は微小な確率的摂動に対して頑健であるべきである。
• 古典的な軌道ベースの手法に代わる、Fokker-Planck方程式を用いて確率的転送作用素をモデル化するPDEベースのアプローチを提示すること。
• スペクトルコロケーションによる直接的な広義方程式の解法により、数値的効率を向上させ、数値的拡散を低減すること。
• ラグランジュ的軌道積分やグリッド外補間を必要とせず、一貫集合を高精度に計算可能にすること。
• 古典的手法であるUlam法と比較して、一貫構造の検出において優れた性能を示すことを検証すること。

提案手法

• 本手法は、確率的摂動を受ける流れにおける確率密度の時間発展をモデル化する確率的転送作用素を、Fokker-Planck方程式を解くことで得る。
• スペクトルコロケーションを用いてFokker-Planck方程式を離散化し、固定された粗いグリッド上で高次の空間的精度を達成する。
• Fokker-Planck方程式に含まれるラプラシアン項に起因する剛性を扱うために、指数時刻差分（etd）スキームを採用する。
• 高広義領域におけるアリーキングエラーを軽減するために、広義項の歪対称形を用いる。
• 得られた転送作用素は、基底関数に対するFokker-Planck解の時間積分により構築され、その特異ベクトルを計算することで一貫集合を特定する。
• 本手法は明示的な軌道積分や補間を回避し、グリッド上のベクトル場のサンプリングにのみ依存する。

実験結果

リサーチクエスチョン

• RQ1 Fokker-Planck方程式は、動的系における一貫集合を直接計算するための有効なツールとして用いることができるか？
• RQ2 Fokker-Planckに基づく転送作用素は、古典的手法のUlam法と比較して、精度と数値的安定性において優れているか？
• RQ3 スペクトルコロケーションと指数時刻積分は、一貫集合検出における軌道ベースの手法に代わる、より効率的かつ高精度な代替手段を提供できるか？
• RQ4 Galerkin法や有限差分法と比較して、Fokker-Planck方程式を直接解くことで、数値的拡散がどの程度低減されるか？
• RQ5 Fokker-Planck方程式の直接解法により得られる一貫集合は、幾何学的に意味を持つものであり、摂動に対して頑健であるか？

主な発見

• Fokker-Planckに基づく転送作用素はコンパクトであり、$L^2$ 上の有界集合を $H^1$ 上の有界集合に写像するため、特異ベクトルの正しく定義された存在と収束が保証される。
• 本手法は粗いグリッドでも高精度を達成でき、ラグランジュ的軌道積分やグリッド外補間の必要がない。
• 数値実験の結果、Fokker-Planckアプローチは特に高広義領域において、古典的手法のUlam法よりも一貫構造をより頑健に検出できる。
• スペクトルコロケーションと指数時刻積分の使用により、標準的な有限差分法やGalerkin法と比較して、数値的拡散とアリーキングエラーが顕著に低減される。
• 四重渦流れの例において、本手法は一貫集合を適切に同定でき、特異ベクトルが既知の輸送障壁に沿った空間的一致性のある領域を示している。
• 理論的解析により、Fokker-Planckに基づく作用素が正の性質や確率的性質といった重要な性質を保持していることが確認され、物理的整合性が保たれている。