Investigation of three-pulse photon echo in thick crystal using finite-difference time-domain method

This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo＇s amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^（3＋）：YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp（α＇L）,which is previously employed to describe the relationship between echo＇s efficiency and thickness,should be modified as 1.3 · 0.09 exp（2.4 ·α＇L） with the propagation of echo considered.

Optical simulation of in-plane-switching blue phase liquid crystal display using the finite-difference time-domain method

The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell＇s equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.

An efficient locally one-dimensional finite-difference time-domain method based on the conformal scheme

An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.

A Revised Piecewise Linear Recursive Convolution FDTD Method for Magnetized Plasmas

The piecewise linear recursive convolution （PLRC） finite-different time-domain （FDTD） method improves accuracy over the original recursive convolution （RC） FDTD approach and current density convolution （JEC） but retains their advantages in speed and efficiency. This paper describes a revised piecewise linear recursive convolution PLRC-FDTD formulation for magnetized plasma which incorporates both anisotropy and frequency dispersion at the same time, enabling the transient analysis of magnetized plasma media. The technique is illustrated by numerical simulations of the reflection and transmission coefficients through a magnetized plasma layer. The results show that the revised PLRC-FDTD method has improved the accuracy over the original RC FDTD method and JEC FDTD method.

Theoretical Proof of Unconditional Stability of the 3-D ADI-FDTD Method

In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is introduced recently.In this paper,a theoretical proof of the stabilityof the three-dimensional(3-D)ADI-FDTD method is presented.It is shown that the 3-D ADI-FDTDmethod is unconditionally stable and free from the CFL condition restraint.

Analysis of Complex Electromagnetic Structures by Hybrid FDTD/WCIP Method

This paper proposes a hybrid full-wave analysis using Finite-Difference Time-Domain (FDTD) and Wave Concept Iterative Process (WCIP) methods, developed to analyze locally arbitrarily shaped microwave structures and Multilayer Planar structure. Using the equivalence principle, the original problem can be decomposed into two sub regions and solve each sub region separately. An interpolation scheme is proposed for communicating between the FDTD fields and WCIP wave, which will not require the effort of fitting the WCIP mesh to the FDTD cells in the interface region. This method is applied to calculate the scattering parameters of arbitrary (3-D) microwave structures. Applying FDTD to 3D discontinuity and WCIP to the remaining region preserves the advantages of both WCIP flexibility and FDTD efficiency. A comparison of the results with the FDTD staircasing data verifies the accuracy of the proposed method.

Unconditionally stable Crank-Nicolson algorithm with enhanced absorption for rotationally symmetric multi-scale problems in anisotropic magnetized plasma

Large calculation error can be formed by directly employing the conventional Yee’s grid to curve surfaces.In order to alleviate such condition,unconditionally stable CrankNicolson Douglas-Gunn(CNDG)algorithm with is proposed for rotationally symmetric multi-scale problems in anisotropic magnetized plasma.Within the CNDG algorithm,an alternative scheme for the simulation of anisotropic plasma is proposed in body-of-revolution domains.Convolutional perfectly matched layer(CPML)formulation is proposed to efficiently solve the open region problems.Numerical example is carried out for the illustration of effectiveness including the efficiency,resources,and absorption.Through the results,it can be concluded that the proposed scheme shows considerable performance during the simulation.

THE STUDY ON MEASURING TECHNIQUE OF PARTIAL DISCHARGE IN GAS INSULATED SWITCHGEAR USING ULTRA HIGH FREQUENCY METHOD WITH EXTERNAL SENSORS

Objective In order to find early latent faults and prevent catastrophic failures, diagnosis of insulation condition by measuring technique of partial discharge(PD) in gas insulated switchgear (GIS) is applied in this paper, which is one of the most basic ways for diagnosis of insulation condition. Methods Ultra high frequency(UHF) PD detection method by using internal sensors has been proved efficient, because it may avoid the disturbance of corona, but the sensor installation of this method will be limited by the structure and operation condition of GIS. There are some of electromagnetic (E-M) waves leak from the place of insulation spacer, therefore, the external sensors UHF measuring PD technique is applied, which isn't limited by the operation condition of GIS. Results This paper analyzes propagated electromagnetic (E-M) waves of partial discharge pulse excited by using the finite-difference time-domain (FDTD) method. The signal collected at the outer point is more complex than that of the inner point, and the signals' amplitude of outer is about half of the inner, because it propagates through spacer and insulation slot. Set up UHF PD measuring system. The typical PD in 252kV GIS bus bar was measured using PD detection UHF technique with external sensors. Finally, compare the results of UHF measuring technique using external sensors with the results of FDTD method simulation and the traditional IEC60270 method detection. Conclusion The results of experiment shows that the UHF technique can realize the diagnosis of insulation condition, the results of FDTD method simulation and the result UHF method detection can demonstrate each other, which gives references to further researches and application for UHF PD measuring technique.

FDTD Simulation on Terahertz Waves Propagation Through a Dusty Plasma

