Finite difference time domain method forward simulation of complex geoelectricity ground penetrating radar model

The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.

Numerical Solutions of Finite Well in Two Dimensions Using the Finite Difference Time Domain Method

The higher excited states for two dimensional finite rectangular well potential are calculated numerically,by solving the Schrödinger equation using the finite difference time domain method.Although,this method is suitable to calculate the ground state of the quantum systems,it has been improved to calculate the higher excited states directly.The improvement is based on modifying the iterative process involved in this method to include two procedures.The first is known as cooling steps and the second is known as a heating step.By determining the required length of the cooling iteration steps using suitable excitation energy estimate,and repeating these two procedures using suitable initial guess function for sufficient times.This modified iteration will lead automatically to the desired excited state.In the two dimensional finite rectangular well potential problem both of the suitable excitation energy and the suitable initial guess wave function are calculated analytically using the separation of variables technique.

Design of the sample cell in near-field surface-enhanced Raman scattering by finite difference time domain method

The rational design of the sample cell may improve the sensitivity of surface-enhanced Raman scattering （SERS） detection in a high degree. Finite difference time domain （FDTD） simulations of the configuration of Ag film-Ag particles illuminated by plane wave and evanescent wave are performed to provide physical insight for design of the sample cell. Numerical solutions indicate that the sample cell can provide more ＂hot spots＂ and the massive field intensity enhancement occurs in these ＂hot spots＂. More information on the nanometer character of the sample can be got because of gradient-field Raman （GFR） of evanescent wave. OCIS codes： 290.5860, 240.0310, 240.6680, 999.9999 （surface-enhanced Raman scattering）.

Finite-difference time-domain studies of low-frequency stop band in superconductor-dielectric superlattice

The transmission coefficients of electromagnetic （EM） waves due to a superconductor-dielectric superlattice are numerically calculated. Shift operator finite difference time domain （SO-FDTD） method is used in the analysis. By using the SO-FDTD method, the transmission spectrum is obtained and its characteristics are investigated for different thicknesses of superconductor layers and dielectric layers, from which a stop band starting from zero frequency can be apparently observed. The relation between this low-frequency stop band and relative temperature, and also the London penetration depth at a superconductor temperature of zero degree are discussed, separately. The low-frequency stop band properties of superconductor-dielectric superlattice thus are well disclosed.

Finite difference time domain analysis of two-dimensional surface acoustic wave piezoelectric phononic crystals at radio frequency

The finite-difference time-domain （FDTD） method is proposed for analyzing the surface acoustic wave （SAW） propagation in two-dimensional （2D） piezoelectric phononic crystals （PCs） at radio frequency （RF）, and also experiments are established to demonstrate its analysis result of the PCs＇ band gaps. The FDTD method takes the piezoelectric effect of PCs into account, in which periodic boundary conditions are used to decrease memory/time consumption and the perfectly matched layer boundary conditions are adopted as the SAW absorbers to attenuate artificial reflections. Two SAW delay lines are established with/without piezoelectric PCs located between interdigital transducers. By removing several echoes with window gating function in time domain, delay lines transmission function is achieved. The PCs＇ transmission functions and band gaps are obtained by comparing them in these two delay lines. When Aluminum/128°YX-LiNbO3 is adopted as scatter and substrate material, the PCs＇ band gap is calculated by this FDTD method and COMSOL respectively. Results show that computational results of FDTD agree well with experimental results and are better than that of COMSOL.

Full-Wave Analysis of Slotline Using Time-Domain Finite-Difference Method

The transmission and dispersive characteristics of slotline are calculated in this paper. The tail of Gaussion pulse is improved because a modified dispersive boundary condition (DBC) is adopted. It leads to a reduction in computer memory requirements and computational time. The computational domain is greatly reduced to enable performance in personal computer. At the same time because edges of a boundary and summits are treated well, the computational results is more accurate and more collector.

Fast division-free parallel structure for convolution perfectly matched layer in finite difference time domain method

Parallel acceleration of convolution perfectly matched layer （CPML） algorithm suffers from massive division operation which is widely accepted as one of the most expensive operations for the equipment such as graphic processing unit （GPU）, field programmable gate array （FPGA） etc. In pursuit of higher efficiency and lower power consumption, this article revisited the CPML theory and proposed a new fast division-free parallel CPML structure. By optimally rearranging the CPML inner iteration process, all the division operators can be eliminated and replaced by recalculating the related field updating coefficients offline. Experiments show that the proposed division-free structure can save more than 50% arithmetic instructions and 25% execution time of the traditional parallel CPML structure without any accuracy loss.

A Novel Absorbing Boundary Condition for the Frequency-Dependent Finite-Difference Time-Domain Method

A new absorbing boundary condition (ABC) for frequency dependent finite difference time domain algorithm for the arbitrary dispersive media is presented. The concepts of the digital systems are introduced to the (FD) 2TD method. On the basis of digital filter designing and vector algebra, the absorbing boundary condition under arbitrary angle of incidence are derived. The transient electromagnetic problems in two dimensions and three dimensions are calculated and the validity of the ABC is verified.

