[論文レビュー] Finite-time Partitions for Lagrangian Structure Identification in Gulf Stream Eddy Transport

Alexandre Fabregat, Igor Mezić+1 more|arXiv (Cornell University)|Jun 23, 2016

Fluid Dynamics and Turbulent Flows被引用数 9

ひとこと要約

本稿では、Koopman作用素理論的拡張と軌道統計のクラスタリングを用いて、非定常な海洋流れにおけるラグランジュ的 coherent 構造（LCS）を特定するための有限時間遍歴的（FiTER）パーティション法を提案する。この手法により、ガルフィールド・ストリーム渦の輸送障壁の3次元再構成と時間発展解析が可能となり、希なまたはノイズの多い速度データに対しても、強固でデータ駆動的な一貫性のある物質領域の同定が達成される。

ABSTRACT

We develop a methodology to identify finite-time Lagrangian structures from data and models using an extension of the Koopman operator-theoretic methods developed for velocity fields with simple (periodic, quasi-periodic) time-dependence. To achieve this, the notion of the Finite Time Ergodic (FiTER) partition is developed and rigorously justified. In combination with a clustering-based approach, the methodology enables identification of the temporal evolution of Lagrangian structures in a classic, benchmark, oceanographic transport problem, namely the cross-stream flux induced by the interaction of a meso- scale Gulf Stream Ring eddy with the main jet. We focus on a single mixing event driven by the interaction between an energetic cold core ring (a cyclone), the strong jet, and a number of smaller scale cyclones and anticyclones. The new methodology enab les reconstruction of Lagrangian structures in three dimensions and analysis of their time-evolution.

研究の動機と目的

• 非定常的かつ周期的でない海洋流れ、特にガルフィールド・ストリーム渦の相互作用の文脈において、有限時間ラグランジュ的 coherent 構造（LCS）を同定するデータ駆動的手法を開発すること。
• 従来のLCSおよびラグランジュ的記述子手法の限界を克服するため、遍歴理論と転送作用素に基づく、幾何学的でないメトリックに基づくアプローチを導入すること。
• 軌道集合の統計的類似性を活用することで、3次元的、希にサンプリングされた、またはノイズの多い速度場においても、一貫性のある物質領域を強固に同定すること。
• 環境輸送イベントのリアルタイムまたはニアリアルタイム解析に適した、数学的に厳密かつ計算的に実行可能なフレームワークを提供すること。
• 測度距離に基づく新しい極値原理を用いて、ほぼ不変集合の理論を有限時間 coherent セットへと拡張すること。

提案手法

• 有限時間遍歴的（FiTER）パーティションを、非定常な非自励的力学系への遍歴パーティション理論の厳密な拡張として提唱する。
• 初期条件分布の進化をモデル化するため、Perron-Frobenius転送作用素フレームワークを用い、位相空間におけるボックス化による近似を実施する。
• 有限時間 coherent セットを、軌道集合間の測度距離の最小化によって定義する：$ A_{ ext{maxcoh}} = \text{argmin}_{A \ni \beta} \frac{1}{2(\nu_A)^2} \bigint_{A \times A} d(\nu_{\bf x}^{t_0,T}, \nu_{\bf y}^{t_0,T}) d\nu({\bf x}) d\nu({\bf y}) $、ここで $ d $ は確率測度上のメトリック（例：Wasserstein）。
• 時間スケール関数 $ \tilde{\rho}_{\tau_s,j} $ と空間基底関数 $ f^m $ を用いたマルチスケール畳み込みフレームワークを採用し、軌道統計を計算することで、転送作用素の有限次元表現を構築する。
• ベクトル化された軌道統計 $ \tilde{f}_{\tau_l,m}^{\tau_s,j}({\bf x}) $ に対してクラスタリングを適用し、類似した統計的挙動を示す領域を同定することで、coherent セットへのパーティション分割を実現する。
• シミュレートされた速度場に埋め込まれた coherent 構造を有するベンチマーク3次元海洋流れ（ガルフィールド・ストリームリング相互作用）を用いて、手法の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 速度勾配やリャプノフ指数の計算に依存せずに、非定常的かつ周期的でない、希にサンプリングされた海洋流れにおいて、ラグランジュ的 coherent 構造をどのように強固に同定できるか？
• RQ2 非自励的系に一般化された遍歴パーティション理論へと拡張できる、有限時間でデータ駆動的なパーティショニング手法を構築できるか？ そして、一貫性のある、一意な coherent 物質領域の同定が可能か？
• RQ3 複雑な3次元海洋渦動力学において、FiTERパーティション手法は、従来のLCSやラグランジュ的記述子手法に比べて、どの程度優れた coherent 構造同定を達成できるか？
• RQ4 境界検出や幾何的仮定を必要とせずに、軌道集合の統計的類似性をどのように用いて coherent セットを定義できるか？
• RQ5 この手法は、地球物理的流れにおける異なる空間的・時間的スケールでの輸送プロセスのマルチスケール解析へと拡張可能か？

主な発見

• FiTERパーティション手法は、3次元ガルフィールド・ストリーム渦相互作用モデルにおいて、高精度でラグランジュ的構造の時間発展を再構成し、coherent 物質領域を同定することに成功した。
• 局所的な変形勾配に依存せず、軌道集合の統計的類似性に依存することで、速度データがノイズ混じりまたは空間的に希であっても、coherent セットの同定が強固に実現された。
• coherence 最小化原理にWassersteinメトリックを用いることで、異なる初期条件パーティションに対して物理的に意味があり、安定的かつ一意な coherent セットの定義が保証された。
• マルチスケール軌道統計のクラスタリングにより、複数の空間的・時間的スケールで coherent 構造を検出でき、流れ内に存在する階層的な混合パターンを明らかにした。
• スカラー平均化や正の関数の選択に起因する一意でない問題を回避する点で、LCSやラグランジュ的記述子の代替として数学的に厳密な代替手段を提供した。
• 実環境モニタリングにおけるリアルタイム応用へも一般化可能であることが、既知の coherent 動的特性を持つベンチマーク海洋輸送問題を用いて実証された。