The frequency dependent permittivity for dusty plasmas is provided by introducing the charging response factor and charge relaxation rate of airborne particles. The field equations that describe the characteristics of Terahertz（THz） waves propagation in a dusty plasma sheath are derived and discretized on the basis of the auxiliary differential equation（ADE） in the finite difference time domain（FDTD） method. Compared with numerical solutions in reference, the accuracy for the ADE FDTD method is validated. The reflection property of the metal Aluminum interlayer of the sheath at THz frequencies is discussed. The effects of the thickness, effective collision frequency, airborne particle density, and charge relaxation rate of airborne particles on the electromagnetic properties of Terahertz waves through a dusty plasma slab are investigated.Finally, some potential applications for Terahertz waves in information and communication are analyzed.

SO-FDTD Analysis of the Plasma Reflectance of Epstein Distribution

The analysis of electromagnetic propagation in a dispersive medium is complicated in the time-domain because its dielectric constant is frequency-dependent. In this paper, the dielectric constant of the dispersive medium is written as a rational polynomial function, and the relationship between D and E is derived in the time-domain. It is referred to as the shift operator finite-different time-domain （SO-FDTD） method. Compared to an analytical solution and a piecewise linear current density recursive convolution （PLJERC） solution, the high accuracy and efl%iency of this method is verified by calculating the reflectance of the electromagnetic wave through a cold plasma slab. As the electron density in plasma is distributed as the Epstein formula, the effect of distribution grads and electron collision frequency on the reflectance is calculated by using the SO-FDTD method. The result shows that the increase in the distribution grads coefficient affects the reflectance sharply. When it comes to a smaller distribution grads coelBcient, the increase of the collision frequency showed a significant effect on the reflectance, but on the contrary, there is actually less and less effect till it disappears.

FDTD Simulations on Plasmonic Properties of End-to-End and Side-by-Side Assembled Au Nanorods

The surface plasmon resonance( SPR) of gold( Au) nanorod can be tuned in a large visible-near infrared( Vis-NIR) region by changing the aspect ratio of nanorod. Compared with the SPR of isolated Au nanorod,assembly of Au nanorods exhibits strong coupling effect in the nanogap and wealthy changes in the optical spectra. The SPR coupling effects and localized electronic fields for end-to-end( E-E) and side-by-side( S-S) assembled Au nanorods dimers are studied through finite-difference time-domain( FDTD) simulation. With decreasing the gap spacing,the longitudinal SPR( SPRL) red-shifts for the E-E orients Au nanorods dimer and blue-shifts for S-S orientes Au nanorods dimer. The transverse SPR( SPRT) has slight red-shifting for S-S assembly and no shifting for the E-E assembly. Moreover,a new coupling SPR appears for the E-E assembly in a long wavelength in the NIR region,blue-shifting and enhancing with decreasing the gap spacing. Based on the spring oscillator model and the polarization of nanoparticles with incident electric field,the SPR shifting and the appearance of new coupling SPR of assembled Au nanorods are proposed.

Optimization of band-notched UWB antenna using micro-genetic algorithm combined with FDTD

The micro-genetic algorithm (MGA) optimization combined with the finite-difference time-domain (FDTD) method is applied to design a band-notched ultra wide-band (UWB) antenna. A U-type slot on a stepped U-type UWB monopole is used to obtain the band-notched characteristic for 5 GHz wireless local area network (WLAN) band. The measured results show that voltage standing wave ration (VSWR) less than 2 covers 3.1-10.6 GHz operating band and VSWR more than 2 is within 5.150-5.825 GHz notched one with the highest value of 5.6. Agreement among the calculated, HFSS simulated and measured results validates the effiectiveness of this MGA-FDTD method, which is efficient for UWB antennas design.

Characteristic-Based Time Domain Method for Cylindrically Conformal Microstrip Patch Antennas

The characteristics of a cylindrical conformal microstrip patch antenna are analyzed by using the characteristic-based time domain （CBTD） method. A governing equation in the cylindrical coordinate system is formulated directly to facilitate the analysis of cylindrically conformal microstrip patch antennas. The algorithm has second-order accuracy both in time and space domain and has the potential to eliminate the spurious wave reflection from the numerical boundaries of the computational domain, Numerical results demonstrate the important merits and accuracy of the proposed technique in computational electromagnetics,

A Study of Properties of the Photonic Band Gap of Unmagnetized Plasma Photonic Crystal

In this study, the propagation of electromagnetic waves in one-dimensional plasma photonic crystals （PPCs）, namely, superlattice structures consisting alternately of a homogeneous unmagnetized plasma and dielectric material, is simulated numerically using the finite-difference time-domain （FDTD） algorithm. A perfectly matched layer （PML） absorbing technique is used in this simulation. The reflection and transmission coefficients of electromagnetic （EM） waves through PPCs are calculated. The characteristics of the photonic band gap （PBG） are discussed in terms of plasma density, dielectric constant ratios, number of periods, and introduced layer defect. These may provide some useful information for designing plasma photonic crystal devices.

Propagation characteristics of oblique incidence terahertz wave through non-uniform plasma

