Distinct transverse-response signatures of retained-spin, eliminated-spin, and polynomial Burnett-type surrogate closures

This paper demonstrates that transverse linear response analysis can dynamically distinguish between retained-spin micropolar dynamics, eliminated-spin effective theories, and polynomial Burnett-type closures by revealing unique spectral signatures, such as phase lags and pole structures, that are validated through many-particle simulations of rough spheres.

The Gist

Imagine you are trying to figure out how a crowd of people moves through a busy hallway. Usually, we assume everyone just walks forward, and if they bump into each other, they just slow down a bit. This is like the standard Navier-Stokes equations used in fluid dynamics: simple, smooth, and predictable.

But what if the people in the crowd aren't just walking? What if they are also spinning like tops? And what if those spins interact with the crowd's movement in complex ways?

This paper is a detective story about how to tell the difference between three different "stories" of how that spinning crowd behaves, especially when the movement gets very fast or very chaotic (high curvature).

Here is the breakdown of the three suspects and how the author, Satori Tsuzuki, caught them in the act.

The Three Suspects (The Three Theories)

Imagine you see a strange ripple in the crowd. You want to know why it happened. There are three possible explanations:

1. The "Spinning Top" Theory (Retained-Spin):

   • The Story: The people are actually spinning, and that spin is a real, independent thing. They have their own energy and momentum. The spin takes a tiny bit of time to react to the crowd's movement.
   • The Signature: If you watch closely, you'll see two distinct "beats" or rhythms in the movement: the main walking rhythm and a fast, jittery spinning rhythm.

2. The "Instant Adjustment" Theory (Eliminated-Spin):

   • The Story: The people are spinning, but they are so good at it that they adjust instantly. They don't have their own independent rhythm; they just instantly copy the crowd's movement.
   • The Signature: You only see one rhythm (the walking one). However, the way the crowd reacts to fast changes isn't a simple curve; it's a specific, complex shape (a "rational kernel") that looks like a smooth hill.

3. The "Approximation" Theory (Polynomial Burnette):

   • The Story: We don't care about the spinning at all. We just try to guess the crowd's behavior by adding simple math terms (like $x^2$, $x^4$, $x^6$) to our equations to make them fit better.
   • The Signature: This looks similar to the "Instant Adjustment" theory at slow speeds, but it starts to break down and behave wildly when things get fast.

The Detective's Toolkit: The "Transverse Response"

How do you tell these stories apart? You can't just look at a snapshot; you have to poke the system and see how it reacts.

The author uses a technique called Transverse Linear Response. Imagine you are pushing the crowd sideways with a rhythmic wave (like a stadium wave). You measure how the crowd wiggles back.

• The Pole Count (The Heartbeat):

  • If you see two heartbeats (poles) in the reaction, it's the "Spinning Top" Theory. The crowd has a hidden, fast internal rhythm.
  • If you see only one heartbeat, it's either the "Instant Adjustment" or the "Approximation" theory. You need a closer look.

• The Shape of the Curve (The Fingerprint):

  • If the reaction curve is a smooth, rational hill (it levels off nicely), it's the "Instant Adjustment" Theory.
  • If the reaction curve is a polynomial (a simple math equation), it's the "Approximation" Theory.
  • The Trap: The "Approximation" theory is a liar. If you push it too hard (high speed/high frequency), it either gets too stiff (over-damped) or it explodes (becomes unstable). The "Instant Adjustment" theory stays stable no matter how hard you push.

The Experiment: The "Rough Sphere" Simulation

To prove this isn't just math on paper, the author ran a massive computer simulation. Imagine a box filled with 8,000 tiny, perfectly rough billiard balls. "Rough" means they don't just bounce; they grip and spin when they hit each other.

1. The Free Fall Test: They let the balls move and then stopped pushing them. The balls slowed down in a way that perfectly matched the "one heartbeat" (Instant Adjustment) model. The fast spinning was too quick to see in the slow motion, so it looked like a single rhythm.
2. The Rhythmic Push Test: They started shaking the box with a specific rhythm.
   • They measured the spin of the balls versus the swirl of the crowd.
   • The Smoking Gun: The spin didn't react instantly. It lagged behind the swirl by a tiny, measurable fraction of a second.
   • The Verdict: This lag proved that the spin was a real, independent thing (The "Spinning Top" Theory). The "Instant Adjustment" theory (which says spin happens instantly) was proven wrong by the data.

