Blow-up of solutions to the Euler-Poisson-Darbox equation with critical power nonlinearity

Aquí presento un resumen técnico detallado del artículo "Blow-up of solutions to the Euler-Poisson-Darboux equation with critical power nonlinearity" (Explosión de soluciones para la ecuación de Euler-Poisson-Darboux con no linealidad de potencia crítica), escrito por Mengting Fan, Ning-An Lai y Hiroyuki Takamura.

---

1. Planteamiento del Problema

El artículo aborda el problema de Cauchy singular para la ecuación semilineal de Euler-Poisson-Darboux (EPD) en $\mathbb{R}^n$:

$$\begin{cases} u_{tt} - \Delta u + \frac{\mu}{t} u_t = t^\alpha |u|^p, & t > 0, \, x \in \mathbb{R}^n, \\ u(0, x) = u_0(x), \\ u_t(0, x) = 0, \end{cases}$$

donde:

• $\mu > 0$ es un parámetro que describe el efecto de la amortiguación (damping).
• $\alpha \ge 0$ y $p > 1$ son exponentes no lineales.
• El término $\frac{\mu}{t} u_t$ representa una amortiguación singular en $t=0$, lo que introduce dificultades analíticas significativas en comparación con problemas regulares donde el tiempo inicial es $t_0 > 0$.

Objetivo Principal:
Establecer resultados de explosión en tiempo finito (blow-up) y obtener una estimación superior para la vida útil (lifespan) de las soluciones de energía cuando el exponente no lineal $p$ es crítico.

El caso crítico se define como $p = p_S(n + \mu, \alpha)$, donde $p_S$ es la raíz positiva de la ecuación cuadrática:
$$\gamma(n, \mu, \alpha, p) := 2 + (n + \mu + 1 + 2\alpha)p - (n + \mu - 1)p^2 = 0.$$

Este caso es particularmente delicado porque representa la frontera entre la existencia global y la explosión. En el caso crítico, la vida útil de la solución se espera que sea de tipo exponencial ($T(\epsilon) \sim \exp(C\epsilon^{-p(p-1)})$), lo cual es mucho más difícil de probar que en el caso subcrítico (donde la vida útil es polinomial).

2. Metodología

Los autores emplean una combinación de técnicas de análisis funcional, teoría de ecuaciones diferenciales ordinarias y métodos iterativos mejorados. Los componentes clave de su metodología son:

A. Construcción de una Función de Prueba Mejorada

A diferencia de trabajos anteriores que usaban funciones de prueba estándar, los autores construyen una función de prueba $b_q(t, x)$ mediante una integración sobre un parámetro de escala $\lambda$:
$$b_q(t, x) = \int_0^1 h(\lambda t) \phi(\lambda x) \lambda^{q-1} d\lambda,$$
donde:

• $h(t)$ se construye utilizando la función de Bessel modificada $K_{\frac{\mu-1}{2}}(t)$. Esta elección es crucial para manejar la singularidad en $t=0$.
• $\phi(x)$ es una función relacionada con la transformada de Fourier de una esfera unitaria.
• El parámetro $q$ se elige específicamente para que $b_q$ satisfaga la ecuación homogénea asociada a la parte lineal del operador EPD.

Ventaja de esta construcción: Permite obtener un comportamiento asintótico óptimo ($b_q \sim t^{-q}$) y eliminar la singularidad en el origen al integrar por partes, algo que no es posible con funciones de prueba convencionales debido al término $\frac{\mu}{t} u_t$.

B. Funcional de Explosión y Acotación Inferior

Definen un funcional $Y(M)$ que incluye la integral espacio-temporal del término no lineal:
$$Y(M) = \int_1^M \left( \int_{t/2}^t \int_{\mathbb{R}^n} \tau^\alpha |u|^p \eta^2 \dots \right) t^{-1} dt.$$
El núcleo de la prueba consiste en demostrar que este funcional crece logarítmicamente:
$$Y(t) \ge C \epsilon^p \ln t.$$
Para lograr esto, los autores logran una cota inferior mejorada para la integral del término no lineal, obteniendo un factor de crecimiento logarítmico adicional. Esto es esencial para el caso crítico, ya que en el caso subcrítico el crecimiento es polinomial y más fácil de controlar.