The propagation characteristics of oblique incidence terahertz(THz) waves through non-uniform plasma are investigated by the shift-operator finite-difference time-domain(SO-FDTD) method combined with the phase matching condition.The electron density distribution of the non-uniform plasma is assumed to be in a Gaussian profile. Validation of the present method is performed by comparing the results with those obtained by an analytical method for a homogeneous plasma slab.Then the effects of parameters of THz wave and plasma layer on the propagation properties are analyzed. It is found that the transmission coefficients greatly depend on the incident angle as well as on the thickness of the plasma, while the polarization of the incident wave has little influence on the propagation process in the range of frequency considered in this paper. The results confirm that the THz wave can pass through the plasma sheath effectively under certain conditions,which makes it a potential candidate to overcome the ionization blackout problem.

A subwavelength metal-grating assisted sensor of Kretschmann style for investigating the sample with high refractive index

In this paper, a subwavelength metal-grating assisted sensor of Kretschmann style that is capable of detecting the sample with a refractive index higher than that of the substrate is proposed. The sensor configuration is similar to the traditional Kretschmann structure, but the metal film is pattered into a grating. As a TM-polarized laser beam impinges from the substrate, a resonant dip point in reflectance curve is produced at a certain incident angle. Our studies indicate that the sensing sensitivity and resolution are affected by the grating＇s gap and period, and after these parameters havelaeen optimized, a sensing sensitivity of 51.484°/RIU is obtained with a slightly changing resolution.

Exponential time differencing based efficient SC-PML for RCS simulation

To efficiently simulate and calculate the radar cross section(RCS) related electromagnetic problems by employing the finite-difference time-domain(FDTD) algorithm, an efficient stretched coordinate perfectly matched layer(ESC-PML) based upon the exponential time differencing(ETD) method is proposed.The proposed implementation can not only reduce the number of auxiliary variables in the SC-PML regions but also maintain the ability of the original SC-PML in terms of the absorbing performance. Compared with the other existed algorithms, the ETDFDTD method shows the least memory consumption resulting in the computational efficiency. The effectiveness and efficiency of the proposed ESC-PML scheme is verified through the RCS relevant problems including the perfect E conductor(PEC) sphere model and the patch antenna model. The results indicate that the proposed scheme has the advantages of the ETD-FDTD method and ESC-PML scheme in terms of high computational efficiency and considerable computational accuracy.

Rear-surface light intensification caused by a Hertzian-conical crack in 355-nm silica optics

Theoretical studies show that a Hertzian-conical crack can be considered to be composed of double cone faces for simplicity. In the present study, the three-dimensional finite-difference time-domain method is employed to quantify the electric-field distribution within the subsurface in the presence of such a defect under normal incidence irradiation. Both impurities （inside the crack） and the chemical etching have been investigated. The results show that the maximum electric field amplitude │E│max is 9.57374 V/m when the relative dielectric constant of transparent impurity equals 8.5. And the near-field modulation will be improved if the crack is filled with the remainder polishing powders or water vapor/drops. Meanwhile, the laser-induced initial damage moves to the glass-air surface. In the etched section, the magnitude of intensification is strongly dependent on the inclination angle θ. There will be a highest modulation when θ is around π/6, and the maximum value of IEImax is 18.57314 V/m. When θ ranges from π/θ to π/4, the light intensity enhancement factor can easily be larger than 100, and the modulation follows a decreasing trend. On the other hand, the modulation curves become smooth when θ〉π/4 or θ 〈 π/8.

A new parallel meshing technique integrated into the conformal FDTD method for solving complex electromagnetic problems

A new efficient parallel finite-difference time-domain(FDTD) meshing algorithm, based on the ray tracing technique, is proposed in this paper. This algorithm can be applied to construct various FDTD meshes, such as regular and conformal ones. The Microsoft F# language is used for the algorithm coding, where all variables are unchangeable with its parallelization advantage being fully exploited. An improved conformal FDTD algorithm, also integrated with an improved surface current algorithm, is presented to simulate some complex 3D models, such as a sphere ball made of eight different materials, a tank, a J-10 aircraft, and an aircraft carrier with 20 aircrafts. Both efficiency and capability of the developed parallel FDTD algorithm are validated. The algorithm is applied to characterize the induced surface current distribution on an aircraft or a warship.

Approximate CN scheme and its open region problems for metamaterial rotational symmetric simulation

In order to simulate metamaterial rotational symmetric open region problems,unconditionally stable perfectly match layer(PML)implementation is proposed in the body of revolution(BOR)finite-difference time-domain(FDTD)lattice.More precisely,the proposed algorithm is implemented by the Crank-Nicolson(CN)Douglas-Gunn(DG)procedure for BOR metamaterial simulation.The constitutive relationship of metamaterial can be expressed by the Drude model and calculated by the piecewise linear recursive convolution(PLRC)approach.The effectiveness including absorption,efficiency,and accuracy is demonstrated through the numerical example.It can be concluded that the proposed implementation is to take the advantages of the CNDG-PML procedure,PLRC approach,and BORFDTD algorithm in terms of considerable accuracy,enhanced absorption and remarkable efficiency.Meanwhile,it can be demonstrated that the proposed scheme can maintain its unconditional stability when the time step exceeds the CourantFriedrichs-Levy(CFL)condition.