Collision Operators in Kinetic Codes
Collision operators model the small-angle Coulomb scattering that relaxes distributions, a subtle numerical component that shapes kinetic-code results.
Why collisions need care
Even in hot, nearly collisionless plasmas, Coulomb collisions matter: they set the neoclassical transport floor, damp fine velocity-space structure, and control current-drive efficiency. The collision operator is the term in a kinetic equation that models this scattering. Getting it right is a recognized source of difference between gyrokinetic codes, so its treatment deserves attention in its own right.
The full Landau collision operator is nonlinear and expensive, so codes use approximations of varying fidelity, and the choice affects both accuracy and cost.
Conservation properties
A physical collision operator conserves particles, momentum, and energy and drives the distribution toward a Maxwellian. Simplified operators, such as a pitch-angle-scattering-only Lorentz operator, violate momentum conservation and can produce spurious transport. Careful codes restore conservation by adding correction terms, which is essential for correct neoclassical and current-drive results.
Velocity-space resolution
Collisions create sharp structure at the trapped-passing boundary and smooth fine-scale velocity structure elsewhere. Resolving the boundary layer while damping under-resolved structure requires adequate velocity-space grids, and under-resolution is a common cause of numerical error in kinetic runs.
Design relevance
For the Hyperion breeder, correct collision operators underpin credible predictions of neoclassical transport, bootstrap current, and current-drive efficiency, all of which feed the design. Because the collision treatment is a known code-to-code difference, cross-verification of these predictions is part of the honest simulation record before construction.
- Models small-angle Coulomb scattering
- Must conserve particles, momentum, energy
- Simplified operators need conservation corrections
- Resolution at the trapped-passing boundary matters