A fast and accurate numerical implementation of the envelope model for laser–plasma dynamics
Year: 2019
Authors: Terzani D., Londrillo P.
Autors Affiliation: INO CNR Sez Pisa, Via Giuseppe Moruzzi 1, I-56124 Pisa, Italy; INAF, Osservatorio Astron Bologna, Via Gobetti 93-3, Bologna, Italy; Ist Nazl Fis Nucl, Sez Napoli, Naples, Italy; Univ Napoli Federico II, Dipartimento Fis, Naples, Italy; INAF Sez Bologna, Bologna, Italy; Univ Bologna, Dipartimento Fis & Astron, Bologna, Italy.
Abstract: In laser-driven, plasma wakefield acceleration regimes (LWFA), when relevant scale lengths of the laser envelope and of the driven plasma waves are well separated from the wavelength and frequency of the laser fast oscillating component, a reduced physical model (usually referred to as the envelope model), has been introduced, allowing to formulate the laser–plasma equations in terms of laser cycle-averaged dynamical variables. As a main consequence, physical regimes where this reduced model applies, can be investigated with significant savings of computational resources still assuring comparable accuracy, with respect to standard Particle-In-Cell (PIC) models where all relevant space–time scales have to be resolved.
Here we propose a computational framework characterized by two previously unexplored numerical implementations of the envelope model. The first one is based on explicit second order leapfrog integration of the exact wave equation for laser pulse propagation in a laboratory coordinate system in 3D cartesian geometry, replacing the usually quoted representation in an Eulerian frame moving at the speed of light. Since the laser and driven wakefield wave equations in a laboratory frame are advection dominated, we introduce a proper modification of finite differences approximating longitudinal space derivatives, to minimize dispersive numerical errors coming from the discretized advection operators. The proposed implementation, avoiding semi-implicit procedures otherwise required when dealing with a comoving frame, assures significant saving in computational time and ease of implementation for parallel platforms. The associated equation of motion for plasma particles has been integrated, as in standard PIC codes, using the Boris pusher, properly extended to take into account the specific form of the Lorentz force in the envelope model.
As a second contribution, a novel numerical implementation of the plasma dynamics equations in the cold-fluid approximation, is presented. The scheme is based on the second-order one-step Adams–Bashforth time integrator coupled to upwind non-oscillatory WENO reconstruction for discretized space derivatives. The proposed integration scheme for the Eulerian fluid equations is equivalent to a leapfrog scheme with an added higher order dissipative truncation errors. It can be used either as a much faster, yet of comparable accuracy, alternative to the PIC representation of plasma particle motion, or even in a hybrid fluid–particle combination when kinetic effects and particle injection and acceleration in a wakefield have to be investigated.
Journal/Review: COMPUTER PHYSICS COMMUNICATIONS
Volume: 242 Pages from: 49 to: 59
More Information: We would like to thank Carlo Benedetti, for the useful discussions and the precious advices he gave us during the proceeding of this work. We also acknowledge the A. Marocchino grant IsB15_SIPWFA and the P. Londrillo grant IsC65_ENV-LWFA under the ISCRA initiative at CINECA, for the availability of high performance computing resources and support. Also, we acknowledge support from the European Union?s Horizon2020 research and innovation programme under grant agreement No.653782 (EuPRAXIA project), and from the MIUR funded Italian research Network ELI-Italy (D.M.n. 631/16).KeyWords: Laser plasma acceleration, Particle in cell, envelope model, simulationDOI: 10.1016/j.cpc.2019.04.007ImpactFactor: 3.627Citations: 12data from “WEB OF SCIENCE” (of Thomson Reuters) are update at: 2024-10-13References taken from IsiWeb of Knowledge: (subscribers only)Connecting to view paper tab on IsiWeb: Click hereConnecting to view citations from IsiWeb: Click here