Modeling of Plasma Formation During High-Power Microwave Breakdown in Air Using the Discontinuous Galerkin Time-Domain Method

Rapid plasma formation and evolution during high-power microwave (HPM) air breakdown in an HPM device produce a macroscopic plasma shield to the microwave transmission, which can severely limit the performance of the device. In this paper, the electromagnetic (EM)–plasma interaction and the HPM breakdown in air are modeled by a nonlinearly coupled full-wave Maxwell and plasma fluid system under conditions of high pressure and high collision frequency. The resulting multiphysics and multiscale system is solved by a nodal discontinuous Galerkin time-domain (DGTD) method, which is uniformly high order in both space and time. To demonstrate the capability of the DGTD method in the modeling of the HPM breakdown problems, the air breakdown and plasma formation in a metallic aperture are simulated, from which the underlining physical process can be interpreted and better understood.

[1]  George Em Karniadakis,et al.  The Development of Discontinuous Galerkin Methods , 2000 .

[2]  Jian-Ming Jin,et al.  Fast and Efficient Algorithms in Computational Electromagnetics , 2001 .

[3]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[4]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[5]  Jian-Ming Jin,et al.  A GPU accelerated dynamic p-adaptation for simulation of EM-plasma interaction , 2016, 2016 IEEE International Symposium on Antennas and Propagation (APSURSI).

[6]  J. Maxwell A Treatise on Electricity and Magnetism , 1873, Nature.

[7]  Li Jun Jiang,et al.  A Hybrid Electromagnetics-Circuit Simulation Method Exploiting Discontinuous Galerkin Finite Element Time Domain Method , 2013, IEEE Microwave and Wireless Components Letters.

[8]  J. Butcher Numerical Methods for Ordinary Differential Equations: Butcher/Numerical Methods , 2005 .

[9]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[10]  Christian A. Ringhofer,et al.  Moment Methods for the Semiconductor Boltzmann Equation on Bounded Position Domains , 2001, SIAM J. Numer. Anal..

[11]  Jian-Ming Jin,et al.  An OpenMP-CUDA Implementation of Multilevel Fast Multipole Algorithm for Electromagnetic Simulation on Multi-GPU Computing Systems , 2013, IEEE Transactions on Antennas and Propagation.

[12]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[13]  K. Blotekjaer Transport equations for electrons in two-valley semiconductors , 1970 .

[14]  Li Jun Jiang,et al.  A Hybrid Time-Domain Discontinuous Galerkin-Boundary Integral Method for Electromagnetic Scattering Analysis , 2014, IEEE Transactions on Antennas and Propagation.

[15]  Carretera de Valencia,et al.  The finite element method in electromagnetics , 2000 .

[16]  François Rogier,et al.  ADI-FDTD modeling of microwave plasma discharges in air towards fully three-dimensional simulations , 2015, Comput. Phys. Commun..

[17]  Jean-Pierre Boeuf,et al.  Computational Studies of Filamentary Pattern Formation in a High Power Microwave Breakdown Generated Air Plasma , 2010, IEEE Transactions on Plasma Science.

[18]  L. Warne,et al.  Ionization coefficient approach to modeling breakdown in nonuniform geometries. , 2003 .

[19]  K. Kunz,et al.  A frequency-dependent finite-difference time-domain formulation for transient propagation in plasma , 1991 .

[20]  Qing Huo Liu,et al.  A Nonspurious 3-D Vector Discontinuous Galerkin Finite-Element Time-Domain Method , 2010, IEEE Microwave and Wireless Components Letters.

[21]  W. Woo,et al.  Microwave absorption and plasma heating due to microwave breakdown in the atmosphere , 1984 .

[22]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[23]  S. Cummer,et al.  An FDTD model for low and high altitude lightning-generated EM fields , 2006, IEEE Transactions on Antennas and Propagation.

[24]  J. Oden,et al.  A discontinuous hp finite element method for convection—diffusion problems , 1999 .

[25]  Paul Bellan,et al.  Fundamentals of Plasma Physics , 2006 .

[26]  Stephen C. Jardin,et al.  Computational Methods in Plasma Physics , 2010 .

[27]  Jian-Ming Jin,et al.  An Accurate and Efficient Finite Element-Boundary Integral Method With GPU Acceleration for 3-D Electromagnetic Analysis , 2014, IEEE Transactions on Antennas and Propagation.

[28]  E. Toro,et al.  An arbitrary high-order Discontinuous Galerkin method for elastic waves on unstructured meshes - V. Local time stepping and p-adaptivity , 2007 .

[29]  Jian-Ming Jin,et al.  A DGTD algorithm with dynamic h-adaptation and local time-stepping for solving Maxwell's equations , 2016, 2016 IEEE International Symposium on Antennas and Propagation (APSURSI).

[30]  Qing Huo Liu,et al.  A new efficient 3D Discontinuous Galerkin Time Domain (DGTD) method for large and multiscale electromagnetic simulations , 2015, J. Comput. Phys..

[31]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems , 1989 .

[32]  Jean-Pierre Boeuf,et al.  Pattern formation and propagation during microwave breakdown , 2010 .

[33]  Pingwen Zhang,et al.  Discontinuous Galerkin methods for dispersive and lossy Maxwell's equations and PML boundary conditions , 2004 .

[34]  B. H. Crichton,et al.  Gas discharge physics , 1996 .

[35]  Baruch Levush,et al.  MAGY: a time-dependent code for simulation of slow and fast microwave sources , 1998 .

[36]  W. V. Roosbroeck Theory of the flow of electrons and holes in germanium and other semiconductors , 1950 .

[37]  David I. Ketcheson,et al.  Strong stability preserving runge-kutta and multistep time discretizations , 2011 .

[38]  Jianming Jin Theory and Computation of Electromagnetic Fields , 2010 .

[39]  Jeffrey L. Young,et al.  A full finite difference time domain implementation for radio wave propagation in a plasma , 1994 .

[40]  J. Callen Fundamentals of Plasma Physics , 2004 .

[41]  Jean-Pierre Boeuf,et al.  Ionization–diffusion plasma front propagation in a microwave field , 2011 .

[42]  Jean-Pierre Boeuf,et al.  Theory and modeling of self-organization and propagation of filamentary plasma arrays in microwave breakdown at atmospheric pressure. , 2010, Physical review letters.

[43]  Siegfried Selberherr,et al.  Using six moments of Boltzmann’s transport equation for device simulation , 2001 .

[44]  J. Maxwell VIII. A dynamical theory of the electromagnetic field , 1865, Philosophical Transactions of the Royal Society of London.

[45]  Chi-Wang Shu Total-variation-diminishing time discretizations , 1988 .

[46]  Jean-Pierre Boeuf,et al.  Finite Volume Time Domain modelling of microwave breakdown and plasma formation in a metallic aperture , 2012, Comput. Phys. Commun..

[47]  S. Selberherr,et al.  A review of hydrodynamic and energy-transport models for semiconductor device simulation , 2003, Proc. IEEE.

[48]  Mengping Zhang,et al.  AN ANALYSIS OF THREE DIFFERENT FORMULATIONS OF THE DISCONTINUOUS GALERKIN METHOD FOR DIFFUSION EQUATIONS , 2003 .

[49]  J. Lambert Numerical Methods for Ordinary Differential Equations , 1991 .

[50]  Qing Huo Liu,et al.  Three‐dimensional unstructured‐grid discontinuous Galerkin method for Maxwell's equations with well‐posed perfectly matched layer , 2005 .