An Advanced EM-Plasma simulator based on the DGTD algorithm with dynamic adaptation and multirate time integration techniques

Su Yan, Jiwei Qian, Jian Ming Jin

Research output: Contribution to journalArticle


In a numerical simulation of the interactions between high-intensity electromagnetic (EM) waves and plasma fluids, a good numerical resolution is required for all the physical quantities of interest, especially the plasma density distribution around the edge of the plasma bulk. Since the plasma formation and evolution is a dynamic process, it is difficult to predetermine the appropriate resolution at a given location. In this paper, we describe a discontinuous Galerkin time-domain (DGTD) based dynamic adaptation algorithm, which is able to adjust either the local mesh size or the local polynomial order in real time of a simulation to provide a sufficient numerical resolution of the physics while keeping the total computational cost low. To alleviate the constraint of the time step size of an explicit time integrator, a multirate time integration method is employed to advance the physics in time, which permits the application of different time step sizes in elements with different sizes or polynomial orders. With these techniques, the DGTD simulation of EM-plasma interactions can achieve a good accuracy and high efficiency.

Original languageEnglish (US)
Article number8651513
Pages (from-to)76-87
Number of pages12
JournalIEEE Journal on Multiscale and Multiphysics Computational Techniques
StatePublished - Jan 1 2019



  • Discontinuous Galerkin time-domain (DGTD) method
  • dynamic adaptation
  • electromagnetic-plasma interaction
  • high-order method
  • high-power microwave (HPM)
  • multirate time integration

ASJC Scopus subject areas

  • Modeling and Simulation
  • Mathematical Physics
  • Physics and Astronomy (miscellaneous)
  • Computational Mathematics

Cite this