Iterative forward and inverse algorithms based on domain integral equations for three-dimensional electric and magnetic objects Aria Abubakar a, * , Peter M. van den Berg b a Schlumberger-Doll Research, 36 Old Quarry Road, Ridgefield, CT 06877-4108, USA b Faculty of Applied Physics, Delft University of Technology, Lorentzweg 1, 2628 CJ, Delft, Netherlands Received 27 November 2002; received in revised form 23 July 2003; accepted 3 October 2003 Abstract An iterative approach to full vector three-dimensional inverse scattering problems, where the unknown objects can have conductivity, permittivity and permeability different from the known background medium, is discussed. Since this problem involves a large number of unknowns, it has to be solved effectively and efficiently so that the results can be obtained in timely manner. The forward modeling is based on a domain integral equation approach formulated in terms of the electric and magnetic contrast sources normalized with the characteristic impedance of the background medium. Our numerical tests indicate that this formulation is prerequisite in order to arrive at a forward solution within an acceptable number of iterations, and hence it is also of significant importance in the optimization process of the inverse problem. The inverse scattering problem is attacked using the Multiplicative Regularized Contrast Source Inversion method as known in the literature. The complexity of this inverse method is approximately equal to the complexity of two equivalent forward algorithms of the conjugate gradient type. Furthermore, this inverse method has been armed with a weighted L 2 -norm regularizer which has been included as a multiplicative constraint. Some representative numerical testings will be presented to illustrate the ability of the our numerical algorithms. Ó 2003 Elsevier Inc. All rights reserved. Keywords: Inverse scattering; Three-dimensional; Full-vectorial; Iterative technique; Nonlinear problem; Integral equation; Electric and magnetic objects 1. Introduction Two- and three-dimensional electromagnetic inverse scattering problem have been studied with various optimization methods, in particular, deterministic gradient type methods. Among the huge amount of literature we only list a few without intention of review: Abubakar and co-authors [2–7,9,30], Barkeshli and Journal of Computational Physics 195 (2004) 236–262 www.elsevier.com/locate/jcp * Corresponding author. Tel.: +1-203-431-5559; fax: +1-203-438-3819. E-mail addresses: aabubakar@ridgefield.oilfield.slb.com, aabubakar@slb.com (A. Abubakar). 0021-9991/$ - see front matter Ó 2003 Elsevier Inc. All rights reserved. doi:10.1016/j.jcp.2003.10.009