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.

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.

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.

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.

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.

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,

Benchmarking of two three-dimensional numerical models in time/space domain to predict railway-induced ground vibrations

In the last 30 years,the scientific community has developed and proposed different models and numerical approaches for the study of vibrations induced by railway traffic.Most of them are formulated in the frequency/wave number domain and with a 2.5D approach.Three-dimensional numerical models formulated in the time/space domain are less frequently used,mainly due to their high computational cost.Notwithstanding,these models present very attractive characteristics,such as the possibility of considering nonlinear behaviors or the modelling of excess pore pressure and non-homogeneous and non-periodic geometries in the longitudinal direction of the track.In this study,two 3D numerical approaches formulated in the time/space domain are compared and experimentally validated.The first one consists of a finite element approach and the second one of a finite difference approach.The experimental validation in an actual case situated in Carregado(Portugal)shows an acceptable fitting between the numerical results and the actual measurements for both models.However,there are some differences among them.This study therefore includes some recommendations for their use in practical soil dynamics and geotechnical engineering.

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

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.

A wavelet multiscale method for inversion of Maxwell equations

This paper is concerned with estimation of electrical conductivity in Maxwell equations. The primary difficulty lies in the presence of numerous local minima in the objective functional. A wavelet multiscale method is introduced and applied to the inversion of Maxwell equations. The inverse problem is decomposed into multiple scales with wavelet transform, and hence the original problem is reformulated to a set of sub-inverse problems corresponding to different scales, which can be solved successively according to the size of scale from the shortest to the longest. The stable and fast regularized Gauss-Newton method is applied to each scale. Numerical results show that the proposed method is effective, especially in terms of wide convergence, computational efficiency and precision.

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.

Spatially dispersive model of wire medium for scattering

The electromagnetic properties of wire medium are studied by analyzing its transmission coefficient, and a homogeneous uniaxial effective permittivity tensor with one spatially dispersive component is obtained. In order to evaluate the validity of the equivalent model for scattering problems, the bistatic and monostatic radar cross sections （RCS） of a block of wire medium are investigated with the finite-difference time-domain （FDTD） method. The difference of RCS from the physical wire medium and the corresponding equivalent model has been compared and analyzed under va- rious parameters of wire medium, which clearly demonstrates the accuracy of the equivalent model of wire medium.

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.

高聚物注浆修复的面板堤坝探地雷达波场特征分析

针对高聚物注浆修复堤坝面板脱空病害的修复效果问题,依据实际堤坝面板层间结构和材料介电参数,基于时域有限差分法和完全匹配层边界条件建立了高聚物注浆修复面板脱空堤坝的探地雷达电磁波计算模型,分析了雷达中心频率、面板脱空修复程度及脱空区域大小、面板厚度及配筋等因素对高聚物注浆修复面板脱空堤坝探地雷达波场特征的影响。结果表明:随着激励源中心频率增大,探地雷达剖面图分辨率逐渐增加;随着脱空区域长度增加,探地雷达剖面图中产生的水平状界面反射波增长;随着脱空区域深度增加,探地雷达剖面图水平状反射线之间时间间隔逐渐增大;随着堤坝面板厚度增加,在面板脱空修复区域绕射波和衍射波幅值减小;探地雷达发射的电磁波遇到钢筋时将产生能量较强的绕射波,强绕射波将脱空区域上下分界面的反射波分割开,难以判断反射波的水平长度。根据高聚物注浆修复面板脱空堤坝的探地雷达波场特征,能有效判断堤坝面板脱空区域的高聚物注浆修复效果,有助于解释实际堤坝面板高聚物注浆修复效果的探地雷达图像。

基于探地雷达的杂散电流环境下混凝土损伤检测研究

因走行轨不能对地完全绝缘,地铁运行时会产生杂散电流,导致混凝土中的钢筋发生电化学腐蚀。由于电化学腐蚀的持续发生,腐蚀产物累积在钢筋表面周围,从而对混凝土产生膨胀应力,导致混凝土开裂脱落,因此有必要对杂散电流影响严重区域的钢筋混凝土结构进行检测。该文首先基于多物理场法对钢筋混凝土的电化学腐蚀进行二维瞬态仿真研究,得到与实际混凝土结构相近似的腐蚀损伤图;进而对得到的损伤图进行探地雷达数值正演模拟检测。研究结果表明,在钢筋混凝土发生电化学腐蚀的过程中,钢筋直径的减小和腐蚀产物厚度的增加对探地雷达检测信号影响较小,因膨胀应力而产生的混凝土裂缝对探地雷达检测信号的影响起主要作用。

基于二维光子晶体全光逻辑门的设计

在二维光子晶体中嵌入了线缺陷,利用线性干涉效应和波导耦合,设计了一种基于二维光子晶体的同或门和与非门结构。主要采用平面波展开法对该二维光子晶体的能带结构进行分析,采用时域有限差分法,结合线性干涉效应,在Rsoft平台对所设计的同或门和与非门进行稳定电场图和归一化功率仿真。仿真结果标明:设计的同或门对比度高达29.5 dB,响应时间为0.073 ps,数据传输速率为13.7 Tbit/s;设计的与非门对比度高达24.15 dB,响应时间为0.08 ps,数据传输速率为12.5 Tbit/s。这些结果表明所设计的结构对比度高、响应时间短和数据传输速率快。

端接复杂电路传输线网络的电磁耦合时域并行计算方法

针对端接复杂电路传输线(TL)网络的电磁耦合问题,仍缺乏高效的场路协同仿真技术。该文将传输线方程与时域有限差分(FDTD)方法、诺顿定理和置换定理以及NGSPICE软件相结合,并引入消息传递接口(MPI)并行技术,提出一种高效的时域混合并行算法(FDTDTL-NGSPICE)。首先,根据诺顿定理和置换定理,将传输线网络分解为传输线子系统和复杂电路子系统,并构建对应的等效电路模型。然后,使用FDTDTL并行算法计算传输线子系统沿线各点的电压和电流,并获取对应诺顿等效电路的电流源和等效导纳大小。最后,使用NGSPICE对复杂电路子系统进行传导干扰分析,获得复杂电路各元件上的瞬态响应,并将端口电压反馈给传输线子系统作为边界,实现传输线网络电磁耦合的场线路联合协同仿真。通过对3类典型场景的计算实例,分别使用时域混合并行算法和电磁仿真软件CST电缆工作室(CS)进行数值模拟并对比,验证所提算法的置信度。