# Implementation of an Approximate Conformal UPML in 2-D DGTD.

1. IntroductionTruncated boundary condition is the key to ensure the accuracy for many electromagnetic numerical methods. In the past decades, different kinds of absorbing boundary conditions (ABC) have been proposed and successfully applied to the finite difference time-domain (FDTD) method and other methods. The perfectly matched layer (PML) ABC is presented by Berenger in 1994 [1]. Later, the coordinate stretched perfectly matched layer (CPML) ABC was proposed by Chew and Weedon [2], and uniaxial anisotropic perfectly matched layer (UPML) ABC was presented by Sacks et al. [3]. The conformal PML was analytically derived and applied to the hyperbolic model by Teixeira et al. in the 1990s [4-6]. The corner-free truncation strategy was proposed and used in FDTD method by Zhang et al. [7]. The discontinuous Galerkin time-domain (DGTD) [8-12] method, which was developed based on the finite volume element time-domain method, has advantages of mesh flexibility of finite element time-domain (FETD) and explicit iterative of FDTD. As a new kind of time-domain algorithm, the study of ABC for DGTD is a hot topic in the computational electromagnetic field.

The first-order Silver-Muller (SM) ABC is widely used owing to its easy realization [13]. However, this kind of ABC only has a better absorption when the wave is perpendicular to the boundary. In order to avoid reflection, it is necessary to increase the distance between the scattering target and the absorption boundary; thus, the calculation amount is increased. Lu et al. introduced nonconforming UPML theory into DG calculations in the case of TM wave and compared the local relative error between PML and first-order SM-ABC [14]. Konig et al. applied the stretched-coordinate PML technique to DGTD [15]. Gedney and Zhao proposed a kind of auxiliary differential equation (ADE) multilevel PML theory by complex frequency domain shift, which is used to cut off the boundary of DGTD region and take good absorption [16]. In the literature [17], the PML theory is used in 3-D DGTD algorithm to realize the absorption of electromagnetic wave through the conformal truncation face by Dosopoulos and Lee. Donderici and Teixeira incorporated conformal PML into the mixed FETD algorithm by utilizing PML constitutive tensors in 2008 [18]. Based on a higher-order curvilinear finite-element and the concept of transformation electromagnetics, Smull et al. proposed anisotropic locally conformal PML; this new PML has favorable accuracy and efficiency [19]. Modave et al. designed PMLs for transient acoustic wave propagation in generally shaped convex truncated domains in DGTD [20]. Yang et al. applied the corner-free truncation strategy to DGTD calculation, which saved the resource effectively [21]. Based on well-posed PML theory, Ren et al. proposed multiaxial full anisotropic media PML for subdomain and applied it to electromagnetic simulation of nonconformal mesh bicrystals; this method avoids the potential late-time instability found in classical PML [22]. However, many physical problems can be solved successfully by two-dimensional models. In two-dimensional models, it has advantages of small memory, fast calculation speed, easy handling the electromagnetic problem of electric large size by computer, and so on in three-dimensional models. Nevertheless, there is only little research focus on the study of DGTD method in two-dimensional problems. Implementation and application of conformal UPML in 2-D are not yet reported. In this paper, the 2-D DGTD iterative formula for elliptic cylindrical UPML is derived and applied to the nested dielectric column scattering.

In order to apply the UPML technique to DGTD in two-dimensional cases, the wave equations of UPML in an orthogonal curvilinear coordinate system are considered with the theory of parameter transformation, the processes of parameters transformation for 2-D UPML, and the specific expressions of transition matrix between the coordinate systems of elliptical and Cartesian are derived. The penalty flux is used for field exchange between units. The 2-D DGTD iterative formula and auxiliary equation in elliptical cylindrical UPML are derived. Numerical results show that very good effectiveness is achieved by using the proposed algorithm.

2. The Expression of UPML Equations from Orthogonal Curvilinear Coordinate System to Cartesian Coordinate System

Suppose that the interface from free space to PML is the isosurface in the orthogonal curve coordinate system (as shown in Figure 1); if the parameters in UPML satisfy the similar condition in Cartesian coordinate system UPML, the incident electromagnetic waves from free space will pass through the interface without reflecting and attenuate rapidly inside the PML.

2.1. UPML Wave Equation in Orthogonal Curvilinear Coordinate System. From the literature [17,23], we know that the form of diagonal matrix of the parameters in the orthogonal curve coordinate system is unchanged. So UPML matching matrix is written as follows:

[mathematical expression not reproducible], (1)

where [u.sub.1], [u.sub.2], [u.sub.3] are unit vectors. And there are

