Efficient parallel FDTD algorithm for modeling infinite graphene sheet simulations

Ground penetrating radar (GPR) forward is one of the geophysical research directions.Through the forward of geological model,the database of radar model can be enriched and the characteristics of typical geological radar echo images can be understood,which in turn can guide the data interpretation of GPR measured profile,thereby improving the GPR data interpretation level.In this article,the interpolating wavelet scale function by using iterative interpolation method is presented,and the derivative of scale function is used in spatial differentiation of discrete Maxwell equations. The forward modeling formula of GPR based on the interpolation wavelet scale method is derived by using fourth-order Runge-Kutta method (RK4) for calculating the higher time derivative.Compared with the conventional finite difference time domain (FDTD) algorithm based on the central difference method,the interpolation wavelet scale algorithm improves the accuracy of GPR wave equation in both space and time discretization.Firstly,the FDTD algorithm and the interpolation wavelet scale method are applied to the forward modeling of a layered model with analytic solution. Single channel radar data and analytical solution fitting indicate that the interpolation wavelet scale method has higher accuracy than FDTD,with the same mesh generation used.Therefore,auxiliary differential equation perfectly matching layer (ADE-PML) boundary condition is used on an interpolation wavelet scale,and the comparisons between reflection errors obtained using CPML (FDTD),RK4ADE-PML (FDTD),and RK4ADE-PML (interpolating wavelet scales) in a homogeneous medium model show that the absorption effect of RK4ADE-PML (interpolating wavelet scales) is better than the other two absorbing boundaries.Finally,interpolation wavelet scale method,with both UPML,FDTD and RK4ADE-PML loaded,is used for two-dimensional GPR forward modeling,showing good absorption effect for evanescent wave.From all the experimental results,the following conclusions are obtained.1) Using the derivative of the interpolating wavelet scale function instead of central difference schemes for the spatial derivative discretization of Maxwell equations and time derivative calculated using the fourth-order Runge Kutta algorithm,the interpolating wavelet scale algorithm has higher accuracy than regular FDTD algorithm due to the improvement in the spatial and time accuracy of GPR wave equation.2) The best absorption layer parameters of interpolating wavelet scale RK4ADE-PML are selected, when the maximum value of the reflection error is the minimum.The maximum reflection error can reach-93 dB,which increases 20 dB compared with that of UMPL boundary in FDTD algorithm.And the higher simulation accuracy of interpolating wavelet scale algorithm than FDTD algorithm is confirmed after calculating single channel radar data.3) Comparing wave field snapshots of GPR forward modeling,radar pictures from wide-angle method and section method indicates that interpolating wavelet scale RK4ADE-PML reduces reflection error of absorption boundary,improves both spatial and time accuracy,is more effective than UPML boundary in eliminating false reflection of large angle incidence, and has better absorption effect for evanescent wave and low-frequency wave.

Discusses a computational procedure, based on the finite difference time domain (FDTD) algorithm, for extremely low-frequency (1 Hz) modeling of electromagnetic scattering by buried objects in lossy sedimentary layers under sea water. The FDTD algorithm is found to be accurate for modeling such objects tens of meters away from a constant current dipole source of 1 A in sea water. For validation of the extremely low-frequency and lossy FDTD modeling, the computed FDTD results are compared with those calculated by using integral equation techniques.

In this paper we discuss a low-frequency computational procedure based on the Finite Difference Time Domain (FDTD) algorithm, for numerical modeling of electromagnetic scattering by buried objects in sediment layers under sea water. The FDTD algorithm is found to be accurate for modeling buried objects in sediment layers tens of meters away from a constant current dipole source of 1A in sea water. For validation of the low-frequency FDTD modeling, the computed FDTD results are compared with those calculated by using analytic expressions and integral,equation techniques. In this paper we also present a technique for detecting conductivity anomalies in sediments, e.g., a buried object in sedimentary layers under sea water, by using the neural network approach. The electric field values are used as the inputs to the neural network and the associated conductivities are used as the targets. The neural network is then trained to associate these conductivities and field values. It is shown in this paper...

