Sobolev-Poincaré inequalities for piecewise $W^{1,p}$ functions over general polytopic meshes

1. Il Problema

Il lavoro si concentra sullo sviluppo e la dimostrazione di disuguaglianze di Sobolev-Poincaré per funzioni a tratti (piecewise) appartenenti allo spazio $W^{1,p}$, definite su famiglie di mesh poligonali/poliopici molto generali in qualsiasi dimensione $d$.

Queste disuguaglianze sono strumenti fondamentali nell'analisi numerica, in particolare per:

• La stabilità e la convergenza di metodi agli elementi finiti non conformi (nonconforming).
• L'analisi di problemi non lineari.
• Metodi basati su spazi polinomiali e non polinomiali (es. metodi Trefftz, estesi).

La sfida principale risiede nel fatto che le disuguaglianze esistenti nella letteratura spesso si limitano a mesh triangolari o cartesiane, richiedono costanti non esplicite, o sono valide solo per spazi polinomiali di ordine finito. Questo lavoro mira a generalizzare tali risultati a mesh poliopiche arbitrarie (con un numero qualsiasi di facce per elemento) e a funzioni generiche $W^{1,p}$, fornendo costanti esplicite che dipendono dalle proprietà geometriche del dominio e della mesh.

2. Metodologia

Gli autori adottano un approccio analitico che non si basa su operatori di arricchimento (come in lavori precedenti di Arnold o Brenner) né su immersioni continue di spazi a variazione limitata. Invece, la dimostrazione poggia sulla generalizzazione di due risultati fondamentali della teoria degli spazi di Sobolev:

1. Disuguaglianze di traccia di Sobolev: Gli autori estendono le disuguaglianze di traccia classiche (validhe per simplex) a famiglie di mesh poliopiche. Viene dimostrato un teorema che lega la norma $L^q$ sulla frontiera di un elemento alla norma $L^s$ e alla seminorma del gradiente all'interno dell'elemento, con costanti che dipendono esplicitamente dal parametro di regolarità della mesh $\gamma$.
2. Inversi a destra stabili dell'operatore divergenza (Disuguaglianze di Babuška-Aziz): Viene costruita un'inversa a destra dell'operatore divergenza che soddisfa condizioni al contorno miste (Dirichlet su una parte del bordo, Neumann sull'altra).
   • Per domini convessi, viene utilizzato un metodo di estensione simmetrica del dominio.
   • Per domini non convessi, viene costruita un'estensione del dominio ($\Omega_{ext}$) tale che l'intersezione con il dominio originale sia una faccia del bordo.
   • Le costanti ottenute per questi inversi sono esplicithe in termini del diametro del dominio, del raggio della sfera interna (star-shapedness) e del rapporto tra le aree dei domini estesi.

Questi due strumenti tecnici vengono combinati tramite argomenti di dualità e disuguaglianze di Hölder per derivare le disuguaglianze di Poincaré richieste.

3. Contributi Chiave e Risultati Principali

Il paper presenta una serie di teoremi e corollari che costituiscono i risultati principali:

• Disuguaglianze di Traccia di Sobolev (Teorema 1.1 e Corollario 1.2):
  Vengono stabilite nuove disuguaglianze di traccia per mesh poliopiche con costanti esplicite. Queste generalizzano i risultati classici a spazi $W^{1,p}$ e norme $L^q$ arbitrarie, gestendo la regolarità geometrica attraverso il parametro $\gamma$.

• Inverso stabile $W^{1,p}$ della divergenza con condizioni miste (Teorema 1.4):
  Viene dimostrato l'esistenza di un inverso a destra dell'operatore divergenza con condizioni al contorno miste. Un contributo significativo è la stima esplicita della costante di stabilità, che distingue tra casi convessi e non convessi, mostrando come la costante dipenda dalla geometria dell'estensione del dominio.

• Disuguaglianze di Sobolev-Poincaré di "Prima Specie" (Teorema 1.5):
  Fornisce una disuguaglianza per funzioni in $W^{1,p}(T_n)$ con esponenti di Lebesgue ottimali ($p^*$). La norma della funzione è controllata dal gradiente a tratti ($\nabla_h v$) e dai termini di salto (jump) sulle facce interne e sul bordo di Dirichlet. Le costanti sono esplicitamente dipendenti dalla geometria.

• Disuguaglianze di Sobolev-Poincaré di "Seconda Specie" (Teorema 1.6):
  Presenta una variante con norme più deboli, più adatta all'analisi di metodi DG (Discontinuous Galerkin) standard e metodi poligonali (polyDG). Questa versione utilizza pesi legati alle dimensioni delle facce ($h_F$) e permette di gestire meglio le mesh con elementi di dimensioni variabili.

• Versioni Mediate (Corollari 1.7 e 1.8):
  Vengono proposte varianti delle disuguaglianze precedenti che utilizzano la proiezione media dei salti sulle facce ($\Pi^0_{F_n} [[v]]$). Queste sono particolarmente utili per l'analisi di metodi specifici come gli elementi Crouzeix-Raviart e il metodo WOPSIP (Weakly Over-Penalty Symmetric Interior Penalty), permettendo di recuperare esponenti ottimali senza richiedere stime inverse aggiuntive che dipendano dal grado polinomiale.

4. Significato e Impatto

L'importanza di questo lavoro risiede in diversi aspetti cruciali per l'analisi numerica moderna:

1. Generalità delle Mesh: I risultati sono validi per mesh poliopiche arbitrarie, inclusi elementi con un numero arbitrario di facce (agglomerati di elementi simpliciali). Questo è essenziale per metodi come il Virtual Element Method (VEM) o il Polytopal Discontinuous Galerkin (PolyDG).
2. Costanti Esplicite: A differenza di molte dimostrazioni esistenti che si basano su argomenti di compattezza (che forniscono costanti non costruttive), questo lavoro fornisce costanti esplicite. Questo è fondamentale per l'analisi di convergenza delle versioni $p$ e $hp$ dei metodi DG, specialmente sotto ipotesi di regolarità minime.
3. Applicabilità ai Problemi Non Lineari: La formulazione in spazi $W^{1,p}$ (e non solo $H^1$) rende questi strumenti direttamente applicabili all'analisi di problemi non lineari, dove la linearità non è garantita.
4. Condizioni al Contorno Miste: La capacità di gestire condizioni al contorno miste (Dirichlet/Neumann) in modo esplicito amplia la portata dell'analisi a problemi fisici realistici con bordi complessi.

In sintesi, il documento fornisce un quadro teorico rigoroso e pratico per l'analisi di metodi numerici avanzati su mesh complesse, colmando il divario tra la teoria classica degli spazi di Sobolev e le esigenze computazionali delle mesh poliopiche moderne.