Résumé : |
This thesis focuses on the Half-Space Matching Method which was developed to treat some scattering problems in complex infinite domains, when usual numerical methods are not applicable. In 2D, it consists in coupling several plane-wave representations in half-spaces surrounding the obstacle(s) with a Finite Element computation of the solution in a bounded domain. To ensure the matching of all these representations, the traces of the solution are linked by Fourier-integral equations set on the infinite boundaries of the half-spaces. In the case of a dissipative medium, this system of integral equations was proved to be coercive plus compact in an $L^2$ framework.
In the present thesis, we derive error estimates with respect to the discretization parameters (both in space and Fourier variables). To handle the non-dissipative case, we propose a modified version of the Half-Space Matching Method, which is obtained by applying a complex-scaling to the unknowns, in order to recover the $L^2$ framework.
We then extend the Half-Space Matching Method to scattering problems in infinite 3D elastic plates for applications to Non-Destructive Testing. The additional complexity compared to the 2D case comes from the decomposition on Lamb modes used in the half-plate representations. Due to the bi-orthogonality relation of Lamb modes, we have to consider as unknowns not only the displacement, but also the normal stress on the infinite bands limiting the half-plates. Some theoretical questions concerning this multi-unknown formulation involving the trace and the normal trace are studied in a 2D scalar case. Connections with integral methods are also addressed in the case where the Green's function is known, at least partially in each subdomain.
The different versions of the method have been implemented in the library XLiFE++ and numerical results are presented for both 2D and 3D cases. |