Adaptive time-stepping Monte Carlo integration of Coulomb collisions

Sarkimaki, Konsta; Hirvijoki, E.; Terava, J.

doi:10.1016/j.cpc.2017.09.025

Title: Adaptive time-stepping Monte Carlo integration of Coulomb collisions

Abstract

Here, we report an accessible and robust tool for evaluating the effects of Coulomb collisions on a test particle in a plasma that obeys Maxwell–Jüttner statistics. The implementation is based on the Beliaev–Budker collision integral which allows both the test particle and the background plasma to be relativistic. The integration method supports adaptive time stepping, which is shown to greatly improve the computational efficiency. The Monte Carlo method is implemented for both the three-dimensional particle momentum space and the five-dimensional guiding center phase space.

Authors:

^[1]; Hirvijoki, E. ^[2]; Terava, J. ^[1]

Aalto Univ., Espoo (Finland)
Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)

Publication Date:: Thu Oct 12 00:00:00 EDT 2017

Research Org.:: Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)

Sponsoring Org.:: USDOE

OSTI Identifier:: 1441004

Alternate Identifier(s):: OSTI ID: 1564327

Grant/Contract Number:: 259675; 298126; AC02–09CH11466; AC02–09-CH11466

Resource Type:: Accepted Manuscript

Journal Name:: Computer Physics Communications

Additional Journal Information:: Journal Volume: 222; Journal Issue: C; Journal ID: ISSN 0010-4655

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; 97 MATHEMATICS AND COMPUTING; Coulomb collision; Monte Carlo; Fokker–Planck equation; Milstein method

Citation Formats


                    Sarkimaki, Konsta, Hirvijoki, E., and Terava, J. Adaptive time-stepping Monte Carlo integration of Coulomb collisions.  United States: N. p., 2017. 
Web.  doi:10.1016/j.cpc.2017.09.025.

Copy to clipboard


                    Sarkimaki, Konsta, Hirvijoki, E., & Terava, J. Adaptive time-stepping Monte Carlo integration of Coulomb collisions.  United States.  https://doi.org/10.1016/j.cpc.2017.09.025

Copy to clipboard


                    Sarkimaki, Konsta, Hirvijoki, E., and Terava, J. Thu .  
"Adaptive time-stepping Monte Carlo integration of Coulomb collisions".  United States.  https://doi.org/10.1016/j.cpc.2017.09.025.  https://www.osti.gov/servlets/purl/1441004.

