Flow instability in Stokes layer of Carreau fluids

This study investigates the instability of Stokes layers in shear-thinning Carreau fluids, revealing that stronger shear-thinning monotonically stabilizes the flow while the fluid's response time has a non-monotonic effect, with instability driven by the phase alignment between perturbations and the oscillatory base flow that enables efficient energy extraction.

The Gist

The Big Picture: Wiggling Fluids

Imagine you have a thick, gooey substance (like honey or ketchup) trapped between two flat plates. Now, imagine you shake these plates back and forth very quickly. This creates a "Stokes layer"—a thin layer of fluid near the plates that wiggles along with them, while the fluid in the middle stays relatively calm.

The researchers wanted to know: If you shake this gooey fluid, will it stay smooth, or will it suddenly become chaotic and turbulent?

Most fluids we know (like water) are "Newtonian," meaning their thickness doesn't change no matter how fast you stir them. But many real-world fluids (like blood, paint, or shampoo) are shear-thinning. This means they get thinner and runnier the faster you move them. The paper investigates how this "getting thinner when shaken" behavior changes the stability of the wiggling fluid.

The Tools: Two Ways to Look at the Fluid

To solve this, the team used two different mathematical "lenses":

1. The Super-Computer Lens (Numerical Method): They used a powerful computer to simulate every tiny detail of the fluid's movement. This is accurate but very slow and difficult, especially when the fluid gets very runny.
2. The "Small Shake" Lens (Expansion Method): They developed a clever math trick. They assumed the fluid's "runniness" change was small and used a series expansion (like adding up terms in a recipe) to predict the flow.
   • The Result: This math trick works perfectly when the fluid isn't changing its thickness too drastically. It's much faster than the computer simulation and gives them a clear formula to understand the physics. However, if the fluid changes thickness too wildly, the math trick breaks down, and they have to rely on the slow computer method.

The Findings: The Goldilocks Zone of Stability

The researchers tested two main knobs on their fluid model:

• Knob A (How much it thins): How drastically the fluid gets runnier when shaken (represented by the power-law index, n).
• Knob B (How fast it reacts): How quickly the fluid's thickness changes in response to the shaking (represented by the time scale, Λ).

Here is what they discovered:

1. The "More Runny" Knob (Decreasing n):
If you make the fluid more shear-thinning (it gets much thinner when shaken), the flow becomes more stable. It's harder to make it go chaotic.

• Analogy: Think of a crowd of people trying to run in place. If everyone is stiff and heavy, they might trip over each other easily. But if everyone is light and fluid, they can move in sync without tripping. Making the fluid "lighter" (more shear-thinning) actually helps it stay organized.

2. The "Reaction Speed" Knob (Increasing Λ):
This is where it gets surprising. The effect of how fast the fluid reacts is not a straight line.

• Slow Reaction: If the fluid reacts slowly to the shaking, it stays stable.
• Medium Reaction: As the reaction speed increases to a medium level, the fluid becomes even more stable. It's like a dancer finding the perfect rhythm.
• Fast Reaction: But if the reaction speed gets too fast (strong shear-thinning), the fluid suddenly becomes unstable and prone to chaos.
• Analogy: Imagine trying to balance a broom on your hand.
  • If you move your hand very slowly, the broom stays up.
  • If you move it at a moderate, rhythmic pace, you can balance it very well.
  • But if you jerk your hand back and forth too frantically, the broom falls over. The fluid behaves similarly: too much "frantic" thinning makes it lose its balance.

The Secret Mechanism: The Dance of Energy

Why does this happen? The team performed an "energy analysis" to see where the chaos comes from.

They found that for the fluid to become unstable, the tiny ripples (disturbances) in the fluid must sync up perfectly with the shaking of the walls to steal energy from them.

• The Stable Phase: When the fluid reacts at a medium speed, the ripples are slightly out of step with the wall's movement. It's like trying to push a swing when the swing is moving away from you; you can't transfer much energy, so the swing (the flow) stays calm.
• The Unstable Phase: When the fluid reacts very quickly (strong shear-thinning), the ripples get back in perfect step with the wall. Now, every time the wall pushes, the ripples push back at the exact right moment, stealing maximum energy. This energy buildup causes the flow to break down into turbulence.

Summary

The paper shows that shear-thinning fluids don't just get "thinner"; they change how they react to being shaken in a complex way.

• Making a fluid more shear-thinning generally helps it stay smooth.
• However, if the fluid's ability to thin out happens too quickly relative to the shaking speed, it can actually trigger chaos.
• The key to stability is timing: if the fluid's internal changes are out of sync with the external shaking, the flow remains calm. If they sync up, the flow explodes into turbulence.

This research helps us understand the fundamental rules of how complex fluids behave when they are oscillated, which is crucial for anything from industrial mixing to understanding blood flow, though the paper itself focuses strictly on the physics of the instability mechanism.

Technical Summary

Technical Summary: Flow Instability in Stokes Layer of Carreau Fluids