The Big Picture: Why Does This Matter?

This paper is like a manual for scientists who study complex fluids (like blood, liquid crystals, or turbulent air).

• The Problem: Sometimes, complex fluids act weirdly. Scientists often try to fix their math by just adding more terms (like $x^4$ or $x^6$) to the equation. This is the "Approximation" theory.
• The Danger: The author shows that these simple fixes are dangerous. They might work for slow speeds, but they will fail catastrophically at high speeds, predicting things that can't happen (like infinite energy or instability).
• The Solution: Instead of guessing with simple math, we should look for the signatures of the hidden physics (like the spin lag). If we see a lag, we know we need a more complex model that treats the spin as a real, moving part of the system.

In short: You can't just "smooth over" the complexity of spinning particles with simple math. If you poke the system hard enough, the hidden spin will reveal itself, and you'll know exactly which story is true.

Technical Summary

1. Problem Statement

In the continuum modeling of complex fluids (e.g., suspensions of rough spheres or polyatomic gases), high-curvature observables (such as $k^4$-weighted energy spectra) can arise from three distinct physical mechanisms that are often indistinguishable at the level of low-order derivative expansions:

1. Explicit Retained-Spin Dynamics: The internal micro-rotation (spin) is treated as an independent, dynamical field coupled to the macroscopic vorticity.
2. Adiabatic Elimination: The fast spin variable is eliminated from the equations, resulting in an effective one-field theory with a rational $k$-dependent kernel.
3. Polynomial Burnett-Type Closures: A conventional higher-gradient closure where the internal field is never introduced, resulting in a polynomial expansion in wavenumber $k$ (e.g., $k^4$, $k^6$ terms).

While the low-wavenumber ($k \to 0$) expansions of mechanisms (2) and (3) appear similar (both generating $k^4$ and $k^6$ terms), their high-wavenumber behavior and pole structures differ fundamentally. The paper addresses the challenge of distinguishing these mechanisms using transverse linear response functions and microscopic simulations.

2. Methodology

The study employs a three-pronged approach combining analytical theory, reduced-model benchmarks, and microscopic simulations:

• Theoretical Framework:

  • The authors utilize a retained-spin micropolar closure derived previously from the Boltzmann–Curtiss equation (companion paper [6]).
  • They derive the fast-spin elimination limit, showing it yields a one-field equation with a rational kernel (non-polynomial in $k$).
  • They construct a polynomial Burnett-type surrogate (Model B) to represent a standard finite-gradient truncation.
  • Four models are compared in the transverse sector:
    • Model A: Standard Navier-Stokes ($k^2$).
    • Model B: Polynomial Burnett surrogate (finite $k^4$ or $k^4+k^6$).
    • Model C: Explicit retained-spin (two-field, two poles).
    • Model D: Eliminated-spin (one-field, rational kernel).

• Reduced-Model Benchmarks (Section IV):

  • The authors perform controlled numerical experiments on single-mode Fourier amplitudes ($\zeta(t; k)$).
  • They analyze dispersion relations (growth/decay rates), free-decay dynamics, and harmonic frequency responses (Bode plots).
  • They define diagnostic metrics (relative error, phase difference, amplitude ratios) to compare models across a $(k, \Omega)$ grid.

• Microscopic Observability (Section V):

  • Event-Driven Molecular Dynamics (EDMD) simulations of perfectly rough hard spheres ($N=8192$) are performed.
  • Two forcing protocols are used:
    1. Free Decay: To isolate late-time hydrodynamic behavior.
    2. Harmonic Forcing: To measure coherent responses in velocity ($u_x$), vorticity ($\zeta_z$), and spin ($\omega_z$).
  • Model Selection: Information Criteria (AIC/BIC) are used to statistically discriminate between the theoretical models based on the noisy microscopic data.

3. Key Contributions

A. Structural Distinction: Rational vs. Polynomial Kernels

