Analysis of a three-dimensional fluid flow in rotating cylinders

Technical Summary: Analysis of a Three-Dimensional Fluid Flow in Rotating Cylinders

Problem Statement
This paper investigates the mathematical modeling and analysis of a capillary-driven, three-dimensional thin-film rimming flow. The physical setting involves a Newtonian, incompressible fluid partially filling a horizontal cylinder of radius $R$ and length $L$, which rotates around its axis with a constant angular velocity $\omega$. The fluid forms a thin coating layer of mean thickness $d \ll R$ on the inner wall of the cylinder. The study focuses on the dynamics of the free surface height $h(t, \theta, z)$, accounting for the competing physical effects of rotation, surface tension, and gravity. The gravitational influence is treated as a small perturbation, measured by a non-dimensional parameter $0 \le \delta \ll 1$.

Methodology
The authors derive a closed evolution equation for the film height starting from the incompressible Navier–Stokes equations with a free boundary. By assuming a small aspect ratio ($\varepsilon = d/R \to 0$), they apply a formal lubrication approximation. This asymptotic limit, combined with specific scaling assumptions for the Weber, Capillary, and Bond numbers, yields a fourth-order quasilinear degenerate-parabolic partial differential equation (PDE) for the rescaled film height $h$:
$$h_t + h_\theta + \gamma \text{div}_\ell \left( h^3 \nabla_\ell (\Delta_\ell h + h) \right) - \delta (h^3 \cos \theta)_\theta = 0$$
defined on a torus $T^2$ (via a reflection principle to handle Neumann boundary conditions at the cylinder ends). Here, $\gamma$ represents the surface tension parameter, and the term $h_\theta$ arises from the rotation of the cylinder.

The analysis proceeds through several stages:

1. Local Well-posedness: Using analytic semigroup theory, the authors establish the local existence and uniqueness of positive maximal solutions for initial data in fractional Sobolev spaces $H^\alpha_{ev}(T^2)$ with $\alpha > 3$.
2. Steady State Analysis: The structure of stationary solutions is characterized based on the ratio $\ell = L/R$. The case $\delta = 0$ (no gravity) is analyzed first, followed by the case $0 < \delta \ll 1$.
3. Stability and Dynamics: Linearization and spectral analysis are employed to determine the stability of steady states. For the critical case $\ell = \pi$ with small gravity, the authors utilize the Poincaré–Lindstedt method of multiple time scales to separate fast dynamics (rotation) from slow dynamics (evolution of the fluid profile shape).

Key Contributions and Results

• Derivation and Well-posedness: The paper provides a rigorous derivation of the 3D rimming-flow equation from the Navier–Stokes system and proves local well-posedness for positive initial data. The solution conserves mass and, in the absence of gravity, dissipates energy.
• Steady States in the Absence of Gravity ($\delta = 0$):
  • If $\ell/\pi \notin \mathbb{Z}$, the only positive steady states are constant profiles $H \equiv m$ (cylindrical films).
  • If $\ell/\pi \in \mathbb{Z}$, a one-dimensional manifold of steady states exists, given by $H(\zeta) = m + a \cos(\frac{\ell}{\pi}\zeta)$, representing deformed cylinders where the cross-section radius varies along the axis.
• Stability Analysis ($\delta = 0$):
  • For short cylinders ($\ell < \pi$), constant steady states are orbitally exponentially stable. Solutions starting near a constant converge exponentially to a traveling wave profile (a shifted circle rotating with the cylinder).
  • For long cylinders ($\ell > \pi$), constant steady states are unstable, with the spectrum containing unstable eigenvalues corresponding to "bubble" formation along the axis.
• Effect of Small Gravity ($0 < \delta \ll 1$):
  • Uniqueness: In the case $\ell/\pi \in \mathbb{Z}$, where steady states were non-unique for $\delta=0$, the introduction of small gravity selects a unique steady state. This is proven using a Lyapunov–Schmidt reduction and a third-order Taylor expansion, contradicting previous suggestions that non-uniqueness might persist.
  • Stability: For $\ell < \pi$, the unique positive steady states $H_\delta$ are exponentially stable. For $\ell > \pi$, instability persists for small $\delta$.
• Critical Case Dynamics ($\ell = \pi$):
  • In the critical case with small gravity, solutions converge exponentially fast to a $\delta$-neighborhood of a three-dimensional manifold $M(m)$ of traveling wave profiles.
  • On a slow time scale $\tau = \delta^2 t$, the dynamics within this neighborhood are approximated by a system of ordinary differential equations (ODEs) governing the coefficients of the manifold.
  • Linearization of this ODE system reveals that the fluid profile's center spirals inward toward the cylinder axis while the cross-section remains circular. The authors provide numerical simulations of this ODE system, showing the spiraling behavior and the interaction between axial and radial perturbations.

Significance and Claims
The paper claims to revisit and extend previous works (specifically [Puk05] and [JLV26]) by providing a comprehensive analysis of the three-dimensional case. Its primary significance lies in:

1. Establishing the local well-posedness of the 3D rimming-flow equation.
2. Characterizing the bifurcation of steady states with respect to the cylinder length ratio $\ell$ and the gravitational parameter $\delta$.
3. Demonstrating that gravity acts as a selection mechanism, breaking the non-uniqueness of steady states in the critical resonant case ($\ell/\pi \in \mathbb{Z}$) to yield a unique solution.
4. Deriving a reduced ODE system that captures the long-time, slow-time-scale dynamics of the fluid film in the critical case $\ell = \pi$, bridging the gap between the full PDE and the dynamics of the fluid profile's shape.

The authors note that while the derivation of the ODE system for the critical case is formal due to challenges with Sobolev embeddings in the quasilinear setting, the resulting system provides a consistent description of the dynamics, supported by numerical simulations. The work does not propose new industrial applications but rather provides a rigorous mathematical foundation for understanding the stability and evolution of thin films in rotating geometries under gravity.