A proper dielectric coating can reduce the electromagnetic scattering of the conducting object significantly so that it cannot be detected by the radar. However, when the object contains a thin coating, fine grid size is needed to discretize the thin coating in the conventional finite-difference time-domain (FDTD) algorithm, which increases the amount of memory and computational time significantly. To overcome this dilemma, we present a transformation optics-based FDTD (TO-FDTD) algorithm to accelerate the solution of electromagnetic scattering from objects with thin dielectric coatings. Two kinds of novel TO-FDTD models are proposed in this paper for a coated cylinder and a coated arbitrary polygonal cylinder, respectively. Through coordinate transformation, the size of the object remains unchanged while its thin coating is enlarged to a thicker one, meaning that it can be simulated by the FDTD algorithm with uniform coarse grids instead of fine grids. The transformed material parameters become inhomogeneous and anisotropic in the transformed region, which can be obtained by solving a Jacobian transformation matrix. We then develop a stable FDTD algorithm for solving anisotropic Maxwell’s equations. Bistatic scatterings of coated cylinders and a coated polygonal cylinder are solved by the TO-FDTD algorithm proposed in this paper, respectively. The result of the TO-FDTD algorithm matches well with the exact value and the result of the commercial software Comsol. The computational efficiency and accuracy of the proposed TO-FDTD algorithm are validated by numerical experiments. Numerical results show that the TO-FDTD algorithm has higher computational accuracy than the conventional FDTD algorithm that fails to simulate the absorbing property of the coating, when the same coarse grid size is used in the simulation. Under the same level of accuracy, the proposed TO-FDTD method can improve the computational efficiency by 62-63 times than the conventional FDTD method with fine grids in the simulations in the paper.

As with the number of cancer increases year by year, so it is important to be found and treated earlier. With biological cells and tissues are sensitive to infrared and visible light, cell morphology and physical structure of the optical properties can easily obtain, we can provide theoretical basis for the early diagnosis of cancer by observing the difference of optical properties between normal and cancerous cells. Compared with Mie scattering theory, finite difference time domain (FDTD) algorithm can analyze any complex structure model. In this paper we use mathematical modeling method to establish the single cell mathematical model and with finite difference time domain algorithm to simulate the propagation and scattering of light in the biological cells, you can calculate the scattering of electromagnetic field distribution at anytime and anywhere. With radar cross section (RCS) to measure the results of the scattering characteristics. Due to the difference between normal cells and cancerous cells are embodied in cell shape, size and the refractive index, through the simulation we can get different cell parameters of light scattering information, Find out the cell parameters change the changing rule of the influence on the scattering characteristics and find out change regularity of scattering characteristics. These data can judge very accurate of the cells is normal or cancerous cells.

In order to study the parameters of cable coupling characteristics, we calculated the voltage and current distribution at each point of the transmission line numerically. An iterative solution method based on Theory of Multiconductor Transmission Line (MTL) and finite difference time domain (FDTD) is used. Firstly, the appropriate FDTD parameters are selected, and in each iteration step, the voltage and current along the line are calculated using the Finite Difference Time Domain (FDTD) algorithm, and the voltage and current distribution is updated accordingly. When the convergence condition is satisfied, the solution of the transmission line equation is obtained. The main advantage of this method is that the inherent parallel algorithm can reduce the computation time and obtain more detailed numerical values within the line.

A three-dimensional (3-D) finite difference time domain (FDTD) algorithm with perfectly matched layer (PML) absorbing boundary condition (ABC) is presented for general inhomogeneous, dispersive, conductive media. The modified time-domain Maxwell's equations for dispersive media are expressed in terms of coordinate-stretching variables. We extend the recursive convolution (RC) and piecewise linear recursive convolution (PLRC) approaches to arbitrary dispersive media in a more general form. The algorithm is tested for homogeneous and inhomogeneous media with three typical kinds of dispersive media, i.e., Lorentz medium, unmagnetized plasma, and Debye medium. Excellent agreement between the FDTD results and analytical solutions is obtained for all testing cases with both RC and PLRC approaches. We demonstrate the applications of the algorithm with several examples in subsurface radar detection of mine-like objects, cylinders, and spheres buried in a dispersive half-space and the mapping of a curved interface. Because of their generality, the algorithm and computer program can be used to model biological materials, artificial dielectrics, optical materials, and other dispersive media.

