Collision operators for partially linearized particle simulation codes
- Lawrence Livermore National Laboratory, University of California, Livermore, California 94550 (United States)
Algorithms are presented for energy- and momentum-conserving like-particle Coulomb collisions in partially linearized ([delta][ital f]) particle simulations. They are developed and implemented in particular for gyrokinetic simulation models of a strongly magnetized plasma. The collision operators include both drag and diffusion terms, are not restricted to a single or few Fourier modes, and approximately conserve both momentum and energy locally in space in a statistical sense. The first algorithm is a many-mode generalization of a test-particle--plus--source algorithm previously proposed. The second is easier to implement and improves upon the first significantly by not requiring many time steps for good conservation. Implementations for the case for ion-ion collisions are given and conservation properties are demonstrated, both directly with non-self-consistent test simulation runs and indirectly with self-consistent runs. The computational cost of particle pushing and solving for fields depends on the relative collisionality and can result in a tripling of the total computational costs if collisions are done at each time step, but typically will be a small fraction of the total simulation cost. It is also shown that binary-collision-based algorithms are unsuitable for partially linearized simulations.
- DOE Contract Number:
- W-7405-ENG-48
- OSTI ID:
- 5282053
- Journal Information:
- Physical Review. E, Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics; (United States), Vol. 49:1; ISSN 1063-651X
- Country of Publication:
- United States
- Language:
- English
Similar Records
Explicit structure-preserving geometric particle-in-cell algorithm in curvilinear orthogonal coordinate systems and its applications to whole-device 6D kinetic simulations of tokamak physics
Center for Gyrokinetic/MHD Hybrid Simulation of Energetic Particle Physics in Toroidal Plasmas (CSEPP). Final report
Related Subjects
PLASMA SIMULATION
ALGORITHMS
CALCULATION METHODS
COMPUTER CODES
CONSERVATION LAWS
COULOMB SCATTERING
DIFFUSION
ION-ION COLLISIONS
MAGNETIC FIELDS
BASIC INTERACTIONS
COLLISIONS
ELASTIC SCATTERING
ELECTROMAGNETIC INTERACTIONS
INTERACTIONS
ION COLLISIONS
MATHEMATICAL LOGIC
SCATTERING
SIMULATION
700380* - Elementary & Classical Processes in Plasmas- (1992-)