Evidence of an inertialess Kapitza instability due to viscosity stratification

This paper demonstrates that continuous viscosity stratification in gravity-driven falling films can induce a surface-mode instability in the zero-inertia (Stokes) limit, driven by a phase shift between perturbation vorticity and interface displacement that reinforces deformation, thereby extending the class of scalar-mediated, inertialess instabilities beyond traditional surfactant-driven Marangoni effects.

The Gist

Imagine a thin sheet of liquid, like honey or paint, flowing down a slanted roof. Usually, if the liquid is thick and heavy (viscous) and moving slowly, it flows smoothly like a calm river. In physics, we call this "Stokes flow," where the liquid is so sluggish that its own momentum (inertia) doesn't matter.

For decades, scientists believed that for waves to form and grow on this flowing sheet (an instability), the liquid needed to be moving fast enough to have "inertia." Without that speed, the flow was thought to be perfectly stable, like a calm pond.

The Big Discovery
This paper reveals a surprising new rule: You don't need speed to create waves. You just need a "gradient."

The authors found that if the liquid's thickness (viscosity) changes gradually from the bottom of the sheet to the top, the flow can become unstable and start wiggling, even if it's moving as slowly as a snail. They call this an "inertialess Kapitza instability."

The Analogy: The "Slippery Slope" and the "Delayed Reaction"

To understand how this happens without speed, let's use a metaphor involving a conga line and a sticky floor.

1. The Setup: The Viscosity Gradient
Imagine a line of people (the fluid) walking down a hallway.

• At the bottom (near the wall), the floor is super sticky (high viscosity).
• At the top (near the ceiling), the floor is slippery (low viscosity).
• Because of this, the people at the top can walk faster than the people at the bottom. This is the "viscosity stratification."

2. The Trigger: A Little Nudge
Now, imagine someone in the middle of the line trips slightly, creating a small bump (a wave) in the line.

• In a normal, uniform fluid, the people would just smooth this out.
• But in our "sticky-to-slippery" fluid, something weird happens.

3. The Mechanism: The "Delayed Reaction" (Phase Lag)
Here is the magic trick that causes the wave to grow:

• The Nudge: The bump pushes a bit of "sticky" fluid from the bottom up toward the "slippery" top.
• The Delay: Because the fluid at the top is moving faster, it drags that sticky patch downstream (forward) before it can settle.
• The Mismatch: By the time the sticky patch arrives at the top, the original bump has moved forward. The "sticky patch" is now sitting behind the bump, or slightly offset from where it should be.
• The Feedback Loop: This misplaced sticky patch creates a weird twist in the flow (vorticity). Because it's in the wrong place (lagging behind), it actually pushes the bump even higher instead of smoothing it out.

It's like a child on a swing. If you push the swing when it's coming toward you, you stop it. But if you push it just a split second after it passes the bottom (a "lag"), you make it go higher. The viscosity gradient creates this perfect "lag" in the fluid's reaction, turning a tiny ripple into a growing wave.

Why Does It Only Happen Sometimes? (The "Goldilocks" Zone)

The paper also found that this instability only happens in a specific "Goldilocks" zone, depending on how fast the "sticky" property spreads (diffusion) versus how fast the fluid moves (advection).

• Too Slow (Too much diffusion): If the "stickiness" spreads out too fast, it smears out the pattern. The fluid can't hold the "lag" long enough to push the wave. The system stays calm.
• Too Fast (Too much advection): If the fluid moves so fast that the "sticky" patches get frozen in place and just get carried away, they never have time to interact with the wave to push it. The system stays calm.
• Just Right (The Window): There is a perfect middle ground where the "sticky" patches move just enough to stay connected to the wave but lag just enough to push it. This is the "finite Péclet number window" mentioned in the paper.

The "Surfactant" Connection

The authors compare this to a phenomenon called Marangoni instability, which happens when soap (surfactant) is added to water. Soap changes the surface tension, creating waves.

• In the soap case, the "soap" is the scalar that moves and causes the instability.
• In this new discovery, the viscosity itself is the "soap." It's not just a property of the fluid; it's a moving variable that creates its own waves.

Why Does This Matter?

This isn't just about honey on a roof. This discovery helps us understand:

1. Industrial Coatings: When painting or coating surfaces with materials that have varying thickness (like thermal sprays), we need to know when they might start rippling uncontrollably.
2. Particle Flows: In fluids filled with particles (like blood or mud), particles often migrate to the surface, changing the viscosity. This paper explains why those flows might become wavy even when they seem too slow to do so.
3. Geology: It might help explain how thick layers of magma or ice flow and deform deep within the Earth or on glaciers, where inertia is negligible but viscosity changes are huge.

In a Nutshell:
The paper proves that you don't need a fast-moving river to create waves. If you have a fluid where the "thickness" changes from bottom to top, the fluid can trick itself into creating waves through a clever timing mismatch, even when it's barely moving at all. It's a new kind of instability driven entirely by the fluid's own internal structure.

Technical Summary

1. Problem Statement

The classical Kapitza instability (surface waves in gravity-driven falling films) is traditionally understood to require finite fluid inertia (Reynolds number $Re > 0$) to operate. In the Stokes flow limit (zero inertia, $Re \to 0$), homogeneous falling films are linearly stable.

However, many practical flows involve viscosity stratification (e.g., particle-laden films with shear-induced migration, thermally stratified coatings, or concentration-graded flows). While viscosity stratification is known to cause interfacial instabilities in two-layer flows (discontinuous viscosity), its effect on the stability of a continuous viscosity profile in a single-phase falling film in the absence of inertia remains an open question. Specifically, it is unclear if viscosity stratification alone, decoupled from complex suspension dynamics or normal stresses, can trigger a surface-mode instability in the Stokes limit.

