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 tra...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.展开更多
Irregular surface flattening,which is based on a boundary conforming grid and the transformation between curvilinear and Cartesian coordinate systems,is a mathematical method that can elegantly handle irregular surfac...Irregular surface flattening,which is based on a boundary conforming grid and the transformation between curvilinear and Cartesian coordinate systems,is a mathematical method that can elegantly handle irregular surfaces,but has been limited to obtaining first arrivals only.By combining a multistage scheme with the fast-sweeping method(FSM,the method to obtain first-arrival traveltime in curvilinear coordinates),the reflected waves from a crustal interface can be traced in a topographic model,in which the reflected wavefront is obtained by reinitializing traveltimes in the interface for upwind branches.A local triangulation is applied to make a connection between velocity and interface nodes.Then a joint inversion of first-arrival and reflection traveltimes for imaging seismic velocity structures in complex terrains is presented.Numerical examples all perform well with different seismic velocity models.The increasing topographic complexity and even use of a high curvature reflector in these models demonstrate the reliability,accuracy and robustness of the new working scheme;checkerboard testing illustrates the method's high resolution.Noise tolerance testing indicates the method's ability to yield practical traveltime tomography.Further development of the multistage scheme will allow other later arrivals to be traced and used in the traveltime inversion.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.61331007 and 61471105)
文摘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.
基金financial support for this work contributed by the National Key Research and Development Program of China(grant nos.2016YFC0600302,2016YFC0600101 and 2016YFC0600201)the National Natural Science Foundation of China(grants 41604075,41430213,41574092 and 41474068)
文摘Irregular surface flattening,which is based on a boundary conforming grid and the transformation between curvilinear and Cartesian coordinate systems,is a mathematical method that can elegantly handle irregular surfaces,but has been limited to obtaining first arrivals only.By combining a multistage scheme with the fast-sweeping method(FSM,the method to obtain first-arrival traveltime in curvilinear coordinates),the reflected waves from a crustal interface can be traced in a topographic model,in which the reflected wavefront is obtained by reinitializing traveltimes in the interface for upwind branches.A local triangulation is applied to make a connection between velocity and interface nodes.Then a joint inversion of first-arrival and reflection traveltimes for imaging seismic velocity structures in complex terrains is presented.Numerical examples all perform well with different seismic velocity models.The increasing topographic complexity and even use of a high curvature reflector in these models demonstrate the reliability,accuracy and robustness of the new working scheme;checkerboard testing illustrates the method's high resolution.Noise tolerance testing indicates the method's ability to yield practical traveltime tomography.Further development of the multistage scheme will allow other later arrivals to be traced and used in the traveltime inversion.
