Open access peer-reviewed chapter

The Orthogonal Expansion in Time-Domain Method for Solving Maxwell Equations Using Paralleling-in-Order Scheme

By Zheng-Yu Huang, Zheng Sun and Wei He

Submitted: July 11th 2018Reviewed: December 7th 2018Published: January 16th 2019

DOI: 10.5772/intechopen.83387

Downloaded: 172

Abstract

The orthogonal expansion in time-domain method is a new kind of unconditionally stable finite-difference time-domain (FDTD) method for solving the Maxwell equation efficiently. Generally, it can be implemented by two schemes: marching-on-in-order and paralleling-in-order, which, respectively, use weighted Laguerre polynomials and associated Hermite functions as temporal expansions and testing functions. This chapter summarized paralleling-in-order-based FDTD method using associated Hermite functions and Legendre polynomials. And a comparison from theoretical analysis to numerical examples is shown. The LD integral transfer matrix can be considered as a “dual” transformation for AH differential matrix, which gives a possible way to find more potential orthogonal basis function to implement a paralleling-in-order scheme. In addition, the differences with these two orthogonal functions are also analyzed. From the numerical results, we can see their agreements in some general cases while differing in some cases such as shielding analysis with the long-time response requirement.

Keywords

  • associated Hermite
  • finite-difference time-domain (FDTD)
  • Legendre polynomials
  • paralleling-in-order
  • unconditionally stable

1. Introduction

To overcome the numerical stability constraints of conventional finite-difference time-domain (FDTD) method [1, 2], many unconditionally stable methods to reduce or eliminate requirements of the stability condition have been proposed and developed, such as alternating-direction implicit method [2, 3] and locally one-dimensional schemes [3], explicit and unconditionally stable FDTD method [4], and orthogonal expansions in time domain [5, 6, 7, 8]. For the orthogonal expansions schemes, field-versus-time variations in the FDTD space lattice are expanded using an appropriate set of orthogonal temporal basis and testing functions, such as weighted Laguerre polynomials (WLP) and associated Hermite (AH) functions, which leads to two different solution schemes: marching-on-in-order and paralleling-in-order, respectively. Both of them appear to be promising according to the reported work where the computational time can be reduced to at least 10% of the conventional FDTD scheme [1]. Recently, the Legendre (LD) polynomials are explored as another possible orthogonal expansion incorporated with FDTD to form a paralleling-in-order-based unconditionally stable FDTD method. Based on it, in this chapter, we made a comparison investigation for these two new methods, which are AH FDTD method and LD FDTD method, especially focused on their differences. Through a numerical example, we validate their effectiveness when compared with the conventional FDTD method and summarized the characteristics of the two methods.

2. Formulation for paralleling-in-order scheme: AH and LD functions

2.1 2D Maxwell’s equations in time domain

The 2D time-domain Maxwell’s equations with the TEz wave case in lossy medium are considered:

εExrtt+σeExrt=HzrtyJxrtE1
μHzrtt+σmHzrt=ExrtyEyrtxMzrtE2
εEyrtt+σeEyrt=HzrtxJyrtE3

where ε, μ, σe, and σmare the permittivity, the permeability, the electric conductivity, and the magnetic loss of the medium, respectively. Eξrtand Jξrt(ξ=x,y) are the electric field component and the electric current densities, respectively. Hzrtand Mzrtare the magnetic field component and magnetic current densities, respectively.

2.2 The differential and integral transfer matrices to deal with the partial differential term in Maxwell’s equations

2.2.1 The associated Hermite function

Associated Hermite function is defined as

ϕnt=2nn!π1/21/2et2/2Hnt,n=01E4

where Hnt=1net2dndtnet2is Hermite polynomials. Although it is not causal, it can be transformed into causal form by virtue of a proper translating and scaling parameters and then used to span the causal electromagnetic responses. The transformed basis function is ϕ¯nt˜=2nn!σπ1/21/2et˜2/2Hnt˜, where transformed time variable t˜=tTf/σ. And Tfis a translating parameter and σis a scaling parameter. By controlling these two parameters, the time-frequency support of the AH functions ϕ¯nt˜space can be changed flexibly. So, arbitrary locally time-supported functions can be spanned by these transformed basis functions, including the causal electromagnetic responses.

