Retained-spin micropolar hydrodynamics from the Boltzmann--Curtiss equation: a generalized Chapman--Enskog construction

This paper derives a retained-spin micropolar hydrodynamic closure from the Boltzmann--Curtiss equation using a generalized Chapman--Enskog construction, explicitly identifying the emergence of standard constitutive coefficients and providing explicit dilute-gas estimates for rotational viscosity and spin-diffusion parameters that are validated by molecular-dynamics simulations.

The Gist

Imagine you are watching a crowded dance floor. In a standard fluid (like water), we usually only care about two things: where the dancers are moving (velocity) and how crowded the room is (density). We assume that if a dancer spins around, it's just a tiny, invisible detail that doesn't change the overall flow of the crowd.

But what if the dancers are spinning tops? What if their individual spin is a major part of how they bump into each other and how the whole crowd moves?

This paper is about building a new set of rules for fluids made of "spinning tops" (like tiny, rough spheres). The author, Satori Tsuzuki, takes a complex mathematical recipe called the Boltzmann–Curtiss equation (which describes how individual particles behave) and simplifies it into a usable "macro" recipe for the whole fluid.

Here is the breakdown of the paper's journey, using everyday analogies:

1. The Problem: The "Spinning" Gap

In standard physics, we have a perfect recipe for how fluids flow. But when particles have their own spin (like a basketball spinning while flying through the air), the old recipes get messy.

• The Old Way: Scientists often treated the spin as something that disappears instantly or is too hard to calculate, so they ignored it or guessed at it.
• The New Way: This paper says, "Let's keep the spin!" It treats the average spin of the crowd as a real, important variable, just like temperature or speed.

2. The Method: The "Generalized Chapman-Enskog" Construction

Think of the Chapman-Enskog method as a way to translate a language.

• The Micro Language: The language of individual particles (colliding, spinning, bouncing).
• The Macro Language: The language of the fluid (pressure, flow, viscosity).

Usually, this translation is strict: you assume the particles settle down instantly. But here, the author uses a "Generalized" translation.

• The Analogy: Imagine trying to describe a busy highway. A strict translator might say, "Cars move at constant speed." But a generalized translator says, "Cars move at constant speed, but we also need to account for the fact that some drivers are slightly distracted and their cars wobble a bit before they settle."
• In this paper, the "wobble" is the spin. The author keeps the spin in the equation as a "quasi-slow variable," meaning it changes slowly enough to matter, but fast enough to need special math to handle.

3. The Big Discovery: Two Types of Friction

One of the most important findings is about friction (viscosity).

• Standard Friction: When you stir honey, it resists because the layers slide past each other. This is "shear viscosity."
• Rotational Friction (The New Guy): When the spinning tops bump into each other, they exchange spin. This creates a new kind of friction called rotational viscosity ($\eta_r$).

The Analogy:
Imagine a room full of people holding spinning fidget spinners.

1. If they just walk past each other, they create standard friction (like honey).
2. But if they bump into each other, their fidget spinners crash. One slows down, the other speeds up. This exchange of spin creates a "twisting" force that resists the flow in a totally different way.
   The paper proves mathematically that this "twisting friction" comes only from the collisions (the bumps), not from the particles just moving through space.

4. The "Rough Sphere" Test

To make sure their new math works, the author tested it on a specific, simple model: Perfectly Rough Elastic Hard Spheres.

• The Analogy: Imagine billiard balls that are covered in sandpaper. When they hit, they don't just bounce; they grip and spin each other.
• The author calculated exactly how much "rotational friction" these sandpaper balls should create based on how dense the crowd is and how rough the balls are.
• The Result: They found a clear rule: The more crowded the room (density), the more the friction grows (specifically, it grows with the square of the density). The rougher the balls, the more friction they create, up to a limit.

5. The Computer Check (The "Reality Check")

Math is great, but does it match reality? The author ran computer simulations (Event-Driven Molecular Dynamics).

