Evolving finite elements for advection diffusion with an evolving interface

Sintesi Tecnica: Elementi Finiti in Evoluzione per Advezione-Diffusione con un'Interfaccia in Evoluzione

Enunciato del Problema
Questo articolo affronta l'approssimazione numerica di equazioni paraboliche definite su un dominio stazionario $\Omega$ contenente un'interfaccia mobile $\Gamma(t)$. L'interfaccia divide il dominio in due sottodomini, $\Omega_1(t)$ e $\Omega_2(t)$. Le equazioni governative modellano il trasporto e la diffusione di campi scalari $u_1$ e $u_2$ all'interno di questi sottodomini, soggetti a:

• Dinamiche di advezione-diffusione-reazione all'interno di ciascun sottodominio con coefficienti dipendenti dal tempo ($A_i, B_i, C_i$).
• Una condizione al contorno di Dirichlet sul bordo esterno $\partial\Omega$.
• Continuità della soluzione attraverso l'interfaccia ($u_1 = u_2$).
• Un salto di flusso prescritto attraverso l'interfaccia ($A_1 \partial_{\nu} u_1 - A_2 \partial_{\nu} u_2 = g$).

Le principali sfide identificate sono l'evoluzione dei sottodomini nel tempo e la presenza del salto di flusso attraverso l'interfaccia mobile. Sebbene esistano applicazioni nella dinamica dei fluidi, nella scienza dei materiali e nella biologia cellulare, l'articolo si concentra sulla formulazione matematica e numerica del problema.

Metodologia
Gli autori propongono un approccio Lagrangiano-Euleriano Arbitrario (ALE) che utilizza Spazi di Elementi Finiti in Evoluzione (ESFEM) basati sul quadro astratto stabilito in [21].

1. Impostazione Analitica Funzionale: Il problema continuo è formulato nel contesto di spazi di Hilbert in evoluzione. Gli autori definiscono coppie compatibili di spazi $(X(t), \phi_t)$ dove $\phi_t$ è una mappa di flusso che trasporta funzioni da un dominio di riferimento al dominio in evoluzione. Ciò permette la definizione di derivate materiali deboli e la derivazione di una formulazione debole ben posta (Teorema 2.16).
2. Costruzione della Mesh:
   • Mesh Iniziale: Il dominio è inizialmente triangolato utilizzando simplex lineari a tratti. Crucialmente, i vertici delle facce dell'interfaccia sono vincolati a giacere esattamente sull'interfaccia iniziale $\Gamma(0)$ (Assunzione M5).
   • Mappatura Isoparametrica: Per ottenere accuratezza di ordine superiore, la mesh lineare è mappata su una mesh isoparametrica utilizzando un diffeomorfismo $\Psi_h$. Questa mappa spinge le facce dell'interfaccia sull'interfaccia curva esatta $\Gamma(0)$ preservando la continuità della mesh.
   • Evoluzione Temporale: La mesh evolve spostando i nodi lagrangiani secondo un campo di velocità globale noto e regolare $w$. Il flusso discreto $\Phi^h$ è costruito come un'interpolazione polinomiale del flusso continuo, garantendo che la mesh rimanga una suddivisione ammissibile conforme in ogni istante.
3. Mappatura di Sollevamento (Lift): Viene introdotto un operatore di sollevamento $\Lambda_h$ per mappare le funzioni dal dominio computazionale discreto (che approssima l'interfaccia curva) indietro al dominio continuo esatto. Ciò è essenziale per confrontare la soluzione discreta con la soluzione esatta nell'analisi dell'errore.
4. Schema Discreto: Il metodo impiega elementi finiti isoparametrici di ordine $k$. Il problema variazionale discreto è formulato utilizzando la derivata materiale discreta e i coefficienti sollevati. Lo schema è risolto utilizzando un metodo a passi temporali (nello specifico, Formule di Differenze Indietro, BDF, negli esperimenti numerici).

Contributi Chiave

• Benessere Posto: L'articolo fornisce un'impostazione analitica funzionale rigorosa che prova l'esistenza e l'unicità della soluzione per il problema continuo (Teorema 2.16), inclusi risultati aggiuntivi di regolarità sotto specifiche assunzioni sui dati.
• Quadro ESFEM Isoparametrico: Gli autori estendono il quadro astratto degli elementi finiti in evoluzione per gestire problemi a interfaccia a due lati con salti di flusso. Dettagliano la costruzione di elementi isoparametrici che approssimano accuratamente il bordo curvo, abilitando una convergenza di ordine superiore.
• Limiti di Errore Ottimali: Il principale risultato teorico (Teorema 5.1) stabilisce limiti di errore di ordine ottimale. Nello specifico, l'errore nella norma $L^2$ è limitato da $O(h^{k+1})$, e l'errore nella semi-norma $H^1$ è limitato da $O(h^k)$, dove $h$ è la dimensione della mesh e $k$ è il grado polinomiale. Questi limiti valgono a meno di una costante e corrispondono all'errore di interpolazione di una soluzione nota regolare.
• Implementazione e Verifica: L'articolo include un'implementazione robusta utilizzando il pacchetto Firedrake e GMSH. Esperimenti numerici in 2D e 3D verificano i tassi di convergenza teorici per vari ordini polinomiali.

Risultati
Gli esperimenti numerici dimostrano l'efficacia del metodo su un problema modello con una soluzione esatta nota che è continua ma non differenziabile attraverso l'interfaccia.

• Convergenza: Per la discretizzazione spaziale di ordine $k$, i risultati numerici confermano una convergenza di ordine $k+1$ nella norma $L^2$ e di ordine $k$ nella semi-norma $H^1$.
• Robustezza: I risultati valgono sia per domini 2D che 3D, anche quando le mesh sono generate indipendentemente (piuttosto che tramite raffinamento uniforme), indicando la stabilità del metodo.
• Disponibilità del Codice: Gli autori forniscono il codice di simulazione per consentire la riproducibilità e un ulteriore utilizzo da parte della comunità.

Significato e Ambito
L'articolo afferma di fornire un quadro teorico e pratico completo per la risoluzione di problemi di advezione-diffusione con interfacce mobili utilizzando elementi finiti isoparametrici di ordine superiore. Il significato risiede nel colmare il divario tra la teoria astratta degli spazi in evoluzione e l'implementazione pratica per problemi di interfaccia con salti di flusso.

Gli autori mantengono un ambito modesto riguardo al campo di velocità: assumono che un campo di velocità globale e regolare sia dato a priori. Affermano esplicitamente che il lavoro non affronta come calcolare tale campo di velocità (ad esempio, tramite re-parametrizzazione per prevenire l'annodamento della mesh), notando che tali metodi esistono in letteratura [8, 18, 19]. Inoltre, sebbene il metodo potrebbe teoricamente applicarsi a sistemi di PDE accoppiati dove l'evoluzione del dominio è parte della soluzione, l'analisi corrente non copre tali problemi accoppiati. Il focus rimane strettamente sull'approssimazione del problema di trasmissione parabolico data un'evoluzione del dominio nota.