From [7], if a causal function urt, such as the electric or magnetic field function, can be expanded by

urt=n=0unrϕ¯nt˜E5

we can deduce the first derivative of uxtwith respect to

turt=1σn=0un+1rn+12un1rn2ϕ¯nt˜E6

Then, the Q-tuple AH domain coefficients for urtand u̇rtfrom (5) and (6) can be obtained as U=U0UQ1Tand U̇=U̇0U̇Q1T. And, we can readily obtain the relationship between Uand U̇as

U̇=αUE7

where

α=22λ1122Q1Q1Q×QE8

By using (8), the partial differential term in Maxwell?s equations can readily be dealt with, and finally, a five-point banded matrix equation for Hz component can be obtained [9].

2.2.2 The associated Legendre polynomial

We expand all the temporal quantities in terms of the associated Legendre polynomial given by [10]:

Pqt=2q+12Lq2tl1,t0lE9

where l is the time support for analyzing a causal response and Lq is the Legendre polynomial with order q, which are orthogonal in the interval [−1,1] satisfying the following recurrence relation:

Lq+1t=2q+1q+1tLqtqq+1Lq1t,E10

and L0t=0, L1t=t. Given a time-support field function urt, it can be expanded by (9) as

urt=q=0uqrPqtE11

where uqris the q-th expanding coefficients, and it can be calculated by

uqr=+urtPqtdtE12

From the intrinsic features of Legendre function, the differential relationship can be described as

Pqt=122q+32q+1Pq+1't122q+12q1Pq1'tE13

If the field derivative of urtto t is expanded as

u'rt=q=0uq1rPqtE14

where uq1ris q-th expanding coefficients for u'rt, then incorporated with (13), it can be deduced as

u'rt=q=0l22q+12q1uq11rl22q+32q+1uq+11rPqt'E15

Connecting (15) and (11), we can get

uqr=l22q+12q1uq11rl22q+32q+1uq+11rE16

When assembling uqrq=0,1Q1as a Q-tuple Uand uq1rq=0,1Q1as U1, a matrix-multiply relationship can be obtained from (16) as the following:

U=αLU1E17

where αLis integral matrix.

αL=l21/131/131/351/351/2Q32Q11/2Q32Q1Q×QE18

Alternatively, Eq. (17) can be rewritten as.

U1=αL1UE19

2.3 From time domain to orthogonal domain and reconstruction

When the differential or integral transfer matrices are obtained, the time-domain Maxwell equation can be transformed directly into AH or LD domain. Here, let us set LD as an example to illustrate the later formulation.

Similar to the paralleling-in-order-based AH FDTD method, we can apply a Q-tuple-domain transformation for LD FDTD method to (1)–(3) and discretize them as the following:

αeijExi,j=Hzi,jHzi,j1/Δy¯jJxi,jE20
αeijEyi,j=Hzi,jHzi1,j/Δx¯iJyi,jE21
αmijHzi,j=Exi,j+1Exi,j/ΔyjEyi+1,jEyi,j/ΔxiMzi,jE22

where

αeij=εi,jαL1+σei,jIE23
αmij=μmi,jαL1+σmi,jIE24

where Exi,j, Eyi,j, Hzi,j, Jxi,j, Jyi,j, and Mzi,jare Q-tuple representations of fields and sources, respectively. And, Iis the Q-dimensional identity matrix. By assembling (20)–(22) and eliminating the electric field components, a five-diagonal banded matrix equation for Hz component can be obtained:

alijHzi1,j+ari+1jHzi+1,j+amijHzi,j+adijHzi,j1+auij+1Hzi,j+1=bi,jE25

where

auij+1=αeij+11/Δy¯j+1/ΔyjE26
adij=αeij1/Δy¯j/ΔyjE27
alij=αeij1/Δx¯i/ΔxiE28
ari+1j=αei+1j1/Δx¯i+1/ΔxiE29
amij=ari+1j+alij+auij+1+adij+αmijE30
bi,j=αmij+11Jxi,j+1αmij1Jxi,j/Δy+αmi+1j1Jyi+1,jαmij1Jyi,j/ΔxMzi,jE31

