A Hybrid Method for Accelerated Simulation of Coulomb Collisions in a Plasma
If the collisional time scale for Coulomb collisions is comparable to the characteristic time scales for a plasma, then simulation of Coulomb collisions may be important for computation of kinetic plasma dynamics. This can be a computational bottleneck because of the large number of simulated particles and collisions (or phase-space resolution requirements in continuum algorithms), as well as the wide range of collision rates over the velocity distribution function. This paper considers Monte Carlo simulation of Coulomb collisions using the binary collision models of Takizuka & Abe and Nanbu. It presents a hybrid method for accelerating the computation of Coulomb collisions. The hybrid method represents the velocity distribution function as a combination of a thermal component (a Maxwellian distribution) and a kinetic component (a set of discrete particles). Collisions between particles from the thermal component preserve the Maxwellian; collisions between particles from the kinetic component are performed using the method of or Nanbu. Collisions between the kinetic and thermal components are performed by sampling a particle from the thermal component and selecting a particle from the kinetic component. Particles are also transferred between the two components according to thermalization and dethermalization probabilities, which are functions of phase space.
- Research Organization:
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- W-7405-ENG-48
- OSTI ID:
- 941391
- Report Number(s):
- UCRL-JRNL-235460; TRN: US0807316
- Journal Information:
- SIAM Journal Multiscale Modeling and Simulation, vol. 7, no. 2, August 20, 2008, pp. 865-887, Vol. 7, Issue 2
- Country of Publication:
- United States
- Language:
- English
Similar Records
ON THE COMPETITION BETWEEN RADIAL EXPANSION AND COULOMB COLLISIONS IN SHAPING THE ELECTRON VELOCITY DISTRIBUTION FUNCTION: KINETIC SIMULATIONS
A semikinetic model for early stage plasmasphere refilling. 1. Effects of Coulomb collisions