Computationally efficient methods for modelling laser wakefield acceleration in the blowout regime

B. M. Cowan, S. Y. Kalmykov, A. Beck, X. Davoine, K. Bunkers, A. F. Lifschitz, E. Lefebvre, D. L. Bruhwiler, B. A. Shadwick, Donald P Umstadter

Research output: Contribution to journalArticle

20 Citations (Scopus)

Abstract

Electron self-injection and acceleration until dephasing in the blowout regime is studied for a set of initial conditions typical of recent experiments with 100-terawatt-class lasers. Two different approaches to computationally efficient, fully explicit, 3D particle-in-cell modelling are examined. First, the Cartesian code vorpal (Nieter, C. and Cary, J. R. 2004 VORPAL: a versatile plasma simulation code. J. Comput. Phys. 196, 538) using a perfect-dispersion electromagnetic solver precisely describes the laser pulse and bubble dynamics, taking advantage of coarser resolution in the propagation direction, with a proportionally larger time step. Using third-order splines for macroparticles helps suppress the sampling noise while keeping the usage of computational resources modest. The second way to reduce the simulation load is using reduced-geometry codes. In our case, the quasi-cylindrical code calder-circ (Lifschitz, A. F. et al. 2009 Particle-in-cell modelling of laser-plasma interaction using Fourier decomposition. J. Comput. Phys. 228(5), 1803-1814) uses decomposition of fields and currents into a set of poloidal modes, while the macroparticles move in the Cartesian 3D space. Cylindrical symmetry of the interaction allows using just two modes, reducing the computational load to roughly that of a planar Cartesian simulation while preserving the 3D nature of the interaction. This significant economy of resources allows using fine resolution in the direction of propagation and a small time step, making numerical dispersion vanishingly small, together with a large number of particles per cell, enabling good particle statistics. Quantitative agreement of two simulations indicates that these are free of numerical artefacts. Both approaches thus retrieve the physically correct evolution of the plasma bubble, recovering the intrinsic connection of electron self-injection to the nonlinear optical evolution of the driver.

Original languageEnglish (US)
Pages (from-to)469-482
Number of pages14
JournalJournal of Plasma Physics
Volume78
Issue number4
DOIs
StatePublished - Aug 1 2012

Fingerprint

lasers
resources
simulation
cells
injection
plasma bubbles
decomposition
laser plasma interactions
propagation
economy
splines
preserving
artifacts
electrons
bubbles
sampling
statistics
interactions
electromagnetism
symmetry

ASJC Scopus subject areas

  • Condensed Matter Physics

Cite this

Computationally efficient methods for modelling laser wakefield acceleration in the blowout regime. / Cowan, B. M.; Kalmykov, S. Y.; Beck, A.; Davoine, X.; Bunkers, K.; Lifschitz, A. F.; Lefebvre, E.; Bruhwiler, D. L.; Shadwick, B. A.; Umstadter, Donald P.

In: Journal of Plasma Physics, Vol. 78, No. 4, 01.08.2012, p. 469-482.

Research output: Contribution to journalArticle

Cowan, BM, Kalmykov, SY, Beck, A, Davoine, X, Bunkers, K, Lifschitz, AF, Lefebvre, E, Bruhwiler, DL, Shadwick, BA & Umstadter, DP 2012, 'Computationally efficient methods for modelling laser wakefield acceleration in the blowout regime', Journal of Plasma Physics, vol. 78, no. 4, pp. 469-482. https://doi.org/10.1017/S0022377812000517
Cowan, B. M. ; Kalmykov, S. Y. ; Beck, A. ; Davoine, X. ; Bunkers, K. ; Lifschitz, A. F. ; Lefebvre, E. ; Bruhwiler, D. L. ; Shadwick, B. A. ; Umstadter, Donald P. / Computationally efficient methods for modelling laser wakefield acceleration in the blowout regime. In: Journal of Plasma Physics. 2012 ; Vol. 78, No. 4. pp. 469-482.
@article{8933a7dcb1884cc58c49a044f4b6978b,
title = "Computationally efficient methods for modelling laser wakefield acceleration in the blowout regime",
abstract = "Electron self-injection and acceleration until dephasing in the blowout regime is studied for a set of initial conditions typical of recent experiments with 100-terawatt-class lasers. Two different approaches to computationally efficient, fully explicit, 3D particle-in-cell modelling are examined. First, the Cartesian code vorpal (Nieter, C. and Cary, J. R. 2004 VORPAL: a versatile plasma simulation code. J. Comput. Phys. 196, 538) using a perfect-dispersion electromagnetic solver precisely describes the laser pulse and bubble dynamics, taking advantage of coarser resolution in the propagation direction, with a proportionally larger time step. Using third-order splines for macroparticles helps suppress the sampling noise while keeping the usage of computational resources modest. The second way to reduce the simulation load is using reduced-geometry codes. In our case, the quasi-cylindrical code calder-circ (Lifschitz, A. F. et al. 2009 Particle-in-cell modelling of laser-plasma interaction using Fourier decomposition. J. Comput. Phys. 228(5), 1803-1814) uses decomposition of fields and currents into a set of poloidal modes, while the macroparticles move in the Cartesian 3D space. Cylindrical symmetry of the interaction allows using just two modes, reducing the computational load to roughly that of a planar Cartesian simulation while preserving the 3D nature of the interaction. This significant economy of resources allows using fine resolution in the direction of propagation and a small time step, making numerical dispersion vanishingly small, together with a large number of particles per cell, enabling good particle statistics. Quantitative agreement of two simulations indicates that these are free of numerical artefacts. Both approaches thus retrieve the physically correct evolution of the plasma bubble, recovering the intrinsic connection of electron self-injection to the nonlinear optical evolution of the driver.",
author = "Cowan, {B. M.} and Kalmykov, {S. Y.} and A. Beck and X. Davoine and K. Bunkers and Lifschitz, {A. F.} and E. Lefebvre and Bruhwiler, {D. L.} and Shadwick, {B. A.} and Umstadter, {Donald P}",
year = "2012",
month = "8",
day = "1",
doi = "10.1017/S0022377812000517",
language = "English (US)",
volume = "78",
pages = "469--482",
journal = "Journal of Plasma Physics",
issn = "0022-3778",
publisher = "Cambridge University Press",
number = "4",

}