2. Methodology

The authors employ a linear stability analysis of a two-dimensional, gravity-driven falling film on an inclined plane.

• Governing Equations: The flow is governed by the incompressible Stokes equations (neglecting inertia) coupled with an advection-diffusion equation for the viscosity field $\mu$.
  • Momentum: $\nabla \cdot \mathbf{T} + \rho \mathbf{g} = 0$
  • Viscosity Transport: $\frac{\partial \mu}{\partial t} + \mathbf{u} \cdot \nabla \mu = D_s \nabla^2 \mu$
• Base State:
  • The velocity profile is fixed as the classical Nusselt parabola ($u_b(y)$), deliberately decoupling the base flow from viscosity variations to isolate the effect of perturbation coupling.
  • The viscosity is prescribed as a linear profile across the film thickness: $\mu_b(y) = 1 + \alpha(y - 0.5)$, where $\alpha$ is the stratification parameter.
• Non-dimensionalization: The system is characterized by the Capillary number ($Ca$) and the Péclet number ($Pe = u_0 h_0 / D_s$), which represents the ratio of advective to diffusive transport of viscosity.
• Analysis Techniques:
  1. Long-wave Asymptotics: An expansion in small wavenumber $k$ ($k \ll 1$) is performed to derive a dispersion relation and identify the instability mechanism analytically.
  2. Chebyshev Spectral Collocation: The full coupled eigenvalue problem for the perturbation streamfunction ($\hat{\psi}$) and viscosity ($\hat{\mu}$) is solved numerically to validate asymptotic results and explore finite wavenumbers.

3. Key Contributions

• Discovery of Inertialess Instability: The paper provides the first evidence that a surface-mode instability can exist in a falling film with zero inertia solely due to continuous viscosity stratification.
• Mechanism Elucidation: The authors identify the instability mechanism as a vorticity-interface phase mismatch driven by the advection-diffusion of viscosity, structurally analogous to surfactant-driven Marangoni instability.
• Finite Pe Window: The study demonstrates that this instability is not monotonic with transport parameters; it exists only within a finite window of Péclet numbers.
• Decoupling of Physics: By fixing the base velocity profile, the authors prove that viscosity stratification is a sufficient condition for instability, independent of the complex couplings (e.g., normal stresses, modified base flow) found in full suspension models.

4. Key Results

A. Stability Characteristics

• Instability Existence: For stratified films ($\alpha > 0$), the surface mode becomes unstable even at $Re=0$.
• Dependence on Parameters:
  • Stratification ($\alpha$): Increasing $\alpha$ lowers the critical Péclet number for onset, broadens the range of unstable wavenumbers, and increases the growth rate.
  • Péclet Number ($Pe$): The instability is confined to a specific window:
    • Low $Pe$ (Diffusion-dominated): Diffusion smooths out viscosity perturbations before they can interact with the shear, stabilizing the flow.
    • Intermediate $Pe$: Optimal balance between advection and diffusion allows viscosity perturbations to persist and feedback into the momentum equation, maximizing growth rates.
    • High $Pe$ (Advection-dominated): Advection dominates, causing viscosity perturbations to become "frozen" and in-phase with the streamfunction. This eliminates the necessary phase lag, restabilizing the flow.
• Neutral Curves: Neutral stability curves in the $(\alpha, Pe)$ and $(k, Pe)$ planes form closed lobes or wedges, confirming the finite window nature of the instability.

B. Instability Mechanism (Hinch's Framework)

The authors trace the causal chain using the vorticity-interface phase framework:

1. Step 1 (Generation of Phase Shift): A sinusoidal interface displacement induces a streamfunction perturbation ($\hat{\psi}_0$). This flow advects the base-state viscosity gradient ($\mu'_b$), generating a viscosity perturbation ($\hat{\mu}_1$). Due to the advection-diffusion balance, $\hat{\mu}_1$ is $90^\circ$ out of phase (purely imaginary relative to $\hat{\psi}_0$) with the interface.
2. Step 2 (Vorticity Feedback): This phase-shifted viscosity perturbation acts as a source term in the momentum equation. It generates a vorticity perturbation that is lagging behind the interface displacement.
3. Result: A lagging vorticity distribution drives fluid upward on the upstream side of wave crests and downward upstream of troughs, reinforcing the interface deformation and causing exponential growth.

C. Structural Analogy

The mechanism is shown to be structurally identical to surfactant-driven Marangoni instability in creeping two-layer flows. In both cases, a transported scalar (viscosity or surfactant concentration) creates a phase-shifted stress that drives the instability, provided the scalar transport is neither too diffusive nor too advective.

5. Significance

• Theoretical Advancement: This work challenges the long-held view that inertia is a prerequisite for Kapitza instability in falling films. It expands the class of known inertialess instabilities to include bulk viscosity stratification.
• Industrial Relevance: The findings are critical for understanding the stability of:
  • Particle-laden films: Where shear-induced migration creates viscosity gradients.
  • Thermally stratified coatings: Where temperature gradients alter viscosity.
  • Concentration-graded flows: Common in polymer processing and chemical engineering.
• Predictive Capability: The identification of a finite $Pe$ window suggests that controlling the transport properties (diffusivity vs. flow speed) can be used to either suppress or enhance wave formation in low-inertia coating processes.

In conclusion, the paper establishes that viscosity stratification alone is a sufficient destabilizing agent for gravity-driven falling films in the Stokes limit, operating through a scalar-mediated phase-shift mechanism that requires an intermediate balance between advection and diffusion.