Calculation Of The Nanbu-Trubnikov Kernel: Implications For Numerical Modeling Of Coulomb Collisions
We investigate the accuracy of and assumptions underlying the numerical binary Monte-Carlo collision operator due to Nanbu [K. Nanbu, Phys. Rev. E 55 (1997)]. The numerical experiments that resulted in Nanbu's parameterized collision kernel are approximate realizations of the Coulomb-Lorentz pitch-angle scattering process, for which an analytical solution is available. It is demonstrated empirically that Nanbu's collision operator quite accurately recovers the effects of Coulomb-Lorentz pitch-angle collisions, or processes that approximate these even for very large values of the collisional time step. An investigation of the analytical solution shows that Nanbu's parameterized kernel is highly accurate for small values of the normalized collision time step, but loses some of its accuracy for larger values of the time step. Finally, a practical collision algorithm is proposed that for small-mass-ratio Coulomb collisions improves on the accuracy of Nanbu's algorithm.
- Research Organization:
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- W-7405-ENG-48
- OSTI ID:
- 962818
- Report Number(s):
- LLNL-CONF-414445; TRN: US0902949
- Resource Relation:
- Conference: Presented at: 2009 SIAM Annual Meeting, Denver, CO, United States, Jul 06 - Jul 10, 2009
- Country of Publication:
- United States
- Language:
- English
Similar Records
Particle simulation of Coulomb collisions: Comparing the methods of Takizuka and Abe and Nanbu
Particle Simulation of Coulomb Collisions: Comparing the Methods of Takizuka & Abe and Nanbu