Quando il flusso compressibile esplode? Scopri il nuovo criterio di blow‑up

Immagina di osservare un gas che scorre dentro un tubo, ma con la particolarità che il tubo può far entrare e uscire il fluido in punti diversi. Questo è il contesto del sistema di Navier‑Stokes compressibile, una delle equazioni più famose per descrivere fluidi come aria o acqua. Un gruppo di ricercatori ha appena proposto un nuovo “criterio di blow‑up”, cioè una regola che indica quando la soluzione dell’equazione resta liscia o, al contrario, si “esplode” in un punto di singolarità.

Il trucco sta nel controllare tre grandezze: la densità inversa (quanto è “leggera” la materia), la velocità del fluido e la variazione spaziale della densità, indicata da ∇ₓρ. Se queste quantità rimangono limitate, la soluzione non va in crisi. Il risultato è importante perché le condizioni al contorno non sono più quelle classiche (velocità normale zero), ma includono flussi in ingresso e in uscita, più realistiche per molte applicazioni ingegneristiche.

Per dimostrare il criterio, gli autori hanno dovuto inventare un nuovo approccio: non hanno potuto usare le stime tradizionali sul derivato materiale, poiché il flusso attraversa il confine. Invece, hanno costruito stime specifiche per ∇ₓρ, riuscendo a tenere sotto controllo l’intera dinamica. Il risultato è una condizione pratica: finché la densità non diventa troppo piccola, la velocità resta moderata e il gradiente di densità è ben comportato, il flusso rimane regolare. Questo apre la porta a simulazioni più affidabili di fenomeni complessi, come l’ingresso di aria in un motore o il flusso di gas in un impianto industriale.

Il sistema di Navier‑Stokes compressibile in tre dimensioni è governato da due equazioni accoppiate: la conservazione della massa, ∂ₜρ + div(ρu) = 0, e la conservazione della quantità di moto, ∂ₜ(ρu) + div(ρu⊗u) + ∇p = μΔu + (λ+μ)∇div u. Qui ρ è la densità, u la velocità, p la pressione, μ e λ i coefficienti di viscosità. Le condizioni al contorno di tipo inflow‑outflow impongono u·n = g(x,t) con g ≠ 0 su ∂Ω, e ρ = ρ_in su Γ_in, dove Γ_in è la parte di frontiera dove g·n < 0. Questo rende il problema non omogeneo e richiede una trattazione più delicata rispetto al caso classico con velocità normale nulla.

Il criterio di blow‑up proposto si basa su una stima a priori in spazi L^p‑L^q per (ρ⁻¹, u) e su un controllo di ∇ₓρ in un Sobolev‑type norm. In pratica, si dimostra che se sup_{0≤t<T}‖ρ⁻¹(t)‖_{L^∞} + sup_{0≤t<T}‖u(t)‖_{L^∞} + ∫₀^T‖∇ₓρ(t)‖_{W^{1,q}}^p dt < ∞, allora la soluzione forte può essere prolungata oltre T. La difficoltà principale è la gestione del termine di convezione ρu·∇u vicino al bordo, dove la componente normale di u è non nulla. Per superare questo ostacolo, gli autori introducono una decomposizione di u in una parte tangenziale e una normale, applicando una stima di Korn modificata che tiene conto del flusso attraverso ∂Ω. Inoltre, sfruttano una versione estesa dell’ineguaglianza di Gagliardo‑Nirenberg per controllare ∇ₓρ in termini di ρ⁻¹ e u.

Il risultato è significativo perché collega la regolarità globale a quantità fisicamente misurabili: la densità non deve avvicinarsi a zero (evitando vacui), la velocità non deve divergere, e il gradiente di densità deve rimanere moderato. Queste condizioni sono realistiche per simulazioni numeriche di flussi compressibili con ingresso/uscita di massa, come nei motori a combustione o nei sistemi di ventilazione industriale. Inoltre, il metodo introdotto apre la strada a ulteriori studi su problemi con condizioni al contorno più complesse, come pareti porose o flussi magnetoidrodinamici.