Parallel computation of unified finite-difference time-domain for underwater sound scattering

In this work, we treat scattering objects, water, surface and bottom in a truly unified manner in a parallel finitedifference time-domain （FDTD） scheme, which is suitable for distributed parallel computing in a message passing interface （MPI） programming environment. The algorithm is implemented on a cluster-based high performance computer system. Parallel computation is performed with different division methods in 2D and 3D situations. Based on analysis of main factors affecting the speedup rate and parallel efficiency, data communication is reduced by selecting a suitable scheme of task division. A desirable scheme is recommended, giving a higher speedup rate and better efficiency. The results indicate that the unified parallel FDTD algorithm provides a solution to the numerical computation of acoustic scattering.

Uniform stable conformal convolutional perfectly matched layer for enlarged cell technique conformal finite-difference time-domain method

Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer（CPML） is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain（ECT-CFDTD） method is used to simulate the wave propagation inside a perfect electric conductor（PEC） waveguide.The algorithm has the same numerical stability as the ECT-CFDTD method.For the long-time propagation problems of an evanescent wave in a waveguide,several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML.Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide.

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.

A spherical higher-order finite-difference time-domain algorithm with perfectly matched layer

A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spherical resonator is modeled in order to demonstrate the advantage of this scheme over the finite-difference time-domain（FDTD） and the multiresolution time-domain（MRTD） schemes with respect to memory requirements and CPU time. Moreover, the Berenger＇s perfectly matched layer（PML） is derived for the spherical HO-FDTD grids, and the numerical results validate the efficiency of the PML.

Adjustable transmission properties through ring-shaped nanotube arrays using finite-difference time-domain method

Metallic ring-shaped nanotube arrays are proposed and its optical transmission properties are studied by using finite-difference time-domain （FDTD） method. Compared with the transmission spectra of conventional circular nanotube arrays, two photonic band gaps are emerged in the transmission spectra offing-shaped nanotube arrays, the two band gaps and transmission spectra are adjusted by the length, inner radius, intertube spacing and the dielectric constants of the core and embedding medium, and magnitude modification, redshift and blueshift of the resonance modes are observed. A metallic ring-shaped nanotube arrays for subwavelength band-stop filter in the range of visible light can be achieved. To understand its physical origin, field-interference mechanism was suggested by the field distributions. The proposed nanostructures and results may have great potential applications in subwavelength near-field optics.

Perfect plane-wave source for a high-order symplectic finite-difference time-domain scheme

The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order symplectic finite- difference time-domain （SFDTD） scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain （FDTD）, the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.

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.

IMPROVED LOCALLY CONFORMAL FINITE-DIFFERENCE TIME-DOMAIN METHOD FOR EDGE INCLINED SLOTS IN A FINITE WALL THICKNESS WAVEGUIDE

An Improved Locally Conformal Finite-Difference Time-Domain (ILC-FDTD) method is presented in this paper, which is used to analyze the edge inclined slots penetrating adjacent broadwalls of a finite wall thickness waveguide. ILC-FDTD not only removes tile instability of the original locally conformal FDTD algorithm, but also improves the computational accuracy by locally modifying magnetic field update equations and the virtual iterative electric fields accordlng to the complexity of tile slot fringe fields. The mutual coupling between two edge inclined slots can also be analyzed by ILC-FDTD effectively.

Finite-difference time-domain modeling of curved material interfaces by using boundary condition equations method

To deal with the staircase approximation problem in the standard finite-difference time-domain（FDTD） simulation,the two-dimensional boundary condition equations（BCE） method is proposed in this paper.In the BCE method,the standard FDTD algorithm can be used as usual,and the curved surface is treated by adding the boundary condition equations.Thus,while maintaining the simplicity and computational efficiency of the standard FDTD algorithm,the BCE method can solve the staircase approximation problem.The BCE method is validated by analyzing near field and far field scattering properties of the PEC and dielectric cylinders.The results show that the BCE method can maintain a second-order accuracy by eliminating the staircase approximation errors.Moreover,the results of the BCE method show good accuracy for cylinder scattering cases with different permittivities.

Numerical modeling of the 2D time-domain transient electromagnetic secondary field of the line source of the current excitation

To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional （2D） time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force （Emf）, and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.

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.

An efficient absorbing boundary for finite-difference time-domain field modelling in acoustics

A highly efficient absorbing boundary condition suitable for use in the finitedtherence time-domain (FDTD) modelling of acoustic fields is presented in this paper. The new method seeks a least square esthoate of a transfer matrix for field components near truncating boundaries by matrir pseud-inversion. The proposed absorbing boundary is considerably more effective than most ekisting ones. The method is also computationally econondcal and robust.The performance of the new method is shown by numerical experiments on a point-source radiation problem, a wedge dimaction problem, and a scattering problem in which a plane wave is scattered by a circular cylinder.

Application of Padé Approximation in Simulating Photonic Crystals

