International Journal of Antennas and Propagation

Volume 2018, Article ID 3046838, 4 pages

https://doi.org/10.1155/2018/3046838

## An Unconditionally Stable Cylindrical FDTD Method to Analyze the EM Ground Wave Propagation

^{1}Key Laboratory of Radar Imaging and Microwave Photonics, Ministry of Education, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China^{2}National Key Laboratory on Electromagnetic Environmental Effects and Electro-optical Engineering, Army Engineering University, Nanjing 210007, China

Correspondence should be addressed to Zheng-Yu Huang; moc.361@jnyzgnauh

Received 24 November 2017; Accepted 26 March 2018; Published 8 May 2018

Academic Editor: Luciano Tarricone

Copyright © 2018 Zheng-Yu Huang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We extended the unconditionally stable associated Hermite (AH) FDTD method to the cylindrical coordinate system for analyzing the electromagnetic (EM) ground wave propagation. With AH domain differentiate operator technology and paralleling-in-order solution scheme, the cylindrical time-domain Maxwell equations with complex frequency shifted perfectly matched layer (CFS-PML) are transformed to AH domain generating a five-point banded equation, which also unified the formula from the central axis of cylinder. The numerical results represent the accuracy and efficiency of the proposed method when comparing for the conventional FDTD method.

#### 1. Introduction

In the finite-difference time-domain (FDTD) method [1], the time step size is limited by the Courant-Friedrich-Levy (CFL) stability condition [2]. So for the problems with fine structures, it should choose a small time step and then a large number of marching-on-in-time steps for simulation, which is a time-consuming process. To eliminate the stability condition for fine structure, many unconditionally stable (US) methods were proposed, such as alternating direction implicit (ADI) method [3] and orthogonal function-based methods [4–9]. The orthogonal function-based methods can be classified into two main categories: the marching-on-in-order scheme and the paralleling-in-order scheme. The former is based on weighted Laguerre polynomials (WLP), and the latter is based on associated Hermite (AH) functions. For the latter, all of the works are based on in Cartesian grid previously.

In many applications, such as EM ground wave propagation from a lightning channel, a cylindrical coordinate is preferred. Here, we extend the AH FDTD method to cylindrical coordinates system to analyze this case. The time-domain Maxwell equations in cylindrical coordinates are transformed to AH domain using a AH differential operator technology [6]. Then a five-point banded equation is derived for calculation of expanding coefficients. And a paralleling-in-order scheme [6] is used for indirectly but efficiently solution. For the CFS-PML is wildly used to effectively attenuate evanescent waves [10, 11], we apply it to AH FDTD method and make a unified formula including a treatment of central axis of cylinder when following the Ampere law. One should note that although CFS-PML is used in [6, 9], it is implemented in Cartesian coordinate system. In this paper, we firstly extended it to cylindrical coordinate system, which might give AH FDTD method much more convenient to analyze the scattering problems for the cylindrical shape-based objective. Finally, a numerical example is used to validate the accuracy and efficiency of the proposed method.

#### 2. Formulation

With the lossy dielectric media, the 2D cylinder Maxwell equations with CFS-PML can be written as [11] where , , , and are the permeability, the magnetic conductivity, the permittivity, and the electric conductivity of the media, respectively. and are the frequency-domain coordinate-stretching variables, which are defined as and , where , , , , , and are the respective PML parameters [11]. Similar to [6], a set of AH basis functions is chosen to expand all of the field components (1), (2), and (3). is the time-scaling parameter, and is the time-translating parameter. Then, a temporal Galerkin testing procedure is used to eliminate the time variables. By applying central difference scheme, we can transform (1), (2), and (3) to AH domain as where

The matrix above is AH domain differential operator [6], which can be considered as the similar time-domain operator or frequency-domain operator and respective identity matrix is for 1. By assembling (4), (5), and (6) and eliminating the electric field components, we can derive a five-point banded matrix equation with as where

When for the boundary condition of central axis, it needs to be treated separately owing to singularities. According to Ampere’s law [12],

Then, the differential form for (14) can be obtained as follows:

Converting it to the AH domain by using the AH domain differential operator and performing discretization for

By applying (16) to (5), the modified magnetic field (12) in central axis can be derived as follows:

After the central axis formula to modify (12) is updated, we can apply the lower-upper (LU) decomposition procedure to calculate this equation. Also, we can perform paralleling-in-order solution scheme [6] for (12). If every magnetic field variable in the entire computational space is calculated, the results of electronic field variables and can be gotten from (4) and (5). In addition, only the lossy medium is considered above, while for the case of dispersive medium, the formulation derivation is almost the same as before, except for a slight modification of (8) with being replaced by its AH domain form from [7].

#### 3. Numerical Verification

To validate the performance of the proposed AH FDTD method, a numerical example with an EM ground wave propagation problem is presented, as shown in Figure 1. The computational domain is discretized into 50 × 50 cells with , including an axial symmetric boundary condition for the left boundary and the other three boundaries terminated by 10 PML layers. The background material is set to be earth with medium parameters: , , , and . The PML parameters are scaled following the expressions in [11]. where represents the interface between FDTD and PML grids, is the thickness of the PML, and is the order of the polynomial. A Gaussian pulse source modulated by sinusoidal signal is located at the grid (0, 27). where , , and . Two observation points (20, 23) and (38, 38) are chosen to record the results.