An adaptive finite-difference method for seismic traveltime modeling based on 3D eikonal equation

3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.

A wavelet finite-difference method for numerical simulation of wave propagation in fluid-saturated porous media

In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.

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.

USE OF FINITE-DIFFERENCE TIME-DOMAIN METHOD FOR CALCULATING EM ABSORPTION IN LOSSY DIELECTRIC SCATTERER

The problem for calculating EM energy absorption by lossy dielectric scatterer ir-radiated by plane wave are discussed.The factors affecting the accuracy of computation arediscussed.The calculated results of EM energy absorption and its distribution in homogeneousand layered homogenous lossy dielectric spheres are presented,and a comparison of these resultswith analytical solution is given.The calculation is carried out for dielectric cylinder on conduct-ing ground as well,and the results are compared with the image theory.All the computationsshew that the finite-difference time-domain method can give satisfactory results.

Numerical Solution of a Problem of Thermal Stresses of a Magnetothermoelastic Cylinder with Rotation by Finite-Difference Method

The present article deals with the investigation thermal stress of a magnetothermoelastic cylinder subjected to rotation,open or closed circuit,thermal and mechanical boundary conditions.The outer and inner surfaces of the cylinder are subjected to both mechanical and thermal boundary conditions.A The transient coupled thermoelasticity in an infinite cylinder with its base abruptly exposed to a heat flux of a decaying exponential function of time is devised solve by the finite-difference method.The fundamental equations’system is solved by utilizing an implicit finite-difference method.This current method is a second-order accurate in time and space;it is also unconditionally stable.To illustrate the present model’s efficiency,we consider a suitable material and acquire the numerical solution of temperature,displacement components,and the components of stresses with time t and through the radial of an infinite cylinder.The results indicate that the effect of coupled thermoelasticity,magnetic field,and rotation on the temperature,stresses,and displacement is quite pronounced.In order to illustrate and verify the analytical developments,the numerical solution of partial differential equations,stress components,displacement components and temperature is carried out and computer simulated results are presented graphically.This study is helpful in the development of piezoelectric devices.

P-and S-wavefield simulations using both the firstand second-order separated wave equations through a high-order staggered grid finite-difference method

In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.

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.

3D finite-difference modeling algorithm and anomaly features of ZTEM

The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the Maxwell＇s equations, and the magnetic components at the center of each edge of the grid cell are evaluated by applying the staggered-grid finite-difference method. The tipper and its divergence are derived to complete the 3D ZTEM forward modeling algorithm. A synthetic model is then used to compare the responses with those of 2D finite-element forward modeling to verify the accuracy of the algorithm. ZTEM offers high horizontal resolution to both simple and complex distributions of conductivity. This work is the theoretical foundation for the interpretation of ZTEM data and the study of 3D ZTEM inversion.

An Investigation of Buckled of Column on Liner Elastic foundation by Using Finite-Difference Mode

The use of columns on elastic foundation is very common in Civil Engineering, like bridge pier, the foundation of the buildings etc. So, it will be useful to find the critical load for the structure, the problem in this paper will be solved by Finite-Difference Mode, that＇ s simple and has an extensive use. The way it works is that by dividing the component into many units. Finite-difference methods （FDM） are numerical methods for anoroximating, the solutions to differential eauations usine finite difference equations to approximate derivatives.

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.

Improved absorbing boundary condition based on linear interpolation for ADI-FDTD method

With the linear interpolation method, an improved absorbing boundary condition（ABC）is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain （ADI-FDTD） method. The reflection of the ABC caused by both the truncated error and the phase velocity error is analyzed. Based on the phase velocity estimation and the nonuniform cell, two methods are studied and then adopted to improve the performance of the ABC. A calculation case of a rectangular waveguide which is a typical dispersive transmission line is carried out using the ADI-FDTD method with the improved ABC for evaluation. According to the calculated case, the comparison is given between the reflection coefficients of the ABC with and without the velocity estimation and also the comparison between the reflection coefficients of the ABC with and without the nonuniform processing. The reflection variation of the ABC under different time steps is also analyzed and the acceptable worsening will not obscure the improvement on the absorption. Numerical results obviously show that efficient improvement on the absorbing performance of the ABC is achieved based on these methods for the ADI-FDTD.

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.

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.

Construction of Conservative Numerical Fluxes for the Entropy Split Method

The entropy split method is based on the physical entropies of the thermally perfect gas Euler equations.The Euler flux derivatives are approximated as a sum of a conservative portion and a non-conservative portion in conjunction with summation-by-parts(SBP)difference boundary closure of(Gerritsen and Olsson in J Comput Phys 129:245-262,1996;Olsson and Oliger in RIACS Tech Rep 94.01,1994;Yee et al.in J Comp Phys 162:33-81,2000).Sj?green and Yee(J Sci Comput)recently proved that the entropy split method is entropy conservative and stable.Stand-ard high-order spatial central differencing as well as high order central spatial dispersion relation preserving(DRP)spatial differencing is part of the entropy stable split methodol-ogy framework.The current work is our first attempt to derive a high order conservative numerical flux for the non-conservative portion of the entropy splitting of the Euler flux derivatives.Due to the construction,this conservative numerical flux requires higher oper-ations count and is less stable than the original semi-conservative split method.However,the Tadmor entropy conservative(EC)method(Tadmor in Acta Numerica 12:451-512,2003)of the same order requires more operations count than the new construction.Since the entropy split method is a semi-conservative skew-symmetric splitting of the Euler flux derivative,a modified nonlinear filter approach of(Yee et al.in J Comput Phys 150:199-238,1999,J Comp Phys 162:3381,2000;Yee and Sj?green in J Comput Phys 225:910934,2007,High Order Filter Methods for Wide Range of Compressible flow Speeds.Proceedings of the ICOSAHOM09,June 22-26,Trondheim,Norway,2009)is proposed in conjunction with the entropy split method as the base method for problems containing shock waves.Long-time integration of 2D and 3D test cases is included to show the com-parison of these new approaches.

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,

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.

电磁计算与电磁工业软件

概要介绍了电磁计算与电磁工业软件的研究进展。1)介绍了电磁计算方法的研究进展,对解析方法、高频近似方法、全波计算方法、快速计算方法,以及近年发展起来的混合方法进行了阐述;2)概述了基于时域有限差分方法(FDTD)、有限元法(FEM)、矩量法(MoM)等全波计算方法的3款国际主流商业电磁计算软件的特点和发展历程;3)介绍了我国具有自主知识产权的“摇光”高频电磁工业软件的算法体系、软件功能以及实际应用情况;4)对电磁计算与电磁工业软件的未来发展进行了总结和展望。