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

Authors: Linqian Li; Bing Wei; Qian Yang; Debiao Ge
Journal: International Journal of Antennas and Propagation
 (2017)

---

## Abstract

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%.


---

## Body

## 1. Introduction

Truncated 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]. König 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ξ3 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.Figure 1
Discrete unit of PML in orthogonal curvilinear coordinate system.


### 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:
(1)Λu1u2u3=diags2s3s1s3s1s2s1s2s3=u1u1s2s3s1+u2u2s3s1s2+u3u3s1s2s3,where u1,u2,u3 are unit vectors. And there are
(2)s1=κ1+σ1jωε0,s2=κ2+σ2jωε0,s3=κ3+σ3jωε0.When the interface is the isosurfaceξ3 (as shown in Figure 1), the no-reflect condition is s1=s2=1, then (1) is rewritten as follows:
(3)Λu1u2u3=diags3s31s3=u1u1s3+u2u2s3+u3u31s3.UPML equations in an orthogonal curvilinear coordinate system are as follows:(4)μ1a⋅∂H∂t+μ1b⋅H+μ1c⋅Ph+∇×E=0,ε1a⋅∂E∂t+ε1b⋅E+ε1c⋅Pe−∇×H=0,where E and H are electromagnetic fields; μ1 and ε1 are permittivity and permeability, respectively. a, b, and c are diagonal tensors. Ph and Pe are auxiliary variables, which are satisfied with the following relationship:
(5)∂Ph∂t+d⋅Ph−κ−1⋅H=0,∂Pe∂t+d⋅Pe−κ−1⋅E=0,where d and κ−1 are the same as a, b, and c. Their expressions are as follows:
(6)a=diagκ2κ3κ1κ3κ1κ2κ1κ2κ3,b=diag1κ1ε0σ2κ3+σ3κ2−a11σ11κ2ε0σ3κ1+σ1κ3−a22σ21κ3ε0σ1κ2+σ2κ1−a33σ3,c=diagσ2σ3ε02−b11σ1ε0σ3σ1ε02−b22σ2ε0σ1σ2ε02−b33σ3ε0,d=diagσ1κ1ε0σ2κ2ε0σ3κ3ε0,κ−1=diag1κ11κ21κ3.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:
(7)T=u1xu1yu1zu2xu2yu2zu3xu3yu3z,where u1x,u2x,u3x,… are components of unit vectors u1, u2, and u3 in x, y, and z direction, respectively. Transformation matrix has a unitary nature [24].Taking the example of tensora, the relationship of a from a curvilinear coordinate system to Cartesian coordinates is as follows:
(8)axyz=T−1⋅a⋅T.It is noted that after the transformation, theaxyz 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 bez-axis, and the elliptic equation is as follows:
(9)x2a2+y2b2=1,where a and b are the semimajor axis and semiminor axis, respectively. The variables w and t are introduced, and let
(10)x=acost,y=bsint,(11)a=ccoshw,b=csinhw,a2−b2=c2cosh2w−sinh2w=c2,where w≥0, 0≤t≤2π, 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 β (as shown in Figure 2), which is defined as included angle between normal vector w^ and the x-axis in the Cartesian coordinate system.Figure 2
Elliptical cylinder coordinate systems.From [25], we know that the tangent line equation of the observation point PxP,yP on the ellipse is as follows:
(12)y=−b2yPxPa2x+b2yP.The slope of the tangent linek is as follows:
(13)k=tanα=−b2xPa2yP,where α is the included angle between tangent line and normal vector w^. Then, the direction cosine and direction sine between tangential unit vectors and normal vector w^ are as follows:
(14)cosα=±11+tan2α,sinα=±tanα1+tan2α.Theβ is as follows:
(15)tanβ=tanα−π2=a2yPb2xP.Then, the direction cosine and direction sine between normal unit vectors andx- and y-axes are as follows:
(16)cosβ=±11+tan2β,sinβ=±tanβ1+tan2β.Thus, the transformation matrix in (7) can be written as follows:
(17)Toval=−sinβcosβ0001cosβsinβ0.Substituting (17) into (8), then
(18)axyz=a11sin2β+a33cos02βa33−a11sinβcosβ0a33−a11sinβcosβa11cos2β+a33sin2β000a22=axxaxy0ayxayy000azz.UPML are distributed between the regionsw0≤w≤wmax. ξ3=w in elliptical cylinder coordinate systems; thus, (2) is rewritten as follows:
(19)s1=sφ=1,s2=sz=1,s3=sw=κw+σjωε0,where
(20)σ=σmaxw−w0wmax−w0m,where w is the isosurface of the ellipse in the barycentric coordinates of the element in UPML region. w0 and wmax are the isosurface of the inside and the outside, respectively. m is taken as an integer.Supposing thatxp, yp, and focal distance c are known, we need to calculate the parameters of a, b, and w; the equation of the ellipse is
(21)xP2a2+yP2b2=1.Thus, the formulas of semimajor axisa and semiminor axis b are
(22)a2=−xP2−yP2−c2±xP2−yP2−c22−4c2xP22,b2=a2−c2,where a2>c2>0.Substituting (22) into (11), we can get
(23)coshw=ac,w=cosh−1ac.Inserting (23) into (20), we can obtain the value of conductivity σ.