TY - JOUR

T1 - Computationally efficient methods for modelling laser wakefield acceleration in the blowout regime

AU - Cowan, B. M.

AU - Kalmykov, S. Y.

AU - Beck, A.

AU - Davoine, X.

AU - Bunkers, K.

AU - Lifschitz, A. F.

AU - Lefebvre, E.

AU - Bruhwiler, D. L.

AU - Shadwick, B. A.

AU - Umstadter, Donald P

PY - 2012/8/1

Y1 - 2012/8/1

N2 - Electron self-injection and acceleration until dephasing in the blowout regime is studied for a set of initial conditions typical of recent experiments with 100-terawatt-class lasers. Two different approaches to computationally efficient, fully explicit, 3D particle-in-cell modelling are examined. First, the Cartesian code vorpal (Nieter, C. and Cary, J. R. 2004 VORPAL: a versatile plasma simulation code. J. Comput. Phys. 196, 538) using a perfect-dispersion electromagnetic solver precisely describes the laser pulse and bubble dynamics, taking advantage of coarser resolution in the propagation direction, with a proportionally larger time step. Using third-order splines for macroparticles helps suppress the sampling noise while keeping the usage of computational resources modest. The second way to reduce the simulation load is using reduced-geometry codes. In our case, the quasi-cylindrical code calder-circ (Lifschitz, A. F. et al. 2009 Particle-in-cell modelling of laser-plasma interaction using Fourier decomposition. J. Comput. Phys. 228(5), 1803-1814) uses decomposition of fields and currents into a set of poloidal modes, while the macroparticles move in the Cartesian 3D space. Cylindrical symmetry of the interaction allows using just two modes, reducing the computational load to roughly that of a planar Cartesian simulation while preserving the 3D nature of the interaction. This significant economy of resources allows using fine resolution in the direction of propagation and a small time step, making numerical dispersion vanishingly small, together with a large number of particles per cell, enabling good particle statistics. Quantitative agreement of two simulations indicates that these are free of numerical artefacts. Both approaches thus retrieve the physically correct evolution of the plasma bubble, recovering the intrinsic connection of electron self-injection to the nonlinear optical evolution of the driver.

AB - Electron self-injection and acceleration until dephasing in the blowout regime is studied for a set of initial conditions typical of recent experiments with 100-terawatt-class lasers. Two different approaches to computationally efficient, fully explicit, 3D particle-in-cell modelling are examined. First, the Cartesian code vorpal (Nieter, C. and Cary, J. R. 2004 VORPAL: a versatile plasma simulation code. J. Comput. Phys. 196, 538) using a perfect-dispersion electromagnetic solver precisely describes the laser pulse and bubble dynamics, taking advantage of coarser resolution in the propagation direction, with a proportionally larger time step. Using third-order splines for macroparticles helps suppress the sampling noise while keeping the usage of computational resources modest. The second way to reduce the simulation load is using reduced-geometry codes. In our case, the quasi-cylindrical code calder-circ (Lifschitz, A. F. et al. 2009 Particle-in-cell modelling of laser-plasma interaction using Fourier decomposition. J. Comput. Phys. 228(5), 1803-1814) uses decomposition of fields and currents into a set of poloidal modes, while the macroparticles move in the Cartesian 3D space. Cylindrical symmetry of the interaction allows using just two modes, reducing the computational load to roughly that of a planar Cartesian simulation while preserving the 3D nature of the interaction. This significant economy of resources allows using fine resolution in the direction of propagation and a small time step, making numerical dispersion vanishingly small, together with a large number of particles per cell, enabling good particle statistics. Quantitative agreement of two simulations indicates that these are free of numerical artefacts. Both approaches thus retrieve the physically correct evolution of the plasma bubble, recovering the intrinsic connection of electron self-injection to the nonlinear optical evolution of the driver.

UR - http://www.scopus.com/inward/record.url?scp=84865679195&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84865679195&partnerID=8YFLogxK

U2 - 10.1017/S0022377812000517

DO - 10.1017/S0022377812000517

M3 - Article

VL - 78

SP - 469

EP - 482

JO - Journal of Plasma Physics

JF - Journal of Plasma Physics

SN - 0022-3778

IS - 4

ER -