Gelfand-type problem for turbulent jets: sharp $L^\infty$ bound on extremal solutions

Sintesi Tecnica: Limiti $L^\infty$ accurati per le soluzioni estreme di problemi di tipo Gelfand in getti turbolenti

Enunciato del problema
Il saggio investiga un problema al contorno che modella l'esplosione termica in getti turbolenti reattivi, originariamente derivato da [GHH]. Il problema è formulato sul disco unitario $B \subset \mathbb{R}^2$ e appartiene alla classe dei problemi di tipo Gelfand. L'equazione governante è:
$$\begin{cases} -\Delta u - \alpha r \phi(r) \frac{\partial u}{\partial r} = \lambda \psi(r) f(u) & \text{in } B, \\ u > 0 & \text{in } B, \\ u = 0 & \text{su } \partial B, \end{cases}$$
dove $u$ rappresenta la temperatura normalizzata, $\lambda$ è il parametro di Frank-Kamenetskii e $\alpha$ è la velocità di iniezione. Le funzioni $\phi$ e $\psi$ descrivono i profili della velocità del flusso e della componente reattiva, rispettivamente, mentre $f$ è un tasso di reazione positivo, crescente e convesso che soddisfa la condizione di integrabilità $\int_0^\infty ds/f(s) < \infty$.

Il lavoro precedente di [GMN] ha stabilito l'esistenza di un parametro critico $\lambda^*(\alpha)$ e di una soluzione estrema $u^*_\alpha$ (il limite delle soluzioni minime per $\lambda \nearrow \lambda^*$). Sebbene [GMN] avesse dimostrato che $u^*_\alpha(0) \to \infty$ per $\alpha \to \infty$, essi potevano stabilire solo un limite superiore su $u^*_\alpha(0)$ sotto ulteriori assunzioni riguardanti la non linearità $f$. L'obiettivo centrale di questo saggio è colmare tale lacuna derivando limiti superiore e inferiore corrispondenti per $u^*_\alpha(0)$ nel limite asintotico di flussi forti ($\alpha \to \infty$) senza imporre condizioni aggiuntive su $f$.

Metodologia
Gli autori sfruttano la simmetria radiale della soluzione estrema per ridurre l'equazione alle derivate parziali a un'equazione differenziale ordinaria. Introducendo la funzione peso $\mu(r) = \exp(\alpha \int_0^r s\phi(s)ds)$ e compiendo un cambio di variabili $z = \alpha r^2$, il problema viene trasformato in un'equazione integrale per la soluzione riscalata $\bar{u}(z)$:
$$\bar{u}'(z) = -\frac{1}{z\bar{\mu}(z)} \int_0^z \bar{\mu}(s)\bar{\psi}(s)F(\bar{u}(s))ds,$$
dove $F = \frac{\lambda^*}{4\alpha}f$.

Il nucleo dell'analisi si basa sulla condizione di semi-stabilità soddisfatta dalla soluzione estrema. Nelle coordinate riscalate, questa condizione assume la forma di una disuguaglianza di tipo Poincaré pesata:
$$\int_0^\alpha z (\bar{\eta}')^2 \bar{\mu} dz \geq \int_0^\alpha \bar{\eta}^2 \bar{\psi} \bar{\mu} F'(\bar{u}) dz,$$
per ogni $\eta \in H^1_0(B)$.

La strategia di dimostrazione prevede due componenti principali:

1. Limite Superiore (Teorema 1.1): Gli autori costruiscono una specifica funzione di test $\bar{\eta}_z$ (costante a tratti/logaritmica) per utilizzare la disuguaglianza di semi-stabilità. Ciò consente di limitare l'integrale della derivata del termine di reazione, $\int F'(\bar{u})$, tramite un termine logaritmico. Combinando ciò con la forma integrale dell'equazione differenziale ordinaria e la convessità di $f$, derivano una disuguaglianza che mette in relazione $f'(u(0))$ con $\log \alpha$.
2. Limite Inferiore (Teorema 1.2): Gli autori impiegano un argomento per contraddizione. Essi assumono che il limite inferiore fallisca, il che implica l'esistenza di una perturbazione $v > 0$ tale che $\bar{u} + v$ rimanga una soluzione dell'equazione. Analizzando l'operatore linearizzato e utilizzando una stima del nucleo (Lemma 4.2) derivata dalle proprietà di $\mu$ e $\psi$, dimostrano che per $\alpha$ sufficientemente grande, tale perturbazione positiva non può esistere senza violare l'estremità di $\lambda^*$. Nello specifico, dimostrano che se il limite fosse falso, si potrebbe costruire una soluzione corrispondente a un parametro $\lambda > \lambda^*$, il che contraddice la definizione di parametro critico.

Contributi Chiave e Risultati
Il saggio stabilisce i seguenti limiti asintotici accurati per la soluzione estrema nell'origine, $u^*_\alpha(0)$, per $\alpha \to \infty$:

• Teorema 1.1 (Limite Superiore): Esistono costanti $\alpha_0, c, C > 0$ tali che per ogni $\alpha \geq \alpha_0$,
  $$u^*_\alpha(0) \leq c A,$$
  dove $A$ è l'unica soluzione di $f'(A) = C \log \alpha$.

• Teorema 1.2 (Limite Inferiore): Esistono costanti $\tilde{\alpha}_0, \tilde{c}, \tilde{C} > 0$ tali che per ogni $\alpha \geq \tilde{\alpha}_0$,
  $$u^*_\alpha(0) \geq \tilde{c} \tilde{A},$$
  dove $\tilde{A}$ è l'unica soluzione di $f'(\tilde{A}) = \tilde{C} \log \alpha$.

Questi risultati confermano che la crescita del picco di temperatura della soluzione estrema è governata dall'inversa della derivata del tasso di reazione valutata su una scala logaritmica della velocità del flusso.

Significatività
Il saggio sostiene di risolvere una questione parzialmente aperta lasciata da [GMN] riguardante il comportamento $L^\infty$ delle soluzioni estreme nei modelli di getti turbolenti. La significatività primaria risiede nella rimozione delle "ulteriori assunzioni sul tasso di reazione non lineare" che erano precedentemente necessarie per ottenere il limite superiore. Dimostrando che i limiti superiore e inferiore corrispondenti valgono per la classe generale di tassi di reazione convessi che soddisfano la standard condizione di integrabilità, gli autori forniscono una caratterizzazione completa del comportamento di blow-up della soluzione estrema nel limite di flusso forte. Questo risultato consolida la comprensione teorica dei limiti di esplosione termica nei getti turbolenti reattivi all'interno del quadro di Gelfand.