A Discontinuous Galerkin Time-Domain Method With Dynamically Adaptive Cartesian Mesh for Computational Electromagnetics

A discontinuous Galerkin time-domain (DGTD) method based on dynamically adaptive Cartesian mesh (ACM) is developed for a full-wave analysis of electromagnetic (EM) fields. The benefits of hierarchical Cartesian grids and adaptive mesh refinement are demonstrated for linear EM propagation problems. The developed DGTD-ACM achieves a desired accuracy by refining nonconformal meshes near material interfaces to reduce stair-casing errors without sacrificing the high efficiency afforded with uniform Cartesian meshes. More importantly, DGTD-ACM can dynamically refine the mesh to resolve the local variation of the fields during propagation of EM pulses. A local time-stepping scheme is adopted to alleviate the constraint on the time-step size due to the stability condition of the explicit time integration. It is shown by numerical examples that the proposed method can achieve a good numerical accuracy and reduce the computational time effectively for linear problems of EM propagation in dispersive media. With further development, the method is expected to provide a powerful tool for solving nonlinear EM problems in plasma physics and electronics.

[1]  Norinobu Yoshida,et al.  A Treatment by the Finite-Difference Time-Domain Method of the Dispersive Characteristics Associated with Orientation Polarization , 1990 .

[2]  Vladimir Kolobov,et al.  Towards adaptive kinetic-fluid simulations of weakly ionized plasmas , 2012, J. Comput. Phys..

[3]  Vladimir I. Kolobov,et al.  Analysis of Compressible Viscous Flow Solvers with Adaptive Cartesian Mesh , 2011 .

[4]  Jan S. Hesthaven,et al.  High-order accurate methods in time-domain computational electromagnetics: A review , 2003 .

[5]  Jianming Jin,et al.  A continuity-preserving and divergence-cleaning algorithm based on purely and damped hyperbolic Maxwell equations in inhomogeneous media , 2017, J. Comput. Phys..

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

[7]  Claus-Dieter Munz,et al.  A Finite-Volume Method for the Maxwell Equations in the Time Domain , 2000, SIAM J. Sci. Comput..

[8]  S. Popinet.  Gerris: a tree-based adaptive solver for the incompressible Euler equations in complex geometries , 2003 .

[9]  J. Hesthaven,et al.  Nodal high-order methods on unstructured grids , 2002 .

[10]  Stephen D. Gedney,et al.  The Discontinuous Galerkin Finite-Element Time-Domain Method Solution of Maxwell ’ s Equations , 2009 .

[11]  Thomas Weiland,et al.  Non-dissipative space-time hp-discontinuous Galerkin method for the time-dependent Maxwell equations , 2013, J. Comput. Phys..

[12]  Jian-Ming Jin,et al.  Modeling of Plasma Formation During High-Power Microwave Breakdown in Air Using the Discontinuous Galerkin Time-Domain Method , 2016, IEEE Journal on Multiscale and Multiphysics Computational Techniques.

[13]  Pedro Velarde,et al.  Development of a Godunov method for Maxwell's equations with Adaptive Mesh Refinement , 2015, J. Comput. Phys..

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

[15]  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..

[16]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[17]  Su Yan,et al.  Waveport modeling for the DGTD simulation of electromagnetic devices , 2018 .

[18]  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.

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

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

[21]  Jianming Jin.  The Finite Element Method , 2010 .

[22]  Pierre Degond,et al.  On a finite-element method for solving the three-dimensional Maxwell equations , 1993 .

[23]  J. Nédélec.  Mixed finite elements in ℝ3 , 1980 .

[24]  J. P. Webb.  Hierarchal vector basis functions of arbitrary order for triangular and tetrahedral finite elements , 1999 .

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

[26]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[27]  Timothy C. Warburton,et al.  Nodal discontinuous Galerkin methods on graphics processors , 2009, J. Comput. Phys..

[28]  Nathan M. Newmark,et al.  A Method of Computation for Structural Dynamics , 1959 .

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

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

[31]  Claus-Dieter Munz,et al.  A three-dimensional finite-volume solver for the Maxwell equations with divergence cleaning on unstructured meshes , 2000 .

[32]  M. F. Pantoja,et al.  Discontinuous Galerkin Time Domain Methods in Computational Electrodynamics : State of the Art , 2015 .

[33]  Claus-Dieter Munz,et al.  Hyperbolic divergence cleaning, the electrostatic limit, and potential boundary conditions for particle-in-cell codes , 2015, J. Comput. Phys..

[34]  Jian-Ming Jin,et al.  A Dynamic $p$ -Adaptive DGTD Algorithm for Electromagnetic and Multiphysics Simulations , 2017, IEEE Transactions on Antennas and Propagation.

[35]  C. Yuan,et al.  FDTD Analysis of Terahertz Wave Propagation in a High-Temperature Unmagnetized Plasma Slab , 2011, IEEE Transactions on Plasma Science.

[36]  Costas D. Sarris,et al.  Adaptive Mesh Refinement for Time-Domain Numerical Electromagnetics , 2006, Adaptive Mesh Refinement for Time-Domain Numerical Electromagnetics.

[37]  Jian-Ming Jin,et al.  Three-dimensional orthogonal vector basis functions for time-domain finite element solution of vector wave equations , 2003 .

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

[39]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[40]  Stéphane Descombes,et al.  Recent advances on a DGTD method for time-domain electromagnetics , 2013 .

[41]  D. Vasileska,et al.  Computational Electronics: Semiclassical and Quantum Device Modeling and Simulation , 2010 .

[42]  Li Jun Jiang,et al.  Integration of Arbitrary Lumped Multiport Circuit Networks Into the Discontinuous Galerkin Time-Domain Analysis , 2013, IEEE Transactions on Microwave Theory and Techniques.

[43]  Li Liu,et al.  Nonuniform time-step Runge-Kutta discontinuous Galerkin method for Computational Aeroacoustics , 2010, J. Comput. Phys..

[44]  U. Navsariwala,et al.  An unconditionally stable finite element time-domain solution of the vector wave equation , 1995 .

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

[46]  Zhiqiang Tan,et al.  CFD-Micromesh: a fast geometric modeling and mesh generation tool for 3D microsystem simulations , 2000, Design, Test, Integration, and Packaging of MEMS/MOEMS.

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

[48]  Vladimir Kolobov,et al.  Adaptive kinetic-fluid solvers for heterogeneous computing architectures , 2015, J. Comput. Phys..

[49]  Claus-Dieter Munz,et al.  Divergence Correction Techniques for Maxwell Solvers Based on a Hyperbolic Model , 2000 .