To save finite-difference time-domain（FDTD） computing time, several methods are proposed to convert the time domain FDTD output into frequency domain. The Pad6 approximation with Baker＇s algorithm and the program are introduced to simulate photonic crystal structures. For a simple pole system with frequency 160THz and quality factor of 5000, the intensity spectrum obtained by the Padé approximation from a 2^8-item sequence output is more exact than that obtained by fast Fourier transformation from a 2^20-item sequence output. The mode frequencies and quality factors are calculated at different wave vectors for the photonic crystal slab from a much shorter FDTD output than that required by the FFT method, and then the band diagrams are obatined. In addition, mode frequencies and Q-factors are calculated for photonic crystal microcavity.

Numerical analysis of multicomponent responses of surface-hole transient electromagnetic method

We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets.

Three-dimensional numerical modeling of fullspace transient electromagnetic responses of water in goaf

The full-space transient electromagnetic response of water-filled goaves in coal mines were numerically modeled. Traditional numerical modeling methods cannot be used to simulate the underground full-space transient electromagnetic field. We used multiple transmitting loops instead of the traditional single transmitting loop to load the transmitting loop into Cartesian grids. We improved the method for calculating the z-component of the magnetic field based on the characteristics of full space. Then, we established the full- space 3D geoelectrical model using geological data for coalmines. In addition, the transient electromagnetic responses of water-filled goaves of variable shape at different locations were simulated by using the finite-difference time-domain （FDTD） method. Moreover, we evaluated the apparent resistivity results. The numerical modeling results suggested that the resistivity differences between the coal seam and its roof and floor greatly affect the distribution of apparent resistivity, resulting in nearly circular contours with the roadway head at the center. The actual distribution of apparent resistivity for different geoelectrical models of water in goaves was consistent with the models. However, when the goal water was located in one side, a false low-resistivity anomaly would appear on the other side owing to the full-space effect but the response was much weaker. Finally, the modeling results were subsequently confirmed by drilling, suggesting that the proposed method was effective.

Borehole-GPR numerical simulation of full wave field based on convolutional perfect matched layer boundary

The absorbing boundary is the key in numerical simulation of borehole radar.Perfect match layer(PML) was chosen as the absorbing boundary in numerical simulation of GPR.But CPML(convolutional perfect match layer) approach that we have chosen has the advantage of being media independent.Beginning with the Maxwell equations in a two-dimensional structure,numerical formulas of finite-difference time-domain(FDTD) method with CPML boundary condition for transverse electric(TE) or transverse magnetic(TM) wave are presented in details.Also,there are three models for borehole-GPR simulation.By analyzing the simulation results,the features of targets in GPR are obtained,which can provide a better interpretation of real radar data.The results show that CPML is well suited for the simulation of borehole-GPR.

MWI in Cylindrical Coordinates and Its Application

Based on FDTD difference expressions and eigenfunctions of Maxwell functions in cylindrical coordinates, mesh wave impedances (MWIs) in 2D and 3D cylindrical coordinates were introduced. Combined with the concept of perfectly matched layer (PML), MWI PML absorbing boundary condition (ABC) algorithm was deduced in 2D cylindrical coordinates. Numerical experiments were done to investigate the validity of MWI and its application in cylindrical coordinates FDTD algorithm. The results showed that MWI in cylindrical coordinates can be used to accurately calculate the numerical reflection error caused by different mesh increments in non uniform FDTD. MWI can also provide theoretical criterion to define the permitted variable range of mesh dimension. MWI PML ABC is easy to be applied and reduces low numerical reflection, which only causes a little higher reflection error compared with Teixeira's PML.

Investigation of composite electromagnetic scattering from ship-like target on the randomly rough sea surface using FDTD method

Composite electromagnetic scattering from a two-dimensional （2D） ship-like target on a one-dimensional sea surface is investigated by using the finite-difference time-domain （FDTD） method. A uniaxial perfectly matched layer is adopted for truncation of FDTD lattices.The FDTD updated equations can be used for the total computation domain by choosing the uniaxial parameters properly. To validate the proposed numerical technique,a 2D infinitely long cylinder over the sea surface is taken into account first.The variation of angular distribution of the scattering changing with incident angle is calculated. The results show good agreement with the conventional moment method. Finally,the influence of the incident angle,the polarization,and the size of the ship-like target on the composite scattering coefficient is discussed in detail.

NUMERICAL STUDIES ON RESONANT AND ENHANCEMENT EFFECTS FOR COUPLING OF MICROWAVE PULSES INTO NARROW SLOTS

In this paper, modifications to the finite-difference time-domain(FD-TD) method for modeling microwave pulse coupling into a slot, which is much narrower than one conventional FD-TD cell, are discussed. The coupling process of microwave pulse into a slot is studied by using the modified FD-TD method, and the dependence of microwave coupling on slot sizes, the carrier frequencies and the polarization directions of the incident waves is analysed. Resonant and enhancement effects which occur in this process are observed. The condition at which the resonant effect takes place is also presented.