Jacobian-Free Newton-Krylov Method for the Simulation of Non-Thermal Plasma Discharges with High-Order Time Integration and Physics-Based Preconditioning: Article No. 112007

Alfredo Gomez, Nicholas Deak, Fabrizio Bisetti

Research output: Contribution to journalArticlepeer-review

6 Scopus Citations

Abstract

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. The plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.
Original languageAmerican English
Number of pages31
JournalJournal of Computational Physics
Volume480
DOIs
StatePublished - 2023

NREL Publication Number

  • NREL/JA-2C00-89515

Keywords

  • backward differentiation formulas
  • Jacobian-free Newton-Krylov
  • non-thermal plasma
  • PETSc
  • physics-based preconditioning
  • streamer discharge

Fingerprint

Dive into the research topics of 'Jacobian-Free Newton-Krylov Method for the Simulation of Non-Thermal Plasma Discharges with High-Order Time Integration and Physics-Based Preconditioning: Article No. 112007'. Together they form a unique fingerprint.

Cite this