## 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:
(1)Λu1u2u3=diags2s3s1s3s1s2s1s2s3=u1u1s2s3s1+u2u2s3s1s2+u3u3s1s2s3,where u1,u2,u3 are unit vectors. And there are
(2)s1=κ1+σ1jωε0,s2=κ2+σ2jωε0,s3=κ3+σ3jωε0.When the interface is the isosurfaceξ3 (as shown in Figure 1), the no-reflect condition is s1=s2=1, then (1) is rewritten as follows:
(3)Λu1u2u3=diags3s31s3=u1u1s3+u2u2s3+u3u31s3.UPML equations in an orthogonal curvilinear coordinate system are as follows:(4)μ1a⋅∂H∂t+μ1b⋅H+μ1c⋅Ph+∇×E=0,ε1a⋅∂E∂t+ε1b⋅E+ε1c⋅Pe−∇×H=0,where E and H are electromagnetic fields; μ1 and ε1 are permittivity and permeability, respectively. a, b, and c are diagonal tensors. Ph and Pe are auxiliary variables, which are satisfied with the following relationship:
(5)∂Ph∂t+d⋅Ph−κ−1⋅H=0,∂Pe∂t+d⋅Pe−κ−1⋅E=0,where d and κ−1 are the same as a, b, and c. Their expressions are as follows:
(6)a=diagκ2κ3κ1κ3κ1κ2κ1κ2κ3,b=diag1κ1ε0σ2κ3+σ3κ2−a11σ11κ2ε0σ3κ1+σ1κ3−a22σ21κ3ε0σ1κ2+σ2κ1−a33σ3,c=diagσ2σ3ε02−b11σ1ε0σ3σ1ε02−b22σ2ε0σ1σ2ε02−b33σ3ε0,d=diagσ1κ1ε0σ2κ2ε0σ3κ3ε0,κ−1=diag1κ11κ21κ3.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:
(7)T=u1xu1yu1zu2xu2yu2zu3xu3yu3z,where u1x,u2x,u3x,… are components of unit vectors u1, u2, and u3 in x, y, and z direction, respectively. Transformation matrix has a unitary nature [24].Taking the example of tensora, the relationship of a from a curvilinear coordinate system to Cartesian coordinates is as follows:
(8)axyz=T−1⋅a⋅T.It is noted that after the transformation, theaxyz 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 bez-axis, and the elliptic equation is as follows:
(9)x2a2+y2b2=1,where a and b are the semimajor axis and semiminor axis, respectively. The variables w and t are introduced, and let
(10)x=acost,y=bsint,(11)a=ccoshw,b=csinhw,a2−b2=c2cosh2w−sinh2w=c2,where w≥0, 0≤t≤2π, 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 β (as shown in Figure 2), which is defined as included angle between normal vector w^ and the x-axis in the Cartesian coordinate system.Figure 2
Elliptical cylinder coordinate systems.From [25], we know that the tangent line equation of the observation point PxP,yP on the ellipse is as follows:
(12)y=−b2yPxPa2x+b2yP.The slope of the tangent linek is as follows:
(13)k=tanα=−b2xPa2yP,where α is the included angle between tangent line and normal vector w^. Then, the direction cosine and direction sine between tangential unit vectors and normal vector w^ are as follows:
(14)cosα=±11+tan2α,sinα=±tanα1+tan2α.Theβ is as follows:
(15)tanβ=tanα−π2=a2yPb2xP.Then, the direction cosine and direction sine between normal unit vectors andx- and y-axes are as follows:
(16)cosβ=±11+tan2β,sinβ=±tanβ1+tan2β.Thus, the transformation matrix in (7) can be written as follows:
(17)Toval=−sinβcosβ0001cosβsinβ0.Substituting (17) into (8), then
(18)axyz=a11sin2β+a33cos02βa33−a11sinβcosβ0a33−a11sinβcosβa11cos2β+a33sin2β000a22=axxaxy0ayxayy000azz.UPML are distributed between the regionsw0≤w≤wmax. ξ3=w in elliptical cylinder coordinate systems; thus, (2) is rewritten as follows:
(19)s1=sφ=1,s2=sz=1,s3=sw=κw+σjωε0,where
(20)σ=σmaxw−w0wmax−w0m,where w is the isosurface of the ellipse in the barycentric coordinates of the element in UPML region. w0 and wmax are the isosurface of the inside and the outside, respectively. m is taken as an integer.Supposing thatxp, yp, and focal distance c are known, we need to calculate the parameters of a, b, and w; the equation of the ellipse is
(21)xP2a2+yP2b2=1.Thus, the formulas of semimajor axisa and semiminor axis b are
(22)a2=−xP2−yP2−c2±xP2−yP2−c22−4c2xP22,b2=a2−c2,where a2>c2>0.Substituting (22) into (11), we can get
(23)coshw=ac,w=cosh−1ac.Inserting (23) into (20), we can obtain the value of conductivity σ.


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