• They created a virtual box of 8,000 spinning, sandpaper-covered balls.
• They spun the whole box and watched how fast the spin died out.
• The Verdict: The computer simulation matched the new math almost perfectly! The "twisting friction" behaved exactly as the new equations predicted.

6. What's Left? (The "Fine Print")

The paper admits that while they nailed the "twisting friction" ($\eta_r$), there are still some tricky parts regarding how spin spreads out sideways (transverse spin diffusion).

• The Analogy: They figured out exactly how much the spinning tops slow each other down when they crash. But figuring out exactly how a "wave" of spinning moves through the crowd from left to right is still a bit fuzzy in the computer simulations. It's a good start, but more work is needed there.

Summary

This paper is a bridge builder.

• It connects the messy, chaotic world of individual spinning particles to the smooth, predictable world of fluid flow.
• It introduces a new, precise way to calculate how "spin" affects fluid friction.
• It proves that if you have a fluid made of rough, spinning particles, you can't ignore their spin—it creates a unique kind of "rotational friction" that standard physics misses.

In short: The author took a complex, scattered set of ideas about spinning fluids, organized them into a single, clear mathematical story, and proved with computer simulations that the story holds up. It's like finally finding the instruction manual for a complex toy that everyone has been guessing how to play with for decades.

Technical Summary

1. Problem Statement

Classical fluid mechanics assumes that local micro-rotation (spin) is slaved to the macroscopic vorticity or relaxes instantaneously. However, for fluids composed of particles with rotational degrees of freedom (e.g., rough spheres, polyatomic molecules), the local mean spin $\omega_0$ can act as an independent hydrodynamic field.

The paper addresses the lack of a unified, self-contained derivation that connects the microscopic Boltzmann–Curtiss kinetic equation to the macroscopic micropolar hydrodynamic equations while explicitly retaining the mean spin as a "quasi-slow" variable. Previous literature often scattered these derivations across continuum mechanics, kinetic theory, and granular flow, leaving ambiguities regarding:

• The exact origin of the antisymmetric stress channel.
• The distinction between kinetic stress and collisional transfer mechanisms.
• The precise conditions under which standard micropolar constitutive laws emerge.

2. Methodology

The author employs a rigorous generalized Chapman–Enskog construction tailored for retained-spin systems. The methodology proceeds through the following stages:

• Exact Balance Laws: Starting from the Boltzmann–Curtiss equation, the author derives exact moment balances for mass, linear momentum, and total angular momentum. By subtracting the orbital angular momentum balance from the total, an exact intrinsic-spin balance is obtained.
• Generalized Expansion: Unlike the strict Chapman–Enskog method where only collision invariants are slow variables, this work treats the local mean spin $\omega_0$ as a quasi-slow variable. This requires an extended thermodynamic ordering where the residual spin-relaxation source on the quasi-equilibrium manifold is treated as $O(\epsilon)$ (same order as gradient corrections).
• Irreducible Decomposition: The first-order source term in the linearized kinetic equation is decomposed into irreducible tensor sectors:
  • Scalar (bulk)
  • Axial (spin-mismatch)
  • Symmetric-traceless (shear/spin-diffusion)
• Constitutive Derivation: The author solves the linearized kinetic equations for the response functions in each sector to derive the constitutive relations for stress ($\sigma_{ij}$) and couple stress ($m_{ij}$).
• Dilute-Gas Estimates: For the specific case of perfectly rough elastic hard spheres, explicit analytical estimates are derived for the rotational viscosity ($\eta_r$) and the transverse spin-diffusion coefficient ($\beta + \gamma$) using collision brackets.
• Numerical Validation: Targeted Event-Driven Molecular Dynamics (EDMD) simulations are performed to verify the theoretical predictions, specifically checking the density scaling ($n^2$) and roughness dependence ($K/(K+1)$) of the rotational viscosity.

3. Key Contributions

A. Theoretical Framework

