An Envelope Tracking Approach for Particle in Cell Simulations

—The state of the art in electromagnetic Finite Ele- ment Particle-in-Cell (EM-FEMPIC) has advanced significantly in the last few years; these have included understanding function spaces that must be used to represent sources and fields consistently, and how currents should be evolved in space and time. In concert, these achieve satisfaction of Gauss’ laws. All of these, were restricted to conditionally stable explicit time stepping. More recently, there has been advances to the state of art: It is now possible to use a implicit EM-FEMPIC method while satisfying Gauss’ law to machine precision. This enables choosing time step sizes dictated by physics as opposed to geometry. In this paper, we take this a step further. For devices characterized by a narrowband high frequency response, choosing a time-step size based on the highest frequency of interest is considerably expensive. In this paper, we use methods derived from envelope tracking to construct an EM-FEMPIC method that analytically provides for the high-frequency oscillations of the system, allow- ing for analysis at considerable coarser time-step sizes even in the presence of non-linear effects from active media such as plasmas. Consequentially, we demonstrate how the pointwise metric used for measuring satisfaction of Gauss’ Laws breaks down when prescribing analytical fast fields and provide a thorough analysis of how charge conservation can be measured. Through a number of examples, we demonstrate that the proposed approach retains the accuracy the regular scheme while requiring far fewer time steps.

[1]  O. H. Ramachandran,et al.  Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell , 2021, Comput. Phys. Commun..

[2]  B. Shanker,et al.  Time integrator agnostic charge conserving finite element PIC , 2021, Physics of Plasmas.

[3]  B. Shanker,et al.  A Set of Benchmark Tests for Validation of 3-D Particle in Cell Methods , 2021, IEEE Transactions on Plasma Science.

[4]  Jianyuan Xiao,et al.  Structure-preserving geometric particle-in-cell methods for Vlasov-Maxwell systems , 2018, Plasma Science and Technology.

[5]  K. Kormann,et al.  GEMPIC: geometric electromagnetic particle-in-cell methods , 2016, Journal of Plasma Physics.

[6]  Martin Campos Pinto,et al.  Charge-conserving FEM–PIC schemes on general grids☆ , 2014 .

[7]  Fernando L. Teixeira,et al.  Exact charge-conserving scatter-gather algorithm for particle-in-cell simulations on unstructured grids: A geometric perspective , 2014, Comput. Phys. Commun..

[8]  Hong Qin,et al.  Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme , 2012, 1401.6723.

[9]  R. Marchand,et al.  PTetra, a Tool to Simulate Low Orbit Satellite–Plasma Interaction , 2012, IEEE Transactions on Plasma Science.

[10]  G. Sangalli,et al.  Isogeometric analysis in electromagnetics: B-splines approximation , 2010 .

[11]  J. Verboncoeur Particle simulation of plasmas: review and advances , 2005 .

[12]  F. Teixeira,et al.  Geometric finite element discretization of Maxwell equations in primal and dual spaces , 2005, physics/0503013.

[13]  B Shahine,et al.  Particle in cell simulation of laser-accelerated proton beams for radiation therapy. , 2002, Medical physics.

[14]  David Chernin,et al.  A simulation study of beam loading in a cavity , 2002, Third IEEE International Vacuum Electronics Conference (IEEE Cat. No.02EX524).

[15]  R. Lemke,et al.  Three-dimensional particle-in-cell simulation study of a relativistic magnetron , 1999 .

[16]  O. Picon,et al.  A finite element method based on Whitney forms to solve Maxwell equations in the time domain , 1995 .

[17]  A. Bossavit Whitney forms: a class of finite elements for three-dimensional computations in electromagnetism , 1988 .

[18]  O. C. Zienkiewicz A new look at the newmark, houbolt and other time stepping formulas. A weighted residual approach , 1977 .

[19]  Albert E. Ruehli,et al.  The modified nodal approach to network analysis , 1975 .

[20]  B. Shanker,et al.  UNCONDITIONALLY STABLE TIME STEPPING METHOD FOR MIXED FINITE ELEMENT MAXWELL SOLVERS , 2020, Progress In Electromagnetics Research C.

[21]  L. Kettunen,et al.  Yee‐like schemes on a tetrahedral mesh, with diagonal lumping , 1999 .