[mathematical expression not reproducible]. (2)

When the interface is the isosurface [[xi].sub.3] (as shown in Figure 1), the no-reflect condition is [s.sub.1] = [s.sub.2] = 1, then (1) is rewritten as follows:

[mathematical expression not reproducible]. (3)

UPML equations in an orthogonal curvilinear coordinate system are as follows:

[mathematical expression not reproducible], (4)

where E and H are electromagnetic fields; [[mu].sub.1] and [[epsilon].sub.1] are permittivity and permeability, respectively. a, b, and c are diagonal tensors. [P.sub.h] and [P.sub.e] are auxiliary variables, which are satisfied with the following relationship:

[mathematical expression not reproducible], (5)

where d and [[kappa].sup.-1] are the same as a, b, and c. Their expressions are as follows:

[mathematical expression not reproducible]. (6)

In practical applications, (4) and (5) often need to be solved in Cartesian coordinates, thus, they are transformed from a curvilinear coordinate system by transformation matrix, which can be written as follows:

[mathematical expression not reproducible], (7)

where [u.sub.1x], [u.sub.2x], [u.sub.3x], ... are components of unit vectors [u.sub.1], [u.sub.2], and [u.sub.3] in x, y, and z direction, respectively. Transformation matrix has a unitary nature [24].

Taking the example of tensor a, the relationship of a from a curvilinear coordinate system to Cartesian coordinates is as follows:

[a.sub.xyz] = [T.sup.-1] x a x T. (8)

It is noted that after the transformation, the [a.sub.xyz] is not the diagonal matrix anymore.

2.2. The Parameter Calculation for Elliptic Cylinder UPML. Since the expression of the transformation matrix depends on the specific coordinate system, the process of calculational parameters in the elliptical cylinder UPML can be described as follows.

Let the ellipse cylindrical axis be z-axis, and the elliptic equation is as follows:

[[x.sup.2]/[a.sup.2]] + [[y.sup.2]/[b.sup.2]] = 1, (9)

where a and b are the semimajor axis and semiminor axis, respectively. The variables w and t are introduced, and let

x = a cos t, y = b sin t, (10)

a = c cosh w, b = c sinh w, [a.sup.2] - [b.sup.2] = [c.sup.2] ([cosh.sup.2] w - [sinh.sup.2] w) = [c.sup.2], (11)

where w [greater than or equal to] 0, 0 [less than or equal to] t [less than or equal to] 2[pi], and c is the focal distance. Its trail is shown in Figure 2. w represents the distance sum between focal points in the ellipse. The transformation matrix between local coordinate of elliptical system and Cartesian coordinate system is decided by rotating angle [beta] (as shown in Figure 2), which is defined as included angle between normal vector [??] and the x-axis in the Cartesian coordinate system.

From [25], we know that the tangent line equation of the observation point P([x.sub.P], [y.sub.P]) on the ellipse is as follows:

y = -[[b.sup.2]/[y.sub.P]][[x.sub.P]/[a.sup.2]]x + [[b.sup.2]/[y.sub.P]]. 12

The slope of the tangent line k is as follows:

k = tan [alpha] = -[[b.sup.2][x.sup.P]/[a.sup.2][y.sub.P]], (13)

where [alpha] is the included angle between tangent line and normal vector [??]. Then, the direction cosine and direction sine between tangential unit vectors and normal vector [??] are as follows:

[mathematical expression not reproducible]. (14)

The [beta] is as follows:

tan [beta] = tan([alpha] - [[pi]/2]) = [[a.sup.2][y.sub.P]/[b.sup.2][x.sub.P]]. (15)

Then, the direction cosine and direction sine between normal unit vectors and x- and y-axes are as follows:

[mathematical expression not reproducible]. (16)

Thus, the transformation matrix in (7) can be written as follows:

[mathematical expression not reproducible]. (17)

Substituting (17) into (8), then

[mathematical expression not reproducible]. (18)

UPML are distributed between the regions [w.sub.0] [less than or equal to] w [less than or equal to] [w.sub.max]. [[xi].sub.3] = w in elliptical cylinder coordinate systems; thus, (2) is rewritten as follows:

[mathematical expression not reproducible], (19)

where

[mathematical expression not reproducible], (20)

where w is the isosurface of the ellipse in the barycentric coordinates of the element in UPML region. [w.sub.0] and [w.sub.max] are the isosurface of the inside and the outside, respectively. m is taken as an integer.

Supposing that [x.sub.p], [y.sub.p], and focal distance c are known, we need to calculate the parameters of a, [b.bar], and w; the equation of the ellipse is

