Structure-preserving dynamical low-rank approximation for parametric elastic guided waves

Sintesi Tecnica: Approssimazione Dinamica a Basso Rango Preservante la Struttura per Onde Guida Elastiche Parametriche

1. Definizione del Problema

Il saggio affronta le sfide computazionali associate alla simulazione di onde guida elastiche per il Monitoraggio della Salute Strutturale (SHM). Nelle applicazioni SHM, come il rilevamento di corrosione o delaminazione, sono richieste simulazioni ad alta fedeltà per modellare la propagazione delle onde, la dispersione e l'interazione con i danni. Tuttavia, queste simulazioni sono computazionalmente costose, in particolare in contesti many-query (ad esempio, quantificazione dell'incertezza, studi parametrici) dove il sistema deve essere risolto per numerose combinazioni di parametri (ad esempio, frequenza e diametro del trasduttore).

Le tecniche standard di Modellazione a Ordine Ridotto (ROM) basate sulla proiezione affrontano due limitazioni primarie in questo contesto:

1. Dinamiche dominate dal trasporto: Le onde guida coinvolgono caratteristiche in movimento con elevati gradienti. Le sottospazi lineari statici (ad esempio, derivati dalla Decomposizione Proporzionale Ortogonale - POD) non riescono a tracciare efficientemente tali fronti d'onda traslanti su finestre temporali arbitrarie, richiedendo spesso ranghi elevati o complessi metodi di registrazione.
2. Preservazione della struttura Hamiltoniana: Le equazioni delle onde elastiche sono sistemi hamiltoniani che conservano l'energia. Le proiezioni ortogonali convenzionali rompono questa struttura simpatica, portando a una deriva energetica e a una perdita di fedeltà fisica nelle simulazioni a lungo termine.

Sebbene l'Approssimazione Dinamica a Basso Rango (DLRA) abbia mostrato potenziale per altre equazioni d'onda adattando dinamicamente la base, la sua applicazione alle onde guida elastiche nell'ambito SHM è rimasta inesplorata. Inoltre, le implementazioni esistenti di DLRA mantengono spesso costi computazionali elevati durante la fase online, rendendole inadatte ai vincoli in tempo reale dello SHM.

2. Metodologia

Gli autori propongono un framework ROM parametrico che preserva la struttura, il quale sfrutta la DLRA all'interno di una strategia offline/online. L'approccio è adattato alle specifiche fisiche delle onde guida elastiche, che possono essere divise in due regimi distinti:

• Regime di Carico ($0 < t < T_L$): Il sistema è eccitato da una forza tempo-dipendente. Il rango della soluzione cresce da zero man mano che l'energia viene iniettata.
• Regime di Propagazione dell'Onda ($t \ge T_L$): L'eccitazione cessa e il sistema evolve in modo conservativo. Il rango della soluzione si stabilizza e la dinamica diventa unitaria (in uno spazio di fase complessificato).

Passaggi Metodologici Chiave:

1. Formulazione Hamiltoniana Complessificata:
   Le equazioni elastodinamiche semi-discretizzate vengono riformulate in un sistema Hamiltoniano del primo ordine. Attraverso una specifica complessificazione lineare delle variabili dello spazio delle fasi ($z = Au + iB\dot{u}$), il flusso simpletico reale viene trasformato in un flusso unitario complesso ($\dot{Z} = -i\Omega Z$). Questa trasformazione semplifica l'algebra e facilita la preservazione della struttura.

