Quadrature error estimates on non-matching grids in a fictitious domain framework for fluid-structure interaction problems

Sintesi Tecnica: Stime dell'Errore di Quadratura su Griglie Non Coincidenti in un Framework di Dominio Fittizio per Problemi di Interazione Fluido-Struttura

Enunciato del Problema
Il documento affronta la discretizzazione numerica di problemi di interazione fluido-struttura (FSI) utilizzando una formulazione a dominio fittizio basata su moltiplicatori di Lagrange distribuiti. In questo framework, il fluido e il solido sono discretizzati utilizzando due mesh completamente indipendenti e non coincidenti. Il solido è descritto in un framework lagrangiano su un dominio di riferimento $B$, mappato nella sua posizione attuale all'interno di un dominio euleriano fisso del fluido $\Omega$. L'accoppiamento cinematico tra la velocità del fluido e la deformazione del solido è imposto tramite un moltiplicatore di Lagrange $\lambda$ definito sul dominio di riferimento del solido.

Una sfida computazionale critica deriva dall'assemblaggio della matrice di accoppiamento (interfaccia). Questa matrice richiede l'integrazione del prodotto tra le funzioni di base del fluido (composte con la mappatura del solido) e le funzioni di base del solido sugli elementi del solido. Poiché queste funzioni sono definite su griglie non coincidenti, l'integrazione esatta richiederebbe il calcolo dell'intersezione geometrica tra il mesh del fluido e il mesh del solido mappato per costruire una regola di quadratura composta. Questo processo è computazionalmente costoso, specialmente poiché il termine di accoppiamento deve essere assemblato ad ogni passo temporale. In alternativa, un approccio approssimato integra direttamente sugli elementi del solido utilizzando regole di quadratura standard senza l'intersezione dei mesh, introducendo un errore di quadratura. Il documento investiga la validità teorica e le stime dell'errore associate a questa strategia di integrazione approssimata.

Metodologia
Gli autori analizzano un modello FSI stazionario basato su Stokes, sebbene notino che i loro risultati siano rappresentativi di casi non lineari più generali. La formulazione variazionale consiste nel trovare la velocità del fluido $u$, la pressione $p$, la mappa di deformazione del solido $X$ e il moltiplicatore di Lagrange $\lambda$ che soddisfano un sistema a punto di sella.

Lo studio si concentra su due scelte specifiche per la forma bilineare $c(\cdot, \cdot)$ utilizzata per definire l'accoppiamento:

1. $c_0$ (Accoppiamento duale): Basata sul prodotto interno $L^2(B)$, dove lo spazio del moltiplicatore è il duale di $H^1(B)$.
2. $c_1$ (Prodotto scalare): Basata sul prodotto interno $H^1(B)$, dove lo spazio del moltiplicatore è $H^1(B)$ stesso.

La metodologia procede attraverso i seguenti passaggi:

• Discretizzazione: Vengono introdotti spazi di elementi finiti per tutte le variabili. Il fluido utilizza l'elemento Bercovier–Pironneau ($P1$-iso-$P2/P1$), mentre la mappatura del solido e il moltiplicatore utilizzano elementi continui a tratti lineari.
• Integrazione Approssimata: Gli autori definiscono forme bilineari discrete in cui i termini di accoppiamento sono calcolati utilizzando regole di quadratura standard sulla mesh del solido, ignorando il fatto che le funzioni di base del fluido composte sono polinomi a tratti con discontinuità rispetto alla mesh del solido.
• Analisi dell'Errore Astratta: Seguendo il framework delle condizioni di inf-sup, gli autori derivano una stima dell'errore generale (Teorema 2) che limita la differenza tra la soluzione discreta esatta e la soluzione ottenuta con l'integrazione imprecisa. Questo limite dipende dai funzionali dell'errore di quadratura.
• Stima dell'Errore di Quadratura: Il contributo tecnico centrale consiste nel dimostrare stime specifiche per l'errore di quadratura dei termini di accoppiamento. Gli autori utilizzano il fatto che la funzione composta $v_h(X)$ è continua ma lineare a tratti sugli elementi del solido, portando a specifiche proprietà di regolarità negli spazi di Sobolev frazionari ($H^{1+s}$). Utilizzano la teoria dell'interpolazione e disuguaglianze inverse per limitare l'errore.
• Dimostrazione di Stabilità: Gli autori dimostrano che il problema discreto rimane ben posto (soddisfa le condizioni di inf-sup discrete) anche quando il termine di accoppiamento è costruito approssimativamente, a patto che vengano soddisfatte specifiche condizioni sulle dimensioni delle mesh.

