Power spectra via the van der Waals effect in the two-dimensional Poiseuille and Couette flow

This paper demonstrates through numerical simulation that the complex dynamics and power decay observed in the Fourier spectra of two-dimensional inertial van der Waals flows near Poiseuille and Couette states are primarily driven by density and velocity divergence variables rather than vorticity.

The Gist

The Big Picture: A "Ghost" Turbulence

Imagine you are watching a river flow smoothly down a channel. To the naked eye, it looks calm and orderly—like a sheet of glass sliding over rocks. This is what scientists call laminar flow.

However, this paper discovers something strange happening under the surface. Even though the river looks calm, if you zoom in on the tiny, invisible ripples of air density and speed, they are going crazy. They are dancing chaotically, breaking down from big waves into tiny, frantic vibrations.

Usually, we think "chaos" and "turbulence" go hand-in-hand with a river turning into white water rapids. But here, the river stays smooth on the surface while the microscopic world underneath is having a wild party. The author calls this "pseudo-laminar" flow.

The Secret Ingredient: The "Van der Waals" Effect

Why is this happening? The paper introduces a specific physical rule called the Van der Waals effect.

Think of air molecules like a crowd of people in a hallway.

• Standard Physics: Usually, we assume these people don't really care about each other until they bump into someone.
• Van der Waals Physics: In this specific scenario, the people have a "personal space" rule. They push away from each other slightly before they even touch, and they pull together slightly when they are far apart.

The author's computer simulation shows that when you add this "personal space" rule to a smooth flow of air, it creates a hidden instability. It's like putting a tiny, invisible spring under a calm table; the table looks still, but the springs underneath are vibrating wildly.

The Experiment: Two Types of Smooth Flows

The author tested this idea on two classic types of smooth flows:

1. Poiseuille Flow: Imagine air flowing through a pipe. It moves fastest in the middle and stops at the walls (like a parabola shape).
2. Couette Flow: Imagine air between two plates where one plate is sliding and dragging the air with it (like a linear ramp).

In both cases, the author added a tiny "nudge" (a small disturbance) to the air density.

• Without the Van der Waals effect: The nudge just drifted along with the flow, like a leaf floating down a calm stream. Nothing exciting happened.
• With the Van der Waals effect: The nudge exploded into complex, chaotic patterns. The big ripples broke down into smaller and smaller ripples, creating a "direct cascade" of energy.

The Magic Trick: The Vorticity "Ghost"

Here is the most surprising part of the paper.

In fluid dynamics, there are two main ways air moves:

1. Squishing/Stretching (Divergence): Air getting compressed or expanding.
2. Spinning (Vorticity): Air swirling like a tornado.

Usually, scientists think the "spinning" (vorticity) is the main cause of turbulence. But the author did a magic trick: They froze the spinning.

They told the computer: "Keep the air spinning exactly as it was in the calm background. Don't let the spinning change at all. Only let the squishing and stretching change."

The Result: The chaos didn't stop! The air still developed the same wild, chaotic patterns and the same "power decay" (a specific mathematical rule about how energy fades out) as before.

The Analogy: Imagine a band playing a song. You might think the drummer (vorticity) is the one making the music complex. But the author turned the drummer's volume down to zero (frozen the spin), and the guitar and bass (density and divergence) kept playing the exact same complex, chaotic song.

Conclusion: The "chaos" isn't coming from the spinning; it's coming from the squishing and stretching of the air, driven by that Van der Waals "personal space" effect.

Why Does This Matter? (The "Power Spectra" Mystery)

Scientists have long observed that in turbulent flows, energy follows a specific mathematical rule (a "power law") as it breaks down from big waves to small waves. This is often called the Kolmogorov slope.

For a long time, the only explanation for this was a "best guess" based on dimensions (like saying, "It must be this way because of how big things are"). No one knew the actual physical machine that made it happen.

This paper suggests:

1. The Machine: The "Van der Waals effect" acting on air density and stretching is the machine that creates this power law.
2. The Location: You don't need a full-blown storm (turbulence) to see this. You can see it in tiny, almost invisible ripples on a calm day.
3. The Future: Because the "spinning" part isn't necessary for this chaos, the math becomes much simpler. Instead of a complicated 3-part puzzle, we might be able to solve a simpler 2-part puzzle to finally understand why the energy fades out the way it does.

Summary in One Sentence

This paper shows that even in a perfectly smooth, calm flow of air, tiny invisible ripples can go chaotic due to the "personal space" of air molecules, and surprisingly, this chaos is driven by the air being squished and stretched, not by it spinning.

Technical Summary

1. Problem Statement

