Numerical approximation of the scattering amplitude in elasticity
We propose a numerical method to approximate the scattering amplitudes for the elasticity system with a non-constant matrix potential in dimensions d=2 and 3. This requires to approximate first the scattering field, for some incident waves, which can be written as the solution of a suitable Lippmann-Schwinger equation. In this work we adapt the method introduced by G. Vainikko in <cit.> to solve such equations when considering the Lamé operator. Convergence is proved for sufficiently smooth potentials. Implementation details and numerical examples are also given.
READ FULL TEXT