• Clarification of Antisymmetric Stress: The paper explicitly demonstrates that the one-particle kinetic stress contributes only to the symmetric part of the stress tensor. The antisymmetric stress (and thus rotational viscosity $\eta_r$) arises entirely from an intrinsic/collisional transfer channel. This resolves long-standing confusion about the origin of rotational viscosity.
• Retained-Spin Closure: It provides a formal derivation of the micropolar equations where $\omega_0$ is an independent field, rather than being eliminated. The resulting equations are:
  $$\rho \frac{D\mathbf{u}}{Dt} = -\nabla P + (\eta + \eta_r)\nabla^2 \mathbf{u} + \left(\frac{\eta}{3} + \xi - \eta_r\right)\nabla(\nabla\cdot\mathbf{u}) + 2\eta_r(\nabla\times\boldsymbol{\omega}_0) + \rho\mathbf{F}$$
  $$\rho J \frac{D\boldsymbol{\omega}_0}{Dt} = (\beta+\gamma)\nabla^2\boldsymbol{\omega}_0 + (\alpha+\beta-\gamma)\nabla(\nabla\cdot\boldsymbol{\omega}_0) + 2\eta_r(\boldsymbol{\zeta} - 2\boldsymbol{\omega}_0) + \rho\mathbf{G}$$
  where $\boldsymbol{\zeta} = \nabla \times \mathbf{u}$ is the vorticity.

B. Analytical Results for Rough Spheres

For perfectly rough elastic hard spheres (diameter $a$, mass $m$, moment of inertia $I$, reduced inertia $K = 4I/ma^2$):

1. Rotational Viscosity ($\eta_r$): Derived via homogeneous spin relaxation.
   $$\eta_r \propto n^2 \frac{K}{K+1} \sqrt{\frac{k_B T}{m}}$$
   This confirms $\eta_r$ is a collisional exchange coefficient (scaling with $n^2$) rather than a kinetic transport coefficient.
2. Transverse Spin-Diffusion ($\beta + \gamma$): Derived via transport-relaxation calculation.
   $$\beta + \gamma \propto n \sqrt{\frac{k_B T}{m}}$$
   This scales linearly with density, characteristic of standard kinetic transport.

C. Numerical Verification

• Homogeneous Relaxation: EDMD simulations confirm the predicted $n^2$ scaling and $K/(K+1)$ dependence for $\eta_r$ across a wide range of densities ($\phi \in [0.005, 0.050]$) and roughness parameters ($K \in [0.05, 1.0]$).
• Finite-$k$ Transverse Mode: Simulations of transverse spin modes qualitatively support the coupled relaxation dynamics but highlight numerical challenges in extracting precise off-diagonal coefficients at higher densities, suggesting the need for refined future studies.

4. Results

• Structural Clarity: The derivation clearly separates exact balance laws, formal first-order constitutive structures, and specific model estimates.
• Coefficient Estimates: Explicit formulas are provided for the rotational viscosity and spin-diffusion coefficients of rough spheres, filling a gap in the literature where such estimates were previously scattered or incomplete.
• Validation: The theoretical scaling laws for $\eta_r$ are robustly supported by molecular dynamics data in the dilute-to-moderate density regime. Deviations at high densities are attributed to the breakdown of simple single-exponential relaxation rather than a failure of the underlying kinetic theory.

5. Significance

This paper serves as a reconstructive and expository cornerstone for micropolar hydrodynamics. Its significance lies in:

1. Unification: It bridges the gap between the Boltzmann–Curtiss kinetic equation and continuum micropolar theory in a single, self-contained derivation.
2. Mechanism Identification: It definitively identifies the collisional transfer mechanism as the sole source of rotational viscosity, correcting misconceptions that it arises from the kinetic stress tensor.
3. Benchmarking: It provides explicit, testable predictions for dilute gases of rough spheres, validated by modern simulation techniques.
4. Foundation for Future Work: By clarifying which parts of the closure are exact and which are estimates, it sets the stage for more complex studies involving dense fluids, inelastic collisions, and higher-order constitutive models (as referenced in the companion manuscript [22]).

In summary, Tsuzuki provides a rigorous, line-by-line derivation that validates the retained-spin micropolar framework, offering precise coefficient estimates and numerical confirmation for the rotational transport properties of rough-particle gases.