By using eigenvalue transformation from αLX=XV, where Xand Vare the eigenvector matrix and diagonal matrix composed of eigenvalues λq, respectively, Eq. (25) can be changed to the paralleling-in-order solution. For the q-th decoupled equation, we have

A1/λqHzq=bqE32

where A·is a banded sparse matrix, with the similar form as from AH FDTD method, and bqis the transformed variables from bi,j=Xbi,j. Finally, we can obtain a paralleling-in-order scheme to calculate all of the expanding coefficients of electromagnetic fields, and then the time-domain responses can be reconstructed from (11).

3. Comparison for the two methods

The above formula can be regarded and classified as a uniform OF differential transfer matrix transformation. Therefore, as long as the LD differential matrix is replaced by the AH domain differential transfer matrix, the FDTD algorithm based on the LD orthogonal basis function, LD FDTD, including the parallel solution AH FDTD algorithm [9], and the alternate direction efficient calculation [11] can be easily realized. The implementation of the program only requires a simple modification.

Table 1 gives a comparison of the relevant properties of the LD FDTD method and the AH FDTD method. It can be seen that the two methods can be considered as a “dual” system, because the AH differential matrix is the basic element of the AH FDTD method and the LD integration matrix is also the basic element of the LD FDTD method. This gives us a revelation that is it possible that any orthogonal basis function can construct a differential or integral transfer matrix and then easily implement a paralleling-in-order scheme similar like AH FDTD algorithm? The answer might be NOT. Such as the Laguerre FDTD method, as introduced before, cannot be calculated in parallel. However, it is undeniable that there may be more basis functions that can implement the paralleling-in-order scheme. If any, we can collectively call these methods as the AH series unconditionally stable FDTD method.

AH FDTDLD FDTD
Differential transfer matrixIntegral transfer matrix
U1=αUU=αLU1
αl=22l1122Q1Q1Q×QαLl=l21/131/131/351/351/2Q32Q11/2Q32Q1Q×Q
tαdtαL1
TQ2lπQ/1.7+1.8FQπQ/1.7+1.82πllQScale factor l = TQ
Finite order of Q
With time-frequency HomomorphismWithout time-frequency
Homomorphism
Antisymmetry
Eigenvalue conjugate symmetry
Antisymmetry
Eigenvalue conjugate symmetry
AλqHq=JqA1/λqHq=Jq

Table 1.

LD comparison of LD FDTD method and AH FDTD method.

4. Numerical verification

4.1 An infinitely large lossy dielectric plate

As AH or LD FDTD method shares with almost the same program, a 1-D program is set for a general verification. Figure 1 shows the simulation results when a uniform plane wave penetrates an infinitely large lossy dielectric plate. The figure includes the electric field waveforms calculated by the AH FDTD method and the LD FDTD method and their relative errors with respect to the conventional FDTD method. It can be seen that the time-domain waveforms of both can be consistent with the results of the FDTD method and the relative errors are basically the same, only differing in the initial part. Therefore, in general, when the order of the two basic functions is the same and the parameters are selected reasonably, the accuracy is basically the same, and the efficiency is almost the same.

Figure 1.

Comparison of calculation results between AH FDTD method and HR FDTD method when simulating an infinitely large lossy dielectric plate. (a) Time-domain waveform. (b) Relative error.

4.2 An nonuniform parallel plate waveguide with a slot

However, the two methods also have the differences when simulating the long-time response applications, such as the example in [12]. The numerical example is set as a TEz wave propagation in a parallel plate waveguide, as shown in Figure 2. It is with a PEC slot of the thickness 0.2 mm and the distance 0.2 mm and a partly filled dielectric material of the thickness 0.8 mm with the dielectric medium parameters given as two cases: case I, ε = 11 ε0, μ = μ0, σe = 0.003 S/m, and σm = 0 Ω/m; case II, ε = 2 ε0, μ = μ0, σe = 30,000 S/m, and σm = 0 Ω/m. There are 140 × 8 uniform cells (Δx=Δy=0.1 mm) in the computational domain. A Gaussian pulse sinusoidally modulated is used as the electric current source profile:

Jyt=expttc/td2sin2πfcttcE33

