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

Resumen Técnico: Cotas $L^\infty$ Agudas para Soluciones Extremales en Problemas de Tipo Gelfand en Jets Turbulentos

Planteamiento del Problema
El artículo investiga un problema de valor en la frontera que modela la explosión térmica en jets turbulentos reactivos, derivado originalmente de [GHH]. El problema se formula en el disco unitario $B \subset \mathbb{R}^2$ y pertenece a la clase de los problemas de tipo Gelfand. La ecuación gobernante es:
$$\begin{cases} -\Delta u - \alpha r \phi(r) \frac{\partial u}{\partial r} = \lambda \psi(r) f(u) & \text{en } B, \\ u > 0 & \text{en } B, \\ u = 0 & \text{en } \partial B, \end{cases}$$
donde $u$ representa la temperatura normalizada, $\lambda$ es el parámetro de Frank-Kamenetskii y $\alpha$ es la velocidad de inyección. Las funciones $\phi$ y $\psi$ describen los perfiles de velocidad de flujo y de componente reactivo, respectivamente, mientras que $f$ es una tasa de reacción positiva, creciente y convexa que satisface la condición de integrabilidad $\int_0^\infty ds/f(s) < \infty$.

Un trabajo previo por [GMN] estableció la existencia de un parámetro crítico $\lambda^*(\alpha)$ y una solución extremal $u^*_\alpha$ (el límite de las soluciones mínimas cuando $\lambda \nearrow \lambda^*$). Aunque [GMN] demostró que $u^*_\alpha(0) \to \infty$ cuando $\alpha \to \infty$, solo pudieron establecer una cota superior para $u^*_\alpha(0)$ bajo supuestos adicionales respecto a la no linealidad $f$. El objetivo central de este artículo es cerrar esta brecha derivando cotas superior e inferior coincidentes para $u^*_\alpha(0)$ en el límite asintótico de flujos fuertes ($\alpha \to \infty$) sin imponer condiciones adicionales sobre $f$.

Metodología
Los autores explotan la simetría radial de la solución extremal para reducir la ecuación en derivadas parciales a una ecuación diferencial ordinaria. Al introducir la función de peso $\mu(r) = \exp(\alpha \int_0^r s\phi(s)ds)$ y realizar un cambio de variables $z = \alpha r^2$, el problema se transforma en una ecuación integral para la solución reescalada $\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,$$
donde $F = \frac{\lambda^*}{4\alpha}f$.

El núcleo del análisis se basa en la condición de semi-estabilidad satisfecha por la solución extremal. En las coordenadas reescaladas, esta condición toma la forma de una desigualdad de tipo Poincaré ponderada:
$$\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,$$
para toda $\eta \in H^1_0(B)$.

La estrategia de demostración consta de dos componentes principales:

1. Cota Superior (Teorema 1.1): Los autores construyen una función de prueba específica $\bar{\eta}_z$ (constante por partes/logarítmica) para utilizar la desigualdad de semi-estabilidad. Esto les permite acotar la integral de la derivada del término de reacción, $\int F'(\bar{u})$, mediante un término logarítmico. Combinando esto con la forma integral de la EDO y la convexidad de $f$, derivan una desigualdad que relaciona $f'(u(0))$ con $\log \alpha$.
2. Cota Inferior (Teorema 1.2): Los autores emplean un argumento de contradicción. Suponen que la cota inferior falla, lo que implica la existencia de una perturbación $v > 0$ tal que $\bar{u} + v$ sigue siendo una solución a la ecuación. Al analizar el operador linealizado y utilizar una estimación de núcleo (Lema 4.2) derivada de las propiedades de $\mu$ y $\psi$, demuestran que para $\alpha$ suficientemente grande, tal perturbación positiva no puede existir sin violar la extremalidad de $\lambda^*$. Específicamente, demuestran que si la cota fuera falsa, se podría construir una solución correspondiente a un parámetro $\lambda > \lambda^*$, lo cual contradice la definición del parámetro crítico.

Contribuciones Clave y Resultados
El artículo establece las siguientes cotas asintóticas agudas para la solución extremal en el origen, $u^*_\alpha(0)$, cuando $\alpha \to \infty$:

• Teorema 1.1 (Cota Superior): Existen constantes $\alpha_0, c, C > 0$ tales que para todo $\alpha \geq \alpha_0$,
  $$u^*_\alpha(0) \leq c A,$$
  donde $A$ es la solución única de $f'(A) = C \log \alpha$.

• Teorema 1.2 (Cota Inferior): Existen constantes $\tilde{\alpha}_0, \tilde{c}, \tilde{C} > 0$ tales que para todo $\alpha \geq \tilde{\alpha}_0$,
  $$u^*_\alpha(0) \geq \tilde{c} \tilde{A},$$
  donde $\tilde{A}$ es la solución única de $f'(\tilde{A}) = \tilde{C} \log \alpha$.

Estos resultados confirman que el crecimiento del pico de temperatura de la solución extremal está gobernado por la inversa de la derivada de la tasa de reacción evaluada en una escala logarítmica de la velocidad de flujo.

Significado
El artículo afirma resolver una pregunta parcialmente abierta dejada por [GMN] con respecto al comportamiento $L^\infty$ de las soluciones extremales en modelos de jets turbulentos. La significancia primaria radica en la eliminación de los "supuestos adicionales sobre la tasa de reacción no lineal" que eran necesarios previamente para obtener la cota superior. Al demostrar que las cotas superior e inferior coincidentes se mantienen para la clase general de tasas de reacción convexas que satisfacen la condición de integrabilidad estándar, los autores proporcionan una caracterización completa del comportamiento de explosión (blow-up) de la solución extremal en el límite de flujo fuerte. Este resultado consolida la comprensión teórica de los límites de explosión térmica en jets turbulentos reactivos dentro del marco de Gelfand.