Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Part 2: From Boltzmann to Fluid Equations

From Particles to Stars | Statistical Thinking Module 2 | ASTR 596

San Diego State University

Learning Objectives

By the end of Part 2, you will be able to:


Part 2: From Statistics to Fundamental Physics

2.1 The Boltzmann Equation: The Master Evolution Equation

Priority: 🟡 Standard Path

The Boltzmann equation is the master equation governing how probability distributions evolve in phase space. Think of it as Newton’s F=maF = ma but for probability clouds rather than individual particles. While Newton tells us how one particle’s position and velocity change over time, Boltzmann tells us how the probability of finding particles at various positions and velocities evolves. It’s the fundamental equation that bridges the microscopic world of individual particles to the macroscopic world of fluid dynamics and thermodynamics.

Imagine tracking not one particle but the probability cloud of where particles might be. This cloud flows through space (particles move), deforms under forces (acceleration changes velocities), and gets scrambled by collisions (randomization). The Boltzmann equation captures all three processes in one elegant framework.

The Boltzmann equation governs how distribution functions evolve in phase space:

ft+vrf+Fmvf=(ft)coll\boxed{\frac{\partial f}{\partial t} + \vec{v} \cdot \nabla_r f + \frac{\vec{F}}{m} \cdot \nabla_v f = \left(\frac{\partial f}{\partial t}\right)_{\text{coll}}}

Let’s understand each term physically:

This equation is exact but unsolvable for 1057 particles. The magic happens when we take moments.

Remember from Part 1 that the collision time between particles is much shorter than the star’s dynamical timescale in stellar interiors (τcollτdyn\tau_\text{coll} \ll \tau_\text{dyn})? That’s why we can often set the collision integral to zero locally — collisions have already done their work establishing the Maxwell-Boltzmann distribution. The distribution has thermalized so thoroughly that it maintains its equilibrium shape even as the star evolves. This is the magic of LTE: the collision term has already won the race, so we can ignore it in our macroscopic equations.

2.3 The Moment-Taking Machine: From Boltzmann to Fluid Equations

Priority: 🔴 Essential

Now comes the magic trick that transforms statistical mechanics into the equations you know and love. We’re going to multiply the unsolvable Boltzmann equation by different powers of velocity and integrate. Each multiplication extracts different physics—like using different filters on the same photograph reveals different features. The blue filter shows the sky, the red filter shows the sunset, the infrared filter shows the heat. Similarly, multiplying by 1 extracts mass flow, by vv extracts momentum flow, by v2v^2 extracts energy flow.

This procedure seems almost too simple to work, yet it transforms an equation tracking 1057 individual particles into the handful of smooth equations that govern stars, galaxies, and gas clouds. Watch carefully—this is where statistics becomes physics.

The Universal Procedure:

  1. Multiply Boltzmann equation by vnv^n

  2. Integrate over all velocities

  3. Get evolution equation for the nn-th moment

Let’s do this explicitly for the first few moments.

Zero-th Moment: Mass Conservation

Multiply the Boltzmann equation by particle mass mm (a constant) and integrate:

m[ft+vrf+Fmvf]d3v=0m \int \left[\frac{\partial f}{\partial t} + \vec{v} \cdot \nabla_r f + \frac{\vec{F}}{m} \cdot \nabla_v f\right] d^3v = 0

(We set the collision integral to zero since it conserves particle number by definition.)

Let’s work through each term carefully:

Term 1:

mftd3v=mtfd3v=ρtm \int \frac{\partial f}{\partial t} d^3v = m \frac{\partial}{\partial t} \int f d^3v = \frac{\partial \rho}{\partial t}

where ρ=nm\rho = n m is the mass density in g/cm³.

Term 2: mvrfd3v=mrvfd3v=(ρu)m \int \vec{v} \cdot \nabla_r f d^3v = m \nabla_r \cdot \int \vec{v} f d^3v = \nabla \cdot (\rho\vec{u})

where u=v\vec{u} = \langle \vec{v} \rangle is the mean velocity in cm/s.

Term 3 (Force term): mFmvfd3v=Fvfd3vm \int \frac{\vec{F}}{m} \cdot \nabla_v f d^3v = \int \vec{F} \cdot \nabla_v f d^3v