2. Fase Offline (Training):

   • Regime di Carico: Una base ridotta simpatica statica ($A_l$) viene costruita utilizzando una SVD Complessa di snapshot di addestramento. Viene utilizzata una proiezione di Galerkin per risolvere il sistema ridotto durante l'eccitazione.
   • Regime di Propagazione: Una base ridotta dinamica ($U(t)$) viene costruita. Fondamentalmente, gli autori derivano una soluzione analitica in forma chiusa per l'equazione di evoluzione non lineare della base sinistra in questo regime. Tale soluzione si basa sulla decomposizione spettrale dell'operatore ridotto e sulla base iniziale al tempo $T_L$.
   • Condizione di Giunzione (Gluing Condition): Per garantire la coerenza tra la base statica ($A_l$) e la base dinamica ($U(t)$), gli autori impongono una condizione di allineamento geometrico (minimizzazione dell'angolo principale tra i sottospazi) all'interfaccia $T_L$.

3. Fase Online (Query):

   • Per nuovi parametri, il sistema viene risolto in modo economico durante il regime di carico utilizzando la base statica pre-calcolata.
   • Al tempo $T_L$, la soluzione viene proiettata sulla base dinamica.
   • Per $t > T_L$, la soluzione viene propagata analiticamente utilizzando la formula in forma chiusa derivata. Ciò elimina la necessità di integrazione temporale durante la fase di propagazione dell'onda.
   • Il campo d'onda completo viene ricostruito combinando la base tempo-dipendente con la soluzione a basso rango.

3. Contributi Chiave

• Prima Applicazione allo SHM: Il saggio introduce la prima applicazione della DLRA Simpletica ai problemi di onde guida elastiche per il Monitoraggio della Salute Strutturale.
• Propagatore Analitico: Gli autori derivano una soluzione in forma chiusa per l'equazione di evoluzione della base sinistra nel regime di propagazione in un contesto parametrico semplificato (dove i parametri influenzano solo la funzione di carico). Ciò consente un propagatore ridotto che preserva l'energia che non richiede integrazione temporale dopo la fase di carico.
• Strategia Ibrida Offline/Online: Il framework concilia efficacemente l'efficienza della logica offline/online con le equazioni d'onda hamiltoniane, trattando diversamente i regimi di carico e di propagazione (Galerkin statica vs. DLRA dinamica).
• Preservazione della Struttura: Il metodo preserva rigorosamente la struttura simpatica, garantendo la conservazione dell'energia a lungo termine senza deriva.

4. Risultati Numerici

Il metodo è stato validato su un problema di elasticità piano-deformata 2D che coinvolge onde di Lamb dispersive che interagiscono con un danno circolare (foro). I parametri variati erano il diametro del trasduttore e la frequenza centrale di carico.

• Compressione: Il metodo ha ottenuto elevati rapporti di compressione con ranghi dinamici ($r$) compresi tra 10 e 30, rispetto alla dimensione del modello a pieno ordine (FOM) ($N \approx 10.800$).
• Accuratezza:
  • Gli errori di ricostruzione del campo completo sono stati bassi, compresi tra $10^{-3}$ e $10^{-2}$ (circa lo 0,8% - 7,5% a seconda dei parametri).
  • Gli errori di ricostruzione delle serie temporali dei sensori sono rimasti inferiori al 4%.
• Conservazione dell'Energia: Il modello ridotto ha mostrato un'eccellente conservazione dell'energia a lungo termine, con una deriva dell'energia relativa che è rimasta al livello della precisione di macchina durante il regime conservativo.
• Accelerazione Computazionale:
  • Rispetto al Modello a Pieno Ordine (FOM), il DLRA-ROM ha ottenuto accelerazioni di due o tre ordini di grandezza (fino a 820$\times$ per i dati dei sensori e 1000$\times$ per la ricostruzione del campo completo).
  • Ha superato significativamente un modello ROM modale di base (che richiedeva circa 500 modi), ottenendo un miglioramento di circa 4$\times$ rispetto all'approccio modale grazie alla dimensione ridotta molto più piccola ($r=20$ contro $r=500$).

5. Significato e Limitazioni

Significato:
Il saggio dimostra che le tecniche di basso rango dinamico sono altamente efficaci per i problemi di propagazione delle onde guida, che sono notoriamente difficili per i ROM statici a causa della loro natura dominata dal trasporto. Derivando un propagatore analitico, gli autori superano l'ostacolo computazionale dell'integrazione temporale online, rendendo l'approccio praticabile per applicazioni SHM in tempo reale o many-query. La preservazione della struttura Hamiltoniana assicura che la fedeltà fisica della simulazione sia mantenuta per lunghe durate.

Limitazioni:
Gli autori notano esplicitamente che l'attuale framework si basa su una dipendenza parametrica semplificata, in cui i parametri influenzano solo la funzione di carico (caratteristiche del trasduttore). La metodologia non gestisce ancora le variazioni parametriche negli operatori del sistema stessi (ad esempio, proprietà dei materiali, geometria o configurazione del danno). L'estensione dell'approccio a questi casi più generali è identificata come una direzione per lavori futuri.