A Hybrid Method for Accelerated Simulation of Coulomb Collisions in a Plasma
Abstract
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 and Abe and of 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 Takizuka and Abe or of 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.
MSC codes
Keywords
Get full access to this article
View all available purchase options and get full access to this article.
References
1.
G. A. Bird, Molecular Gas Dynamics and the Direct Simulation of Gas Flows, Oxford University Press, Oxford, UK, 1998.
2.
A. V. Bobylev and K. Nanbu, Theory of collision algorithms for gases and plasmas based on the Boltzmann equation and the Landau-Fokker-Planck equation, Phys. Rev. E (3), 61 (2000), pp. 4576–4582.
3.
A. M. Dimits and W. W. Lee, Partially linearized algorithms in gyrokinetic particle simulation, J. Comput. Phys., 107 (1993), pp. 309–323.
4.
A. M. Dimits, T. J. Williams, J. A. Byers, and B. I. Cohen, Scalings of ion-temperature-gradient-driven anomalous transport in tokamaks, Phys. Rev. Lett., 77 (1996), pp. 71–74.
5.
G. Hu and J. A. Krommes, Generalized weighting scheme for f particle-simulation method, Phys. Plasmas, 1 (1994), pp. 863–874.
6.
M. Kotschenreuther, Numerical simulation, Bull. Am. Phys. Soc., 33 (1988), pp. 2107–2108.
7.
N. A. Krall and A. W. Trivelpiece, Principles of Plasma Physics, McGraw–Hill, New York, 1973.
8.
W. M. Mannheimer, M. Lampe, and G. Joyce, Langevin representation of Coulomb collisions in PIC simulations, J. Comput. Phys., 138 (1997), pp. 563–584.
9.
K. Nanbu, Theory of cumulative small-angle collisions in plasmas, Phys. Rev. E (3), 55 (1997), pp. 4642–4652.
10.
L. Pareschi and R. E. Caflisch, An implicit Monte Carlo method for rarefied gas dynamics, J. Comput. Phys., 154 (1999), pp. 90–116.
11.
M. N. Rosenbluth, W. M. MacDonald, and D. L. Judd, Fokker-Planck equation for an inverse-square force, Phys. Rev. (2), 107 (1957), pp. 1–6.
12.
M. Sherlock, A Monte-Carlo method for Coulomb collisions in hybrid plasma models, J. Comput. Phys., 227 (2008), pp. 2286–2292.
13.
L. Spitzer, Jr., Physics of Fully Ionized Gases, 2nd ed., Interscience, New York, 1967.
14.
T. Takizuka and H. Abe, A binary collision model for plasma simulation with a particle code, J. Comput. Phys., 25 (1977), pp. 205–219.
15.
B. A. Trubnikov, Particle interactions in a fully ionized plasma, in Reviews of Plasma Physics, Vol. 1, Consultant Bureau, New York, 1965, pp. 105–204.
16.
C. Wang, T. Lin, R. E. Caflisch, B. Cohen, and A. Dimits, Particle simulation of Coulomb collisions: Comparing the methods of Takizuka & Abe and Nanbu, J. Comput. Phys., 227 (2008), pp. 4308–4329.
Information & Authors
Information
Published In
Copyright
Copyright © 2008 Society for Industrial and Applied Mathematics.
History
Submitted: 9 October 2007
Accepted: 5 May 2008
Published online: 20 August 2008
MSC codes
Keywords
Authors
Metrics & Citations
Metrics
Citations
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.
Cited By
- A Coulomb collision model based on collision cross section and its effects on kinetic characteristics of atmospheric thermal microplasmasPhysics of Fluids, Vol. 36, No. 10 | 3 October 2024
- The collisional particle-in-cell method for the Vlasov–Maxwell–Landau equationsJournal of Plasma Physics, Vol. 90, No. 4 | 18 October 2024
- Particle simulation methods for the Landau-Fokker-Planck equation with uncertain dataJournal of Computational Physics, Vol. 503 | 1 Apr 2024
- A tutorial overview of the angular scattering models of electron–neutral, ion–neutral, neutral–neutral, and Coulomb collisions in Monte Carlo collision modeling on low-temperature plasmaPlasma Sources Science and Technology, Vol. 33, No. 2 | 16 February 2024
- Stochastic Galerkin particle methods for kinetic equations of plasmas with uncertaintiesJournal of Computational Physics, Vol. 479 | 1 Apr 2023
- Stochastic Galerkin Particle Methods for Kinetic Equations of Plasmas with UncertaintiesSSRN Electronic Journal, Vol. 106 | 1 Jan 2022
- Accurate particle time integration for solving Vlasov-Fokker-Planck equations with specified electromagnetic fieldsJournal of Computational Physics, Vol. 387 | 1 Jun 2019
- Hybrid fluid–kinetic neutral model based on a kinetic correction on a fluid solutionContributions to Plasma Physics, Vol. 58, No. 6-8 | 2 March 2018
- Asymptotic preserving and time diminishing schemes for rarefied gas dynamicKinetic & Related Models, Vol. 10, No. 3 | 1 Jan 2017
- Multi-scale kinetic scheme for the collisional Vlasov-Poisson systemComputers & Fluids, Vol. 140 | 1 Nov 2016
- A hybrid method with deviational particles for spatial inhomogeneous plasmaJournal of Computational Physics, Vol. 309 | 1 Mar 2016
- A Monte Carlo method with negative particles for Coulomb collisionsJournal of Computational Physics, Vol. 298 | 1 Oct 2015
- Monte Carlo SimulationEncyclopedia of Applied and Computational Mathematics | 21 November 2015
- Numerical methods for plasma physics in collisional regimesJournal of Plasma Physics, Vol. 81, No. 1 | 10 October 2014
- An entropy based thermalization scheme for hybrid simulations of Coulomb collisionsJournal of Computational Physics, Vol. 273 | 1 Sep 2014
- Numerical methods for kinetic equationsActa Numerica, Vol. 23 | 12 May 2014
- First principles calculation of the effect of Coulomb collisions in partially ionized gasesPhysics of Plasmas, Vol. 21, No. 4 | 16 April 2014
- Monte Carlo methods and their analysis for Coulomb collisions in multicomponent plasmasJournal of Computational Physics, Vol. 246 | 1 Aug 2013
- On efficient simulations of multiscale kinetic transportPhilosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, Vol. 371, No. 1982 | 13 January 2013
- Fast elliptic solvers in cylindrical coordinates and the Coulomb collision operatorJournal of Computational Physics, Vol. 230, No. 21 | 1 Sep 2011
- Understanding the accuracy of Nanbu’s numerical Coulomb collision operatorJournal of Computational Physics, Vol. 228, No. 13 | 1 Jul 2009
- Small-angle Coulomb collision model for particle-in-cell simulationsJournal of Computational Physics, Vol. 228, No. 5 | 1 Mar 2009
View Options
- Access via your Institution
- Questions about how to access this content? Contact SIAM at [email protected].