Brinkman's term for the steady-state Navier-Stokes equations with prescribed flux rate or pressure drop in perforated pipes

Sintesi Tecnica: Omogeneizzazione delle Equazioni di Navier-Stokes in Regime Stazionario in Tubi Perforati con Condizioni al Contorno Miste

Enunciato del Problema
Questo lavoro investiga il comportamento asintotico del moto stazionario di un fluido viscoso e incomprimibile all'interno di un tubo distorto contenente un gran numero di piccole particelle solide. Il dominio $\Omega_\varepsilon$ è un sottoinsieme perforato di un dominio limitato e ammissibile $\Omega \subset \mathbb{R}^3$, ottenuto rimuovendo $N(\varepsilon)$ particelle di diametro $\varepsilon^3$ separate da una distanza reciproca $\varepsilon$. Questa configurazione corrisponde al regime critico ($\alpha = 3$) nella teoria dell'omogeneizzazione, dove la dimensione delle particelle scala con il cubo della distanza inter-particellare.

Il flusso del fluido è governato dalle equazioni di Navier-Stokes stazionarie con condizioni al contorno miste. Nello specifico:

1. Condizioni di non scivolamento (no-slip) sono imposte sulle pareti fisiche del tubo ($\Gamma_W$) e sulle superfici delle particelle ($\partial K_n^\varepsilon$), dove le particelle possono muoversi con velocità costanti prescritte $\mu_n^\varepsilon$.
2. Condizioni di pressione di Bernoulli sono applicate all'ingresso ($\Gamma_I$) e all'uscita ($\Gamma_O$). Queste coinvolgono la derivata normale della velocità e la pressione di Bernoulli ($\Phi = p + \frac{1}{2}|u|^2$).
3. Vengono analizzati due scenari distinti:
   • Flusso Netto Prescritto: Il flusso totale $F$ attraverso l'uscita è fissato, mentre la caduta di pressione è una costante incognita determinata dalla soluzione.
   • Caduta di Pressione Prescritta: La caduta di pressione tra ingresso e uscita è fissata, mentre il flusso risultante è una quantità incognita.

La principale sfida matematica consiste nel derivare le equazioni efficaci (omogeneizzate) per il limite per $\varepsilon \to 0$ senza imporre assunzioni di piccolezza sulla forza esterna, sul tasso di flusso o sulla caduta di pressione. Ciò è complicato dal fatto che le standard stime di energia per le equazioni di Navier-Stokes spesso non forniscono limiti uniformi in questi problemi con valori al contorno misti senza ulteriori vincoli.

Metodologia
Gli autori utilizzano il metodo dell'energia all'interno del framework della teoria dell'omogeneizzazione, seguendo la strategia generale di Tartar. La dimostrazione procede in tre fasi principali:

1. Limiti Uniformi (stime $\varepsilon$-indipendenti):
   Stabilire i limiti uniformi per le soluzioni $(u_\varepsilon, \Phi_\varepsilon)$ è il passo più critico e innovativo. A differenza dei casi standard in cui bastano semplici stime di energia, gli autori non possono fare affidamento esclusivamente sul test dell'equazione del momento con la soluzione stessa a causa delle condizioni al contorno miste e della caduta di pressione incognita.

   • Invece, utilizzano un argomento per contraddizione basato sulla legge di Bernoulli per le soluzioni delle equazioni di Euler stazionarie.
   • Assumono che le norme delle soluzioni esplodano quando $\varepsilon \to 0$, normalizzano la sequenza e mostrano che il limite soddisfa le equazioni di Euler.
   • Analizzando il comportamento al contorno della pressione di Bernoulli limite sulle pareti del tubo (utilizzando risultati di regolarità locale e la geometria specifica di ingresso/uscita), dimostrano che la pressione limite deve essere costante sulle pareti. Ciò porta a una contraddizione con le condizioni di flusso o di caduta di pressione, provando così l'uniformità dei limiti della sequenza.
   • Per il problema del flusso prescritto, viene fatta un'ulteriore assunzione geometrica secondo cui le particelle non si accumulano vicino ai cilindri di ingresso e di uscita per facilitare l'analisi del contorno della pressione limite.

