Direct and Inverse Scattering Problems
Work Group Inverse Problems
Prof. Dr. A. Kirsch, Dr. T. Arens, PD. Dr. F. Hettlich
Wave phenomena occur in various situations for instance as acustic waves at sound sources, as electromagnetic waves in tele communications or as elastic waves propagating from earthquakes. A common feature of all of these waves is, that they will be scattered, absorbed and/or transmitted by different media, by scattering objects'.
The mathematical modelling is done by partial differential equations like the wave equation
with spatial variable and time variable . A plane monochromatic wave with frequency and incident direction has the form
where denotes the wave number and the velocity in the medium.
The scattering of leads to a scattered field and the total field as a superposition of the incident and the scattered wave. If the scattering obstacle is constant in time, it holds
where and satisfy the reduced wave equation, or Helmholtz equation , in the exterior of the scattering obstacle. Assuming a bounded scattering object the scattered field will be of the asymptotic form
in large distance to the obstacle. Here and denote the polar coordinates of .
The picture shows level lines of far field pattern and with respect to the angles given by and in to (for two dimensional examples).
Which object belongs to which far field pattern?
- The direct scattering problem consists in showing existence, uniqueness, stability, and the numerical computation (and its graphical illustration) of the fields, if the scattering object is known.
- The invers scattering problem consists in the identification of the scattering object, if far field patterns for directions are known by measurements (i.e. if such plots are given).
Research of the work group
The direct problems can be written as integral equations using fundamental solutions, which can be helpful for theoretical questions on existence and stability as well as for numerical approximations. Currently in the working group algorithms for numerical solution of three dimensional scattering problems in case of biperiodic gratings are developed and investigated.
The fast and efficient evaluation of solutions of the direct problems is essential for using iterative methods (see below) in solving the corresponding inverse problems. Additionally, by the direct solver we obtain simulated data, which can be used to verify and compare various inversion schemes. The influence of parameter, which describe the shape or the material, can be analysed systematically by such data.
Iterative Methods: Already in case of the scattering of only one incident field at an obstacle the inverse problem can be considered. It consists in solving a non-linear and ill posed equation of the form
Iterative Newton type methods like
will require the domain derivative of the operator. Additionally a regularization is necessary, since otherwise the linearised equation has no stable inversion. By stopping the iteration early enough a good reconstruction from a far field pattern is possible. The picture shows the best, a mean, and the worst reconstruction from experiments with random noise in the data and with computed by a modified regularization scheme, which uses also the second domain derivative.
The factorization method: Iterative algorithm are all-purpose methods and lead still to the best reconstruction results. But they cost a high computational effort (in any step a direct boundary value problem has to be solved) and admit at most local convergence. Therefore, recently different approaches were established, which avoid these drawbacks. Various reconstruction algorithms based on the factorization method have been developed in the work group.
The far field pattern gives the kernel of the far field operator , defined by
with the unit sphere im . The far field operator is compact and normal in case of non-absorbing scattering objects (i.e. ) and can be factorised in the form
with certain operators and , where is compact and satisfies a Garding inequality. This is considerably used in representing the characteristic function of the set in the form
where , , and a spectral system of .
The following reconstructions are comupted by the factorization method from data: (The picture with two obstacles corresponds to the lower one of the two presented far field pattern above.)
Department of Mathematics,
Institut for Algebra and Geometry,
Work Group Inverse Problems