High-order spline interpolations in the particle simulation
Journal of Computational Physics
An electromagnetic field algorithm for 2d implicit plasma simulation
Journal of Computational Physics
Electromagnetic direct implicit plasma simulation
Journal of Computational Physics
Performance and optimization of direct implicit particle simulation
Journal of Computational Physics
A second-order implicit particle mover with adjustable damping
Journal of Computational Physics
Computations of the complex error function
SIAM Journal on Numerical Analysis
Numerical recipes in Fortran 90 (2nd ed.): the art of parallel scientific computing
Numerical recipes in Fortran 90 (2nd ed.): the art of parallel scientific computing
Three-dimensional perfectly matched layer for the absorption of electromagnetic waves
Journal of Computational Physics
Plasma Physics Via Computer
An Eulerian gyrokinetic-Maxwell solver
Journal of Computational Physics
STRSCNE: A Scaled Trust-Region Solver for Constrained Nonlinear Equations
Computational Optimization and Applications
On the elimination of numerical Cerenkov radiation in PIC simulations
Journal of Computational Physics
Particle simulations of space weather
Journal of Computational Physics
Hi-index | 31.45 |
Implicit particle-in-cell codes offer advantages over their explicit counterparts in that they suffer weaker stability constraints on the need to resolve the higher frequency modes of the system. This feature may prove particularly valuable for modeling the interaction of high-intensity laser pulses with overcritical plasmas, in the case where the electrostatic modes in the denser regions are of negligible influence on the physical processes under study. To this goal, we have developed the new two-dimensional electromagnetic code ELIXIRS (standing for ELectromagnetic Implicit X-dimensional Iterative Relativistic Solver) based on the relativistic extension of the so-called Direct Implicit Method [D. Hewett, A.B. Langdon, Electromagnetic direct implicit plasma simulation, J. Comput. Phys. 72 (1987) 121-155]. Dissipation-free propagation of light waves into vacuum is achieved by an adjustable-damping electromagnetic solver. In the high-density case where the Debye length is not resolved, satisfactory energy conservation is ensured by the use of high-order weight factors. In this paper, we first derive the electromagnetic direct implicit method as a simplified Newton scheme. Its linear properties are then investigated through numerically solving the relation dispersions obtained for both light and plasma waves, accounting for finite space and time steps. Finally, our code is successfully benchmarked against explicit particle-in-cell simulations for two kinds of physical problems: plasma expansion into vacuum and relativistic laser-plasma interaction. In both cases, we will demonstrate the robustness of the implicit solver for crude discretizations, as well as the gains in efficiency which can be realized over standard explicit simulations.