Copy to clipboard


                    
@article{osti_1441004,

  title        = {Adaptive time-stepping Monte Carlo integration of Coulomb collisions},

  author       = {Sarkimaki, Konsta and Hirvijoki, E. and Terava, J.},

  abstractNote = {Here, we report an accessible and robust tool for evaluating the effects of Coulomb collisions on a test particle in a plasma that obeys Maxwell–Jüttner statistics. The implementation is based on the Beliaev–Budker collision integral which allows both the test particle and the background plasma to be relativistic. The integration method supports adaptive time stepping, which is shown to greatly improve the computational efficiency. The Monte Carlo method is implemented for both the three-dimensional particle momentum space and the five-dimensional guiding center phase space.},

  doi          = {10.1016/j.cpc.2017.09.025},

  journal      = {Computer Physics Communications},

  number       = C,

  volume       = 222,

  place        = {United States},

  year         = {Thu Oct 12 00:00:00 EDT 2017},

  month        = {Thu Oct 12 00:00:00 EDT 2017}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (Publisher)

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.1016/j.cpc.2017.09.025

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 6 works

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Works referenced in this record:

Variable Step Size Control in the Numerical Solution of Stochastic Differential Equations
journal, October 1997

Gaines, J. G.; Lyons, T. J.
SIAM Journal on Applied Mathematics, Vol. 57, Issue 5
DOI: 10.1137/S0036139995286515

Multilevel Monte Carlo simulation of Coulomb collisions
journal, October 2014

Rosin, M. S.; Ricketson, L. F.; Dimits, A. M.
Journal of Computational Physics, Vol. 274
DOI: 10.1016/j.jcp.2014.05.030

Higher-order time integration of Coulomb collisions in a plasma using Langevin equations
journal, June 2013

Dimits, A. M.; Cohen, B. I.; Caflisch, R. E.
Journal of Computational Physics, Vol. 242
DOI: 10.1016/j.jcp.2013.01.038

Dynamical friction in a relativistic plasma
journal, May 2014

Pike, O. J.; Rose, S. J.
Physical Review E, Vol. 89, Issue 5
DOI: 10.1103/PhysRevE.89.053107

A guiding-center Fokker–Planck collision operator for nonuniform magnetic fields
journal, September 2004

Brizard, Alain J.
Physics of Plasmas, Vol. 11, Issue 9
DOI: 10.1063/1.1780532

Monte Carlo implementation of a guiding-center Fokker-Planck kinetic equation
journal, September 2013

Hirvijoki, E.; Brizard, A.; Snicker, A.
Physics of Plasmas, Vol. 20, Issue 9
DOI: 10.1063/1.4820951

An adaptive timestepping algorithm for stochastic differential equations
journal, December 2003

Lamba, H.
Journal of Computational and Applied Mathematics, Vol. 161, Issue 2
DOI: 10.1016/j.cam.2003.05.001

Monte Carlo evaluation of transport coefficients
journal, January 1981

Boozer, Allen H.
Physics of Fluids, Vol. 24, Issue 5
DOI: 10.1063/1.863445

Volume-preserving algorithm for secular relativistic dynamics of charged particles
journal, April 2015

Zhang, Ruili; Liu, Jian; Qin, Hong
Physics of Plasmas, Vol. 22, Issue 4
DOI: 10.1063/1.4916570

Slowing-Down Dynamics of Fast Particles in Plasmas via the Fokker-Planck Equation
journal, January 2004

Anderson, D.; Lisak, M.; Andersson, F.
Nuclear Science and Engineering, Vol. 146, Issue 1
DOI: 10.13182/NSE04-A2396

An adaptive timestepping algorithm for stochastic differential equations
journal, December 2003

Lamba, H.
Journal of Computational and Applied Mathematics, Vol. 161, Issue 2
DOI: 10.1016/j.cam.2003.05.001

On the simulation of iterated Itô integrals
journal, January 2001

Rydén, Tobias; Wiktorsson, Magnus
Stochastic Processes and their Applications, Vol. 91, Issue 1
DOI: 10.1016/s0304-4149(00)00053-3

Monte Carlo implementation of a guiding-center Fokker-Planck kinetic equation
journal, September 2013

Hirvijoki, E.; Brizard, A.; Snicker, A.
Physics of Plasmas, Vol. 20, Issue 9
DOI: 10.1063/1.4820951

Differential form of the collision integral for a relativistic plasma
journal, October 1987

Braams, Bastiaan J.; Karney, Charles F. F.
Physical Review Letters, Vol. 59, Issue 16
DOI: 10.1103/physrevlett.59.1817

A guiding-center Fokker-Planck collision operator for nonuniform magnetic fields
text, January 2004

Brizard, Alain J.
arXiv
DOI: 10.48550/arxiv.physics/0405052

Works referencing / citing this record:

Electron acceleration in a JET disruption simulation
journal, August 2018

Sommariva, C.; Nardon, E.; Beyer, P.
Nuclear Fusion, Vol. 58, Issue 10
DOI: 10.1088/1741-4326/aad47d

Electron acceleration in a JET disruption simulation
text, January 2018

Sommariva, Cristian; Nardon, Eric; Beyer, Peter
arXiv
DOI: 10.48550/arxiv.1805.05655

Similar Records in DOE PAGES and OSTI.GOV collections:

Kinetics of electron-positron pair plasmas using an adaptive Monte Carlo method

Journal Article Pilla, R P ; Shaham, J - Astrophysical Journal

A new algorithm for implementing the adaptive Monte Carlo method is given. It is used to solve the Boltzmann equations that describe the time evolution of a nonequilibrium electron-positron pair plasma containing high-energy photons. These are coupled nonlinear integro-differential equations. The collision kernels for the photons as well as pairs are evaluated for Compton scattering, pair annihilation and creation, bremsstrahlung, and Coulomb collisions. They are given as multidimensional integrals which are valid for all energies. For an homogeneous and isotropic plasma with no particle escape, the equilibrium solution is expressed analytically in terms of the initial conditions. For two specificmore »« less
https://doi.org/10.1086/304534
Linearized Coulomb Collision Operator for Simulation of Interpenetrating Plasma Streams

Journal Article Dimits, A. M. ; Banks, J. W. ; Berger, R. L. ; ... - IEEE Transactions on Plasma Science

We present an extension of a linearized Coulomb collision operator, previously used in several Eulerian kinetic codes for like-species collisions and unlike-species collisions in the case where the backgrounds about which the linearization is made all are in collisional equilibrium, to the situation of interpenetrating plasma streams. In the latter case, the backgrounds can not be taken to be in equilibrium and a significant generalization is required. Our development is targeted towards the Eulerian kinetic plasma code LOKI, which evolves the Vlasov- Poisson or Vlasov-Maxwell system in a Cartesian “2+2-dimensional” phase space. The extended operator has been implemented in amore »« less
Cited by 3
https://doi.org/10.1109/TPS.2019.2897790

Full Text Available
Vector implementation of nonlinear Monte Carlo Coulomb collisions

Journal Article Wang, W X ; Okamoto, M ; Nakajima, N ; ... - Journal of Computational Physics

A nonlinear Monte Carlo collision model based on Coulomb binary collisions is presented with emphasis on its efficient implementation. The presented Monte Carlo collision operator, which has a simple form, fulfills particle number, momentum, and energy conservation laws quasi-locally and is equivalent to the nonlinear Fokker-Planck operator in the limit of infinitesimally small time interval of the binary collisions. Two vectorizable algorithms are designed for its fast implementation. Drastic speedup is obtained by the algorithms on vector computers. Various test simulations regarding relaxation processes, electrical conductivity, etc. are carried out in the velocity space. The test results, which are inmore »« less
https://doi.org/10.1006/jcph.1996.0204
Differential form of the collision integral for a relativistic plasma

Technical Report Braams, B J ; Karney, C F.F.

A differential formulation for the Beliaev and Budker relativistic collision integral is presented. This permits the rapid numerical evaluation of the collision integral. The decomposition into spherical harmonics allows the collision operator to be expressed in terms of one-dimensional integrals for simple background distributions. This is useful in carrying out analytical work. It also provides a convenient method for calculating the boundary conditions for the potentials. 6 refs.
https://doi.org/10.2172/6268312

Full Text Available
Self-consistent inclusion of classical large-angle Coulomb collisions in plasma Monte Carlo simulations

Journal Article Turrell, A.E., E-mail: a.turrell09@imperial.ac.uk ; Sherlock, M. ; Rose, S. J. - Journal of Computational Physics

Large-angle Coulomb collisions allow for the exchange of a significant proportion of the energy of a particle in a single collision, but are not included in models of plasmas based on fluids, the Vlasov–Fokker–Planck equation, or currently available plasma Monte Carlo techniques. Their unique effects include the creation of fast ‘knock-on’ ions, which may be more likely to undergo certain reactions, and distortions to ion distribution functions relative to what is predicted by small-angle collision only theories. We present a computational method which uses Monte Carlo techniques to include the effects of large-angle Coulomb collisions in plasmas and which self-consistentlymore »« less
https://doi.org/10.1016/J.JCP.2015.06.034

Similar Records

Title: Adaptive time-stepping Monte Carlo integration of Coulomb collisions

Abstract

Citation Formats

Variable Step Size Control in the Numerical Solution of Stochastic Differential Equations journal, October 1997

Multilevel Monte Carlo simulation of Coulomb collisions journal, October 2014

Higher-order time integration of Coulomb collisions in a plasma using Langevin equations journal, June 2013

Dynamical friction in a relativistic plasma journal, May 2014

A guiding-center Fokker–Planck collision operator for nonuniform magnetic fields journal, September 2004

Monte Carlo implementation of a guiding-center Fokker-Planck kinetic equation journal, September 2013

An adaptive timestepping algorithm for stochastic differential equations journal, December 2003

Monte Carlo evaluation of transport coefficients journal, January 1981

Volume-preserving algorithm for secular relativistic dynamics of charged particles journal, April 2015

Slowing-Down Dynamics of Fast Particles in Plasmas via the Fokker-Planck Equation journal, January 2004

An adaptive timestepping algorithm for stochastic differential equations journal, December 2003

On the simulation of iterated Itô integrals journal, January 2001

Monte Carlo implementation of a guiding-center Fokker-Planck kinetic equation journal, September 2013

Differential form of the collision integral for a relativistic plasma journal, October 1987

A guiding-center Fokker-Planck collision operator for nonuniform magnetic fields text, January 2004

Electron acceleration in a JET disruption simulation journal, August 2018

Electron acceleration in a JET disruption simulation text, January 2018

Variable Step Size Control in the Numerical Solution of Stochastic Differential Equations
journal, October 1997

Multilevel Monte Carlo simulation of Coulomb collisions
journal, October 2014

Higher-order time integration of Coulomb collisions in a plasma using Langevin equations
journal, June 2013

Dynamical friction in a relativistic plasma
journal, May 2014

A guiding-center Fokker–Planck collision operator for nonuniform magnetic fields
journal, September 2004

Monte Carlo implementation of a guiding-center Fokker-Planck kinetic equation
journal, September 2013

An adaptive timestepping algorithm for stochastic differential equations
journal, December 2003

Monte Carlo evaluation of transport coefficients
journal, January 1981

Volume-preserving algorithm for secular relativistic dynamics of charged particles
journal, April 2015

Slowing-Down Dynamics of Fast Particles in Plasmas via the Fokker-Planck Equation
journal, January 2004

An adaptive timestepping algorithm for stochastic differential equations
journal, December 2003

On the simulation of iterated Itô integrals
journal, January 2001

Monte Carlo implementation of a guiding-center Fokker-Planck kinetic equation
journal, September 2013

Differential form of the collision integral for a relativistic plasma
journal, October 1987

A guiding-center Fokker-Planck collision operator for nonuniform magnetic fields
text, January 2004

Electron acceleration in a JET disruption simulation
journal, August 2018

Electron acceleration in a JET disruption simulation
text, January 2018