An Improved Subdomain Level Nonconformal Discontinuous Galerkin Time Domain (DGTD) Method for Materials With Full-Tensor Constitutive Parameters

Ultrawideband simulation of negative refraction in bicrystals is important for the design optimization of devices involving such anisotropic media, but it is a nontrivial task, especially when low-order methods are utilized. This work proposes an improved discontinuous Galerkin time domain (DGTD) me...

Full description

Saved in:
Bibliographic Details
Main Authors: Qiang Ren, Qiwei Zhan, Qing Huo Liu
Format: Article
Language:English
Published: IEEE 2017-01-01
Series:IEEE Photonics Journal
Subjects:
Online Access:https://ieeexplore.ieee.org/document/7862232/
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Ultrawideband simulation of negative refraction in bicrystals is important for the design optimization of devices involving such anisotropic media, but it is a nontrivial task, especially when low-order methods are utilized. This work proposes an improved discontinuous Galerkin time domain (DGTD) method for simulating time-dependent electromagnetic fields for inhomogeneous media with full anisotropic constitutive parameters (full anisotropic media). It employs the electric field intensity <inline-formula><tex-math notation="LaTeX">${\mathbf E}$</tex-math> </inline-formula> and magnetic flux density <inline-formula><tex-math notation="LaTeX">${\mathbf B}$</tex-math> </inline-formula> to solve Maxwell&#x0027;s equations. The <inline-formula><tex-math notation="LaTeX">${\mathbf EB}$ </tex-math></inline-formula>-scheme-based anisotropic Riemann solver and nonconformal mesh are employed for domain decomposition to allow efficient spatial discretization. An unsplit-field Maxwellian multiaxial perfectly matched layer for full anisotropic media is derived and shown to be effective to absorb outgoing waves and suppress the potential late-time instability found in classical PML. In addition, the total-field/scattered-field technique is further studied to allow a nonconformal mesh, vector basis functions, and half-space situation. This newly improved DGTD method is validated with test cases and applied to the negative reflection in YVO<inline-formula> <tex-math notation="LaTeX">$_{4}$</tex-math></inline-formula> bicrystal.
ISSN:1943-0655