Abstract

The FDFD electromagnetic model computes wave scattering by directly discretizing Maxwell's equations along with specifying the material characteristics in the scattering volume. No boundary conditions are need except for the outer grid termination absorbing boundary. We use a sparse matrix Matlab code with generalized minimum residue (GMRES) Krylov subspace iterative method to solve the large sparse matrix equation, along with the Perfectly Matched Layer (PML) absorbing boundary condition. The PML conductivity profile employs the empirical optimal value from[4-6]. The sparse Matlab-based model is about 100 times faster than a previous Fortran-based code implemented on the same Alpha-class supercomputer. The 3D FDFD model is easily manipulated, it can handle all types of layer-based geometries if the target region is less than 25% of the total computational space. Several cases have been investigated. The scattered electromagnetic fields due to spherical and elliptic mine-like TNT targets buried in simulated Bosnian soil are computed and compared to reference solutions. The field distribution of dipole in the half space is computed and compared to the half-space Born approximation forward Modeling[8].

Notes

Poster presented at the 2006 Thrust R1B Effective Forward Models Conference

Keywords

half-space Born approximation, Bosnian soil, Matlab-based model, PML, GMRES, FDFD

Subject Categories

Computer simulation--Technological innovations

Disciplines

Computer Engineering

Publisher

Bernard M. Gordon Center for Subsurface Sensing and Imaging Systems (Gordon-CenSSIS)

Publication Date

2006

Rights Holder

Bernard M. Gordon Center for Subsurface Sensing and Imaging Systems (Gordon-CenSSIS)



Click button above to open, or right-click to save.

Share

COinS