[논문 리뷰] Exact coherent states underlying chaotic falling-film dynamics

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

Fluid Dynamics and Thin Films인용 수 0

한 줄 요약

이 논문은 낙하 필름에 대한 이차원 장파(long-wave) 방정식을 도출하고, 그 거동 영역을 매핑하며, DNS를 저차원 관성 매니폴드에 투영하고 DNS와 축소 모델에서 Newton-Krylov 방법을 사용해 혼돈의 계면 역학 속에 내재된 정확한 일관 상태를 식별한다.

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, 데이터 기반 매니폴드 학습, Newton-Krylov 보정을 결합한 워크플로우를 사용해 불변 해를 찾는 방법을 제시한다.

제안 방법

• 필름 두께 H(x,y,t)에 대한 이차원 장파 방정식을 분산 매개변수 delta와 함께 도출하고 사용한다.
• 스펙트럴 DNS 해석기으로 이중 주기 도메인에서 방정식을 해결하고 거동을 전력 스펙트럼과 에너지 측정으로 분석한다.
• 고차원 DNS 데이터를 POD로 축소한 뒤 암시적 랭크 최소화 오토인코더(IRMAE-WD)로 관성 매니폴드를 차원 d_M으로 식별한다.
• 매니폴드 위의 역학에 대한 뉴럴 ODE를 학습해 잠재 좌표에서의 시간 진행을 근사적으로 얻는다.
• 대칭 축소 및 Jacobian-free Newton-Krylov(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)에서의 거동 지도는 이동 파동, 버스트 이동 파동, 그리고 완전히 혼돈스러운 계면 역학을 나타낸다.
• 혼돈 영역에서 attractor는 유한 차원 관성 매니폴드 위에 놓이며 그 차원 d_M은 도메인 크기에 따라 대략 증가한다.
• IRMAE-WD는 명확한 스펙트럴 간극을 가지는 비선형 매니폴드를 식별하여 지배적 에너지 보유 구조를 포착하는 compact한 d_M 표현을 가능하게 한다.
• POD만으로는 매니폴드 차원을 과대평가하는 반면 IRMAE-WD는 더 낮고 내재적 차원을 보여주며 저위-중 주파수 스펙트럼의 재구성도 정확하다.
• 이동 파동, 상대 주기궤도, 평형 상태 등 정확한 일관 상태가 혼돈 attractor 내부에 내재되어 있으며 이 ECS 근처로 궤도가 반복적으로 접근한다.
• 축소 차원 동역학을 전체 Newton-Krylov 보정의 프리컨디셔너로 활용한 워크플로우는 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