[論文レビュー] Exact coherent states underlying chaotic falling-film dynamics

Isaac J. G. Lewis, C. Ricardo Constante-Amores|arXiv (Cornell University)|Mar 7, 2026

Fluid Dynamics and Thin Films被引用数 0

ひとこと要約

この論文は、落下膜の二次元長波動方程式を導出し、その領域をマッピングし、DNSを低次元慣性多様体に射影し、DNSと縮約モデルの双方でニュートン・クライロフ法を用いて、混沌とした界面力学に埋め込まれた正確なコヒーレント状態を同定する。

ABSTRACT

Dynamical-systems approaches to spatiotemporal chaos have been developed primarily for single-phase flows, where the system state is defined by bulk velocity fields. Extending these ideas to two-phase flows remains challenging because the dynamics are intrinsically coupled to the evolution of a deforming interface. Here, we address this challenge for a two-dimensional vertical falling film by formulating the dynamics in terms of the interface evolution. Starting from the Navier--Stokes equations, we recover a classical long-wave interface evolution equation, originally derived by Topper & Kawahara (1978). Using this formulation, we perform an extensive parametric study to construct a regime map in the space of domain size and dispersion parameter. The resulting map reveals a rich range of interfacial behaviors, including travelling waves, bursting travelling waves, and fully chaotic regimes. In the chaotic falling film regime, we exploit the dissipative nature of the governing equation, which suggests that the long-time dynamics evolve onto an inertial manifold. Using a data-driven approach, we parameterize this inertial manifold and estimate its intrinsic dimension, suggesting approximately linear growth with domain size. We then construct low-dimensional models in manifold coordinates to facilitate the search for exact coherent states of the full system. Using this approach, we identify travelling waves, relative periodic orbits and equilibria embedded within the chaotic attractor. Chaotic trajectories repeatedly approach the neighbourhoods of these invariant solutions, indicating that the recurrent interfacial patterns observed in the dynamics correspond to visits to these coherent states. To the best of our knowledge, this constitutes the first identification of exact coherent structures embedded in chaotic falling-film dynamics.

研究の動機と目的

• 二相落下膜流の時空的混沌をダイナミカルシステムの視点から理解する動機付け。
• 界面力学の低次元表現を慣性多様体上に構築する。
• 混沌とした落下膜力学に埋め込まれた正確なコヒーレント状態を同定・分類する。
• DNS、データ駆動型多様体学習、ニュートン-クライロフの改良を組み合わせるワークフローを示す。

提案手法

• 厚さ H(x,y,t) の二次元長波方程式を分散パラメータ delta を用いて導出・適用する。
• スペクトルDNSソルバーを用いて二重周期境域で方程式を解き、スペクトル解析とエネルギー指標で領域を分析する。
• POD に続く暗黙的ランク最小化オートエンコーダー(IRMAE-WD)で高次元DNSデータを縮約し、次元 d_M の慣性多様体を同定する。
• 多様体上のダイナミクスを学習するニューラルODEを用いて潜在座標での近似時間発展を得る。
• 対称性の削減とヤコビ行列を介さないニュートン-クライロフ法(JFNK)を用いて、DNSと縮約座標の丛において正確なコヒーレント状態（平衡解、移動波、周期的軌道、相対周期的軌道）を同定し、最終収束のためDNSで精密化する。

Figure 1 : Streamwise and spanwise power spectra $P(k_{x})$ and $P(k_{y})$ for cases $(L=22,\delta=0.002)$ and $(L=30,\delta=0.65)$ corresponding to panels (a) and (b), respectively. Spectra are shown for Fourier discretizations with $N=32,64,$ and $128$ .

実験結果

リサーチクエスチョン

• RQ1 二次元落下膜方程式の領域サイズと分散パラメータを変えると、界面挙動はどのような乱れの域で現れるのか？
• RQ2 正確なコヒーレント状態は混沌とした落下膜力学に埋め込まれ、観測軌道を組織化するのに用いられるか？
• RQ3 縮約モデリングとデータ駆動的多様体は、界面カオスにおけるECSの同定と収束を助けるか？
• RQ4 混沌とした落下膜力学の慣性多様体の本質的次元は、領域サイズを超えてどのようになるのか？
• RQ5 縮約表現はDNSのスペクトルをどれだけ再現し、ECS計算をどれだけ効率化できるか？

主な発見

• (L, delta) のレジームマップは、移動波、暴走移動波、そして完全にカオスな界面力学を明らかにする。
• カオス領域では、アトラクターは有限次元の慣性多様体上にあり、その次元 d_M は領域サイズとともに概ね増加する。
• IRMAE-WD は明確なスペクトルギャップを持つ非線形多様体を同定し、支配的エネルギー含有構造を捉える compact な d_M 表現を可能にする。
• POD のみは多様体の次元を過大評価するのに対し、IRMAE-WD は低〜中波数スペクトルの再構成を正確に行い、 intrinsic な次元を低く保つ。
• 正確なコヒーレント状態（移動波、相対周期的軌道、平衡解）はカオス的アトラクター内に埋め込まれており、これらのECSの近傍に繰り返し近づく。
• 縮約系を前処理として全ニュートン-クライロフ精緻化を用いるワークフローは、DNSで機械精度のECSを得つつ計算コストを大幅に削減する。

Figure 2 : (a) Dependence of the Reynolds number as a function of the dispersion parameter in the range $\delta\in\left[10^{-3},5\right]$ for a falling film with $h_{0}$ =0.2 mm and $\theta=1.5608$ . The blue star and green point correspond to $\delta$ values of 0.002 and 0.65373, respectively, wher