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

Resumo Técnico: Limites $L^\infty$ Agudos para Soluções Extremais em Problemas do Tipo Gelfand em Jatos Turbulentos

Enunciado do Problema
O artigo investiga um problema de valor de contorno que modela a explosão térmica em jatos turbulentos reativos, originalmente derivado em [GHH]. O problema é formulado no disco unitário $B \subset \mathbb{R}^2$ e pertence à classe dos problemas do tipo Gelfand. A equação governante é:
$$\begin{cases} -\Delta u - \alpha r \phi(r) \frac{\partial u}{\partial r} = \lambda \psi(r) f(u) & \text{em } B, \\ u > 0 & \text{em } B, \\ u = 0 & \text{em } \partial B, \end{cases}$$
onde $u$ representa a temperatura normalizada, $\lambda$ é o parâmetro de Frank-Kamenetskii e $\alpha$ é a velocidade de injeção. As funções $\phi$ e $\psi$ descrevem os perfis de velocidade de fluxo e de componente reativo, respectivamente, enquanto $f$ é uma taxa de reação positiva, crescente e convexa que satisfaz a condição de integrabilidade $\int_0^\infty ds/f(s) < \infty$.

O trabalho anterior por [GMN] estabeleceu a existência de um parâmetro crítico $\lambda^*(\alpha)$ e uma solução extremal $u^*_\alpha$ (o limite das soluções mínimas conforme $\lambda \nearrow \lambda^*$). Embora [GMN] tenha provado que $u^*_\alpha(0) \to \infty$ quando $\alpha \to \infty$, eles puderam apenas estabelecer um limite superior para $u^*_\alpha(0)$ sob suposições adicionais relativas à não linearidade $f$. O objetivo central deste artigo é fechar essa lacuna ao derivar limites superior e inferior correspondentes para $u^*_\alpha(0)$ no limite assintótico de fluxos fortes ($\alpha \to \infty$) sem impor condições extras sobre $f$.

Metodologia
Os autores exploram a simetria radial da solução extremal para reduzir a equação diferencial parcial a uma equação diferencial ordinária. Ao introduzir a função de peso $\mu(r) = \exp(\alpha \int_0^r s\phi(s)ds)$ e realizar uma mudança de variáveis $z = \alpha r^2$, o problema é transformado em uma equação integral para a solução reescalonada $\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,$$
onde $F = \frac{\lambda^*}{4\alpha}f$.

O núcleo da análise baseia-se na condição de semi-estabilidade satisfeita pela solução extremal. Nas coordenadas reescalonadas, esta condição assume a forma de uma desigualdade do 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 todo $\eta \in H^1_0(B)$.

A estratégia de prova envolve dois componentes principais:

1. Limite Superior (Teorema 1.1): Os autores constroem uma função de teste específica $\bar{\eta}_z$ (constante por partes/logarítmica) para utilizar a desigualdade de semi-estabilidade. Isso permite que eles limitem a integral da derivada do termo de reação, $\int F'(\bar{u})$, por um termo logarítmico. Combinando isso com a forma integral da EDO e a convexidade de $f$, eles derivam uma desigualdade relacionando $f'(u(0))$ com $\log \alpha$.
2. Limite Inferior (Teorema 1.2): Os autores empregam um argumento de contradição. Eles assumem que o limite inferior falha, o que implica a existência de uma perturbação $v > 0$ tal que $\bar{u} + v$ permanece uma solução para a equação. Ao analisar o operador linearizado e utilizar uma estimativa de núcleo (Lema 4.2) derivada das propriedades de $\mu$ e $\psi$, eles mostram que, para $\alpha$ suficientemente grande, tal perturbação positiva não pode existir sem violar a extremalidade de $\lambda^*$. Especificamente, eles demonstram que, se o limite fosse falso, seria possível construir uma solução correspondente a um parâmetro $\lambda > \lambda^*$, o que contradiz a definição do parâmetro crítico.

Contribuições Principais e Resultados
O artigo estabelece os seguintes limites assintóticos agudos para a solução extremal na origem, $u^*_\alpha(0)$, conforme $\alpha \to \infty$:

• Teorema 1.1 (Limite Superior): Existem constantes $\alpha_0, c, C > 0$ tais que para todo $\alpha \geq \alpha_0$,
  $$u^*_\alpha(0) \leq c A,$$
  onde $A$ é a solução única de $f'(A) = C \log \alpha$.

• Teorema 1.2 (Limite Inferior): Existem constantes $\tilde{\alpha}_0, \tilde{c}, \tilde{C} > 0$ tais que para todo $\alpha \geq \tilde{\alpha}_0$,
  $$u^*_\alpha(0) \geq \tilde{c} \tilde{A},$$
  onde $\tilde{A}$ é a solução única de $f'(\tilde{A}) = \tilde{C} \log \alpha$.

Estes resultados confirmam que o crescimento do pico de temperatura da solução extremal é governado pela inversa da derivada da taxa de reação avaliada em uma escala logarítmica da velocidade de fluxo.

Significância
O artigo afirma resolver uma questão parcialmente aberta deixada por [GMN] a respeito do comportamento $L^\infty$ de soluções extremais em modelos de jatos turbulentos. A principal significância reside na remoção das "suposições adicionais sobre a taxa de reação não linear" que eram anteriormente necessárias para obter o limite superior. Ao provar que os limites superior e inferior correspondentes se mantêm para a classe geral de taxas de reação convexas que satisfazem a condição de integrabilidade padrão, os autores fornecem uma caracterização completa do comportamento de blow-up da solução extremal no limite de fluxo forte. Este resultado consolida a compreensão teórica dos limites de explosão térmica em jatos turbulentos reativos dentro do framework de Gelfand.