[[x.sup.2.sub.P]/[a.sup.2]] + [[y.sup.2.sub.P]/[b.sup.2]] = 1. (21)

Thus, the formulas of semimajor axis a and semiminor axis b are

[mathematical expression not reproducible], (22)

where [a.sup.2] > [c.sup.2] > 0.

Substituting (22) into (11), we can get

cosh w = [a/c], w = [cosh.sup.-1] ([a/c]). (23)

Inserting (23) into (20), we can obtain the value of conductivity [sigma].

3. Iterative Formulas of 2-D DGTD in Elliptical Cylinder UPML

In the case of two-dimensional TM wave, [E.sub.x] = [E.sub.y] = 0 and [H.sub.z] = 0. Substituting (18) into (4), we can get

[mathematical expression not reproducible]. (24)

The auxiliary equations in a two-dimensional case are

[mathematical expression not reproducible]. (25)

There are nondiagonal elements in (24) and (25), such as [a.sub.xy], [a.sub.yx], [b.sub.xy], [b.sub.yx], and so on. This will lead to a complex form of matrix equations. In order to facilitate the acquisition of the time-domain iterative formula, [partial derivative][H.sub.y]/[partial derivative]t and [partial derivative][H.sub.x]/[partial derivative]t in (24) need to be eliminated. Suppose that

[mathematical expression not reproducible]. (26)

Thus

[mathematical expression not reproducible]. (27)

After using the penalty flux [11] as well as expending it by basic functions in (27), we can get

[mathematical expression not reproducible]. (28)

After (28) is discrete in time axis, then there are

[mathematical expression not reproducible]. (29)

The iterative formulas of auxiliary equations are

[mathematical expression not reproducible], (30)

where [[alpha].sub.hx], [[beta].sub.hx], [[alpha].sub.hy], [[beta].sub.hy], [[alpha].sub.ez], [c.sub.h], [c.sub.h], [c.sub.j], and [c.sub.k] are iterative coefficients.

4. Numerical Examples

In the elliptical cylindrical computing region, the semimajor and semishort axes of the outer PML boundary are 4 m and 3 m, respectively. And the inner PML boundary has a semimajor axis of 3.5 m. Both the inner and outer PML boundaries have the same focus length of 2.29 m. The region is divided into 46,122 triangular units and 23,298 nodes with the discrete scale of 0.05 m. The time discrete interval is taken as dt = 0.16 x [10.sup.-10] s. The radiant source of Gaussian impulse which has the width and peak of [tau] = 200 dt and [t.sub.0] = 100 dt is set at the center of the region (0 m, 0 m). Time-domain waveform and spectrum at the monitoring point of (0, 2.5 m) are shown in Figure 3(a). The solid and dotted line are the impulse waveform in time-domain and the spectrum, respectively. And "o" is the analytic solution obtained by the second kind of zero order Hankel function. Refer to [18]; the resulting reflection error is shown in Figure 3(b). It can be seen from the figure that the reflection error is below -37 dB.

The PEC elliptic column has a semimajor axis and a semishort axis of 1.5 mm and 0.5 mm, respectively. Incident plane wave has a wavelength of 2 mm. The computational region is truncated with conformal PML layer. The result of bistatic scattering for metal elliptical column is shown in Figure 4 ("o"). For easy comparison, the result of method of moment (solid line) is also shown in figure. As illustrated, the results are consistent.

Two nested dielectric elliptical columns target truncated by conformal PML. Sparse discrete conformal model is shown in Figure 5(a). The inner elliptical target has a relative dielectric coefficient of 4, and the semimajor axis and semiminor axis are 5 mm and 2 mm, respectively. The outer has a relative dielectric coefficient of 2, which are 6 mm and 3 mm. The rest of conformal boundaries from inside to outside are connected boundary, extrapolated boundary, and PML boundary. Both the inner and outer PML boundaries have the same focus distance. The calculating region is divided into 94,634 triangular units and 47,662 nodes with the discrete scale of 0.08 mm. The incident plane wave has a wavelength of 3 mm. Bistatic RCS for coated dielectric elliptical column target is shown in Figure 5(b) ("o"). The results of FDTD (solid line) and DGTD ("[??]") (its region is truncated with rectangular PML) are also shows in the figure. As can be seen, the calculation results of the three coincide.

