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.

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.

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.

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）.

A novel method for calculating broadband electrical performance of high-speed aircraft radome under thermo-mechanical-electrical coupling

The electrical performance of radomes on high-speed aircraft can be influenced by the thermal and mechanical loads produced during high-speed flight,which can affect the detection dis-tance and accuracy of the guidance system.This paper presents a new method that uses the Finite Difference Time Domain(FDTD)method to calculate the electrical performance of radomes under Thermo-Mechanical-Electrical(TME)coupling.This method can accurately characterize the effects of material dielectric temperature drift and structural deformation on the electrical performance of the radome under flight conditions,enabling high-precision full-wave calculations of the broadband electrical performance of the radome.The method initiates by utilizing a Finite Element Grid Model(FE-GM)of the radome to sequentially acquire the radome's response temperature field and structural deformation field through thermal and mechanical simulations.Subsequently,spatial mapping techniques are developed to accurately incorporate the dielectric temperature drift and structural deformation of the radome into its Yee grid Electromagnetic(EM)simulation model.A verification case was designed to test the proposed method,and the results confirmed its high computational accuracy.Additionally,the effectiveness and necessity of the method were further demonstrated by analyzing the electrical performance of a fused silica ceramic radome used on a high-speed aircraft.

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.

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.

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.

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 infrared radiative properties of nanoscale metallic complex slit arrays

The radiative properties(absorptance, reflectance, and transmittance) of deep slits with five nanoscale slit profile variations at the transverse magnetic wave incidence were numerically investigated by employing the finite difference time domain method. For slits with attached features, their radiative properties can be much different due to the modified cavity geometry and dangled structures, even at wavelengths between 3 and 15 μm. The shifts of cavity resonance excitation result in higher transmittance through narrower slits at specific wavelengths and resonance modes are confirmed with the electromagnetic fields. Opposite roles possibly played by features in increasing or decreasing absorptance are determined by the feature position and demonstrated by Poynting vectors. Correlations among all properties of a representative slit array and the slit density are also comprehensively studied. When multiple slit types coexist in an array(complex slits), a wide-band transmittance or absorptance enhancement is feasible by merging spectral peaks contributed from each type of slits distinctively. Discrepancy among infrared properties of four selected slit combinations is explained while effects of slit density are also discussed.

Modeling bidirectional reflection distribution function of microscale random rough surfaces

The radiative properties of three different materials surfaces with one-dimensional microscale random roughness were obtained with the finite difference time domain method(FDTD) and near-to-far-field transformation.The surface height conforms to the Gaussian probability density function distribution.Various computational modeling issues that affect the accuracy of the predicted properties were discussed.The results show that,for perfect electric conductor(PEC) surfaces,as the surface roughness increases,the magnitude of the spike reduces and eventually the spike disappears,and also as the ratio of root mean square roughness to the surface correlation distance increases,the retroreflection becomes evident.The predicted values of FDTD solutions are in good agreement with the ray tracing and integral equation solutions.The overall trend of bidirectional reflection distribution function(BRDF) of PEC surfaces and silicon surfaces is the same,but the silicon's is much less than the former's.The BRDF difference from two polarization modes for the gold surfaces is little for smaller wavelength,but it is much larger for the longer wavelength and the FDTD simulation results agree well with the measured data.In terms of PEC surfaces,as the incident angle increases,the reflectivity becomes more specular.

Optimizing bowtie structure parameters for specific incident light

We investigate optical properties of a bowtie-shaped aperture using the finite difference time domain method to optimize its geometric parameters for specific incident lights. The influence of the parameters on local field enhancement and resonant wavelength in the visible frequency range is numerically analysed. It is found that the major resonance of the spectrum is exponentially depended on the bowtie angle but independent of the whole aperture size. The simulation also demonstrates that increasing the aperture size raises the local field intensity on the exit plane due to an enlarged interaction area between the light and the metal medium. And the near-field spot size is closely related to the gap. Based on these results, the design rules of the bowtie structure can be optimized for specific wavelengths excited.

Optimization of Short-Pulse GPR Transmit Antenna

By analyzing the current distribution of Bow-Tie antenna used in short-pulse ground penetrating radar, the methods of antenna load and driving are presented in this paper to reduce strength of reflective wave both at antenna end and excitation point. The numerical simulation results show the strength of reflective wave is smaller than ?55 dB comparing with the driving wave when the methods are adopted. Key words short-pulse GPR - resistive loading - finite difference time domain method CLC number TN 820 Foundation item: Supported by the National Natural Science Foundation of China (49984001)Biography: LI Tai-quan (1961-), male, Ph. D candidate, research direction: antenna model and design