Problem Statement
This study investigates the linear stability of the Stokes layer—a prototype time-periodic flow occurring over an oscillating flat plate—when the fluid exhibits shear-thinning behavior. While the stability of Newtonian Stokes layers is well-documented, the influence of non-Newtonian rheology, specifically shear-thinning, on time-periodic flow instability remains largely unexplored. The authors focus on fluids modeled by the Carreau constitutive equation, which captures the transition from zero-shear-rate viscosity to infinite-shear-rate viscosity. A primary challenge addressed is the absence of an analytical solution for the laminar base flow in shear-thinning fluids, necessitating the development of robust numerical and semi-analytical methods to characterize the time-dependent viscosity and velocity profiles before performing stability analysis.

Methodology
The research employs a Floquet analysis to examine the linear stability of the time-periodic base flow. The methodology is structured around three core components:

1. Base Flow Calculation: Due to the nonlinearity of the Carreau model, an analytical solution for the base flow is unavailable. The authors utilize two distinct methods:

   • Numerical Method: A spectral collocation method using Chebyshev polynomials in space and trigonometric polynomials in time is used to solve the full nonlinear governing equations. This approach handles arbitrary values of the Carreau number ($\Lambda$) and power-law index ($n$).
   • Expansion Method: For small $\Lambda$ (representing weak shear-thinning), a binomial expansion of the viscosity term is performed. The solution is expanded in powers of $\Lambda$ up to the sixth order ($\Lambda^6$), yielding semi-analytical solutions for the base flow harmonics. This method is validated against the numerical solution and found to be accurate for $\Lambda \lesssim 0.2$.

2. Linear Stability Analysis: The linearized perturbation equations are derived using Reynolds decomposition. The perturbation field is assumed to follow a Floquet-Fourier-Hill (FFH) form, where the solution is a product of a periodic shape function and an exponential growth/decay term. The resulting eigenvalue problem is solved numerically to determine the linear growth rate ($\omega_i$) and frequency ($\omega_r$) of the most unstable modes.

3. Energy Analysis: To elucidate the physical mechanisms driving instability, an energy budget analysis is conducted. This involves integrating the perturbation kinetic energy equation over the channel height and time to quantify contributions from energy production, viscous dissipation, pressure work, and additional terms arising from viscosity stratification and the difference between tangent and bulk viscosity.

Key Contributions and Results

• Validation of Methods: The study establishes that the binomial expansion method accurately reproduces the numerical base flow and linear stability results for weak shear-thinning regimes ($\Lambda \leq 0.2$), providing a computationally efficient alternative for theoretical analysis in this limit.
• Non-Monotonic Effect of $\Lambda$: The investigation reveals a non-monotonic relationship between the Carreau number ($\Lambda$, the ratio of viscosity relaxation time to wall-oscillation penetration time) and flow stability.
  • Increasing $\Lambda$ from zero initially stabilizes the flow, raising the critical Reynolds number ($Re_c$).
  • Beyond an intermediate threshold (approximately $\Lambda \approx 3$), further increasing $\Lambda$ destabilizes the flow, lowering $Re_c$ below the Newtonian value. This destabilization is observed in both weakly and strongly shear-thinning regimes.
• Monotonic Stabilizing Effect of $n$: In contrast to $\Lambda$, decreasing the power-law index $n$ (which signifies stronger shear-thinning) exerts a monotonic stabilizing effect on the flow within the investigated parameter space.
• Frequency Dependence: The study demonstrates that varying the dimensional oscillation frequency ($\Omega$) also produces a non-monotonic effect on stability. The flow is stable at very low and very high frequencies but becomes unstable within an intermediate frequency range.
• Instability Mechanism: The energy analysis identifies the phase relationship between the perturbation field and the oscillatory base shear as the dominant mechanism.
  • Stabilization: Occurs when there is a significant phase mismatch between the disturbance and the base flow shear, suppressing efficient energy extraction.
  • Destabilization: Arises when the perturbation field aligns (is in phase) with the time-dependent shear, allowing for efficient energy extraction from the base flow. This mechanism parallels the energy production process in steady shear flows but is identified here for the first time in a time-periodic shear flow context.

Significance and Claims
The authors claim that this work provides the first Floquet analysis of the Stokes layer in Carreau fluids, filling a gap in the literature regarding shear-thinning effects on time-periodic flows. The study highlights that shear-thinning does not have a uniform effect on stability; rather, it can either delay or promote transition to turbulence depending on the specific rheological parameters ($\Lambda$ and $n$).

The identification of the phase-matching mechanism as the driver for instability in time-periodic shear-thinning flows offers a new physical perspective distinct from steady flow analyses. The authors note that these findings have potential implications for applications involving oscillatory flows of complex fluids, such as in microfluidic mixers or biological systems (e.g., blood flow in arteries), where controlling flow instability could enhance mixing or prevent turbulence. However, the paper modestly frames these as potential implications, noting that nonlinear instabilities and experimental verification are required for definitive conclusions. The work also serves as a reference for future studies on viscoelastic fluids, as the current model isolates shear-thinning effects without elastic contributions.