This paper introduces a direct parallel partial FFT-type algorithm for the numerical solutions of the two- and three-dimensional Helmholtz equations. The governing equations are discretized by high-order compact finite difference methods. The resulting discretized system is indefinite, making the convergence of most iterative methods deteriorate as frequency increases. In this situation, the parallel direct approaches are a better alternative, especially for the systems with discontinuous and singular right-hand sides. The research focuses on the efficient parallel implementation of the proposed algorithm in shared memory environments (OpenMP). The complexity and scalability of the direct parallel method are investigated on scattering problems with realistic ranges of parameters in soil and mine-like targets.
The goal of this paper is to describe a novel parallel high-resolution 3D numerical method for the solution of high-frequency electromagnetic wave propagation. The sequential numerical method was developed by the first author in 2014. The discussed parallel algorithm will be used later by the authors to computationally simulate data for the solution of the inverse problem of imaging mine-like targets. Thus the solution of the forward problem presented in this paper is a necessary prelude to the future solution of a related inverse problem. In this paper, land mines are modeled as small abnormalities embedded in an otherwise uniform media with an air-ground interface. These abnormalities are characterized by the electrical permittivity and the conductivity, whose values differ from those of the host media. The main challenge in the calculation of the scattered electromagnetic signal in these settings is the requirement of solving the Helmholtz equation for high frequencies. This is excessively time-consuming using standard direct solution techniques. A high-resolution and scalable numerical procedure for the solution of this equation is described in this paper. The kernel of this algorithm is a combination of a second, fourth or sixth order compact finite-difference scheme and a preconditioned Krylov subspace approach. Both fourth and sixth order compact approximations for the Helmholtz equation are considered to reduce approximation and pollution errors, thereby softening the point-per-wavelength constraint. The coefficient matrix of the resulting system is not Hermitian and possesses positive as well as negative eigenvalues. This represents a significant challenge for constructing an efficient iterative solver. In our approach, this system is solved by a combination of Krylov subspace-type method with a direct parallel FFT-type preconditioner. The resulting numerical method allows a natural and efficient implementation on parallel computers. Numerical results for realistic ranges of parameters in soil and mine-like targets confirm the high efficiency of the proposed parallel iterative algorithm.
Access to the requested content is limited to institutions that have purchased or subscribe to SPIE eBooks.
You are receiving this notice because your organization may not have SPIE eBooks access.*
*Shibboleth/Open Athens users─please
sign in
to access your institution's subscriptions.
To obtain this item, you may purchase the complete book in print or electronic format on
SPIE.org.
INSTITUTIONAL Select your institution to access the SPIE Digital Library.
PERSONAL Sign in with your SPIE account to access your personal subscriptions or to use specific features such as save to my library, sign up for alerts, save searches, etc.