A domain hemivariational inequality for 2D and 3D convective Brinkman-Forchheimer extended Darcy equations

Titolo: Un'Ineguaglianza Emivariazionale di Dominio per le Equazioni di Brinkman-Forchheimer Estese con Convezione (CBFeD) in 2D e 3D

Autore: Jyoti Jindal, Sagar Gautam, Manil T. Mohan.

1. Il Problema

Il lavoro si concentra sullo studio delle equazioni di Brinkman-Forchheimer estese con convezione (CBFeD) in regime non stazionario, in domini limitati bidimensionali (2D) e tridimensionali (3D). Queste equazioni descrivono il flusso di fluidi viscosi incomprimibili attraverso mezzi porosi saturi.

Il sistema matematico (1.1) è dato da:
$$\begin{cases} \partial_t y - \mu \Delta y + (y \cdot \nabla)y + \alpha|y|^{q-1}y + \beta|y|^{r-1}y + \nabla \pi = f + g, & \text{in } D \times (0, T) \\ \nabla \cdot y = 0, & \text{in } D \times [0, T) \\ y = 0, & \text{su } \partial D \times [0, T] \\ y(0) = y_0, & \text{in } D \end{cases}$$
Dove:

• $y$ è il campo di velocità, $\pi$ la pressione.
• $\mu > 0$ è il coefficiente di viscosità efficace (Brinkman).
• $\beta > 0$ è il coefficiente di Forchheimer (damping/assorbimento).
• $\alpha < 0$ rappresenta un meccanismo di pompaggio (termine di "pumping").
• $r \ge 1$ è l'esponente di assorbimento (critico se $r=3$).
• Novità del modello: Il termine di forzante $g$ non è una funzione regolare, ma soddisfa un'ineguaglianza emivariazionale di dominio: $-g(x,t) \in \partial_C j(x,t,y(x,t))$, dove $\partial_C j$ è il sottodifferenziale di Clarke di un superpotenziale $j$ localmente lipschitziano. Questo introduce non-smoothness e non-convessità direttamente all'interno del dominio (non solo al bordo), modellando fenomeni fisici complessi come attriti non lineari o leggi di contatto adesivo.

2. Metodologia

Gli autori adottano un approccio analitico rigoroso basato su diverse fasi chiave:

• Formulazione Variazionale: Viene definita una formulazione debole del problema nello spazio di Banach appropriato $W$, che combina spazi di Sobolev ($L^2(0,T;V)$, $L^{r+1}(0,T;\tilde{L}^{r+1})$) e la derivata temporale.
• Regolarizzazione del Sottodifferenziale di Clarke: Poiché il termine non-smooth $\partial_C j$ rende difficile l'applicazione diretta dei metodi classici, gli autori introducono una regolarizzazione. Il sottodifferenziale viene approssimato tramite una famiglia di funzioni monodrome $\theta_{i,\varepsilon}$ ottenute mediante convoluzione con un mollificatore. Questo trasforma il problema in una successione di problemi di dimensione finita ben posti.
• Metodo di Galerkin: Viene utilizzato uno schema di discretizzazione di Galerkin. Si proietta il problema regolarizzato su spazi di dimensione finita generati dalle prime $n$ autofunzioni dell'operatore di Stokes.
• Stime A Priori: Vengono derivate stime energetiche uniformi per le soluzioni approssimate. Un punto cruciale è il trattamento del termine di regolarizzazione $\int \theta_{i,\varepsilon} y_i dx$, gestito tramite il Lemma 3.6, che permette di ottenere stime senza imporre restrizioni severe sul parametro di viscosità $\mu$.
• Passaggio al Limite: Utilizzando il teorema di Banach-Alaoglu per la convergenza debole e il lemma di Aubin-Lions per la convergenza forte, si estrae una sottosuccessione convergente. La proprietà di monotonicità e la densità permettono di identificare il limite come soluzione del problema originale, dimostrando che il limite del termine regolarizzato appartiene al sottodifferenziale di Clarke.

3. Risultati Chiave

A. Esistenza di Soluzioni Deboli