High Precision Time Domain Forward Modeling for Crosshole Electromagnetic Tomography

To improve the resolution of crosshole electromagnetic tomography, high precision of forward modeling is necessary. A pseudo-spectral time domain （PSTD） forward modeling was used to simulate electromagnetic wave propagation between two boreholes. The PSTD algorithm is based on the finite difference time domain （FDTD） method and uses the fast Fourier transform （FFT） algorithm for spatial derivatives in Maxwell＇s equations. Besides having the strongpoint of the FDTD method, the calculation precision of the PSTD algorithm is higher than that of the FDTD method under the same calculation condition. The forward modeling using the PSTD method will play an important role in enhancing the resolution of crosshole electromagnetic tomography.

Orthogonal experiment design of EMI of security monitoring system in coal mines

Security monitoring system of coal mines is indispensable to ensure the safe and efficient production of colliery. Due to the special and narrow underground field of the coal mine, the electromagnetic interference can cause a series of misstatements and false positives on the monitoring system, which will severely hamper the safe production of coal industry. In this paper, first, the frequency characteristics of the interference source on the power line are extracted when equipment runs normally. Then the finite difference time domain method is introduced to analyze the effects of the electromagnetic interference parameters on the security monitoring signal line. And the interference voltage of the two terminal sides on the single line is taken as evaluating indexes. Finally, the electromagnetic interference parameters are optimized by orthogonal experimental design based on the MATLAB simulation on the normal operation of equipment.

ANALYSIS OF FDTD TO UPML FOR MAXWELL EQUATIONS IN POLAR COORDINATES

An FDTD system associated with uniaxial perfectly matched layer（UPML） for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-boundary value problems of the transverse magnetic（TM） mode to Maxwell＇s equations is obtained by Yee＇s algorithm,and the open domain of the scattering problem is truncated by a circle with a UPML.Besides,an artificial boundary condition is imposed on the outer boundary of the UPML.Afterwards,stability of the FDTD system on the truncated domain is established through energy estimates by the Gronwall inequality.Numerical experiments are designed to approve the theoretical analysis.

Photocurrent improvement of an ultra-thin silicon solar cell using the localized surface plasmonic effect of clustering nanoparticles

The cluster-shaped plasmonic nanostructures are used to manage the incident light inside an ultra-thin silicon solar cell.Here we simulate spherical,conical,pyramidal,and cylindrical nanoparticles in a form of a cluster at the rear side of a thin silicon cell,using the finite difference time domain(FDTD)method.By calculating the optical absorption and hence the photocurrent,it is shown that the clustering of nanoparticles significantly improves them.The photocurrent enhancement is the result of the plasmonic effects of clustering the nanoparticles.For comparison,first a cell with a single nanoparticle at the rear side is evaluated.Then four smaller nanoparticles are put around it to make a cluster.The photocurrents of 20.478 mA/cm2,23.186 mA/cm2,21.427 mA/cm2,and 21.243 mA/cm2 are obtained for the cells using clustering conical,spherical,pyramidal,cylindrical NPs at the backside,respectively.These values are 13.987 mA/cm2,16.901 mA/cm2,16.507 mA/cm2,17.926 mA/cm2 for the cell with one conical,spherical,pyramidal,cylindrical NPs at the backside,respectively.Therefore,clustering can significantly improve the photocurrents.Finally,the distribution of the electric field and the generation rate for the proposed structures are calculated.

Absorption enhancement and sensing properties of Ag diamond nanoantenna arrays

Noble metal nanoantenna could effectively enhance light absorption and increase detection sensitivity. In this paper,we propose a periodic Ag diamond nanoantenna array to increase the absorption of thin-film solar cells and to improve the detection sensitivity via localized surface plasmon resonance. The effect of nanoantenna arrays on the absorption enhancement is theoretically investigated using the finite difference time domain（FDTD） method with manipulating the spectral response by geometrical parameters of nanoantennas. A maximum absorption enhancement factor of 1.51 has been achieved in this study. In addition, the relation between resonant wavelength（intensity reflectivity） and refractive index is discussed in detail. When detecting the environmental index using resonant wavelengths, a maximum detection sensitivity of about 837 nm/RIU（refractive index unit） and a resolution of about 10-3RIU can be achieved. Moreover, when using the reflectivity, the sensitivity can be as high as 0.93 AU/RIU. Furthermore, we also have theoretically studied the effectiveness of nanoantennas in distinguishing chemical reagents, solution concentrations, and solution allocation ratios by detecting refractive index. From the results presented in this paper, we conclude that this work might be useful for biosensor detection and other types of detections.

Perfectly matched layer implementation for ADI-FDTD in dispersive media

