A New Form of the Numerical Dispersion Relation of ADI-FDTD for the 3D Maxwell's Equations

Author(s):  
Liping Gao
2013 ◽  
Vol 13 (4) ◽  
pp. 1107-1133 ◽  
Author(s):  
Tony W. H. Sheu ◽  
L. Y. Liang ◽  
J. H. Li

AbstractIn this paper an explicit finite-difference time-domain scheme for solving the Maxwell’s equations in non-staggered grids is presented. The proposed scheme for solving the Faraday’s and Ampere’s equations in a theoretical manner is aimed to preserve discrete zero-divergence for the electric and magnetic fields. The inherent local conservation laws in Maxwell’s equations are also preserved discretely all the time using the explicit second-order accurate symplectic partitioned Runge-Kutta scheme. The remaining spatial derivative terms in the semi-discretized Faraday’s and Ampere’s equations are then discretized to provide an accurate mathematical dispersion relation equation that governs the numerical angular frequency and the wavenumbers in two space dimensions. To achieve the goal of getting the best dispersive characteristics, we propose a fourth-order accurate space centered scheme which minimizes the difference between the exact and numerical dispersion relation equations. Through the computational exercises, the proposed dual-preserving solver is computationally demonstrated to be efficient for use to predict the long-term accurate Maxwell’s solutions.


1976 ◽  
Vol 16 (1) ◽  
pp. 47-55 ◽  
Author(s):  
V. Atanassov ◽  
I. Zhelyazkov ◽  
A. Shivarova ◽  
Zh. Genchev

In this paper we propose an exact solution of Vlasov and Maxwell's equations for a bounded hot plasma in order to derive the dispersion relation of the axially-symmetric surface waves propagating along a plasma column. Assuming specular reflexion of plasma particles from the boundary, expressions for the components of the electric displacement vector are obtained on the basis of the Vlasov equation. Their substitution in Maxwell's equations, neglecting the spatial dispersion in the transverse plasma dielectric function, allows us to determine the plasma impedance. The equating of plasma and dielectric impedances gives the wave dispersion relation which, in different limiting cases, coincides with the well-known results.


Geophysics ◽  
1999 ◽  
Vol 64 (5) ◽  
pp. 1369-1377 ◽  
Author(s):  
Tim Bergmann ◽  
Joakim O. Blanch ◽  
Johan O. A. Robertsson ◽  
Klaus Holliger

The Lax‐Wendroff correction is an elegant method for increasing the accuracy and computational efficiency of finite‐difference time‐domain (FDTD) solutions of hyperbolic problems. However, the conventional approach leads to implicit solutions for staggered‐grid FDTD approximations of Maxwell’s equations with frequency‐dependent constitutive parameters. To overcome this problem, we propose an approximation that only retains the purely acoustic, i.e., lossless, terms of the Lax‐Wendroff correction. This modified Lax‐Wendroff correction is applied to an O(2, 4) accurate staggered‐grid FDTD approximation of Maxwell’s equations in the radar frequency range (≈10 MHz–10 GHz). The resulting pseudo-O(4, 4) scheme is explicit and computationally efficient and exhibits all the major numerical characteristics of an O(4, 4) accurate FDTD scheme, even for strongly attenuating and dispersive media. The numerical properties of our approach are constrained by classical numerical dispersion and von Neumann‐Routh stability analyses, verified by comparisons with pertinent 1-D analytical solutions and illustrated through 2-D simulations in a variety of surficial materials. Compared to the O(2, 4) scheme, the pseudo-O(4, 4) scheme requires 64% fewer grid points in two dimensions and 78% in three dimensions to achieve the same level of numerical accuracy, which results in large savings in core memory.


2011 ◽  
Vol 3 (6) ◽  
pp. 663-688 ◽  
Author(s):  
Huajun Zhu ◽  
Songhe Song ◽  
Yaming Chen

AbstractIn this paper, we develop a multi-symplectic wavelet collocation method for three-dimensional (3-D) Maxwell’s equations. For the multi-symplectic formulation of the equations, wavelet collocation method based on autocorrelation functions is applied for spatial discretization and appropriate symplectic scheme is employed for time integration. Theoretical analysis shows that the proposed method is multi-symplectic, unconditionally stable and energy-preserving under periodic boundary conditions. The numerical dispersion relation is investigated. Combined with splitting scheme, an explicit splitting symplectic wavelet collocation method is also constructed. Numerical experiments illustrate that the proposed methods are efficient, have high spatial accuracy and can preserve energy conservation laws exactly.


2017 ◽  
Vol 2017 ◽  
pp. 1-11
Author(s):  
Liping Gao ◽  
Shouhui Zhai

In this paper, we develop a new method to reduce the error in the splitting finite-difference method of Maxwell’s equations. By this method two modified splitting FDTD methods (MS-FDTDI, MS-FDTDII) for the two-dimensional Maxwell equations are proposed. It is shown that the two methods are second-order accurate in time and space and unconditionally stable by Fourier methods. By energy method, it is proved that MS-FDTDI is second-order convergent. By deriving the numerical dispersion (ND) relations, we prove rigorously that MS-FDTDI has less ND errors than the ADI-FDTD method and the ND errors of ADI-FDTD are less than those of MS-FDTDII. Numerical experiments for computing ND errors and simulating a wave guide problem and a scattering problem are carried out and the efficiency of the MS-FDTDI and MS-FDTDII methods is confirmed.


Sign in / Sign up

Export Citation Format

Share Document