Recent Achievements on a DGTD Method for Time-Domain Electromagnetics (original) (raw)
Related papers
2009
This work is concerned with the design of a hp-like discontinuous Galerkin (DG) method for solving the two-dimensional time-domain Maxwell equations on non-conforming locally refined triangular meshes. The proposed DG method allows non-conforming meshes with arbitrary-level hanging nodes. This method combines a centered approximation for the evaluation of fluxes at the interface between neighboring elements of the mesh, with a leapfrog time integration scheme. It is an extension of the DG formulation recently studied in [13]. Several numerical results are presented to illustrate the efficiency and the accuracy of the method, but also to discuss its limitations, through a set of 2D propagation problems in homogeneous and heterogeneous media.
Mathematics and Computers in Simulation, 2015
We present a high-order discontinuous Galerkin (DG) method for solving the time-dependent Maxwell equations on non-conforming hybrid meshes. The hybrid mesh combines unstructured tetrahedra for the discretization of irregularly shaped objects with a hexahedral mesh for the rest of the computational domain. The transition between tetrahedra and hexahedra is completely non-conform, that is, no pyramidal or prismatic elements are introduced to link these elements. Within each mesh element, the electromagnetic field components are approximated by a arbitrary order nodal polynomial and a centered approximation is used for the evaluation of numerical fluxes at inter-element boundaries. The time integration of the associated semi-discrete equations is achieved by a fourth-order leapfrog scheme. The method is described and discussed, including algorithm formulation, stability, and practical implementation issues such as the hybrid mesh generation and the computation of flux matrices with cubature rules. We illustrate the performance of the proposed method on several two-and three-dimensional examples involving comparisons with DG methods on single element-type meshes. The results show that the use of non-conforming hybrid meshes in DG methods allows for a notable reduction in computing time without sacrificing accuracy.
Numerical Mathematics: Theory, Methods and Applications, 2009
A high-order leapfrog based non-dissipative discontinuous Galerkin timedomain method for solving Maxwell's equations is introduced and analyzed. The proposed method combines a centered approximation for the evaluation of fluxes at the interface between neighboring elements, with a N th-order leapfrog time scheme. Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable under some CFL-like condition on the time step. The convergence of the semi-discrete approximation to Maxwell's equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with high-order elements show the potential of the method.
IEEE Transactions on Antennas and Propagation
A discontinuous Galerkin time-domain (DGTD) method based on dynamically adaptive Cartesian meshes (ACM) is developed for a full-wave analysis of electromagnetic fields in dispersive media. Hierarchical Cartesian grids offer simplicity close to that of structured grids and the flexibility of unstructured grids while being highly suited for adaptive mesh refinement (AMR). The developed DGTD-ACM achieves a desired accuracy by refining non-conformal meshes near material interfaces to reduce stair-casing errors without sacrificing the high efficiency afforded with uniform Cartesian meshes. Moreover, DGTD-ACM can dynamically refine the mesh to resolve the local variation of the fields during propagation of electromagnetic 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. Simulations of electromagnetic wave diffraction over conducting and dielectric cylinders and spheres demonstrate that the proposed method can achieve a good numerical accuracy at a reduced computational cost compared with uniform meshes. For simulations of dispersive media, the auxiliary differential equation (ADE) and recursive convolution (RC) methods are implemented for a local Drude model and tested for a cold plasma slab and a plasmonic rod. With further advances of the charge transport models, the DGTD-ACM method is expected to provide a powerful tool for computations of electromagnetic fields in complex geometries for applications to high-frequency electronic devices, plasmonic THz technologies, as well as laser-induced and microwave plasmas.
An implicit hybridized discontinuous Galerkin method for the 3D time-domain Maxwell equations
Applied Mathematics and Computation, 2017
We present a time-implicit hybridizable discontinuous Galerkin (HDG) method for numerically solving the system of three-dimensional (3D) time-domain Maxwell equations. This method can be seen as a fully implicit variant of classical so-called DGTD (Discontinuous Galerkin Time-Domain) methods that have been extensively studied during the last 10 years for the simulation of time-domain electromagnetic wave propagation. The proposed method has been implemented for dealing with general 3D problems discretized using unstructured tetrahedral meshes. We provide numerical results aiming at assessing its numerical convergence properties by considering a model problem on one hand, and its performance when applied to more realistic problems. We also include some performance comparisons with a centered flux time-implicit DGTD method.
A high-order non-conforming discontinuous Galerkin method for time-domain electromagnetics
Journal of Computational and Applied Mathematics, 2010
In this paper, we discuss the formulation, stability and validation of a high-order non-dissipative discontinuous Galerkin (DG) method for solving Maxwell's equations on non-conforming simplex meshes. The proposed method combines a centered approximation for the numerical fluxes at inter element boundaries, with either a second-order or a fourth-order leapfrog time integration scheme. Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary-level hanging nodes. The method is proved to be stable and conserves a discrete counterpart of the electromagnetic energy for metallic cavities. Numerical experiments with high-order elements show the potential of the method.
2007
This work is concerned with the design of a hp-like discontinuous Galerkin (DG) method for solving the 2D time-domain Maxwell's equations on non-conforming locally refined triangular meshes. The proposed DG method allows non-conforming meshes with arbitrary-level hanging nodes. This method combines a centered approximation for the evaluation of fluxes at the interface between neighboring elements of the mesh, with a leap-frog time integration scheme. It is an extension of the DG formulation recently studied in . After reviewing the stability properties of the DG method introduced in [18], we present a new implementation of this method which allows local h-refinement and p-enrichment and which is based on a numerical quadrature formula for the computation of the flux matrices associated to non-conforming interfaces. Numerical experiments are presented which both validate the theoretical results of and provide further insights regarding the practical performance of the resulting hp-like DG method, particulary when non-conforming locally refined meshes are employed.
IEEE Transactions on Magnetics, 2000
This paper is concerned with the design of a high-order discontinuous Galerkin (DG) method for solving the 2-D time-domain Maxwell equations on nonconforming triangular meshes. The proposed DG method allows for using nonconforming meshes with arbitrary-level hanging nodes. This method combines a centered approximation for the evaluation of fluxes at the interface between neighboring elements of the mesh, with a leap-frog time integration scheme. Numerical experiments are presented which both validate the theoretical results and provide further insights regarding to the practical performance of the proposed DG method, particulary when nonconforming meshes are employed.
International Journal of Computer Mathematics, 2011
The paper discusses high-order geometrical mapping for handling curvilinear geometries in high accuracy discontinuous Galerkin simulations for time-domain Maxwell problems. The proposed geometrical mapping is based on a quadratic representation of the curved boundary and on the adaptation of the nodal points inside each curved element. With high-order mapping, numerical fluxes along curved boundaries are computed much more accurately due to the accurate representation of the computational domain. Numerical experiments for 2D and 3D propagation problems demonstrate the applicability and benefits of the proposed high-order geometrical mapping for simulations involving curved domains.