C. Manejo de la Singularidad en $t=0$

Un desafío técnico mayor es la integración por partes en la definición de solución débil, donde el término $\frac{\mu}{t} u_t$ diverge. Los autores utilizan las propiedades de recurrencia y el comportamiento asintótico de las funciones de Bessel modificadas $K_\nu$ para demostrar que el límite $\lim_{t \to 0} (-\partial_t b_q + \frac{\mu}{t} b_q)$ existe y es positivo, permitiendo que la integración por partes sea válida y no genere términos divergentes.

3. Contribuciones Clave

1. Resolución del Caso Crítico para EPD Singular: Extienden los resultados de explosión conocidos para el caso subcrítico ($p < p_S$) al caso crítico ($p = p_S$) para la ecuación EPD con amortiguación singular.
2. Mejora de la Función de Prueba: Introducen una función de prueba basada en una integral de escala de funciones de Bessel modificadas, lo que permite un control preciso del comportamiento asintótico y la eliminación de singularidades en $t=0$.
3. Respuesta Parcial a un Problema Abierto: El trabajo proporciona una respuesta parcial al Problema Abierto 1 planteado por D'Abbicco (2021), quien señaló la falta de conocimiento completo sobre los resultados de explosión para la ecuación EPD semilineal.
4. Estimación de Vida Útil Óptima: Demuestran que para datos iniciales pequeños con soporte compacto, la vida útil $T(\epsilon)$ satisface:
   $$T(\epsilon) \le \exp\left( C \epsilon^{-p(p-1)} \right).$$
   Esta estimación confirma que el comportamiento en el caso crítico es de tipo exponencial, consistente con la teoría de la ecuación de onda semilineal clásica.

4. Resultados Principales

El Teorema 1.2 establece que bajo las siguientes condiciones:

• $\alpha \ge 0$, $n \ge 3$ (o $n=2$ con $\alpha > 0$).
• $0 < \mu \le \mu^*(n, \alpha)$, donde $\mu^*$ es un umbral específico derivado de la comparación entre los exponentes de Fujita y Strauss.
• $p = p_S(n + \mu, \alpha)$.
• Datos iniciales $u_0 \in H^1(\mathbb{R}^n)$ positivos y con soporte compacto.

Entonces, existe un $\epsilon_0 > 0$ tal que si la amplitud de los datos iniciales es pequeña ($\epsilon$), la solución de energía explota en tiempo finito y su vida útil está acotada superiormente por una función exponencial inversa de $\epsilon$.

5. Significado e Impacto

• Avance Teórico: Este trabajo cierra una brecha importante en la teoría de ecuaciones de onda con amortiguación dependiente del tiempo. Muestra que, a pesar de la singularidad en el coeficiente de amortiguación, el comportamiento crítico sigue siendo gobernado por el exponente de Strauss desplazado.
• Herramientas Nuevas: La técnica de construir funciones de prueba mediante integración sobre un parámetro de escala utilizando funciones de Bessel modificadas es un aporte metodológico valioso que puede aplicarse a otros problemas de ondas con coeficientes singulares o no regulares.
• Validación de Conjeturas: El resultado valida la conjetura de que $p_S(n+\mu, \alpha)$ actúa como el exponente crítico real para la existencia global frente a la explosión en este rango de parámetros, alineándose con la intuición física de que la amortiguación débil no cambia la naturaleza crítica del problema, solo desplaza el umbral.

En resumen, el artículo demuestra rigurosamente que las soluciones a la ecuación de Euler-Poisson-Darboux con amortiguación singular y no linealidad crítica no pueden existir globalmente para datos iniciales pequeños, estableciendo una cota superior precisa para su tiempo de vida mediante un análisis refinado de funciones de prueba y estimaciones logarítmicas.