2. Costruzione di Correttori e Operatori di Restrizione:
   Per gestire il dominio perforato, gli autori costruiscono correttori di Stokes $(\omega_{n,1}^\varepsilon, \omega_{n,2}^\varepsilon)$ all'esterno delle particelle e nelle regioni anulari che le circondano. Questi correttori sono utilizzati per definire un operatore di restrizione $R_\varepsilon$.

   • A differenza dei lavori precedenti nel regime subcritico dove bastavano correttori capacitari scalari, il regime critico richiede questi correttori vettoriali di Stokes per tenere conto della capacità relativa non nulla del banco di particelle.
   • L'operatore di restrizione $R_\varepsilon(\phi)$ mappa le funzioni di test dal dominio omogeneizzato $\Omega$ al dominio perforato $\Omega_\varepsilon$ modificandole vicino alle particelle per soddisfare le condizioni al contorno nulle, preservando al contempo la proprietà di divergenza nulla.

3. Passaggio al Limite:
   Utilizzando i limiti uniformi e l'operatore di restrizione, gli autori passano al limite nella formulazione debole delle equazioni di Navier-Stokes.

   • Il termine convettivo non lineare converge tramite convergenza forte in $L^4$ (garantita dai limiti uniformi e dagli embedding compatti).
   • L'interazione tra il fluido e le particelle genera un termine aggiuntivo nell'equazione limite. Questo termine deriva dal comportamento limite della somma delle forze di trascinamento esercitate dalle particelle, caratterizzato dalle matrici di resistenza di Stokes $G_n$.

Contributi Chiave e Risultati
L'articolo stabilisce due teoremi principali (Teoremi 2.4 e 2.5) che descrivono il limite omogeneizzato del sistema:

• Emergenza del Termine di Brinkman: Le equazioni omogeneizzate nel dominio limite $\Omega$ sono le equazioni di Navier-Stokes aumentate da un termine di Brinkman. L'equazione efficace assume la forma:
  $$-\Delta u + (\nabla \times u) \times u + \nabla \Phi + Gu - J = f$$
  dove:

  • $Gu$ rappresenta un termine di attrito proporzionale alla velocità del fluido $u$, derivante dal trascinamento collettivo delle particelle (con $G$ che è il limite della somma scalata delle matrici di resistenza di Stokes).
  • $J$ è un termine sorgente derivante dal flusso limite delle particelle stesse in movimento.
  • Le condizioni al contorno sull'ingresso e sull'uscita rimangono di tipo misto coinvolgendo la pressione di Bernoulli.

• Estensione al Regime Critico con Dati Non Nulli: I risultati estendono i precedenti risultati di omogeneizzazione (specificamente in [59] per il regime subcritico) al regime critico ($\alpha=3$). Fondamentalmente, questi risultati valgono senza alcuna assunzione di piccolezza sulla forza esterna, sul tasso di flusso o sulla caduta di pressione. Questo rappresenta un avanzamento significativo, poiché molti risultati di omogeneizzazione per Navier-Stokes si basano su assunzioni di piccoli dati per garantire l'unicità o l'esistenza di limiti.

• Gestione delle Particelle in Movimento: Il framework accoglie particelle con velocità costanti non nulle $\mu_n^\varepsilon$, a condizione che la loro energia cinetica sia uniformemente limitata (Assunzione 2.11). Ciò consente una modellazione più realistica delle interazioni fluido-particella rispetto ai modelli di particelle fisse.

Significato
L'articolo sostiene che il suo contributo più significativo sia la derivazione rigorosa della legge di Brinkman per le equazioni di Navier-Stokes stazionarie nel regime critico sotto dati non ristretti. Gli autori sottolineano che la derivazione dei limiti uniformi tramite l'argomento di contraddizione che coinvolge la legge di Bernoulli è una questione delicata, in particolare perché la complessa struttura delle perforazioni deve essere tenuta conto per garantire che l'argomento sia valido al livello $\varepsilon$.

Costruendo con successo un vero operatore di restrizione che non modifica le funzioni di test che si annullano sulle particelle e gestendo la convergenza del termine di pressione in un contesto polidisperso con velocità delle particelle non nulle, il lavoro fornisce una solida base matematica per modellare il flusso di fluidi in tubi perforati dove la densità delle particelle è sufficientemente alta da indurre un trascinamento di tipo Brinkman ma non così alta da risultare nella legge di Darcy. I risultati sono direttamente applicabili a campi come l'emodinamica, l'ingegneria petrolifera e l'aerodinamica, dove il flusso di fluido attraverso array di ostacoli solidi è un fenomeno fondamentale.