Figure 2.

The geometry configuration for a 2D parallel plate waveguide with a PEC slot and a partly filled dielectric medium [12].

where td= 1/2fc, tc= 4td, and fc= 12 GHz. And the total simulation time is set as l= 1.28 ns for case I and l= 12.8 ns for case II; then it leads to the marching-in-on-time steps for N = 6000 and N = 60,000, respectively. And the number of orders for LD functions is chosen as 80 and 300, respectively, to obtain a good approximation of field components.

The Ey electric field responses at measurement point p1 and p2, located at the center of the slot and behind the medium, respectively, are calculated, which are both in agreement with the conventional FDTD method as shown in Figures 2 and 3. For comparison, the AH FDTD method is also used in these two cases. One can find the good results in Figure 3, but the errors come out in Figure 4 for AH FDTD method when the same number of orthogonal functions (Q = 80 for case I or 300 for case II) is used as LD FDTD method. However, when Q reaches 800, the results from AH FDTD method can achieve a comparable accuracy with the ones from LD FDTD method. One should note that for case II the waveform at point p2 has larger amplitude attenuation and longer delay than the result at point p1 due to the high dielectric medium located between them.

Figure 3.

The calculated results of transient electric field Ey for the case of I [12].

Figure 4.

The calculated results of transient electric field Ey for the case of II [12].

Tables 2 and 3 show the comparison of the computational resources. We can see that the simulation takes much more time for the FDTD method compared with proposed method, especially for the case of II, while the trade-off for the proposed method is that it consumes more memory than conventional FDTD method, which is similar to the AH FDTD method. In addition, from Table 3, we can find the advantages compared with AH FDTD method that the proposed method can use relative smaller memory storage and slightly fewer CPU times to get a readily results.

Δt (ps)Memory (MB)CPU time (s)
FDTD (N = 6000)0.211.82.97
AH FDTD (Q = 80)212.91.32
LD FDTD (Q = 80)212.91.32

Table 2.

The comparison of computational resources for the case of I [12].

Δt (ps)Memory (MB)CPU time (s)
FDTD (N = 60,000)0.211.830.8
AH FDTD (Q = 300)2111.81.55
AH FDTD (Q = 800)2128.91.95
LD FDTD (Q = 300)2111.81.55

Table 3.

The comparison of computational resources for the case of II [12].

5. Conclusions and future developments

The paralleling-in-order-based unconditionally stable FDTD methods are introduced using associated Hermite and Legendre polynomials in this chapter. The direct Q-tuple-domain transformation for time-domain Maxwell equation is guaranteed by using the integral matrix and differential matrix for Legendre function and associated Hermite functions that are introduced from the intrinsic integral or differential features for these orthogonal functions. Normally, the integral matrix of Legendre function can be considered as an inverse relationship from the differential operator, similar to the AH differential matrix. From this view, we can consider them as a uniform algorithm organized from the paralleling-in-order solution scheme. In addition, this chapter also detailed the different properties and the formula with these two methods theoretically and tested by numerical examples. Numerical examples for 1D and 2D cases validate their effectiveness and show LD FDTD with a better performance than AH FDTD method, in long-time simulation applications. In the next step, the more general paralleling-in-order scheme should be summarized, and then find or construct other possible orthogonal functions for their specific applications.

Acknowledgments

This work is supported by the National Natural Science Foundation of China under Grants 61801217 and 51477183 and Natural Science Foundation of Jiangsu Province under Grant BK20180422. This support is gratefully acknowledged.

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Zheng-Yu Huang, Zheng Sun and Wei He (January 16th 2019). The Orthogonal Expansion in Time-Domain Method for Solving Maxwell Equations Using Paralleling-in-Order Scheme, Polynomials - Theory and Application, Cheon Seoung Ryoo, IntechOpen, DOI: 10.5772/intechopen.83387. Available from:

chapter statistics

172total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Cyclotomic and Littlewood Polynomials Associated to Algebras

By José-Antonio de la Peña

Related Book

Frontiers in Guided Wave Optics and Optoelectronics

Edited by Bishnu Pal

First chapter

Frontiers in Guided Wave Optics and Optoelectronics

By Bishnu Pal

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us