A nonconforming method for a generalized Darcy-Forchheimer model

Titolo

Un metodo non conforme per un modello generalizzato di Darcy–Forchheimer

1. Il Problema Matematico

Il lavoro si concentra sulla discretizzazione numerica della legge di Darcy–Forchheimer generalizzata, un modello non lineare che descrive il flusso di fluidi in mezzi porosi, particolarmente rilevante per flussi ad alta velocità dove la legge di Darcy lineare non è più sufficiente.

Il problema continuo (equazione 7) è formulato come segue: trovare il flusso $\mathbf{u}$ e il potenziale $p$ tali che:
$$\begin{cases} \nabla p + \frac{\mu}{\rho} K^{-1}\mathbf{u} + \frac{\beta}{\rho} |\mathbf{u}|^{\alpha-2}\mathbf{u} = 0 & \text{in } \Omega \\ \nabla \cdot \mathbf{u} = b & \text{in } \Omega \\ \mathbf{u} \cdot \mathbf{n} = g_N & \text{su } \Gamma_N \\ p = g_D & \text{su } \Gamma_D \end{cases}$$
Dove:

• $\alpha > 2$ è l'esponente della non linearità (il caso classico quadratico corrisponde a $\alpha=3$).
• $K^{-1}$ è il tensore di permeabilità inversa, assunto con regolarità di Lebesgue inferiore rispetto ai lavori precedenti.
• Sono ammesse condizioni al contorno miste e non omogenee ($\Gamma_D$ e $\Gamma_N$).

2. Metodologia e Formulazione

Gli autori propongono una discretizzazione duale mista non conforme di ordine generale. La metodologia si basa sui seguenti pilastri:

• Formulazione Duale Mista: Il problema è formulato in termini di flusso e potenziale, mantenendo la struttura duale.
• Elementi Finiti Non Conformi:
  • Il flusso $\mathbf{u}_h$ è approssimato con spazi di polinomi a tratti di ordine $k-1$ ($V_h^{k-1}$).
  • Il potenziale $p_h$ è approssimato utilizzando elementi di Crouzeix-Raviart (CR) di ordine $k$ ($Q_h^k$). Gli elementi CR sono non conformi, ovvero la continuità tra gli elementi è richiesta solo in media (o su momenti specifici) e non punto per punto.
• Generalità: A differenza dello schema precedente di Girault e Wheeler [23] (che usava solo ordine minimo e non linearità quadratiche), questo lavoro:
  • Gestisce non linearità del tipo $(\alpha-2)$ per $\alpha > 2$.
  • Permette tensori di permeabilità con bassa regolarità ($L^{\frac{\alpha}{\alpha-2}}$).
  • Supporta ordini di approssimazione arbitrari ($k \ge 1$).

3. Contributi Chiave Teorici

Il lavoro presenta diversi risultati teorici fondamentali:

1. Ben-posedness (Benessere del problema): Viene dimostrata l'esistenza e l'unicità della soluzione sia per il problema continuo che per quello discreto, sotto ipotesi di regolarità basse sui dati e sul dominio.
2. Nuove Disuguaglianze di Traccia di Sobolev: Un contributo teorico significativo è la dimostrazione di nuove disuguaglianze di traccia per spazi "broken" (a tratti) su spazi di Sobolev. Un aspetto cruciale è che le costanti in queste disuguaglianze sono indipendenti dal grado polinomiale della discretizzazione. Questo è essenziale per l'analisi di schemi di alto ordine e per strategie $hp$-adattive.
3. Convergenza: Viene provata la convergenza forte delle soluzioni discrete alla soluzione esatta in norme naturali ($L^\alpha$ per il flusso, $L^{\alpha'}$ per il gradiente del potenziale), anche sotto ipotesi di regolarità basse.
4. Stime di Errore: Vengono derivate stime di errore di ordine generale ($O(h^s)$), assumendo una regolarità aggiuntiva per la soluzione esatta e i dati. Le stime mostrano che l'errore sul flusso e sul potenziale dipende dalla regolarità dei dati e dall'ordine $k$ scelto.

4. Risultati Numerici

Gli autori hanno implementato lo schema e testato le prestazioni attraverso esperimenti numerici su due casi di test, utilizzando due schemi iterativi per gestire la non linearità:

1. Schemi Iterativi:
   • Punto fisso standard: Funziona bene per $\alpha \le 3$ e $\alpha$ leggermente superiori, ma fallisce la convergenza per $\alpha$ grandi (es. $\alpha \ge 5.2$) entro il numero massimo di iterazioni.
   • Punto fisso rilassato (damped): Utilizzando un parametro di rilassamento $\omega$, il metodo riesce a convergere anche per valori di $\alpha$ più elevati, sebbene il numero di iterazioni aumenti con $\alpha$.
2. Osservazioni sulle Prestazioni:
   • Ordini di accuratezza: Gli schemi di ordine superiore ($k=2,3,4$) mostrano tassi di convergenza ottimali previsti dalla teoria.
   • Effetto di $\alpha$: Per $\alpha$ vicini a 2 o molto grandi, la convergenza diventa più difficile, richiedendo un'attenta regolazione del parametro di rilassamento.
   • Disaccoppiamento (Caso $k=1$): Come previsto dalla teoria (Remark 4 e Fig. 6), per l'ordine minimo ($k=1$), l'errore del flusso è decouplato dal potenziale e può essere zero fino alla precisione della macchina in certi casi, mentre per ordini superiori ($k>1$) questo disaccoppiamento non si verifica e l'errore del flusso dipende dall'errore del potenziale.

5. Significato e Implicazioni

Questo lavoro è significativo per diversi motivi:

• Estensione del modello: Supera i limiti dei modelli quadratici standard, permettendo di trattare una classe più ampia di fenomeni fisici in mezzi porosi.
• Flessibilità geometrica e materiale: La capacità di gestire tensori di permeabilità poco regolari e condizioni al contorno miste rende il metodo applicabile a problemi reali più complessi.
• Fondamento per strategie $hp$-adattive: L'indipendenza delle costanti dalle disuguaglianze di traccia rispetto al grado polinomiale è un prerequisito fondamentale per lo sviluppo futuro di strategie adattive $hp$ (che variano sia la dimensione della mesh $h$ che l'ordine polinomiale $p$) per risolvere efficientemente questi problemi non lineari.
• Validazione di metodi non conformi: Dimostra che l'uso di elementi non conformi (Crouzeix-Raviart) è non solo fattibile ma vantaggioso anche per problemi non lineari complessi, dove le discretizzazioni conformi standard potrebbero non convergere.

In sintesi, il documento fornisce un quadro matematico rigoroso e un'analisi numerica completa per una nuova classe di metodi agli elementi finiti applicati alla dinamica dei fluidi nei mezzi porosi, ponendo le basi per simulazioni più accurate ed efficienti.