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

Technische Samenvatting: Scherpe $L^\infty$-bounds op extremale oplossingen voor Gelfand-type problemen in turbulente jets

Probleemstelling
Het artikel onderzoekt een randwaardeprobleem dat thermische explosie in reactieve turbulente jets modelleert, oorspronkelijk afgeleid uit [GHH]. Het probleem is geformuleerd op de eenheidsschijf $B \subset \mathbb{R}^2$ en behoort tot de klasse van Gelfand-type problemen. De leidende vergelijking is:
$$\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{op } \partial B, \end{cases}$$
waarbij $u$ de genormaliseerde temperatuur voorstelt, $\lambda$ de Frank-Kamenetskii-parameter is en $\alpha$ de injectiesnelheid is. De functies $\phi$ en $\psi$ beschrijven respectievelijk de stroomsnelheid en de reactieve componentprofielen, terwijl $f$ een positieve, toenemende en convexe reactiesnelheid is die voldoet aan de integrabiliteitsvoorwaarde $\int_0^\infty ds/f(s) < \infty$.

Vorig werk door [GMN] stelde de existentie vast van een kritische parameter $\lambda^*(\alpha)$ en een extremale oplossing $u^*_\alpha$ (de limiet van minimale oplossingen voor $\lambda \nearrow \lambda^*$). Hoewel [GMN] bewees dat $u^*_\alpha(0) \to \infty$ als $\alpha \to \infty$, konden zij slechts een bovengrens vaststellen voor $u^*_\alpha(0)$ onder aanvullende aannames met betrekking tot de niet-lineariteit $f$. Het centrale doel van dit artikel is om dit gat te dichten door het afleiden van evenwichtige (matching) boven- en ondergrenzen voor $u^*_\alpha(0)$ in het asymptotische limiet van sterke stromingen ($\alpha \to \infty$) zonder extra voorwaarden op $f$ op te leggen.

Methodologie
De auteurs maken gebruik van de radiale symmetrie van de extremale oplossing om de partiële differentiaalvergelijking te reduceren tot een gewone differentiaalvergelijking. Door een gewichtsfunctie $\mu(r) = \exp(\alpha \int_0^r s\phi(s)ds)$ te introduceren en een verandering van variabelen $z = \alpha r^2$ uit te voeren, wordt het probleem getransformeerd naar een integraalvergelijking voor de geschaalde oplossing $\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,$$
waarbij $F = \frac{\lambda^*}{4\alpha}f$.

De kern van de analyse berust op de semi-stabiliteitsconditie die door de extremale oplossing wordt vervuld. In de geschaalde coördinaten neemt deze conditie de vorm aan van een gewogen Poincaré-type ongelijkheid:
$$\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,$$
voor alle $\eta \in H^1_0(B)$.

De bewijsstrategie bestaat uit twee componenten:

1. Bovengrens (Theorem 1.1): De auteurs construeren een specifieke testfunctie $\bar{\eta}_z$ (stuksgewijs constante/logaritmische) om de semi-stabiliteitsongelijkheid te benutten. Dit stelt hen in staat om de integraal van de afgeleide van de reactieterm, $\int F'(\bar{u})$, te begrenzen door een logaritmische term. Door de integraalvorm van de ODE te combineren met de convexiteit van $f$, leiden zij een ongelijkheid af die $f'(u(0))$ relateert aan $\log \alpha$.
2. Ondergrens (Theorem 1.2): De auteurs gebruiken een tegensprekingsargument. Zij nemen aan dat de ondergrens faalt, hetgeen impliceert dat er een perturbatie $v > 0$ bestaat zodanig dat $\bar{u} + v$ een oplossing blijft van de vergelijking. Door de geliniariseerde operator te analyseren en gebruik te maken van een kernel-schatting (Lemma 4.2) afgeleid van de eigenschappen van $\mu$ en $\psi$, tonen zij aan dat voor voldoende grote $\alpha$ een dergelijke positieve perturbatie niet kan bestaan zonder de extremaliteit van $\lambda^*$ te schenden. Specifiek demonstreren zij dat indien de grens onjuist zou zijn, men een oplossing zou kunnen construeren die overeenkomt met een parameter $\lambda > \lambda^*$, wat in strijd is met de definitie van de kritische parameter.

Belangrijkste bijdragen en resultaten
Het artikel stelt de volgende scherpe asymptotische grenzen vast voor de extremale oplossing in het nulpunt, $u^*_\alpha(0)$, als $\alpha \to \infty$:

• Theorem 1.1 (Bovengrens): Er bestaan constanten $\alpha_0, c, C > 0$ zodat voor alle $\alpha \geq \alpha_0$,
  $$u^*_\alpha(0) \leq c A,$$
  waarbij $A$ de unieke oplossing is van $f'(A) = C \log \alpha$.

• Theorem 1.2 (Ondergrens): Er bestaan constanten $\tilde{\alpha}_0, \tilde{c}, \tilde{C} > 0$ zodat voor alle $\alpha \geq \tilde{\alpha}_0$,
  $$u^*_\alpha(0) \geq \tilde{c} \tilde{A},$$
  waarbij $\tilde{A}$ de unieke oplossing is van $f'(\tilde{A}) = \tilde{C} \log \alpha$.

Deze resultaten bevestigen dat de groei van de piektemperatuur van de extremale oplossing wordt bepaald door de inverse van de afgeleide van de reactiesnelheid, geëvalueerd op een logaritmische schaal van de stroomsnelheid.

Betekenis
Het artikel claimt een gedeeltelijk openstaande vraag van [GMN] op te lossen met betrekking tot de $L^\infty$-gedraging van extremale oplossingen in turbulente jet-modellen. De primaire significantie ligt in het wegnemen van de "aanvullende aannames op de niet-lineaire reactiesnelheid" die voorheen noodzakelijk waren om de bovengrens te verkrijgen. Door te bewijzen dat de evenwichtige boven- en ondergrenzen gelden voor de algemene klasse van convexe reactiesnelheden die voldoen aan de standaard integrabiliteitsvoorwaarde, bieden de auteurs een volledige karakterisering van het blow-up gedrag van de extremale oplossing in het limiet van sterke stroming. Dit resultaat verstevigt het theoretische begrip van thermische explosielimieten in reactieve turbulente jets binnen het Gelfand-kader.