Using integration by parts: =f(vF)d3v+[Ff]v=v=+=0= -\int f (\nabla_v \cdot \vec{F}) d^3v + \underbrace{\left[\vec{F} \cdot f\right]_{v=-\infty}^{v=+\infty}}_{=0}

The boundary term vanishes because ff must decay faster than any polynomial as v|v| \to \infty (required for finite mass and energy). For uniform forces (like gravity), vF=0\nabla_v \cdot \vec{F} = 0, so the entire term vanishes.

Mathematical note: This is where the assumption of “no particles at infinite velocity” enters. It’s physically reasonable but mathematically crucial.

Result - The Continuity Equation:

ρt+(ρu)=0\boxed{\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho\vec{u}) = 0}

This is mass conservation! The density changes only due to flow divergence.

First Moment: Momentum Conservation

Multiply by mvim v_i (ii-th component of momentum) and integrate. After working through the algebra:

(ρui)t+xj(ρuiuj+Pij)=ρFi\boxed{\frac{\partial (\rho u_i)}{\partial t} + \frac{\partial}{\partial x_j}(\rho u_i u_j + P_{ij}) = \rho F_i}

where

Pij=ρ(viui)(vjuj)P_{ij} = \rho \langle (v_i - u_i)(v_j - u_j) \rangle

is the pressure tensor in dyne/cm².

For isotropic pressure (same in all directions): Pij=PδijP_{ij} = P\delta_{ij} where δij\delta_{ij} is the Kronecker delta.

This simplifies to:

ρDuDt=P+ρF\boxed{\rho \frac{D\vec{u}}{Dt} = -\nabla P + \rho \vec{F}}

This is the Euler equationNewton’s second law for fluids!

Second Moment: Energy Conservation

Multiply by 12mv2\frac{1}{2}mv^2 (kinetic energy) and integrate to get:

Et+[(E+P)u]=ρFu\boxed{\frac{\partial E}{\partial t} + \nabla \cdot [(E + P)\vec{u}] = \rho \vec{F} \cdot \vec{u}}

where EE is the energy density in erg/cm³. The pressure PP appears naturally in the energy flux — pressure does work on flowing fluid!

The Beautiful Pattern

Let’s step back and see what we’ve accomplished:

Moment OperationMultiply Boltzmann byIntegrate to getPhysical MeaningConservation Law
0th momentmmmfd3v=ρm \int f d^3v = \rhoMass densityMass conservation
1st momentmvmvmvfd3v=ρum \int vf d^3v = \rho uMomentum densityMomentum conservation
2nd momentmv2m v^2mv2fd3vEm \int v^2f d^3v \propto EEnergy densityEnergy conservation

Key Takeaway: Each moment extracts a different conservation law.

The procedure is universal — it works for any system where particles follow the Boltzmann equation!

Part 2 Synthesis: The Moment-Taking Framework

Priority: 🔴 Essential

You’ve just learned one of the most powerful techniques in physics: transforming unsolvable microscopic equations into tractable macroscopic ones through taking moments. This brief synthesis will cement your understanding before we see it in action.

The Universal Recipe

Building on the statistical foundation from Module 1 — where we learned that temperature is a distribution parameter and pressure emerges from ensemble averages — we’ve discovered a recipe that works for any system of particles:

Step 1: Start with the distribution f(r,v,t)f(r,v,t)

Step 2: Write the Boltzmann equation

Step 3: Take moments (multiply by vnv^n and integrate)

Step 4: Get conservation laws

Why This Works: Information Compression

Taking moments is fundamentally about information compression. Consider what we’re doing:

We’ve compressed 1058 numbers down to 103\sim 10^3 — a reduction factor of 1055!

The “lost” information? The precise velocity of particle number 8,745,293,048,571,293 at this exact instant. We don’t care, and neither does nature at macroscopic scales.

Key Insight: Pressure IS Variance

The most profound realization from our moment-taking:

P=ρ(vu)2=nmVar(v)P = \rho \langle(v - u)^2\rangle = nm \cdot \text{Var}(v)

Pressure isn’t just “related to” velocity spread — it IS mass density times velocity variance. This identity (not approximation!) means:

With this framework firmly in mind, you’re ready to see it applied to real stellar physics.