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

Résumé Technique : Bornes $L^\infty$ Précises sur les Solutions Extremales pour les Problèmes de Type Gelfand dans les Jets Turbulents

Énoncé du Problème
L'article étudie un problème aux limites modélisant l'explosion thermique dans les jets turbulents réactifs, initialement dérivé de [GHH]. Le problème est formulé sur le disque unité $B \subset \mathbb{R}^2$ et appartient à la classe des problèmes de type Gelfand. L'équation directrice est :
$$\begin{cases} -\Delta u - \alpha r \phi(r) \frac{\partial u}{\partial r} = \lambda \psi(r) f(u) & \text{dans } B, \\ u > 0 & \text{dans } B, \\ u = 0 & \text{sur } \partial B, \end{cases}$$
où $u$ représente la température normalisée, $\lambda$ est le paramètre de Frank-Kamenetskii et $\alpha$ est la vitesse d'injection. Les fonctions $\phi$ et $\psi$ décrivent respectivement les profils de vitesse d'écoulement et de composant réactif, tandis que $f$ est un taux de réaction positif, croissant et convexe satisfaisant la condition d'intégrabilité $\int_0^\infty ds/f(s) < \infty$.

Les travaux précédents de [GMN] ont établi l'existence d'un paramètre critique $\lambda^*(\alpha)$ et d'une solution extrémale $u^*_\alpha$ (la limite des solutions minimales quand $\lambda \nearrow \lambda^*$). Bien que [GMN] ait prouvé que $u^*_\alpha(0) \to \infty$ quand $\alpha \to \infty$, ils n'avaient pu établir qu'une borne supérieure sur $u^*_\alpha(0)$ sous des hypothèses supplémentaires concernant la non-linéarité $f$. L'objectif central de cet article est de combler cette lacune en dérivant des bornes supérieure et inférieure concordantes pour $u^*_\alpha(0)$ dans la limite asymptotique des écoulements forts ($\alpha \to \infty$) sans imposer de conditions supplémentaires sur $f$.

Méthodologie
Les auteurs exploitent la symétrie radiale de la solution extrémale pour réduire l'équation aux dérivées partielles en une équation différentielle ordinaire. En introduisant la fonction de poids $\mu(r) = \exp(\alpha \int_0^r s\phi(s)ds)$ et en effectuant un changement de variable $z = \alpha r^2$, le problème est transformé en une équation intégrale pour la solution rescalée $\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,$$
où $F = \frac{\lambda^*}{4\alpha}f$.

Le cœur de l'analyse repose sur la condition de semi-stabilité satisfaite par la solution extrémale. Dans les coordonnées rescalées, cette condition prend la forme d'une inégalité de Poincaré pondérée :
$$\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,$$
pour toute fonction $\eta \in H^1_0(B)$.

La stratégie de preuve comprend deux composantes :

1. Borne Supérieure (Théorème 1.1) : Les auteurs construisent une fonction de test spécifique $\bar{\eta}_z$ (constante par morceaux/logarithmique) pour utiliser l'inégalité de semi-stabilité. Cela leur permet de borner l'intégrale de la dérivée du terme de réaction, $\int F'(\bar{u})$, par un terme logarithmique. En combinant cela avec la forme intégrale de l'EDO et la convexité de $f$, ils dérivent une inégalité reliant $f'(u(0))$ à $\log \alpha$.
2. Borne Inférieure (Théorème 1.2) : Les auteurs emploient un argument par l'absurde. Ils supposent que la borne inférieure échoue, ce qui implique l'existence d'une perturbation $v > 0$ telle que $\bar{u} + v$ reste une solution de l'équation. En analysant l'opérateur linéarisé et en utilisant une estimation de noyau (Lemme 4.2) dérivée des propriétés de $\mu$ et $\psi$, ils démontrent que pour $\alpha$ suffisamment grand, une telle perturbation positive ne peut exister sans violer l'extrémalité de $\lambda^*$. Plus précisément, ils démontrent que si la borne était fausse, on pourrait construire une solution correspondant à un paramètre $\lambda > \lambda^*$, ce qui contredit la définition du paramètre critique.

Contributions Clés et Résultats
L'article établit les bornes asymptotiques précises suivantes pour la solution extrémale à l'origine, $u^*_\alpha(0)$, lorsque $\alpha \to \infty$ :

• Théorème 1.1 (Borne Supérieure) : Il existe des constantes $\alpha_0, c, C > 0$ telles que pour tout $\alpha \geq \alpha_0$,
  $$u^*_\alpha(0) \leq c A,$$
  où $A$ est l'unique solution de $f'(A) = C \log \alpha$.

• Théorème 1.2 (Borne Inférieure) : Il existe des constantes $\tilde{\alpha}_0, \tilde{c}, \tilde{C} > 0$ telles que pour tout $\alpha \geq \tilde{\alpha}_0$,
  $$u^*_\alpha(0) \geq \tilde{c} \tilde{A},$$
  où $\tilde{A}$ est l'unique solution de $f'(\tilde{A}) = \tilde{C} \log \alpha$.

Ces résultats confirment que la croissance du pic de température de la solution extrémale est régie par l'inverse de la dérivée du taux de réaction évaluée à une échelle logarithmique de la vitesse d'écoulement.

Signification
L'article affirme résoudre une question partiellement ouverte laissée par [GMN] concernant le comportement $L^\infty$ des solutions extrémales dans les modèles de jets turbulents. La principale signification réside dans la suppression des « hypothèses supplémentaires sur le taux de réaction non linéaire » qui étaient auparavant nécessaires pour obtenir la borne supérieure. En prouvant que les bornes supérieure et inférieure concordantes tiennent pour la classe générale de taux de réaction convexes satisfaisant la condition d'intégrabilité standard, les auteurs fournissent une caractérisation complète du comportement de blow-up de la solution extrémale dans la limite d'écoulement fort. Ce résultat consolide la compréhension théorique des limites d'explosion thermique dans les jets turbulents réactifs au sein du cadre de Gelfand.