• È stata dimostrata l'esistenza di soluzioni deboli per il sistema CBFeD non stazionario in 2D e 3D con ineguaglianze emivariazionali di dominio.
• Condizioni: I risultati valgono per tutti gli esponenti di assorbimento $r \ge 1$ e viscosità $\mu > 0$.
• Vantaggio rispetto alla letteratura: A differenza di lavori precedenti (es. [10]) che richiedevano condizioni restrittive su $\mu$ (es. $\mu > C \|l\|$), questo lavoro dimostra l'esistenza senza alcuna restrizione aggiuntiva sul coefficiente di viscosità.

B. Validità dell'Equazione Energetica

• È stato provato che le soluzioni deboli soddisfano l'uguaglianza energetica (non solo la disuguaglianza), un risultato non banale per sistemi non lineari in 3D.
• Condizioni:
  • Dimensione 2D: per ogni $r \ge 1$.
  • Dimensione 3D: per $r \ge 3$.
• Questo risultato è ottenuto approssimando la soluzione tramite spazi finiti generati dalle autofunzioni dell'operatore di Stokes, evitando l'assunzione aggiuntiva $2\beta\mu \ge 1$ spesso richiesta in letteratura per il caso critico $r=3$.

C. Unicità delle Soluzioni

• È stata dimostrata l'unicità della soluzione debole sotto condizioni specifiche:
  • 2D: Per ogni $r \ge 1$.
  • 3D: Per $r > 3$.
  • Caso Critico 3D ($r=3$): L'unicità vale se $2\beta\mu > 1$.
• L'unicità è garantita dall'effetto dissipativo del termine di assorbimento $\beta|y|^{r-1}y$, che controlla la non linearità convettiva. Il risultato non richiede restrizioni su $\mu$ diverse da quelle necessarie per la dissipazione nel caso critico.

4. Contributi e Novità

1. Estensione alle Ineguaglianze di Dominio: Mentre la maggior parte della letteratura sulle emivariational inequalities (HVI) si concentra su condizioni al contorno, questo lavoro affronta HVIs di dominio, dove la non-smoothness influenza l'intero volume del fluido.
2. Rimozione di Restrizioni su $\mu$: L'approccio analitico sviluppato (basato sul Lemma 3.6 e stime energetiche raffinate) elimina la necessità di vincoli artificiali sulla viscosità, rendendo i risultati più generali e fisicamente rilevanti.
3. Generalizzazione ai Sistemi CBFeD: Estende la teoria esistente dalle classiche equazioni di Navier-Stokes (NSE) ai sistemi CBFeD, che includono meccanismi di smorzamento e pompaggio, cruciali per la modellazione dei flussi in mezzi porosi ad alta velocità.
4. Trattamento del Caso Critico 3D: Fornisce una prova di unicità per $r=3$ in 3D con la condizione $2\beta\mu > 1$, migliorando i risultati noti che spesso richiedevano condizioni più stringenti o non trattavano l'unicità in questo regime.

5. Significato e Impatto

Questo lavoro colma un vuoto teorico significativo nell'analisi matematica dei flussi in mezzi porosi non lineari.

• Fisica: Offre un modello matematico più robusto per fluidi in mezzi porosi dove le forze di attrito o di contatto sono descritte da leggi non-smooth (es. adesione, usura) distribuite nel volume.
• Matematica: Dimostra la fattibilità di applicare il metodo di Galerkin a problemi emivariazionali di dominio complessi, fornendo un quadro metodologico che può essere esteso ad altri sistemi di equazioni alle derivate parziali non lineari.
• Applicazioni: I risultati sono rilevanti per l'ingegneria petrolifera, la geotermia e lo studio dei flussi biologici, dove la modellazione accurata della dissipazione e delle interazioni non lineari è fondamentale.

In sintesi, il paper stabilisce un nuovo standard per l'analisi di esistenza, unicità e proprietà energetiche dei flussi viscosi in mezzi porosi soggetti a condizioni di non-smoothness interne, superando limitazioni tecniche presenti negli studi precedenti.