The finite difference time domain (FDTD) algorithm has become a popular-technique with which to solve many different electromagnetic problems. However, for real problems the mesh size becomes large and as a result the computation time on a workstation is in the region of tens of hours. In fact, in some cases there will be insufficient memory available on a workstation to solve large problems. We describe a general method based on distributed computing to reduce both the memory requirements and the computation time of the FDTD algorithm. An implementation of the FDTD algorithm on a network of workstations is then evaluated. The implementation is shown to give a peak speed up of 7.5 on a network of 8 workstations over the sequential computation time. >

In this paper, the well known Finite Difference Time Domain (FDTD) algorithm is utilized to design and analyze a number of planar printed microwave integrated circuits (MICs). Such designs and analyses are based on a developed FDTD computer solver. Prior to the designs, the solver has been validated through comparisons of its computed results with the ones which came from some analyzed structures previously published in other literature, and a good agreement has been observed. The results synthesized by the proposed designs lead to a strong belief that the FDTD is indeed a robust, versatile and powerful algorithm in solving various electromagnetic problems.

This paper introduces a parallel FDTD (Finite Difference Time Domain) algorithm based on MPI (Message Passing Interface) parallel environment and NUMA (Non-Uniform Memory Access) architecture workstation. The FDTD computation is carried out independently in local meshes in each process. The data are exchanged by communication between adjacent subdomains to achieve the FDTD parallel method. The results show the consistency between serial and parallel algorithms, and the computing efficiency is improved effectively.

Highly integrated RF-frontends are a key component for K/Ka-band satellite communications on-the-move, multimedia entertainment systems, 5G mobile communication as well as 60 GHz broadband home access equipment and backbone networks. Designing such complex modules poses real challenges. These RF-frontends are not only composed of a large number of antenna elements but also of highly integrated complex multilayer feeding structures, RF-transitions, RF-components and power combiners. This paper presents how the XPU technology, which is a smart implementation of the Finite Difference Time Domain (FDTD) algorithm on modern CPU architectures, enables a time-efficient design of frontend modules for Ka-band and WiGig communication, speeding up the whole development process.

A new broadband planar offset antenna is described based on the T-shaped microstrip feeding structure in this paper. This antenna is optimized by the finite difference time domain (FDTD) and genetic algorithms (GA). The results show the impedance bandwidth of this antenna is from 1.864 to 6.643 GHz for S11 les -10 dB (VSWRles2), much higher than that of usual slot patch antenna, the far-filed pattern almost remain the same as that of the T-shaped microstrip feeding antenna

New formulations based on the state-space approach are presented for implementing the anisotropic perfectly matched layer into the finite difference time domain (FDTD) algorithm. The formulations are simple and independent from the material properties of the FDTD computational domain. Numerical tests are included to validate the proposed formulations.

A novel implementation of the finite difference time domain (FDTD) algorithm to analyze periodic structures over a wideband with oblique incident waves is presented. Based on the field transformation treatment, the method splits each component of the transformed fields of the electric displacement vector and the magnetic induction intensity vector into three parts in a simple manner. This treatment is independent of material and easily implemented, so that it could be utilized in a wide range of periodic structure analysis. The advantage of the method is demonstrated through two numerical examples. One is Jerusalem cross frequency selective surfaces (JCFSS) with anisotropic lossy media in the millimeter wave band, and excited by the TM wave at θ = 20° and the TE wave at θ = 45°, respectively. Another is a photonic periodic structure with rectangular holes for the visible and infrared bands, where it acts as a dispersive medium, with the incident TE plane wave at θ = 30°. The results are in good agreement with those of CST Microwave Studio, and thereby the validity of the method is verified.

SUMMARY A novel stable anisotropic finite-difference time-domain (FDTD) algorithm based on the overlapping cells is developed for solving Maxwell's equations of electrodynamics in anisotropic media with interfaces between different types of materials, such as the interface between anisotropic dielectrics and dispersive medium or perfect electric conductor (PEC). The previous proposed conventional anisotropic FDTD methods suffer from the late-time instability due to the extrapolation of the field components near the material interface. The proposed anisotropic overlapping Yee FDTD method is stable, as it relies on the overlapping cells to provide the collocated field values without any interpolation or extrapolation. Our method has been applied to simulate electromagnetic invisibility cloaking devices with both anisotropic dielectrics and PEC included in the computational domain. Numerical results and eigenvalue analysis confirm that the conventional anisotropic FDTD method is weakly unstable, whereas our method is stable. Copyright © 2013 John Wiley & Sons, Ltd.