In the case of two-dimensional TM wave,Ex=Ey=0 and Hz=0. Substituting (18) into (4), we can get
(24)μ1axx∂Hx∂t+axy∂Hy∂t+μ1bxxHx+bxyHy+μ1cxxPh,x+cxyPh,y+∂Ez∂y=0,μ1ayx∂Hx∂t+ayy∂Hy∂t+μ1byxHx+byyHy+μ1cyxPh,x+cyyPh,y−∂Ez∂x=0,ε1azz∂Ez∂t+ε1bzzEz+ε1czzPe,z−∂Hy∂x−∂Hx∂y=0.The auxiliary equations in a two-dimensional case are(25)∂Ph,x∂t+dxxPh,x+dxyPh,y−κxx−1Hx+κxy−1Hy=0,∂Ph,y∂t+dyyPh,y+dyxPh,x−κyx−1Hx+κyy−1Hy=0,∂Pe,z∂t+dzzPe,z−κzz−1Ez=0⇒∂Pe,z∂t−Ez=0.There are nondiagonal elements in (24) and (25), such as axy,ayx,bxy,byx, 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, ∂Hy/∂t and ∂Hx/∂t in (24) need to be eliminated. Suppose that
(26)A1=axxayy−ayxaxy,A2=ayybxx−axybyx,A3=ayybxy−axybyy,A4=ayycxx−axycyx,A5=ayycxy−axycyy,A6=ayy,A7=axy,B2=axxbyy−ayxbxy,B3=axxbyx−ayxbxxB4=axxcyx−ayxcxx,B5=axxcyy−ayxcxy,B6=axx,B7=ayx,C1=czz,D1=dxy,D2=dyx.Thus
(27)μ1A1∂Hx∂t+μ1A2Hx+μ1A3Hy+μ1A4Ph,x+μ1A5Ph,y+A6∂Ez∂y+A7∂Ez∂x=0,μ1A1∂Hy∂t+μ1B2Hy+μ1B3Hx+μ1B4Ph,x+μ1B5Ph,y−B6∂Ez∂x−B7∂Ez∂y=0.After using the penalty flux [11] as well as expending it by basic functions in (27), we can get
(28)μ1A1Mm∂Hx∂t+μ1A2MmHx+μ1A3MmHy+μ1A4MmPh,x+μ1A5MmPh,y+A6SyEz+A7SxEz+A6nyκemFe,z+A6nyνhmnxFh,y−A6nyνhmnyFh,x+A7nxκemFe,z+A7nxνhmnxFh,y−A7nxνhmnyFh,x=0,μ1A1Mm∂Hy∂t+μ1B2MmHy+μ1B3MmHx+μ1B4MmPh,x+μ1B5MmPh,y−B6SxEz−B7SyEz−B6nxκemFe,z+B6nxνhmnxFh,y−B6nxνhmnyFh,x+B7nyκemFe,z+B7nyνhmnxFh,y−B7nyνhmnyFh,x=0,ε1azzMm∂Ez∂t+ε1bzzMmEz+ε1C1MmPe,z−SxHy+SyHx+κhmnyFh,x−nxFh,y−νemFe,z=0.After (28) is discrete in time axis, then there are
(29)Hx,n+1/2=αhxHx,n−1/2−βhxA3Hy,n−1/2−βhxA4Ph,x,n−βhxA5Ph,y,n−βhx1μ1Mm−1A6Sy+A7SxEz,n+A6κemnyFe,z,n+A6νhmnxnyFh,y,n−1/2−A6νhmnynyFh,x,n−1/2+A7κemnxFe,z,n+A7νhmnxnxFh,y,n−1/2−A7νhmnxnyFh,x,n−1/2.Hy,n+1/2=αhyHy,n−1/2−βhyB3Hx,n−1/2−βhyB4Ph,x,n−βhyB5Ph,y,n+βhy1μ1Mm−1B6Sx+B7SyEz,n+B6κemnxFe,z,n+B6νhmnxnxFh,y,n−1/2−B6νhmnxnyFh,x,n−1/2+B7κemnyFe,z,n+B7νhmnynxFh,y,n−1/2−B7νhmnynyFh,x,n−1/2,Ez,n+1=αezEz,n−C1βezPe,z,n+1/2+βez1ε1Mm−1SxHy,n+1/2−SyHx,n+1/2−κhmnyFh,x,n+1/2−nxFh,y,n+1/2+νemFe,z,n.The iterative formulas of auxiliary equations are(30)Ph,x,n+1=ch⋅Ph,x,n−ci⋅D1Ph,y,n+ci⋅κxx−1Hx,n+1/2+κxy−1Hy,n+1/2,Ph,y,n+1=cj⋅Ph,y,n−ck⋅D2Ph,x,n+ck⋅κyx−1Hx,n+1/2+κyy−1Hy,n+1/2,Pe,z,n+1/2=Pe,z,n−1/2+ΔtEz,n,where αhx, βhx, αhy, βhy, αez, βez, ch, ch, cj, and ck 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 asdt=0.16×10−10s. The radiant source of Gaussian impulse which has the width and peak of τ=200dt and t0=100dt 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.Figure 3
Radiant field and reflection error.
(a)
Time-domain waveform and spectrum
(b)
Reflection errorThe 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 Figure4 (“o”). For easy comparison, the result of method of moment (solid line) is also shown in figure. As illustrated, the results are consistent.Figure 4
Bistatic RCS of PEC elliptic column.Two nested dielectric elliptical columns target truncated by conformal PML. Sparse discrete conformal model is shown in Figure5(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.Figure 5
Bistatic scattering for target.
(a)
Sparse discrete model
(b)
Bistatic RCS for two nested elliptical columnsIn order to verify the benefits of conformal UPML, the numbers of unit and node in two different regions by three discrete scales of 0.1 mm, 0.08 mm, and 0.04 mm are given in Table1. 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.Table 1
Comparison of the number of unit and node in different calculated regions.
Shape
Discrete scale
0.1 mm
0.04 mm
0.08 mm
Number of unit
Number of node
Number of unit
Number of node
Number of unit
Number of node
Storages (MB)
Time (mins)
Nonconformal
77,670
39,184
489,534
245,638
121,474
61,175
116.5
134.7
Conformal
61,088
30,821
383,690
192,534
94,634
47,662
92.8
101.6


## 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.



---

*Source: 6358561-2017-12-17.xml*