Alternating direction implicit finite difference time domain （ADI-FDTD） method is unconditionally stable and the maximum time step is not limited by the Courant stability condition, but rather by numerical error. Compared with the conventional FDTD method, the time step of ADI-FDTD can be enlarged arbitrarily and the CPU cost can be reduced. 2D perfectly matched layer （PML） absorbing boundary condition is proposed to truncate computation space for ADI-FDTD in dispersive media using recursive convolution（RC） method and the 2D PML formulations for dispersive media are derived. ADI-FDTD formulations for dispersive media can be obtained from the simplified PML formulations. The scattering of target in dispersive soil is simulated under sine wave and Gaussian pulse excitations and numerical results of ADI-FDTD with PML are compared with FDTD. Good agreement is observed. At the same time the CPU cost for ADI-FDTD is obviously reduced.

Design and analysis of multi-channel drop filter using dual L defected hexagonal photonic crystal ring resonator

In this paper,dual L defected hexagonal Photonic Crystal Ring Resonator(PCRR)using Channel Drop Filter(CDF)is designed for Coarse Wavelength Division Multiplexing(CWDM)systems.In this structure,the external rods of the ring resonator are arranged in a hexagon and the internal rods are removed in L arrangement for introducing defects.Scatter rods are used to prevent leakage.By using the L defected hexagonal resonator,a multi-channel CDF is designed,which exhibits multiple wavelengths of CWDM(1500 nm–1600 nm)region.In addition,the selection of rod size and the position of rods in the proposed multi-channel CDF are validated by varying the radius of coupling and scattering rods,as well as the position of resonators,respectively.By using plane wave expansion and opti Finite Difference Time Domain(FDTD)method,the electromagnetic wave propagation and the photonic band gap are obtained.

Selection of DBO measurement wavelength for bottom mark asymmetry based on FDTD method

A physical model for simulating overlay metrology employing diffraction based overlay(DBO)principles is built.It can help to optimize the metrology wavelength selection in DBO.Simulation result of DBO metrology with a model based on the finite-difference time-domain(FDTD)method is presented.A common case(bottom mark asymmetry)in which error signals are always induced in DBO measurement due to the process imperfection were discussed.The overlay sensitivity of the DBO measurement across the visible illumination spectrum has been performed and compared.After adjusting the model parameters compatible with the actual measurement conditions,the metrology wavelengths which provide the accuracy and robustness of DBO measurement can be optimized.

Temperature field simulation of polyolefin-absorber mixture by FDTD-FDM model during microwave heating

This work has performed a numerical simulation of the temperature field during microwave heating of polyolefin-absorber mixture by means of a combined electric and thermal model.A finite difference time domain was used to model the electric field distribution within the cavity,while the finite difference method was used to calculate the temperature field distribution in different reactors.This study has focused only on the process from room temperature to 500 K for reducing heating time and energy consumption.This temperature range is a process with high energy consumption,difficult to control and great influence on the follow-up reaction.Temperature dependence of dielectric properties and thermal properties of heated materials are fully considered and simulated through an iterative process.The simulation results show that input power,the size and location of the heated materials,the position of the waveguide,and the kinds of microwave absorbers are important factors affecting the heating process.As a result,the uniform temperature distribution(the temperature difference Tdb10 K)can be achieved by choosing the appropriate input power(500–2000 W),the appropriate proportion of microwave absorber(the volume ratio of Si C to HDPE is 30:70),and combining with the moving and rotating of the heated materials.The uniform temperature field obtained without mechanical stirring is very important for reducing energy consumption and subsequent reactions.

Main Factors for Affecting Photonic Bandgap of Photonic Crystals

The factors affecting one dimensional （1D） and two dimensional （2D） photonic crystals （PhCs） are systemically analyzed in this paper by numerical simulation. Transfer matrix method （TMM） is employed for 1D PCs, both finite difference time domain method （FDTD） and plane wave expansion method （PWE） are employed for 2D PCs. The result shows that the photonic bandgaps （PBG） are directly affected by crystal type, crystal lattice constant, modulation of refractive index and periodicity, and it is should be useful for design of different type photonic crystals with the required PBG and functional devices. Finally, as an example, a near-IR 1D PCs narrow filter was designed.

Modeling transmittance through submicron silver slit arrays

Mid-infrared transmittance of submicron silver slit arrays was numerically studied with the finite difference time domain method. The slit width varies from 50 nm to 300 nm and a square feature may attach at either or both slit sides. Although the side length of features is one or two orders of magnitude shorter than the wavelength, the attached nanoscale features can modify the transmittance significantly. The transmittance was also further investigated in detail by looking into the electromagnetic fields and Poynting vectors of selected slit geometries. The investigation results show that such change can be attributed to the cavity resonance effect inside the slit arrays. The work is of great importance to the wavelength-selective devices design in optical devices and thermal application fields.

OPTIMAL DESIGN OF SMART ANTENNA ARRAY

This letter investigates an efficient design procedure integrating the Genetic Algorithm (GA) with the Finite Difference Time Domain (FDTD) for the fast optimal design of Smart Antenna Arrays (SAA). The FDTD is used to analyze SAA with mutual coupling. Then,on the basis of the Maximal Signal to Noise Ratio (MSNR) criteria, the GA is applied to the optimization of weighting elements and structure of SAA. Finally, the effectiveness of the analysis is evaluated by experimental antenna arrays.