In order to verify the benefits of conformal UPML, the numbers of unit and node in two different regions by three discrete scales of 0.1mm, 0.08 mm, and 0.04 mm are given in Table 1. Compared to the elliptical region, their numbers in rectangular computational region (the region surrounded by rectangular boundary which is shown in Figure 5(a)) increase by about 27%. In order to further compare the time and memory consumption, the following calculation is given. Take the case of discrete scale of 0.08 mm example (shown in Table 1), the numbers of unit and node in rectangular region are 121,474 and 61,175, respectively, which increase by 28.3%. And the memory and time increase by 25.5% and 32.5%, respectively. It is clear that conformal UPML has obvious advantages in terms of computational resources and computational time savings.

5. Conclusions

According to the DGTD equations and the parameter transformation theory in orthogonal curvilinear coordinate systems, the processes of parameter transformation for 2-D UPML between the coordinate systems of elliptical and Cartesian are presented, and expressions of transition matrix are derived. The 2-D iterative formulas of DGTD and that of auxiliary equation in elliptical cylinder UPML are derived, and the conformal UPML calculation in DGTD are realized. Numerical results show that good absorption for the outgoing wave can be realized by using conformal UPML. Because of the discrete characteristics of the unstructured grids, the absorption layer can be set up flexibly according to the shape of the target. Compared with the computational cost in rectangular PML layer, the memory and computation time for conformal structure decrease by 20.3% and 24.5%, respectively. Thus, the method presented in this paper will be very beneficial to the engineering computation since lots of computing resources can be saved.

https://doi.org/10.1155/2017/6358561

Conflicts of Interest

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research is supported in part by the National Natural Science Foundation of China under Grant nos. 61571348, 61231003, and 61401344.

References

[1] J. P. Berenger, "A perfectly matched layer for the absorption of electromagnetic waves," Journal of Computational Physics, vol. 114, no. 2, pp. 185-200, 1994.

[2] W. C. Chew and W. H. Weedon, "A 3D perfectly matched medium from modified Maxwell's equations with stretched coordinates," Microwave and Optical Technology Letters, vol. 7, no. 13, pp. 599-604, 1994.

[3] Z. S. Sacks, D. M. Kingsland, R. Lee, and J.-F. Lee, "A perfectly matched anisotropic absorber for use as an absorbing boundary condition," IEEE Transactions on Antennas and Propagation, vol. 43, no. 12, pp. 1460-1463, 1995.

[4] F. L. Teixeira and W. C. Chew, "Analytical derivation of a conformal perfectly matched absorber for electromagnetic waves," Microwave and Optical Technology Letters, vol. 17, no. 4, pp. 231-236, 1998.

[5] K.-P. Hwang and J.-M. Jin, "Application of a hyperbolic grid generation technique to a conformal PML implementation," IEEE Microwave and Guided Wave Letters, vol. 9, no. 4, pp. 137-139, 1999.

[6] F. L. Teixeira, K.-P. Hwang, W. C. Chew, and J.-M. Jin, "Conformal PML-FDTD schemes for electromagnetic field simulations: a dynamic stability study," IEEE Transactions on Antennas and Propagation, vol. 49, no. 6, pp. 902-907, 2001.

[7] Y. Zhang, N. Feng, H. Zheng, and H. Liu, "A corner-free truncation strategy in three-dimensional FDTD computation," IEEE Transactions on Electromagnetic Compatibility, vol. 58, no. 2, pp. 512-522, 2016.

[8] I. Perugia and D. Schotzau, "The hp-local discontinuous Galerkin method for low-frequency time-harmonic Maxwell equations," Mathematics of Computation, vol. 72, no. 243, pp. 1179-1214, 2003.

[9] B. Cockburn, F. Li, and C. W. Shu, "Locally divergence-free discontinuous Galerkin methods for the Maxwell equations," Journal of Computational Physics, vol. 194, no. 2, pp. 588-610, 2004.

[10] J. Chen and Q. H. Liu, "Discontinuous Galerkin time-domain methods for multiscale electromagnetic simulations: a review," Proceedings of the IEEE, vol. 101, no. 2, pp. 242-254, 2013.

[11] J. Alvarez, A Discontinuous Galerkin Finite Element Method for the Time-Domain Solution of Maxwell Equations, University of Granada, Spain, 2013.

[12] S. Yan and J. M. Jin, "A dynamic p-adaptive DGTD algorithm for electromagnetic and multiphysics simulations," IEEE Transactions on Antennas and Propagation, vol. 65, no. 5, pp. 2446-2459, 2017.

[13] P. Bonnet, X. Ferrieres, B. L. Michielsen, P. Klotz, and J. L. Roumiguieres, "Finite-volume time domain method," in Time Domain Electromagnetics, pp. 307-368, Acade Press, San Diego, CA, USA, 1999.

