Improved hybrid implicit–explicit finite-difference time-domain method with higher accuracy and user-defined stability condition

A new explicit and nondissipative finite-difference time-domain (FDTD) method in two and three dimensions is proposed for the relaxation of the Courant condition. The third-degree spatial difference terms with second- and fourth-order accuracies are added with coefficients to the time-development equations of FDTD(2,4). Optimal coefficients are obtained by a brute-force search of the dispersion relations, which reduces phase velocity errors but satisfies the numerical stabilities as well. The new method is stable with large Courant numbers, whereas the conventional FDTD methods are unstable. The new method also has smaller numerical errors in the phase velocity than conventional FDTD methods with small Courant numbers.

This paper develops a perfectly matched layer (PML) for the hybrid implicit explicit finite-difference time-domain (HIE-FDTD) method. It is demonstrated that the inclusion of the PML is a simple extension of the original HIE-FDTD method. It doesn't affect the weakly condition stability and computational accuracy. Simulation results show that when the time step size satisfies with the time stability condition, the reflection error of the HIE-FDTD method is below −90dB, which demonstrates good absorbing performance of the PML method.

In this paper, a hybrid algorithm of explicit finite-element time-domain (FETD) and finite-difference time-domain (FDTD) methods with a large time step is proposed. The FETD method based on the unstable modes filtering (UMF-FETD) is adopted in the region involving fine structures while the FDTD method is employed in the other region. The proposed hybrid algorithm reduces the scale of the system matrix in UMF-FETD and guarantees that the whole calculation region adopts the uniform, large time step. The result proves the good ability of the proposed hybrid algorithm.

In this work, they propose a one-step leapfrog hybrid implicit-explicit finite-difference time-domain (HIE-FDTD) method for body-of-revolution (BOR). Meanwhile, its Convolutional Perfect Matched Layer (CPML) absorbing boundary condition is implemented. In this method, the implicit difference is applied in the angular direction. All the resultant updating equations are still explicit. However, the stability condition of the proposed method is relaxed. The analytical analysis shows that its time step is only determined by the smaller one of spatial increments Δρ and Δz. A scattering example is provided to demonstrate the new algorithm. At the same time, the relative of reflection error of the CPML is given with comparisons of Mur.

This paper presents a new 3-D hybrid implicit explicit (HIE) finite-difference time-domain (FDTD) method with a stability criterion relaxed from the discretization space steps in one direction. The proposed method is particularly useful and efficient for electromagnetic (EM) simulation of structures having fine details in one Cartesian direction such as thin shields and thin slots on a metallic enclosure whose EM analysis are popular and important in electromagnetic compatibility. In this algorithm, unlike common FDTD methods in which both electric and magnetic fields are updated, only the electric field needs to be updated during the time-marching process of the algorithm; hence, it is called as a “single-field” FDTD method. This could be achieved by applying the Crank-Nicolson (CN) scheme only to the electric field, while in earlier HIE-FDTD (and also unconditionally stable FDTD) methods, the CN scheme is applied to both electric and magnetic fields. The proposed algorithm consists of two implicit and three explicit updating equations that can be solved simultaneously (in parallel) by applying multithreading. Thanks to few and simple updating equations of the proposed algorithm, it is computationally more efficient than other HIE-FDTD methods in terms of runtime. The high computational efficiency and accuracy of the proposed method are demonstrated by producing several numerical examples and providing comparison to the results obtained using methods available in the literature.

New version of FDTD with GFM program in Schwarzschild space–time

The nearly perfectly matched layer (NPML) is incorporated into a novel four-steps weakly conditionally stable hybrid implicit explicit finite difference time domain (HIE-FDTD) method in three-dimensional (3-D) domains in this paper. By dividing the matrix which obtained by the Maxwell's equation into four sub-matrices with the symmetric operation and uniform division and incorporating NPML into it, the time step of the proposed algorithm is divided into four sub-steps, then the formulation of the proposed method is generated. Moreover, numerical results are carried out for different Courant-Friedrichs-Lewy (CFL) numbers at different positions with different NPML layers in 3-D domains, and it has been obtained that the presented formulation is accurate as well as weakly conditionally stable.

In this letter, the current source implementations are presented for the hybrid implicit–explicit finite-difference time-domain (HIE-FDTD), the improved HIE-FDTD, and the leapfrog HIE-FDTD methods. The asymmetry and the field errors of these HIE-FDTD methods are investigated for all possible values of the time index parameters. The proposed implementations give very low asymmetry and field errors for all the three HIE–FDTD methods. The field errors given by these HIE-FDTD methods are lower than those given by the alternating direction implicit finite-difference time-domain (ADI-FDTD) method.