The paper addresses the long-standing challenge of explaining the power-law decay of Fourier spectra observed in gas flows (both in atmospheric observations and laboratory experiments). While empirical dimensional hypotheses (like Kolmogorov's $k^{-5/3}$) have historically described these spectra, they lack a rigorous physical mechanism derived from first principles.

Specifically, the author investigates:

• Whether a two-dimensional (2D) flow, which is often considered insufficient for modeling real-world 3D turbulence, can generate the complex dynamics and power spectra associated with turbulence.
• The role of the Van der Waals effect (arising from pressure equilibration at low Mach numbers) in driving instabilities and spectral decay.
• The relative importance of vorticity versus density and velocity divergence in generating these spectra.
• The nature of the transition between laminar and turbulent states, specifically looking for "pseudo-laminar" regimes where macroscopic flow remains ordered while microscopic fluctuations exhibit chaotic spectral decay.

2. Methodology

The study employs numerical simulations of compressible inertial gas flow in a straight, periodic 2D channel.

• Governing Equations: The author uses a modified set of inertial flow equations where the pressure is constant (inertial flow), causing the Van der Waals effect (the second virial coefficient) to become the leading-order term in the momentum transport equation.
• Formulation: The equations are recast into divergence-vorticity variables ($\chi = \nabla \cdot \mathbf{u}$ and $\omega = \nabla^\perp \cdot \mathbf{u}$). This separates the compressible (divergence) and rotational (vorticity) effects.
  • The system includes a coupling term $-2\det(\nabla \mathbf{u})$ and the Van der Waals term $\frac{4p_0}{\rho_{HS}} \nabla \cdot (\frac{\nabla \rho}{\rho})$ in the divergence equation.
• Numerical Setup:
  • Solver: OpenFOAM with implicit Euler time discretization and second-order finite volume schemes (van Leer flux limiter).
  • Domain: A rectangular channel ($40 \text{ cm} \times 25 \text{ cm}$) with periodic boundary conditions in the streamwise direction and no-slip walls.
  • Background States: Simulations are run around two stationary states: Poiseuille flow (parabolic velocity profile) and Couette flow (linear velocity profile).
  • Initial Conditions: The density is perturbed by a small (1%), large-scale sinusoidal deviation, while divergence and vorticity start at their background values.
• Experimental Variations:
  1. Full System: Simulating the complete coupled system ($\rho, \chi, \omega$).
  2. Reduced System: "Pinning" the vorticity $\omega$ to its constant background state, leaving only density $\rho$ and divergence $\chi$ as dynamic variables.
  3. Control: Simulations without the Van der Waals effect ($p_0=0$) to verify code stability and decoupling.

3. Key Contributions

• 2D Turbulence Mechanism: The paper demonstrates that 2D flows can spontaneously generate chaotic dynamics and power-law spectral decay without breaking down into macroscopic turbulence, challenging the notion that 3D is strictly required for such phenomena.
• Decoupling Vorticity: A major theoretical contribution is the finding that vorticity fluctuations are not essential for the generation of power spectra. By fixing vorticity to its background state, the reduced system ($\rho, \chi$) reproduces the qualitative dynamics and spectral slopes of the full system.
• Identification of the Driver: The study identifies the Van der Waals effect coupled with velocity divergence as the primary physical mechanism driving the instability and spectral decay, rather than rotational dynamics.
• Pseudo-Laminar Regime: The paper defines and demonstrates a "pseudo-laminar" regime where tracer streaks remain intact (macroscopically laminar) while the Fourier spectra of fluctuations exhibit complex, turbulent-like power decay.

4. Key Results

A. Macroscopic vs. Microscopic Behavior

• Macroscopic: In all simulations, the flow remains macroscopically laminar. Tracer streaks seeded in the flow do not mix or break, remaining distinct even after 0.05 seconds, despite minor smudging.
• Microscopic: Small fluctuations in density and velocity divergence develop chaotic, wave-like patterns that transition from large scales to small scales (a direct cascade).

B. Power Spectra Decay Rates

The study computed time-averaged Fourier spectra for density ($\rho$), divergence ($\chi$), vorticity ($\omega$), and kinetic energy components. The decay rates ($k^{-\alpha}$) were found to be:

Variable	Poiseuille Flow Slope	Couette Flow Slope
Density ($\rho^2$)	$k^{-7/3}$	$k^{-7/3}$
Divergence ($\chi^2$)	$k^{-5/3}$	$k^{-4/3}$
Vorticity ($\omega^2$)	$k^{-8/3}$	$k^{-7/3}$
Kinetic Energy Components	Varied ($k^{-3}$ to $k^{-11/3}$)	Varied ($k^{-3}$ to $k^{-10/3}$)

• Observation: The decay rates are distinct for different variables and depend on the background flow profile (Poiseuille vs. Couette).
• Reduced System: When vorticity was fixed, the reduced system produced identical power slopes to the full system. This confirms that the spectral physics resides in the density-divergence coupling.

C. Comparison with Kolmogorov Theory

• The results show that different components of kinetic energy (which share the same physical units and origin) exhibit different power slopes. This contradicts the standard Kolmogorov dimensional hypothesis, which predicts a universal slope (typically $k^{-5/3}$) based solely on dimensional analysis.
• The paper argues that the background flow profile significantly influences the spectral slope, suggesting that dimensional hypotheses fail to capture the specific physical mechanisms (like the Van der Waals coupling) at play.

5. Significance and Implications

• Redefining Turbulence: The work suggests that turbulence is not a binary state (laminar vs. turbulent) but may consist of distinct parts. A flow can be macroscopically laminar yet possess the spectral signatures of turbulence in its small-scale fluctuations.
• Simplification of Analysis: By proving that vorticity can be pinned without altering the spectral dynamics, the complexity of linearizing the system for theoretical analysis is reduced from a $3 \times 3$ system to a $2 \times 2$ system. This makes the "inertial" regime (the intermediate stage between instability onset and asymptotic decay) more accessible to analytical solution.
• Critique of Dimensional Analysis: The paper challenges the reliance on Buckingham $\pi$ theorem and Kolmogorov's dimensional hypotheses, arguing that they fail to account for specific physical mechanisms (like the Van der Waals effect) that dictate the actual spectral slopes in real gas flows.
• Future Directions: The findings pave the way for explicit analytical solutions of the reduced 2D system, potentially providing the first rigorous physical derivation of power-law spectra in gas flows, moving beyond empirical observation.

In conclusion, Abramov demonstrates that the Van der Waals effect acting on density and velocity divergence is sufficient to generate complex, chaotic spectral decay in 2D flows, even when the flow remains macroscopically laminar and vorticity is static. This provides a new physical mechanism for power spectra that challenges traditional turbulence theories.