A novel discontinuous Galerkin time-domain (DGTD) method called the shift operator DGTD (SO-DGTD) method in the 3-D case is proposed for the simulation of electromagnetic (EM) waves in a magnetized plasma. In this paper, the electric permittivity tensor of the magnetized plasma is written as a rational matrix polynomial function; a time-domain iteration formula between the electric displacement field D and the electric field E is obtained by the shift operator (SO) method. The proposed algorithm is validated by a simulation of a magnetoplasma-coated sphere and a cavity filled with a plasma. Lastly, the transmission coefficients of magnetized plasma periodic slabs are analyzed.

A discontinuous Galerkin time-domain (DGTD) method based on the current density equation in a 3-D case is proposed for the simulation of electromagnetic (EM) waves in magnetized plasma. In this article, two paths are given for the derivation of time-domain iteration formulas, auxiliary differential equation method for the current density equation (ADE-JE) and LT-JE. The implementation is explained in detail. The proposed algorithm avoids additional data storage and has a relatively high calculation speed, which makes it a powerful tool for magnetized plasma. The algorithm is validated by numerical examples, such as the Faraday rotation effect test and time-varying example.

The plasma sheath is known as a popular topic of computational electromagnetics, and the plasma case is more resource-intensive than the non-plasma case. In this paper, a parallel shift-operator discontinuous Galerkin time-domain method using the MPI (Message Passing Interface) library is proposed to solve the large-scale plasma problems. To demonstrate our algorithm, a plasma sheath model of the high-speed blunt cone was established based on the results of the multiphysics software, and our algorithm was used to extract the radar cross-section (RCS) versus different incident angles of the model.

Using the numerical discrete technique with unstructured grids, conformal perfectly matched layer (PML) absorbing boundary in the discontinuous Galerkin time-domain (DGTD) can be set flexibly so as to save lots of computing resources. Based on the DGTD equations in an orthogonal curvilinear coordinate system, the processes of parameter transformation for 2-D UPML between the coordinate systems of elliptical and Cartesian are given; and the expressions of transition matrix are derived. The calculation scheme of conductivity distribution in elliptic cylinder absorbing layer is given, and the calculation coefficient of DGTD in elliptic UPML is calculated. Furthermore, the 2-D iterative formulas of DGTD and that of auxiliary equation in the elliptical cylinder UPML are derived; the conformal UPML calculation in DGTD is realized. Numerical results show that very good accuracy and computational efficiency are achieved by using the method in this paper. Compared to the rectangular computational region, both the memory and computation time of conformal UPML absorbing boundary are reduced by more than 20%.