A direct imaging method for inverse scattering problem of biharmonic wave with phased and phaseless data

Résumé Technique : Une Méthode d'Imagerie Directe pour les Problèmes de Diffusion Inverse Biharmonique

Énoncé du Problème
Ce papier aborde le problème de diffusion inverse biharmonique, qui vise à déterminer la forme et la localisation d'un obstacle étendu $D$ dans $\mathbb{R}^2$ à l'aide de données de mesure. Contrairement aux problèmes de diffusion inverse plus courants associés aux opérateurs différentiels d'ordre deux (acoustique, élastique, électromagnétique), ce travail se concentre sur l'équation d'onde biharmonique, une équation aux dérivées partielles d'ordre quatre. L'étude considère à la fois des données de phase (où l'information de phase est disponible) et des données sans phase (où seul le module du champ est mesuré). L'obstacle est soumis à l'une des quatre conditions aux limites enracinées dans la théorie des plaques : encastrée, simplement appuyée, appuyée sur rouleaux ou libre.

Méthodologie
Les auteurs proposent une méthode d'imagerie directe basée sur la technique de Migration en Temps Inversé (RTM). Contrairement aux méthodes itératives, cette approche ne nécessite ni des hypothèses géométriques préalables sur le diffuseur, ni la résolution numérique du problème de diffusion directe à chaque itération.

Le cœur de la méthodologie consiste à construire des fonctions d'imagerie spécifiques ($I_j(z)$) pour les points d'échantillonnage $z$ à l'intérieur d'un domaine $\Omega$. Ces fonctions sont conçues pour présenter des caractéristiques distinctes selon que le point d'échantillonnage se trouve près ou loin de la frontière de l'obstacle.

1. Données de Phase (IP-1 à IP-10) : Les auteurs définissent dix fonctions d'imagerie utilisant diverses combinaisons de données de mesure, notamment :

   • Le champ diffusé ($u^{sc}$) et sa dérivée normale ($\partial_n u^{sc}$).
   • Le moment de flexion ($M u^{sc}$) et la force transversale ($N u^{sc}$).
   • Les diagrammes de champ lointain ($u^\infty$).
   • Des données générées par des sources ponctuelles et des ondes planes.

   Les fonctions d'imagerie sont construites en utilisant la solution fondamentale de l'équation de Helmholtz ($\Phi_\kappa$) plutôt que la solution fondamentale de l'équation d'onde biharmonique. Cette formulation vise à réduire la complexité computationnelle. Les fonctions impliquent des intégrales sur les frontières de mesure ($\Gamma_r, \Gamma_s$) ou le cercle unité ($S$), combinant les données mesurées avec la solution fondamentale de Helmholtz et ses dérivées.

2. Données Sans Phase (IP-11) : Pour le cas où seul le module du champ total $|u|$ est disponible, une fonction d'imagerie spécifique $I_{11}(z)$ est construite. Cette fonction repose sur la différence entre le carré du module du champ total et le carré du module du champ incident, normalisé par le champ incident.

Analyse Théorique et Contributions Clés
Le papier fournit une analyse mathématique rigoureuse des fonctions d'imagerie proposées :

• Analyse de Résolution : Les auteurs prouvent que pour tout point d'échantillonnage $z$, la fonction d'imagerie $I_j(z)$ peut être décomposée en un terme principal et un terme de bruit/erreur. Le terme principal est proportionnel à l'intégrale du carré du module du diagramme de champ lointain d'une solution auxiliaire $\psi(\cdot, z)$, qui satisfait l'équation d'onde biharmonique avec des données aux limites dérivées de la partie imaginaire de la solution fondamentale de Helmholtz.
• Mécanisme de Contraste : L'analyse démontre que les fonctions d'imagerie présentent un fort contraste lorsque le point d'échantillonnage $z$ est proche de la frontière de l'obstacle et décroissent lorsque $z$ s'éloigne. Ce comportement est attribué aux propriétés de la partie imaginaire de la solution fondamentale de l'équation de Helmholtz (une fonction de Bessel), qui atteint un pic lorsque la source et le point d'observation coïncident.
• Nouveauté dans la Formulation : Une contribution significative est la dérivation de fonctions d'imagerie qui ne dépendent que de la solution fondamentale de l'équation de Helmholtz (ou d'une onde plane acoustique) plutôt que de la solution fondamentale biharmonique. Les auteurs notent que les études RTM précédentes pour les ondes acoustiques, élastiques et électromagnétiques ne peuvent pas être directement étendues au cas biharmonique en raison de l'ordre supérieur de l'opérateur différentiel.
• Robustesse aux Conditions aux Limites : Le cadre théorique est établi comme applicable aux quatre types de conditions aux limites (encastrée, simplement appuyée, appuyée sur rouleaux et libre), avec des comportements asymptotiques spécifiques dérivés pour chacun.

Résultats Numériques
Le papier présente des expériences numériques pour valider les résultats théoriques :

• Reconstruction d'un Obstacle Unique : L'algorithme reconstruit avec succès un cercle unité en utilisant toutes les 11 fonctions d'imagerie (10 avec phase, 1 sans phase) sans bruit.
• Robustesse au Bruit : Des expériences avec un obstacle en forme de cerf-volant démontrent la robustesse de l'algorithme face à des niveaux de bruit de 5 % et 10 %. Les formes reconstruites restent précises et distinctes.
• Obstacles Multiples : La méthode s'avère efficace pour localiser et façonner simultanément deux obstacles distincts (un cercle et un cerf-volant), même en présence de bruit.

Signification et Revendications
Le papier revendique que la méthode d'imagerie directe proposée offre deux avantages principaux par rapport aux méthodes itératives :

1. Elle élimine le besoin d'hypothèses géométriques préalables sur le diffuseur.
2. Elle évite la charge computationnelle liée à la résolution répétée du problème de diffusion directe.

Les auteurs soulignent que leur approche est la première à appliquer l'imagerie directe basée sur la RTM aux problèmes de diffusion inverse biharmonique en utilisant un seul type de données de mesure (allant des champs diffusés aux champs totaux sans phase) tout en s'appuyant sur la solution fondamentale plus simple de Helmholtz. L'analyse de résolution confirme que la méthode fournit un contraste clair pour identifier la frontière de l'obstacle, et les expériences numériques vérifient son efficacité dans la reconstruction de géométries complexes et sa stabilité dans des environnements bruyants. Ce travail étend l'applicabilité des méthodes d'imagerie directe des équations d'onde d'ordre deux aux équations d'onde biharmoniques d'ordre quatre.