A Nonlocal $p$-Laplacian Interface Model with Sharp Interface

Sintesi Tecnica: Un Modello di Interfaccia Nonlocale del $p$-Laplaciano con Interfaccia Netta

Definizione del Problema
Il saggio affronta la sfida di formulare problemi di interfaccia all'interno di un quadro nonlocale. I classici problemi di interfaccia coinvolgono equazioni alle derivate parziali (PDE) definite su un dominio $\Omega$ diviso da un'interfaccia $\Gamma$, dove la soluzione o le sue derivate presentano discontinuità. Nello specifico, gli autori considerano l'estensione del problema di interfaccia di Poisson e del problema dell'interfaccia del $p$-Laplaciano con condizioni di membrana a contesti nonlocali. Una difficoltà primaria nelle formulazioni nonlocali è che le classiche condizioni di interfaccia locali (come la continuità della soluzione e i salti prescritti nel flusso normale) non possono essere applicate direttamente. Gli approcci esistenti spesso approssimano l'interfaccia netta $\Gamma$ con uno strato volumetrico di spessore $\delta$, il che può complicare la modellazione della diffusione indipendente nei sottodomini e la gestione dei dati al contorno.

Metodologia
Gli autori propongono una nuova approssimazione nonlocale basata sull'energia che mantiene un'interfaccia netta $\Gamma$. Invece di modificare la geometria del dominio, il modello utilizza specifici funzionali di energia in cui le condizioni di interfaccia sono codificate implicitamente.

1. Formulazione dell'Energia:

   • Per il caso lineare ($p=2$), il funzionale di energia locale $E(u)$ è approssimato da un funzionale nonlocale $E_\delta(u)$. Il termine di diffusione nonlocale è costruito utilizzando un kernel $\Lambda(x,y)$ che restringe le interazioni all'interno dei sottodomini $\Omega_1$ e $\Omega_2$, garantendo una diffusione nonlocale indipendente in ciascuna regione.
   • La condizione di interfaccia di Dirichlet (salto nel flusso) è gestita tramite un termine di penalità che coinvolge medie nonlocali ($K_{\delta,i} * u$) lungo l'interfaccia $\Gamma$.
   • La condizione al contorno di Neumann è incorporata naturalmente attraverso la struttura variazionale.
   • Il quadro viene esteso al caso nonlineare del $p$-Laplaciano con condizioni di membrana, dove il salto della soluzione attraverso l'interfaccia viene penalizzato. Il risultante funzionale nonlocale $F_\delta(u)$ condivide la stessa forma strutturale di $E_\delta(u)$.

2. Analisi Teorica:

   • Esistenza e Unicità: Gli autori stabiliscono l'esistenza e l'unicità dei minimizzatori per i funzionali nonlocali $E_\delta$ e $F_\delta$ utilizzando il metodo diretto del calcolo delle variazioni. Ciò si basa su una disuguaglianza di Poincaré nonlocale appena derivata che tiene conto della penalità del salto all'interfaccia.
   • $\Gamma$-Convergenza: Il nucleo del contributo teorico è la dimostrazione della $\Gamma$-convergenza. Gli autori dimostrano che, al tendere al vuoto dell'orizzonte nonlocale $\delta$ ($\delta \to 0$), i funzionali nonlocali $\Gamma$-convergono ai loro controparti locali ($E$ e $F$).
   • Compattezza: Combinata con i risultati di compattezza per gli operatori del $p$-Laplaciano nonlocale, la $\Gamma$-convergenza implica che i minimizzatori dei problemi nonlocali convergano ai minimizzatori dei problemi di interfaccia locali nella norma $L^p$.

3. Implementazione Numerica:

   • Gli autori impiegano un metodo degli elementi finiti efficiente basato sulla struttura tensoriale dei kernel gaussiani e su basi polinomiali multi-cubiche.
   • Questo approccio decompone gli integrali ad alta dimensione in prodotti di integrali monodimensionali, riducendo significativamente il costo computazionale.
   • Regole di quadratura speciali (Gauss-Legendre) sono adattate per gestire gli strati limite derivanti dalle funzioni base cubiche vicino alle interfacce delle celle.

Risultati Chiave

• Preservazione dell'Interfaccia Netta: A differenza dei precedenti approcci a vincolo volumetrico, il modello proposto mantiene un'interfaccia netta, consentendo una diffusione nonlocale indipendente in sottodomini disgiunti senza estendere i dati di interfaccia nel bulk.
• Dimostrazione di Convergenza: Prove rigorose confermano che i minimizzatori dei funzionali nonlocali convergono ai minimizzatori locali quando $\delta \to 0$. Il tasso di convergenza osservato è del primo ordine.
• Validazione Numerica: Esperimenti numerici in una e due dimensioni verificano la convergenza teorica. L'errore $L^2$ tra la soluzione nonlocale e la soluzione locale di riferimento decade linearmente con $\delta$, e i risultati sono robusti rispetto alle diverse risoluzioni della mesh.

Significatività e Rivendicazioni
Il saggio sostiene che il modello nonlocale a interfaccia netta proposto offre un quadro più naturale e flessibile per i problemi di interfaccia rispetto agli approoli di approssimazione volumetrica esistenti. Formulando il problema tramite funzionali di energia, le difficili condizioni di interfaccia del flusso normale sono soddisfatte implicitamente. Gli autori evidenziano che questo approccio:

1. Facilita l'estensione a problemi di interfaccia più complessi, come dimostrato dalla riuscita applicazione al $p$-Laplaciano nonlineare con condizioni di membrana.
2. Evita la necessità di estendere i dati di interfaccia (ad esempio, $g$) da $L^2(\Gamma)$ a $H^1(\Omega)$.
3. Fornisce un legame matematicamente rigoroso tra i modelli di interfaccia non locali e locali attraverso la $\Gamma$-convergenza, garantendo che il modello nonlocale sia un'approssimazione valida del classico problema fisico al diminuire dell'orizzonte.

Il lavoro non propone nuove applicazioni fisiche, ma fornisce uno strumento matematico e numerico robusto per la modellazione di fenomeni di interfaccia esistenti in campi come la fluidodinamica e la scienza dei materiali, dove i fenomeni non locali sono rilevanti.