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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

A 1D time-domain method for in-plane wave motions in a layered half-space

A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the layered half-space is discretized on the basis of the propagation characteristic of elastic wave according to the Snell law. Then, the finite element method with lumped mass and the central difference method are incorporated to establish 2D wave motion equations, which can be transformed into 1D equations by discretization principle and explicit finite element method. By solving the 1D equations, the displacements of nodes in any vertical line can be obtained, and the wave motions in layered half-space are finally determined based on the characteristic of traveling wave. Both the theoretical analysis and the numerical results demonstrate that the proposed method has high accuracy and good stability.

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.

Second-Order Wave Diffraction Around 3-D Bodies by A Time-Domain Method

A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.

Terahertz Time-Domain Spectroscopy method for optical parameter extraction of plastic materials

A method for extracting optical parameters of plastics materials based on terahertz time domain spectroscopy is presented. The transmission-type Terahertz Time-Domain Spectroscopy(THz TDS) system is adopted to detect the refractive index and extinction coefficient on different plastic materials. Then the corresponding spectral information is obtained by Fourier transform of the terahertz time domain waveform of the sampling points, including the corresponding amplitude and phase information of the waveform. The optical parameter extraction model is built. By using the simplex optimization method, the curves of the refractive index and extinction coefficient for the plastic material are obtained. The experimental samples are made of different plastic parallel plate materials. The experimental results show that the optimization of optical parameters can improve their extraction accuracy, and the error of refractive index is ±0.005. Extraction technology with the simplex optimization method of optical parameter based on THz TDS can help to extract the optical parameters of engineering plastics. It is of great significance for the research of terahertz nondestructive testing.

Using a Time-domain Higher-order Boundary Element Method to Simulate Wave and Current Diffraction from a 3-D Body

To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.