Viscoacoustic prestack reverse time migration based onthe optimal time-space domain high-order finite-difference method

Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.

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.

Efficient Finite Difference/Spectral Method for the Time Fractional Ito Equation Using Fast Fourier Transform Technic

A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.

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.

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.

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

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.

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.

Modeling of Surface Waves in a Fluid Saturated Poro-Elastic Medium under Initial Stress Using Time-Space Domain Higher Order Finite Difference Method

In this present context, mathematical modeling of the propagation of surface waves in a fluid saturated poro-elastic medium under the influence of initial stress has been considered using time dependent higher order finite difference method (FDM). We have proved that the accuracy of this finite-difference scheme is 2M when we use 2nd order time domain finite-difference and 2M-th order space domain finite-difference. It also has been shown that the dispersion curves of Love waves are less dispersed for higher order FDM than of lower order FDM. The effect of initial stress, porosity and anisotropy of the layer in the propagation of Love waves has been studied here. The numerical results have been shown graphically. As a particular case, the phase velocity in a non porous elastic solid layer derived in this paper is in perfect agreement with that of Liu et al. (2009).

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.

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 Methods for the Time Fractional Advection-diffusion Equation

In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald-Letnikov formula of order α ∈(0, 1). We investigate the stability analysis by using von Neumann method with mathematical induction and prove that these three proposed methods are unconditionally stable. Numerical results are presented to demonstrate the effectiveness of the schemes mentioned in this paper.

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.

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.

高阶SF-SFDTD方法在含时薛定谔方程求解中的应用研究

时域有限差分方法(finite-difference time-domain,FDTD(2,2))被广泛用于量子力学中薛定谔方程的求解,然而受Courant-Friedrichs-Lewy(CFL)条件的影响,计算空间中的网格尺寸会限制时间步长的取值范围,极大降低了FDTD(2,2)方法的数值计算效率.另外,FDTD(2,2)方法在时间域和空间域只具有二阶数值精度,在计算中往往会导致较大的误差累计,影响仿真结果的正确性.为了克服这些问题,结合空间滤波方法(spatial filtering,SF)和高阶辛时域有限差分(symplectic finite-difference time-domain,SFDTD(3,4))方法(3和4分别表示时间和空间数值精度),提出了一种时间稳定性条件可扩展的SF-SFDTD(3,4)方法用于求解含时薛定谔方程.SF-SFDTD(3,4)方法无需对传统SFDTD(3,4)方法的迭代公式进行进一步的推导,只需要在每一次的数值迭代过程中加入空间滤波操作,滤除因采用不满足CFL条件的时间步长而产生的不稳定空间域高频分量,保证数值方法的稳定性,因此所提方法与传统SFDTD(3,4)方法具有较高的兼容性.同时,理论分析了SF-SFDTD(3,4)方法的数值色散误差.最后,通过数值算例验证了本文所提方法的正确性和有效性.

基于稀疏采样的FDTD/TDPO混合优化算法

基于时域有限差分(finite-difference time-domain,FDTD)和时域物理光学(time-domain physical optics,TDPO)的混合算法提出用稀疏采样优化近场外推的方法,解决了大口径反射面天线电磁计算用时过长、计算效率低的问题。传统混合算法中馈源的近场外推到反射面天线的步骤计算量较大,稀疏采样法从空间和时间两个维度进行间隔采样,减少了此步骤的计算量,从而节省了用时。给出了馈源为喇叭天线、抛物面口径为100个波长的抛物面天线算例,所得计算结果与传统方法相比最大相对误差仅有-36.03 dB,用时最快可缩短为原来的2.51%,与CST软件计算结果对比吻合度较高,用时缩减了93.18%,验证了本文优化算法可在不影响计算精度的前提下提高计算效率。

Finite Difference Solution of Response Time Delay of Magneto-rhelological Damper

Magneto-rhelological（MR） dampers are devices that employ rheological fluids to modify their mechanical properties. Their mechanical simplicity, high dynamic range, lower power requirements, large force capacity, robustness and safe maimer of operation in cases of failure have made them attractive devices for semi-active real-time control in civil, aerospace and automotive applications. Time response characteristic is one of the most important technical performances of MR dampers, and response time directly affects the control frequency, application range and the actual effect of MR dampers. In this study, one kind of finite difference solution for predicting the response time ofmagneto-rheological dampers from ＂off-state＂ to ＂on-state＂ is put forward. A laminar flow model is used to describe the flow in the MR valve, and a bi-viscous fluid flow model is utilized to describe the relationship of shear stress and shear rate of MR fluid. An explicit difference format is used to diseretize the Novior-Stoks equation, and stability condition of this algorithm is built by Von-Neumann stability criterion. The pressure gradient along the flow duct is solved by a dichotomy algorithm with iteration, and the changing curve of the damping force versus time of MR damper is obtained as well. According to the abovementioned numerical algorithm, the damping forces versus time curves from ＂off-state＂ to ＂on-state＂ of a cylindrical piston type MR damper are computed. Moreover, the MR damper is installed in a material test system（MTS）, the magnetic field in the wire circles of the MR damper is ＂triggered＂ when the MR damper is imposed to do a constant speed motion, and the damping force curves are recorded. The comparison between numerical results and experimental results indicates that this finite difference algorithm can be used to estimate the response time delay of MR dampers.