The paper establishes that eliminating a fast internal variable produces a rational kernel in wavenumber space, whereas a Burnett truncation produces a polynomial.

• Rational Kernel (Model D): As $k \to \infty$, the damping scales effectively as $k^2$ (with a constant shift), preserving stability.
• Polynomial Kernel (Model B):
  • Strict $k^4$ truncation: Remains stable but becomes over-damped at large $k$ (decay rate grows too fast).
  • Matched $k^6$ truncation: To match the low-$k$ expansion of the rational kernel, a negative $k^6$ coefficient is required. This leads to finite-$k$ instability (the damping coefficient becomes negative) beyond a critical wavenumber $k_{crit}$.
• Conclusion: No finite polynomial truncation can simultaneously reproduce the low-$k$ expansion and the large-$k$ asymptotic structure of the eliminated-spin theory.

B. Pole Counting and Transients

• Model C (Explicit Spin): Possesses two poles (a slow hydrodynamic mode and a fast spin-relaxation mode).
• Models A, B, D: Possess only one pole.
• In the fast-spin regime, the fast pole in Model C manifests as a short transient in free decay and a negligible correction to the steady-state response, making Model C and Model D appear similar in coarse measurements. However, the presence of the second pole is the definitive signature of retained spin.

C. Microscopic Discrimination via Spin-Vorticity Lag

The paper demonstrates that phase-locked spin response is a critical discriminator:

• Adiabatic Elimination (Model D): Predicts a frequency-independent, real spin-to-vorticity ratio.
• Retained Spin (Model C): Predicts a complex ratio with a frequency-dependent phase lag ($R \propto 1/(1+i\Omega\tau)$).
• EDMD simulations successfully resolve this finite phase lag, strongly favoring the retained-spin model over instantaneous elimination.

4. Key Results

• Stability Analysis: The matched polynomial surrogate (Model B with $k^6$) exhibits a "near-critical" instability at finite $k$, whereas the rational eliminated-spin model (Model D) remains stable for all $k$. This provides a clear diagnostic: if a system is stable at high $k$ but shows high-curvature effects, a polynomial closure is likely incorrect.
• EDMD Validation:
  • Free Decay: Confirms a clean late-time one-pole hydrodynamic sector.
  • Harmonic Forcing: The spin channel ($\omega_z$) is measurable and exhibits a systematic phase lag relative to vorticity ($\zeta_z$).
  • Model Selection:
    • The spin-to-vorticity ratio strongly rejects adiabatic elimination in favor of retained spin ($\Delta$AIC/BIC $\gg 0$).
    • The multi-k vorticity response rejects a pure $k^2$ (Navier-Stokes) closure and, with sufficient statistics, favors the rational eliminated-spin kernel over the polynomial surrogate.
• Parameter Sensitivity: The distinction between retained spin and elimination becomes more pronounced as the spin relaxation time ($\tau_s$) increases (i.e., as the micro-inertia $J$ increases).

5. Significance

• Diagnostic Framework: The paper provides a practical "response-theoretic" language to distinguish between microscopic rotational physics, effective eliminated dynamics, and phenomenological higher-gradient closures. It moves beyond coefficient-matching to structural identification (pole counting and kernel shape).
• Resolution of Ambiguity: It resolves the ambiguity of high-curvature observables (like $k^4$ spectra), showing that their presence alone does not prove the existence of a specific closure type; the response function and stability properties are required for identification.
• Validation of Micropolar Theory: By using EDMD on rough spheres, the authors provide the first direct microscopic evidence that the "retained-spin" dynamics are physically real and distinguishable from adiabatic elimination in a many-particle system.
• Implications for Turbulence and Complex Fluids: The findings suggest that standard polynomial Burnett closures may be fundamentally flawed for predicting high-wavenumber behavior (due to instability or over-damping), and that rational kernels or explicit spin fields are necessary for accurate modeling of complex fluids with internal rotation.

In summary, Tsuzuki demonstrates that transverse linear response is a powerful tool for model discrimination, capable of separating retained micro-rotational physics from effective eliminated dynamics and identifying the structural limitations of polynomial higher-gradient approximations.