Contributi Chiave e Risultati

1. Ben Posto con Accoppiamento Approssimato: Il documento dimostra che il problema discreto è stabile e ben posto anche quando la matrice di accoppiamento viene assemblata utilizzando la quadratura approssimata, senza richiedere l'intersezione delle mesh.
2. Stime dell'Errore di Quadratura:
   • Per l'accoppiamento basato su $L^2$ ($c_0$), l'errore di quadratura è limitato da $O(h_B^{3/2} |\log h_B|)$. Gli autori mostrano che il metodo approssimato fornisce tassi di convergenza ottimali a condizione che $h_B \to 0$.
   • Per l'accoppiamento basato su $H^1$ ($c_1$), l'analisi dell'errore è più complessa a causa del termine del gradiente. Il limite dell'errore include termini dipendenti dal rapporto tra le dimensioni delle mesh. Nello specifico, la convergenza ottimale richiede non solo che $h_B \to 0$, ma anche che il rapporto $h_B/h_\Omega \to 0$ (dove $h_\Omega$ è la dimensione della mesh del fluido). Se il rapporto $h_B/h_\Omega$ viene mantenuto costante, il metodo può perdere ottimalità.
3. Disuguaglianze Inverse: Gli autori derivano e utilizzano specifiche disuguaglianze inverse per gestire le norme del moltiplicatore di Lagrange nelle stime dell'errore, particolarmente nel caso della norma duale.
4. Validazione Numerica: Sono presentati quattro test numerici per validare le stime teoriche:
   • Test 1 & 2: Confrontano i tassi di convergenza per $c_0$ e $c_1$ sotto diverse strategie di raffinamento della mesh (rapporto costante vs. rapporto decrescente). I risultati confermano che $c_0$ è robusto con l'integrazione approssimata, mentre $c_1$ richiede che $h_B/h_\Omega \to 0$ per mantenere l'accuratezza.
   • Test 3: Utilizza una soluzione con un salto nel gradiente di velocità attraverso l'interfaccia (tipico della FSI), confermando il tasso di convergenza subottimale atteso di $O(h^{1/2})$ dovuto alle singolarità della soluzione.
   • Test 4: Indaga elementi di ordine superiore ($P2/P1$), mostrando che mentre $c_0$ rimane robusto, $c_1$ richiede comunque l'integrazione esatta o rapporti di mesh specifici per raggiungere tassi ottimali.

Significatività e Rivendicazioni
Il documento rivendica di fornire una rigorosa giustificazione teorica per l'uso dell'integrazione approssimata nei metodi FSI a dominio fittizio, una pratica spesso adottata per l'efficienza computazionale. Gli autori dimostrano che:

• L'approccio "impreciso" non compromette la stabilità dello schema numerico.
• L'accoppiamento basato su $L^2$ ($c_0$) è particolarmente adatto all'integrazione approssimata, poiché non richiede che le mesh del fluido e del solido vengano raffinate in un rapporto specifico per mantenere l'accuratezza.
• L'accoppiamento basato su $H^1$ ($c_1$), pur offrendo potenzialmente diverse proprietà di stabilità, impone requisiti di mesh più severi ($h_B/h_\Omega \to 0$) per evitare che gli errori di quadratura dominino l'errore della soluzione.

Gli autori concludono che i loro risultati teorici confermano le osservazioni sperimentali del loro lavoro precedente [12, 13], fornendo una solida base matematica per la scelta tra tecniche di assemblaggio esatte (con intersezione di mesh) e approssimate (senza intersezione) in base alla specifica forma di accoppiamento e configurazione della mesh.