On shifting the thermal explosion threshold by a vortical flow in dimension two

Résumé technique : Sur le déplacement du seuil d'explosion thermique par un écoulement vortical en dimension deux

Énoncé du problème
Cet article étudie une généralisation naturelle du modèle classique de Frank-Kamenetskii pour l'explosion thermique, spécifiquement au sein d'un domaine borné bidimensionnel $\Omega \subset \mathbb{R}^2$. L'étude se concentre sur la distribution de température stationnaire $u$ dans un récipient de combustion à frontière froide, régie par l'équation elliptique semi-linéaire avec un terme d'advection :
$$\begin{cases} -\Delta u + A \mathbf{v} \cdot \nabla u = \lambda f(u) & \text{dans } \Omega, \\ u = 0 & \text{sur } \partial\Omega, \end{cases}$$
où $\lambda > 0$ est le paramètre de Frank-Kamenetskii (mesurant l'intensité de la réaction), $A \in \mathbb{R}$ est l'amplitude de l'écoulement, et $\mathbf{v}$ est un écoulement vortical incompressible prescrit satisfaisant $\mathbf{v} = \nabla^\perp \psi$ et la condition de non-pénétration $\mathbf{v} \cdot \nu = 0$. Le terme de réaction $f$ est supposé être une fonction $C^2$ positive, convexe, croissante, avec une croissance suffisamment rapide à l'infini.

Le phénomène physique central est l'existence d'un seuil critique $\lambda^*(A)$. Pour $\lambda < \lambda^*(A)$, le système admet une solution minimale classique représentant un équilibre thermique stable. Pour $\lambda > \lambda^*(A)$, aucune solution classique n'existe, ce qui correspond à une explosion thermique (emballement thermique). L'article examine si un écoulement de faible amplitude peut ajuster ce seuil d'explosion. Plus précisément, il s'agit de savoir si l'on peut construire un écoulement tel que $\lambda^*(A)$ varie linéairement avec $A$ (c'est-à-dire $\lambda^*(A) = \lambda^*(0) + \theta A + o(A)$ avec $\theta \neq 0$) en inversant la direction de l'écoulement, à condition que le domaine ne soit pas un disque.

Méthodologie
Les auteurs emploient une combinaison d'analyse qualitative de problèmes aux limites elliptiques, de théorie de la perturbation et d'analyse spectrale. La méthodologie se déroule selon plusieurs étapes clés :

1. Existence et régularité des solutions extrémales : Les auteurs établissent d'abord que, sous les hypothèses données, la solution extrémale $u^*$ (la limite ponctuelle des solutions minimales lorsque $\lambda \nearrow \lambda^*(A)$) est une solution classique ($C^{2,\alpha}$). Cela étend les résultats connus pour le cas $A=0$ au cas avec advection. La preuve repose sur la construction d'inégalités intégrales imitant la condition de semi-stabilité, en utilisant la fonction propre principale de l'opérateur linéarisé et en appliquant des estimations spécifiques sur le terme de réaction $f$.
2. Unicité et dégénérescence : Il est prouvé que $u^*$ est l'unique solution classique à $\lambda = \lambda^*(A)$ et que la linéarisation du problème autour de $u^*$ est dégénérée (la valeur propre principale est nulle).
3. Différentiabilité du seuil : En utilisant le théorème des fonctions implicites et des arguments de l'alternative de Fredholm, les auteurs démontrent que le paramètre critique $\lambda^*(A)$ est une fonction dérivable de l'amplitude de l'écoulement $A$. Une formule explicite pour la dérivée $\dot{\lambda}^*(A)$ est dérivée, impliquant la solution extrémale $u^*$, la fonction propre adjacente $\tilde{\phi}_1$ et le champ d'écoulement.
4. Analyse géométrique de la dérivée : Le cœur du résultat principal consiste à analyser la dérivée à $A=0$. L'expression pour $\dot{\lambda}^*(0)$ est montrée comme dépendant de l'intégrale du produit de la fonction de courant $\psi$ et du jacobien de la solution extrémale $u$ et de la fonction propre principale $\phi_1$ du problème non perturbé ($A=0$).
5. Brisure de symétrie : Les auteurs prouvent que si le jacobien $\nabla u \cdot \nabla^\perp \phi_1$ s'annule identiquement, le domaine $\Omega$ doit être un disque (via un argument de problème aux limites surdéterminé). Par conséquent, si $\Omega$ n'est pas un disque, il existe une région où ce jacobien est non nul, permettant la construction d'une fonction de courant $\psi$ qui produit une dérivée non nulle.

Contributions et résultats clés

• Régularité des solutions extrémales : L'article établit que, sous l'hypothèse d'une croissance suffisamment rapide de $f$, la solution extrémale pour le problème de Gelfand avec écoulement vortical en 2D est toujours classique ($C^{2,\alpha}$). Cela confirme que les solutions « faibles » souvent discutées dans des dimensions supérieures ne surviennent pas dans ce cadre 2D avec ces conditions d'écoulement spécifiques.
• Différentiabilité du seuil d'explosion : Les auteurs prouvent que le seuil d'explosion $\lambda^*(A)$ est dérivable par rapport à l'amplitude de l'écoulement $A$. Ils fournissent la formule explicite :
  $$\dot{\lambda}^*(A) = \frac{\int_\Omega (\mathbf{v} \cdot \nabla u^*) \tilde{\phi}_1 \, dx}{\int_\Omega f(u^*) \tilde{\phi}_1 \, dx}$$
  où $\tilde{\phi}_1$ est la fonction propre positive de l'opérateur linéarisé adjoint.
• Théorème principal (Déplacement du seuil) : Le résultat principal (Théorème 1.1) stipule que si le domaine $\Omega$ n'est pas un disque, il existe un écoulement vortical incompressible régulier tel que le seuil d'explosion se déplace linéairement avec l'amplitude de l'écoulement. Spécifiquement, $\lambda^*(A) = \lambda^*(0) + \theta A + o(A)$ avec $\theta \neq 0$. Cela implique qu'en inversant la direction de l'écoulement, on peut soit augmenter, soit diminuer le seuil d'explosion.
• Contrainte géométrique : L'article démontre rigoureusement que cette capacité à déplacer le seuil via un faible écoulement est impossible si le domaine est un disque, car le terme jacobien pertinent s'annule identiquement dans ce cas de symétrie.

Signification
L'article revendique une importance particulière en fournissant un cadre théorique constructif pour contrôler les limites d'explosion thermique dans les récipients de combustion par la dynamique des fluides. Il répond à une question naturelle concernant la dépendance non monotone du paramètre critique sur l'amplitude de l'écoulement observée dans les études numériques. En prouvant que le seuil peut être ajusté en inversant la direction de l'écoulement (sous réserve que la géométrie ne soit pas un disque parfait), ce travail offre une base mathématique pour la conception de dispositifs de combustion où la manipulation du flux est utilisée pour améliorer la sécurité (augmenter le seuil) ou contrôler les taux de réaction. De plus, la caractérisation détaillée de la régularité et de l'unicité des solutions extrémales en présence d'advection contribue à la compréhension plus large des problèmes de type Gelfand.

Les auteurs notent que la condition de « croissance rapide » de $f$ est principalement utilisée pour assurer la régularité de la solution extrémale et peut être assouplie pour inclure des non-linéarités typiques telles que $f(u) = e^u$ ou $f(u) = (1+u)^p$ ($p>1$). Les résultats sont strictement théoriques, reposant sur les propriétés des opérateurs elliptiques et de la théorie spectrale, sans proposer de nouveaux dispositifs expérimentaux.