EPJ Appl. Metamat.
Volume 6, 2019
Metamaterials Research and Development in Japan
|Number of page(s)||7|
|Published online||30 July 2019|
FDTD modeling of nonperiodic antenna located above metasurface using surface impedance boundary condition
Department of Electrical and Electronic Engineering, Tokyo University of Agriculture and Technology, Koganei, Tokyo 184-8588, Japan
* e-mail: email@example.com
Accepted: 12 July 2019
Published online: 30 July 2019
This paper investigates an FDTD modeling method for precisely calculating the characteristics of a single, that is, a nonperiodic antenna located above a metasurface that consists of an infinite periodic conducting element on a flat dielectric substrate. The original FDTD method requires enormous computational resources to analyze such structures because an appropriate periodic boundary condition (PBC) is not supported, and a brute force approach has to be used for this reason. Another option is to use the array scanning method in which a single source is synthesized from a superposition of infinite phased array of point sources. In this method, some problems such as a mutual coupling between the single antenna and the metasurface, a computational error contained in a numerical integration over the Brillouin zone and so on have not been resolved yet. In order to resolve these difficulties and to reduce computational resources, a surface impedance boundary condition (SIBC) is incorporated into the FDTD method in this paper. The validity of the method is numerically confirmed by calculating an input impedance and a radiation pattern of a horizontal dipole antenna located above the metasurface.
© T. Uno et al., published by EDP Sciences, 2019
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The use of periodic artificial materials or metamaterials has been investigated for antennas, microwave devices and many other applications in the field of optics due to their own high ability to enhance the original characteristics and to control electromagnetic wave propagation [1–3]. An example is a periodic leaky-wave antenna that a direction of antenna beam is changed by an exciting frequency keeping at an enhanced directivity. One of other applications is an artificial magnetic conductor (AMC) that can be used to create low profile antennas . The AMC is often realized by a mushroom type electromagnetic band gap (EBG) structure that is used to suppress surface-wave propagation as well . The purpose of this paper is to provide a computational method of electromagnetic fields excited by a single(nonperiodic) antenna that is placed near a two dimensional periodic structure. This type of antenna is sometimes referred to as a metamaterial-inspired antenna in the field of antenna engineering.
The finite difference time domain (FDTD) method [6–9] has been widely recognized as one of the most powerful electromagnetic computational methods. Since the computational algorithm is simple and straightforward, and furthermore a practically acceptable accuracy can be easily obtained, the FDTD method has been applied to analyze the metamaterial antennas as well. For a perfectly periodic geometry such as a plane wave scattering by the metasurface, calculation costs can drastically be reduced using a periodic boundary condition (PBC), that is, Floquet's condition or Bloch's condition because the field computations are completed within a unit cell only [8–10]. On the contrary, the PBC can not be used when the single, that is, the nonperiodic antenna is included in the periodic structure. Therefore, a sufficiently large computational space that contains both quite a lot of periodic elements and the nonperiodic antenna is required to precisely analyze the electromagnetic fields when the original, that is, the brute force FDTD technique is used. Accordingly the computational resources increase remarkably in this brute-force FDTD(BF-FDTD) scheme. In order to overcome this difficulty, the array scanning method  has been incorporated into the FDTD calculation in which the single source is synthesized from a superposition of infinite phased array of point sources [12,13]. Although this method is theoretically attractive, an integral representation over the Brillouin zone are included in an electric field expression and its numerical calculation can not be accurately performed because an integrand exhibits a severe singularity. Furthermore, it may be extremely difficult to extend it to the method that can precisely calculate a mutual coupling between the single antenna and the metasurface. On the other hand, accurate surface impedance expressions for various metasurface have already been given, or can be easily calculated by the FDTD method or other numerical methods as well. Considering these situations, a computational scheme that the metasurface is replaced by the surface impedance and then a surface impedance boundary condition (SIBC) is incorporated into the FDTD calculation [14–16], is effectively available. In this paper, the validity and the accuracy of the SIBC-based FDTD method are numerically tested in calculations of an input impedance and a radiation pattern of a horizontal dipole antenna located above the metasurface, because this method is fundamentally an approximation method.
In the next section, a transmission line approximation of the surface impedance of the metasurface is briefly reviewed. These results are interpolated by a rational function of an angular frequency in order to successfully introduce it into the FDTD algorithm for analyzing a scattering problem due to a frequency dispersive medium. The metasurface treated here consists of the infinite periodic rectangular patch which is placed on the dielectric slab backed by the perfect conductor. The FDTD formulation of the electromagnetic fields for the frequency dispersive surface impedance is discussed in Section 3. In the FDTD method, a far-field is usually calculated using a far-zone transformation formula in time domain, but this paper uses the equivalence theorem in frequency domain directly. In Section 4, the computational scheme of the radiation pattern of the dipole antenna located above the metasurface is indicated after an exact electric field representation in far-zone is derived. It is numerically demonstrated using some numerical examples that the SIBC-FDTD is extremely effective for the far-field calculations.
An impedance boundary condition on a curved surface S illustrated in Figure 1 is expressed as [17,18](1)or simply , where is a unit normal to the boundary surface S, and are an electric current and a magnetic current on S, respectively. is a dyadic surface impedance and satisfies a condition .
Electric and magnetic currents on curved surface S.
The metasurface considered in this paper is illustrated in Figure 2. A W x × W y PEC patch is periodically placed on the dielectric slab backed by the PEC. The periods in x and y directions are p x and p y , respectively. The surface impedance Z s of the flat structure as shown in Figure 2 can be represented by a scalar quantity in accordance with a polarization of an incident plane wave and has already been derived using the transmission line approximation , and is given by where Z g (θ, ω) corresponds to the impedance for periodic patch and Z d (θ , ω) the impedance of a shorted transmission line of a length d. θ is an incident angle as shown in Figure 2 and ω is the angular frequency. When μ r = 1, Z g (θ, ω) and Z d (θ , ω) are expressed by(2) (3)where p xy = p x /p y , ε eff = (ε r + 1)/2, ,(4) , and . A refraction angle θ t is determined by Snell's law .
When the approximate expression of the surface impedance is not explicitly given, the equivalent numerical data can be available by calculating the frequency and angular dependencies of the reflection coefficients for the plane wave incidence using the FDTD method or other computational techniques.
Rectangular patch on dielectric slab backed by PEC.
A time domain surface impedance is needed in order to incorporate the SIBC into the FDTD computation, because the electric field is related to the magnetic field through the frequency dependent surface impedance as discussed in the above section. However, it is quite difficult to transform the surface impedance Z s (θ, ω) to the time domain impedance Z s (θ, t) using equations (2)–(4), because Z s indicates a complicated frequency dependence. On the other hand, the electromagnetic scattering investigated in this paper is considered as a kind of passive system such as a passive electric circuit. Therefore, it is reasonable that the impedance is expressed by a positive-real rational function with respect to the angular frequency ω like a Foster function for a loss-less electrical circuit.
In order to obtain the time domain surface impedance expression that is easily applicable to the FDTD method, we assume that the surface impedance is expressed as(5)where s = jω, and A n and p n , (n = 1, 2, … , N) are unknown coefficients to be decided from the given data of Z s (s) for a fixed incident angle θ . The time domain surface impedance can be easily obtained from eq. (5), and is given by(6)
In order to confirm the validity and the accuracy of the interpolation method for the surface impedance given in Section 2.1, the coefficients A n and p n were decided by solving the least square equations obtained from given data of Z s (jω) from 10 GHz to 20 GHz. In this calculation, we set , d = 2.3 mm and ε r = 4.0, and then 3 appropriate poles were found, that is N = 3 the approximated surface impedance and the original surface impedance are shown in Figure 3. It is found that the approximated surface impedance agrees with the original surface impedance perfectly.
Rational function interpolation of surface reactance.
In order to make full use of the capabilities of metasurface, the antenna is naturally placed very near the surface. Hence a height h shown in Figure 4 is small compared with an operating wavelength in practice. In this section, the FDTD method using the SIBC is studied for such situation. On the other hand, the radiated fields can be considered as the superposition from various angular spectrums and are incident to the metasurface with various angular spectrums. The reflect-backed field components to the antenna may be confined within the partial components incident to an area just under the antenna, so that the surface impedance at normal incident dominates the antenna impedance. Consequently, only the surface impedance at normal incident is incorporated into the FDTD formulation.
Transforming the SIBC to the time domain, next two convolutional representations are obtained at the boundary surface shown in Figure 4 (7)where Z s (t) is the Fourier transformed time domain surface impedance and is given by equation (6) in the proposed SIBC-FDTD calculation. Since the Z s (t) is expressed by an exponential function, the convolution integrals (7) can easily be evaluated by the recursive convolution (RC) method  or the piecewise linear RC (PLRC) method . In this paper, the RC method was used for simplicity. The electromagnetic fields in other computational region are calculated by the conventional Yee's algorithm.
Antenna above periodic metasurface or surface impedance boundary.
In this section, the validity of the impedance calculation using the SIBC-FDTD is confirmed. The overall structure of the metasurface and its unit cell are the same as Figure 2, and the same parameters given in the above section were also used for all numerical examples. Above the interface at z = 0, the very thin horizontal dipole antenna whose length is 2l = 10 mm, is placed at a position of z = h as shown in Figure 4. In all FDTD calculations shown below, the cell size was set . A PML (perfectly matched layer) absorbing boundary [6–9] was used for truncating the computational space.
The SIBC-FDTD is compared with two different kinds of the brute-force FDTD methods. One of which is an infinite slab model shown in Figure 5a in which the dielectric slab is equivalent to infinite but the patches outside the PML are not included in the calculation. Other the finite slab model is shown in Figure 5b. This model is often used in the original FDTD and/or the commercially-available simulators that are not provided with an absorbing boundary capable to a piecewise homogeneous medium as treated in this paper.
Since the computational space is truncated by the PML, we checked how large computational space is required for obtaining the accurate input impedance. The input impedance for the number of patches including the computational space is shown in Figure 6 where the total number of patches is given by N x × N x and . ‘BF∞' denotes the brute-force FDTD using the infinite slab model shown in Figure 5a, and ‘BFf' is the one shown in Figure 5b. The operating frequency was set 15 GHz at which the surface impedance resonates as shown in Figure 3 and hence the metasurface acts as a perfect magnetic conductor (PMC) for a normal incident plane wave. It is found that ‘BF∞' converges approximately to the SIBC at N x = 30 which corresponds 3 wavelengths at 15GHz. On the other hand, the ‘BFf' slowly converges to the other value near ‘BF∞' in less than about 10 Ω. This relatively small difference may be caused by the reflected surface wave at the boundary between the dielectric slab and the air.
The frequency characteristics of the input impedance are shown in Figure 7 for h = 10 mm, 7.5 mm and 5 mm. In all BF-FDTD calculations, 30 × 30 unit cells are included in the computational space. It is found that both methods agree fairly well with each other and that the impedance calculation by the SIBC method is effective even if the antenna is very close to the metasurface.
Two BF-FDTD models.
Input impedance of dipole antenna vs. Nx at resonant frequency 15 GHz.
Input impedance of dipole antenna for h = 10 mm, 7.5 mm and 5 mm.
In order to investigate the scattering field from the metasurface, let us start with the exact field representation. The dyadic Green's functions for the dielectric slab backed by the PEC as shown in Figure 8 can be easily derived by the same way described in  or , and some formulas expressed by vector wave functions in cylindrical coordinate can be obtained. However, these expressions are inconvenient to use the scattering problem due to the rectangular patches shown in Figure 2. In this case, the rectangular coordinate representations of the dyadic Green's functions are available . Using the Green's functions, the electric field in the region of z > 0 is expressed by (8)where and are equivalent electric and magnetic surface currents on the surface S A enclosing the antenna, is the surface current on the conducting patch S P , as shown in Figure 8. is an electric type dyadic Green's function that contains a free space component and a scattered field component. is a magnetic type and is obtained by exchanging ε and μ contained in due to the duality of electromagnetic fields. Noted here that is different from the definition of the magnetic type Green's function in .
The far-field electric field is obtained by applying a saddle point method to equation (8), and is expressed as follows(9)where D 0 (θ, ϕ) is a vector directivity function in free space and can be calculated in a familiar way. On the other hand, is the directivity function of the scattered field component, and given as follows.(10)where R (TE) and R (TM) are reflection coefficients of the dielectric slab backed by the PEC for TE and TM plane wave incidences, respectively, and , . The vector functions and are given by(11)where .
Although the far field is usually calculated using the far-zone transformation formula in the time domain, a diffetent manner was adopted for efficient calculations in this paper. First, the antenna is excited by a sinusoidal voltage and then the amplitude and the phase of the equivalent currents , and are calculated at all points on the surfaces. The directivity function for an observation angle (θ, ϕ) is obtained by calculating the integrals (11) using the complex representations of equivalent currents.
When the SIBC is used, R (TE) and R (TM) are replaced by the next equations and the second term of in equation (11) is omitted.(12)where , , and the negative sign is adopted for the TM incidence.
In order to confirm the validity of the SIBC-FDTD method for the far-field calculation, the radiation patterns for the same structure described above, were calculated using the BF∞-FDTD method and SIBC-FDTD method.
The radiation pattern D θ (θ) in z − x plane is shown in Figure 9 for the case that the length of the antenna 2l is half-wavelength at 10 GHz and the height h is quarter-wavelength. A phase of the reflection coefficient is roughly estimated as 75 degrees from Figure 3, and hence the pattern is closely resemble the half-wavelength dipole antenna in shape even though h is is quarter-wavelength. The number of patches including in the computational space of the BF∞-FDTD method was set 60 × 60 which is 4 times larger than the impedance calculations. Both patterns are normalized by the maximum value of the SIBC-FDTD. It is found that the SIBC-FDTD agrees completely with BF∞-FDTD, and that the validity of the SIBC is numerically demonstrated.
Figure 10 shows the radiation pattern at 15GHz where 2l and h are half-wavelength and quarter-wavelength at this frequency respectively. In this geometry, the metasurface acts as the PMC and then a direct field from the antenna is reduced by the reflected wave from the metasurface near θ = 0. It is found that the noticeable discrepancy between the BF∞-FDTD and the SIBC-FDTD is observed. In order to investigate the cause, the surface current density along x axis on the interface at z = 0 is shown in Figure 11. The discrete peaks correspond to the current density on the conducting patch surfaces. Thus the current on the patch attenuates very slowly with respect to the distance from the antenna at the resonant frequency (15 GHz) but relatively rapidly decays at non-resonant frequency (10 GHz), and hence larger computational space is needed in the BF∞-FDTD calculation at 15GHz. Similar numerical error characteristics are hold even when the antennais extremely close to the interface as shown in Figure 12.
It is often that the BFf-FDTD method is applied to the far-field calculation in expectation of getting practically acceptable results as is the same with the input impedance. In order to confirm it, the radiation pattern of the same antenna as Figure 10 is shown in Figure 13. It is found that BFf-FDTD method is obviously incorrect in spite of the huge computational space. This failure is caused by a large reflection from the edge of the dielectric slab. It is also suggested that the surface wave is effectively excited in the dielectric slab and has an extremely large effect to the far-field.
Radiation patterns of half-wavelength dipole antenna at 10 GHz in z–x plane (h = λ/4).
Radiation patterns of half-wavelength dipole antenna at 15 GHz in z–x plane (h = λ/4).
Surface current density on interface along x axis.
Radiation patterns of half-wavelength dipole antenna in z–x plane (h = λ/10).
Radiation patterns of half-wavelength dipole antenna calculated by BFf-FDTD method.
The SIBC-based FDTD method has been investigated to calculate the properties of the dipole antenna located above the metasurface that consists of the infinite periodic conducting element on the flat dielectric substrate. First, the rational function interpolation of the given surface impedance data was introduced in order to obtain the expression that can be successfully incorporated into the FDTD method. Next, the input impedance of the horizontal dipole antenna located above the metasurface which consists of the rectangular PEC patch element, was calculated using the SIBC-FDTD method and two kinds of BF-FDTD methods, that is BF∞-FDTD and BFf-FDTD. It has been shown that the SIBC-FDTD agrees fairly well with the BF∞-FDTD even if the computational space is relatively small and that the numerical accuracy of the BFf-FDTD method is interior to the BF∞-FDTD method but is acceptable. Finally, it has been shown that this method can be applied to the radiation pattern calculation successfully by enlarging the computational space. It has also been shown that BFf-FDTD is inappropriate because the the unecessary large reflection occurs in the dielectric slab.
Since the computational results of the BF-FDTD method depend on the number of the unit cell contained in the computational space, it must be checked for every metasurface structure how large computational space is required for obtaining the accurate results, especially the radiation pattern. On the contrary, the scattering fields by the surface impedance sheet may hardly be affected by the dimension of the computational space even if the appropriate absorbing boundary such as the PML is used. This is an additional advantage of the SIBC-FDTD method. However, only the surface impedance at normal incidence is incorporated to the FDTD impedance calculation and the allocated positions of the magnetic field are different from the position of the electric field. Hence, the numerical error resulted from these inherent features are contained in the proposed SIBC-FDTD method. It would be the remaining problems to be solved for more accurate calculation. The paper has been treated the single dipole antenna, but is easily applicable to any other finite size antennas.
- Special Issues on electromagnetic application of photonic bandgap structure, progress in electromagnetics, PIER 41 (2003) [Google Scholar]
- Special Issue on metamaterial structures phenomena and applications, IEEE Trans. Microw. Theory Technol. 52 , 2 (2004) [CrossRef] [Google Scholar]
- F. Capolino, ed., Metamaterials handbook, Part I: Theory and phenomena of metamaterials, Pert II: Applications of metamaterials (CRC Press, Boca Raton, 2009) [Google Scholar]
- D. Sievenpiper, L. Zhang, R.F. Broas, N. Alexopolous, E. Yablonovitch, High-impedance electromagnetic surfaces with a forbidden frequency band, IEEE Trans. Microw. Theory Technol. 47 , 7–17 (1995) [Google Scholar]
- F. Yang, Y. Rahmat-Samii, Microstrip antenna integrated with electromagnetic band-gap (EBG) structure: A low mutual coupling design for array application, IEEE Trans. Antennas Propagat. 51, 2936–2946 (2003) [Google Scholar]
- A. Taflove, S.C. Hagness, Computational electrodynamics the Finite-Difference Time-Domain method , 3rd edn. (Artech House, Boston, 2005) [Google Scholar]
- K.S. Kunz, R.J. Luebbers, The finite difference time domain method for electromagnetics (CRC Press, Boca Raton, 1993) [Google Scholar]
- H. Yang, R. Mittra, FDTD modeling of metamaterials: theory and applications (Artech House, Boston, 2008) [Google Scholar]
- T. Uno, T. Arima, Y. He, FDTD method for computational electromagnetics − Fundamentals and practical application (Corona Publ. Co, Tokyo, 2016), in Japanese [Google Scholar]
- P. Harms, R. Mittra, W. Ko, Implementation of the periodic boundary-condition in the Finite-Difference Time-Domain algorithm for FSS structures, IEEE Trans. Antennas Propagat. 42 , 1317–1324 (1994) [CrossRef] [Google Scholar]
- R. Sigelmann, A. Ishimaru, Radiation from periodic structures excited by an aperiodic source, IEEE Trans. Antennas Propagat. 13 , 354–364 (1965) [CrossRef] [Google Scholar]
- F. Capolino, D.R. Jackson, D.R. Wilton, Fundamental properties of the field at the interface between air and a periodic artificial material excited by a line source, IEEE Trans. Antennas Propagat. 53 , 91–99 (2005) [CrossRef] [Google Scholar]
- F. Capolino, D.R. Jackson, D.R. Wilson, L.B. Felsen, Comparison of methods for calculating the field excited by a dipole near a 2-D periodic material, IEEE Trans. Antennas Propagat. 55 , 1644–1655 (2007) [CrossRef] [Google Scholar]
- J.G. Maloney, G.S. Smith, The use of surface impedance concepts in the finite-difference time-domain method, IEEE Trans. Antennas Propagat. 40 , 38–48 (1992) [CrossRef] [Google Scholar]
- J.H. Beggs, R.J. Luebbers, K.S. Yee, K.S. Kunz, Finite-Difference Time-Domain implementation of surface impedance boundary-conditions, IEEE Trans. Antennas Propagat. 40 , 49–56 (1992) [CrossRef] [Google Scholar]
- T. Arima, T. Uno, Meta-surface modeling method in FDTD analysis by introducing surface impedance boundary conditions, in 2017 IEEE Int. Conf. Computational Electromag. (ICCEM) (2017), 261–262 [CrossRef] [Google Scholar]
- T.B.A. Senior, J.L. Volakis, Approximate boundary conditions in electromagnetics (IEE Press, London, 1995) [CrossRef] [Google Scholar]
- I.V. Lindell, Methods for electromagnetic field analysis , 2nd edn. (IEEE Press, New York, 1995) [Google Scholar]
- S. Tretyakov, Analytical modeling in applied electromagnetics (Artech House, Boston, 2000) [Google Scholar]
- S.A. Tretyakov, C.R. Simovski, Dynamic model of artificial reactive impedance surfaces, J. Electromagn. Waves Appl. 17, 131–145 (2003) [CrossRef] [Google Scholar]
- D.B. Kuznetsov, J.E. SchuttAine, Optimal transient simulation of transmission lines, IEEE Trans. Circ. Syst. I-Fundam. Theory Appl. 43 , 110–121 (1996) [CrossRef] [Google Scholar]
- R.J. Luebbers, F. Hunsberger, K.S. Kunz, R.B. Standler, M. Schneider, A frequency-dependent finite-difference time-domain formulation for dispersive materials, IEEE Trans. Electromagn. Compat. 32 , 222–227 (1990) [Google Scholar]
- D.F. Kelly, R.J. Luebbers, Piecewise linear recursive convolution for dispersive medium using FDTD, IEEE Trans. Antennas Propagat. 44 , 792–797 (1996) [CrossRef] [Google Scholar]
- C.-T. Tai, Dyadic Green functions in electromagnetic theory , 2nd edn. (IEEE Press, New York, 1994) [Google Scholar]
- L.B. Felsen, N. Marcuvitz, Radiation and scattering of waves (Prentice-Hall Inc., Englewood Cliffs, NJ, 1973) [Google Scholar]
- T. Uno, Simplification of dyadic Green's function for plane stratified media and its application to lossless DNG slab, IEICE Trans. Commun. (Japanese edition) J89-B, 1661–1671 (2006) [Google Scholar]
Cite this article as: Toru Uno, Takuji Arima, Akihide Kurahara, FDTD modeling of nonperiodic antenna located above metasurface using surface impedance boundary condition, EPJ Appl. Metamat. 6, 17 (2019)
Radiation patterns of half-wavelength dipole antenna at 10 GHz in z–x plane (h = λ/4).
|In the text|
Radiation patterns of half-wavelength dipole antenna at 15 GHz in z–x plane (h = λ/4).
|In the text|
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.