[14] T. Lu, P. Zhang, and W. Cai, "Discontinuous Galerkin methods for dispersive and lossy Maxwell's equations and PML boundary conditions," Journal of Computational Physics, vol. 200, no. 2, pp. 549-580, 2004.

[15] M. Konig, C. Prohm, K. Busch, and J. Niegemann, "Stretched-coordinate PMLs for Maxwell's equations in the discontinuous Galerkin time-domain method," Optics Express, vol. 19, no. 5, pp. 4618-4631, 2011.

[16] S. D. Gedney and B. Zhao, "An auxiliary differential equation formulation for the complex-frequency shifted PML," IEEE Transactions on Antennas and Propagation, vol. 58, no. 3, pp. 838-847, 2010.

[17] S. Dosopoulos and J. F. Lee, "Interior penalty discontinuous Galerkin finite element method for the time-dependent first order Maxwell's equations," IEEE Transactions on Antennas and Propagation, vol. 58, no. 12, pp. 4085-4090, 2010.

[18] B. Donderici and F. L. Teixeira, "Conformal perfectly matched layer for the mixed finite element time-domain method," IEEE Transactions on Antennas and Propagation, vol. 56, no. 4, pp. 1017-1026, 2008.

[19] A. P. Smull, A. B. Manic, S. B. Manic, and B. M. Notaros, "Anisotropic locally-conformal perfectly matched layer for higher order curvilinear finite-element modeling," IEEE Transactions on Antennas and Propagation, vol. 65, no. 12, pp. 7157-7165, 2017.

[20] A. Modave, J. Lambrechts, and C. Geuzaine, "Perfectly matched layers for convex truncated domains with discontinuous Galerkin time domain simulations," Computers & Mathematics with Applications, vol. 73, no. 4, pp. 684-700, 2017.

[21] Q. Yang, B. Wei, L. Li, and D. Ge, "Implementation of corner-free truncation strategy in DGTD method," Waves in Random and Complex Media, vol. 27, no. 2, pp. 367-380, 2017.

[22] Q. Ren, Q. Zhan, and Q. H. Liu, "An improved subdomain level non-conformal discontinuous Galerkin time domain (DGTD) method for materials with full-tensor constitutive parameters," IEEE Photonics Journal, vol. 9, no. 2, article 2600113, 2017.

[23] F. L. Teixeira and W. C. Chew, "Systematic derivation of anisotropic PML absorbing media in cylindrical and spherical coordinates," IEEE Microwave and Guided Wave Letters, vol. 7, no. 11, pp. 371-373, 1997.

[24] D. B. Ge and B. Wei, Electromagnetic Wave Theory, Science Press, Beijing, 2011.

[25] L. X. Wang, D. Z. Fang, and M. Y. Zhang, Mathematics Manual, People's Education Press, Beijing, 1979.

Linqian Li, (1,2) Bing Wei, (1,2) Qian Yang, (1,2) and Debiao Ge (1,2)

(1) School of Physics and Optoelectronic Engineering, Xidian University, Xi'an 710071, China

(2) Collaborative Innovation Center of Information Sensing and Understanding, Xidian University, Xi'an 710071, China

Correspondence should be addressed to Linqian Li; llq11123@163.com

Received 26 May 2017; Revised 26 October 2017; Accepted 12 November 2017; Published 17 December 2017

Academic Editor: Ikmo Park

Caption: Figure 1: Discrete unit of PML in orthogonal curvilinear coordinate system.

Caption: Figure 2: Elliptical cylinder coordinate systems.

Caption: Figure 3: Radiant field and reflection error.

Caption: Figure 4: Bistatic RCS of PEC elliptic column.

Caption: Figure 5: Bistatic scattering for target.

Table 1: Comparison of the number of unit and node in different calculated regions. Discrete scale 0.1 mm 0.04 mm Shape Number Number Number Number Number of unit of node of unit of node of unit Nonconformal 77,670 39,184 489,534 245,638 121,474 Conformal 61,088 30,821 383,690 192,534 94,634 Discrete scale 0.08 mm Shape Number Storages Time of node (MB) (mins) Nonconformal 61,175 116.5 134.7 Conformal 47,662 92.8 101.6

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Li, Linqian; Wei, Bing; Yang, Qian; Ge, Debiao |

Publication: | International Journal of Antennas and Propagation |

Date: | Jan 1, 2018 |

Words: | 3596 |

Previous Article: | Chiral All-Dielectric Metasurface Based on Elliptic Resonators with Circular Dichroism Behavior. |

Next Article: | Flexible and Conformal Antennas and Applications. |