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

Technische Samenvatting: Een Directe Afbeeldingsmethode voor Inverse Biharmoonse Verspreidingsproblemen

Probleemstelling
Dit artikel behandelt het inverse biharmonische verspreidingsprobleem, dat erop gericht is de vorm en locatie van een uitgebreid obstakel $D$ in $\mathbb{R}^2$ te bepalen aan de hand van meetgegevens. In tegenstelling tot de gebruikelijkere inverse verspreidingsproblemen die geassocieerd zijn met differentiaaloperatoren van de tweede orde (akoestisch, elastisch, elektromagnetisch), richt deze studie zich op de biharmonische golfvergelijking, een partiële differentiaalvergelijking van de vierde orde. De studie beschouwt zowel gefaseerde data (waarbij fase-informatie beschikbaar is) als fasevrije data (waarbij slechts de modulus van het veld wordt gemeten). Het obstakel is onderworpen aan een van de vier randvoorwaarden die zijn geworteld in plaattheorie: ingeklemde, simpel ondersteunde, rollend ondersteunde of vrije randvoorwaarden.

Methodologie
De auteurs stellen een directe afbeeldingsmethode voor die gebaseerd is op de techniek van Reverse Time Migration (RTM). In tegenstelling tot iteratieve methoden vereist deze aanpak geen voorafgaande geometrische aannames over de verstrooier, noch is het noodzakelijk om bij elke iteratie het directe verspreidingsprobleem numeriek op te lossen.

De kern van de methodologie bestaat uit het construeren van specifieke afbeeldingsfuncties ($I_j(z)$) voor de bemonsteringspunten $z$ binnen een domein $\Omega$. Deze functies zijn ontworpen om onderscheidende kenmerken te vertonen, afhankelijk van of het bemonsteringspunt zich dicht bij of ver van de rand van het obstakel bevindt.

1. Gefaseerde Data (IP-1 tot IP-10): De auteurs definiëren tien afbeeldingsfuncties die gebruikmaken van verschillende combinaties van meetgegevens, waaronder:

   • Het verstrooide veld ($u^{sc}$) en zijn normale afgeleide ($\partial_n u^{sc}$).
   • Het buigend moment ($M u^{sc}$) en de dwarskracht ($N u^{sc}$).
   • Ver-veldpatronen ($u^\infty$).
   • Data gegenereerd door puntbronnen en vlakke golven.

   De afbeeldingsfuncties worden geconstrueerd met behulp van de fundamentele oplossing van de Helmholtz-vergelijking ($\Phi_\kappa$) in plaats van de fundamentele oplossing van de biharmonische golfvergelijking. Deze formulering is bedoeld om de rekencomplexiteit te verminderen. De functies omvatten integralen over de meetranden ($\Gamma_r, \Gamma_s$) of de eenheidscirkel ($S$), waarbij de gemeten data wordt gecombineerd met de fundamentele oplossing van Helmholtz en zijn afgeleiden.

2. Fasevrije Data (IP-11): Voor het geval dat slechts de modulus van het totale veld $|u|$ beschikbaar is, wordt een specifieke afbeeldingsfunctie $I_{11}(z)$ geconstrueerd. Deze functie is gebaseerd op het verschil tussen het kwadraat van de modulus van het totale veld en het kwadraat van de modulus van het invallende veld, genormaliseerd door het invallende veld.

Theoretische Analyse en Belangrijkste Bijdragen
Het artikel biedt een rigoureuze wiskundige analyse van de voorgestelde afbeeldingsfuncties:

• Resolutieanalyse: De auteurs bewijzen dat voor elk bemonsteringspunt $z$ de afbeeldingsfunctie $I_j(z)$ kan worden ontbonden in een hoofdterm en een ruis/foutterm. De hoofdterm is evenredig met de integraal van het kwadraat van de modulus van het ver-veldpatroon van een hulpoplossing $\psi(\cdot, z)$, die voldoet aan de biharmonische golfvergelijking met randgegevens die zijn afgeleid uit het imaginaire deel van de fundamentele oplossing van Helmholtz.
• Contrastmechanisme: De analyse toont aan dat de afbeeldingsfuncties een hoog contrast vertonen wanneer het bemonsteringspunt $z$ zich dicht bij de rand van het obstakel bevindt, en afnemen naarmate $z$ verder weg beweegt. Dit gedrag wordt toegeschreven aan de eigenschappen van het imaginaire deel van de fundamentele oplossing van de Helmholtz-vergelijking (een Besselfunctie), die piekt wanneer de bron- en observatiepunten samenvallen.
• Nieuwheid in Formulering: Een significante bijdrage is de afleiding van afbeeldingsfuncties die uitsluitend afhankelijk zijn van de fundamentele oplossing van de Helmholtz-vergelijking (of een akoestische vlakke golf) in plaats van de biharmonische fundamentele oplossing. De auteurs merken op dat eerdere RTM-studies voor akoestische, elastische en elektromagnetische golven niet direct kunnen worden uitgebreid naar het biharmonische geval vanwege de hogere orde van de differentiaaloperator.
• Robuustheid ten opzichte van Randvoorwaarden: Het theoretische kader is zodanig opgezet dat het toepasbaar is op alle vier de typen randvoorwaarden (ingeklemd, simpel ondersteund, rollend ondersteund en vrij), met specifieke asymptotisch gedrag afgeleid voor elk.

Numerieke Resultaten
Het artikel presenteert numerieke experimenten ter validatie van de theoretische bevindingen:

• Reconstructie van een Enkel Obstakel: Het algoritme reconstrueert succesvol een eenheidscirkel met alle 11 afbeeldingsfuncties (10 gefaseerde, 1 fasevrije) zonder ruis.
• Robuustheid tegen Ruis: Experimenten met een vlieger-vormig obstakel tonen de robuustheid van het algoritme tegen ruissniveaus van 5% en 10%. De gereconstrueerde vormen blijven nauwkeurig en duidelijk onderscheidbaar.
• Meerdere Obstakels: De methode blijkt effectief twee verschillende obstakels (een cirkel en een vlieger) gelijktijdig te kunnen lokaliseren en van vorm te kunnen voorzien, zelfs in aanwezigheid van ruis.

Betekenis en Claims
Het artikel claimt dat de voorgestelde directe afbeeldingsmethode twee primaire voordelen biedt ten opzichte van iteratieve methoden:

1. Het elimineert de noodzaak van voorafgaande geometrische aannames over de verstrooier.
2. Het vermijdt de rekenkundige last van het herhaaldelijk oplossen van het directe verspreidingsprobleem.

De auteurs benadrukken dat hun aanpak de eerste is die RTM-gebaseerde directe afbeelding toepast op biharmonische inverse verspreidingsproblemen met behulp van één type meetdata (variërend van verstrooide velden tot fasevrije totale velden), terwijl het vertrouwt op de eenvoudigere fundamentele oplossing van Helmholtz. De resolutieanalyse bevestigt dat de methode een duidelijk contrast biedt voor het identificeren van de rand van het obstakel, en de numerieke experimenten verifiëren de efficiëntie bij het reconstrueren van complexe geometrieën en de stabiliteit in omgevings met ruis. Het werk breidt de toepasbaarheid van directe afbeeldingsmethoden uit van golfvergelijkingen van de tweede orde naar biharmonische golfvergelijkingen van de vierde orde.