The finite difference time domain (FDTD) method has been proposed and used for sound field simulation. To reproduce actual sound wave propagation in sound field simulations, it is necessary to apply the radiation characteristics. With the FDTD method, radiation characteristics can be applied by setting sound pressure in a dense grid arrangement. However, conventional techniques for capturing radiation characteristics use a sparse array of microphones and are considered insufficient for the FDTD simulation. Furthermore, the technique required to apply captured acoustic signals in a dense grid arrangement with the FDTD method has not been considered. In this paper, we propose a novel hardware and software system that captures the radiation characteristics for a dense grid arrangement and applies them to the FDTD method, while controlling the sound wave propagation with the non-propagation region. The proposed system produces the average differences from measured values of sound pressure, propagation time, center frequency, and log-spectral distortion of 1.8 dB, 0.04 ms, 700 Hz, and 3.5 dB, respectively, which is more accurate than the conventional techniques. The result shows that this system is useful for improving the accuracy of sound wave propagation reproduction with the sound field simulation.

Various methods are available at the disposal for computing the solution of electromagnetic problems expeditiously. Electromagnetic scattering problems can be classified into two broader categories named as surface and volume problems. These problems can be solved using differential and integral equations. Differential equations based techniques use Finite Element Method (FEM) and Finite Difference Time Domain (FDTD) methods and major integral formulation based technique is Method of Moment (MOM). Any set of linear equations for a given system can be solved using direct or iterative methods. Direct solvers confront many computational complexities and use large amount of memory while iterative solvers reduce computational complexities and use less amount of memory. In this research work some integral equations based solver schemes such ass Finite Element Method (FEM), Finite Difference Time Domain Method (FDTD) and Method of Moments (MOM) for surface and volume scatters are discussed and results are compared. The comparison is with respect to complexity and memory requirements for different solution schemes.

The finite-difference time-domain (FDTD) method as realized in the Yee algorithm is very useful for solving the analytically intractable diffraction problems that arise in the design of diffractive optical elements (DOE). The error, e, of the Yee algorithm is rather large. We introduce a new high accuracy version of the Yee algorithm based on nonstandard finite differences. We have verified its performance by comparing with known analytical solutions. We have used this algorithm to study the near field (0 - 10 wavelengths) diffraction patterns of three binary phase gratings. Grating (1) has a constant period, while grating (2) has a variable one. Grating (3) is generated by interlacing grating (1) with (2). The period of grating (1) is on the order of a wavelength. Some of these computed fields display interesting and possibly useful Fourier spectra.

PV systems are considered one of the most important and popular forms of renewable energy sources worldwide. Due to the installation of PV modules on rooftops and outdoor areas, they are frequently subjected to lightning strikes, leading to a deterioration in PV systems. To provide an efficient lightning protection system, the accurate modeling of PV systems during lightning strikes is required. The main objective of this paper is to assess the transient behavior of PV systems when struck by lightning, where the PV system is modeled using the partial element equivalent circuit (PEEC) and finite difference time-domain (FDTD) methods. The PEEC method is one of the accurate methods used to represent the PV system by its equivalent RLC circuit to study lightning transients. Moreover, the FDTD method provides high accuracy when studying lightning transients by analyzing the electric and magnetic fields inside the PV system. The PEEC and FDTD methods are used to calculate the transient overvoltages at different points of the PV system. Then, a small-scale PV system is tested in the laboratory using an impulse generator representing the lightning strike to verify the correctness of the results obtained by simulation. Most of the existing research works introduce low accuracy in PV system modeling during lightning transients by neglecting the effect of the mounting system, metal frame, and mutual coupling, which are considered in this work. Moreover, the high-frequency nature of lightning current is considered when modeling the elements of the grounding system. Verifying the results by laboratory experiments can guarantee the application of the simulation model on large-scale PV systems.

We have applied the finite difference time domain (FDTD) method with the supersymmetric quantum mechanics (SUSY-QM) procedure to determine excited energies of one dimensional quantum systems. The theoretical basis of FDTD, SUSY-QM, a numerical algorithm and an illustrative example for a particle in a one dimensional square-well potential were given in this paper. It was shown that the numerical results were in excellent agreement with theoretical results. Numerical errors produced by the SUSY-QM procedure was due to errors in estimations of superpotentials and supersymmetric partner potentials.

A finite difference time domain (FDTD) method in combination with the perfectly matched layer (PML) is used for dispersive media. The FDTD method is derived using the convolution relationship between the current density J and the electric field E. The high accuracy of the FDTD method is confirmed by computing the reflection and transmission coefficients for a one dimensional nonmagnetized plasma slab.

The finite difference time domain method and finite element methods have proven to be highly effective techniques for the full wave analysis of microwave circuits and resonators. The finite difference time domain (FDTD) method is highly efficient, explicit and parallelizable. Finite element methods have been used in the frequency domain and in the time domain using implicit methods. Implicit methods require a solution of a linear system of equations at each time step. For the finite element time domain (FETD) algorithm to be competitive with the FDTD, the number of time iterations required for a simulation must be significantly less. An unconditionally stable FETD algorithm is used. The variational formulation is based on the vector wave equation. Unconditional stability is achieved by using a modified Newmark-Beta time integration technique. In addition, this algorithm is implemented on parallel message passing computers using a finite element tearing and interconnecting scheme. Some numerical examples based on cavity resonator problems are presented to demonstrate the efficiency and accuracy of the method and the parallel efficiency of the algorithm.