Efficient combination of a 3D Quasi-Newton inversion algorithm and a vector dual-primal finite element tearing and interconnecting method

• Published in: Inverse problems Vol. 31; no. 8; pp. 85005 - 85025
• Main Authors: Voznyuk, I, Litman, A, Tortel, H
• Format: Journal Article
• Language: English
• Published: IOP Publishing 01.08.2015
• Subjects: adjoint field; Algorithms; Computation; domain decomposition; Electromagnetism; Engineering Sciences; FETI-DPEM; finite element; Finite element method; inverse scattering; Inversions; Mathematical analysis; Mathematics; Optimization and Control; Quasi-newton; Tearing; Three dimensional; Inverse scattering; Inverse problems
• ISSN: 0266-5611; 1361-6420
• DOI: 10.1088/0266-5611/31/8/085005

A Quasi-Newton method for reconstructing the constitutive parameters of three-dimensional (3D) penetrable scatterers from scattered field measurements is presented. This method is adapted for handling large-scale electromagnetic problems while keeping the memory requirement and the time flexibility as low as possible. The forward scattering problem is solved by applying the finite-element tearing and interconnecting full-dual-primal (FETI-FDP2) method which shares the same spirit as the domain decomposition methods for finite element methods. The idea is to split the computational domain into smaller non-overlapping sub-domains in order to simultaneously solve local sub-problems. Various strategies are proposed in order to efficiently couple the inversion algorithm with the FETI-FDP2 method: a separation into permanent and non-permanent subdomains is performed, iterative solvers are favorized for resolving the interface problem and a marching-on-in-anything initial guess selection further accelerates the process. The computational burden is also reduced by applying the adjoint state vector methodology. Finally, the inversion algorithm is confronted to measurements extracted from the 3D Fresnel database.