Numerical investigation of the enhanced unidirectional surface plasmon polaritons generator

A unidirectional surface plasmon polaritons(SPPs) generator with greatly enhanced generation efficiency is proposed. The SPPs generator consists of an asymmetric single nanoslit coated with a polyviny alcohol(PVA) film and a silver rectangle block. The generation efficiency of this SPPs generator is investigated using the finite difference time domain method. Due to the presence of the silver rectangle block, the SPPs generation efficiency of the asymmetric single nanoslit with PVA film can be greatly enhanced and the corresponding wavelength with the maximum enhancement factor can be tuned flexibly. The influence of the structural parameters on the generation efficiency is also investigated for the enhanced unidirectional SPPs generator.

Polarization insensitivity electromagnetically induced reflection in graphene metasurface

The electromagnetically induced reflection(EIR)effect of graphene metamaterials has been investigated by finite difference time domain(FDTD)method.In this study,a metamaterial sandwich structure composed of silica(SiO2),gold and graphene on terahertz band is designed.By changing the width of the two ribbons of graphene length and the incident angle of electromagnetic wave,the EIR effect of the structure is discussed,and it can be found that SiO2 is a kind of excellent dielectric material.The simulation results show that graphene metamaterial is not sensitive to polarized incident electromagnetic wave.Therefore,such EIR phenomena as insensitive polarization and large incident angle can be applied to optical communication filters and terahertz devices.

A functional probe with bowtie aperture and bull's eye structure for nanolithograph

The bowtie aperture surrounded by concentric gratings(the bull’s eye structure) integrated on the near-field scanning optical microscopy(NSOM) probe(aluminum coated fiber tip) for nanolithography has been investigated using the finite-difference time domain(FDTD) method.By modifying the parameters of the bowtie aperture and the concentric gratings,a maximal field enhancement factor of 391.69 has been achieved,which is 18 times larger than that obtained from the single bowtie aperture.Additionally,the light spot depends on the gap size of the bowtie aperture and can be confined to sub-wavelength.The superiority of the combination of the bowtie aperture and the bull’s eye structure is confirmed,and the mechanism for the electric field enhancement in this derived structure is analyzed.

Solution of the time-dependent Schrdinger equation with absorbing boundary conditions

The performances of absorbing boundary conditions （ABCs） in four widely used finite difference time domain （FDTD） methods, i.e. explicit, implicit, explicit staggered-time, and Chebyshev methods, for solving the time-dependent Schrodinger equation are assessed and compared. The computation efficiency for each approach is also evaluated. A typical evolution problem of a single Gaussian wave packet is chosen to demonstrate the performances of the four methods combined with ABCs. It is found that ABCs perfectly eliminate reflection in implicit and explicit staggered-time methods. However, small reflection still exists in explicit and Chebyshev methods even though ABCs are applied.

Analysis of directional emission in square resonator lasers with an output waveguide

Square microcavity laser with an output waveguide is proposed and analyzed by the finite-difference time-domain （FDTD） technique. For a square resonator with refractive index of 3.2, side length of 4μm, and output waveguide of 0.4μm width, we have got the quality factors （Q factors） of 6.7 ×10^2 and 7.3 × 10^3 for the fundamental and first-order transverse magnetic （TM） mode near the wavelength of 1.5μm, respectively. The simulated intensity distribution for the first-order TM mode shows that the coupling efficiency in the waveguide reaches 53%. The numerical simulation shows that the first-order transverse modes have fairly high Q factor and high coupling efficiency to the output waveguide. Therefore the square resonator with an output waveguide is a promising candidate to realize single-mode directional emission microcavity lasers.

Investigation of mode radiation loss for microdisk resonators with pedestals by FDTD technique

Mode radiation loss for microdisk resonators with pedestals is investigated by three-dimensional（3D） finite-difference time-domain（FDTD） technique.For the microdisk with a radius of 1μm,a thickness of 0.2μm,and a refractive index of 3.4,on a pedestal with a refractive index of 3.17,the mode quality（Q） factor of the whispering-gallery mode（WGM） quasi-TE_（7.1） first increases with the increase of the radius of the pedestal,and then quickly decreases as the radius is larger than 0.75μm.The mode radiation loss is mainly the vertical radiation loss induced by the mode coupling between the WGM and vertical radiation mode in the pedestal,instead of the scattering loss around the perimeter of the round pedestal.The WG M can keep the high Q factor when the mode coupling is forbidden.