On the Fujita Phenomenon for a Forced Spatio-Temporal Fractional Diffusion Equation

Resumen Técnico: Sobre el Fenómeno de Fujita para una Ecuación de Difusión Fraccionaria Espacio-Temporal Forzada

Planteamiento del Problema
Este artículo investiga el problema de Cauchy para una ecuación de difusión espacio-temporal fraccionaria semilineal sujeta a un término de forzamiento externo dependiente del tiempo. La ecuación viene dada por:
$$\partial_t^\alpha u + (-\Delta)^s u = |u|^p + t^\sigma w(x), \quad (t, x) \in (0, \infty) \times \mathbb{R}^N,$$
con la condición inicial $u(0, x) = u_0(x) \geq 0$. Aquí, $\partial_t^\alpha$ denota la derivada de Caputo de orden fraccionario $\alpha \in (0, 1)$, $(-\Delta)^s$ es el Laplaciano fraccionario de orden $s \in (0, 1)$, $p > 1$ es el exponente de la no linealidad y $\sigma > -\alpha$ caracteriza el crecimiento o decaimiento temporal del término de forzamiento $w(x)$. El objetivo principal es caracterizar el "fenómeno de Fujita" —específicamente, identificar el exponente crítico $p_F$ que separa el régimen de explosión en tiempo finito del régimen de existencia global—.

Metodología
Los autores emplean una combinación de análisis funcional, cálculo fraccionario y argumentos de escalamiento para analizar el problema:

1. Marco de Solución Débil: El estudio utiliza la fórmula de representación de Duhamel, expresando la solución $u$ como la suma de tres componentes: la evolución lineal de los datos iniciales ($Lu_0$), el término de reacción no lineal ($Nu$) y el término forzado ($Sw$). Esto se basa en los núcleos de la solución fundamental $Z(t, x)$ e $Y(t, x)$ asociados al operador lineal $\partial_t^\alpha + (-\Delta)^s$.
2. Estimaciones de Núcleos: Se lleva a cabo un análisis riguroso de las propiedades de escalamiento, integrabilidad y regularidad de los núcleos $Z$ e $Y$. Estas estimaciones, derivadas de la estructura autosimilar de las soluciones fundamentales, son esenciales para establecer mapeos de contracción en espacios de Lebesgue apropiados.
3. Argumentos de Punto Fijo: Los resultados de existencia local y global se prueban utilizando el teorema del punto fijo de Banach en espacios de funciones adecuados (específicamente $L^\infty(0, T; L^q(\mathbb{R}^N))$ y espacios ponderados con decaimiento temporal).
4. Método de la Función de Prueba (No Existencia): Para probar los resultados de explosión, los autores adaptan el método de la función de prueba. Construyen una función de prueba $\psi_R(t, x)$ cuidadosamente reescalada y aplican fórmulas de integración por partes fraccionaria (que involucran la derivada de Riemann-Liouville de lado derecho) para derivar desigualdades integrales. Al analizar el comportamiento asintótico de estas desigualdades cuando el parámetro de escalamiento $R \to \infty$, derivan contradicciones para exponentes subcríticos.
5. Criterios de Promedio Local: Para la existencia global, el artículo introduce un criterio basado en los promedios locales de los datos iniciales y del término de forzamiento sobre bolas, en lugar de depender únicamente de la pequeñez de la norma global.

Contribuciones Clave y Resultados

• Bien Poseimiento Local: El artículo establece la existencia y unicidad local en el tiempo de soluciones débiles en $L^q(\mathbb{R}^N)$ para datos iniciales $u_0$ y forzamiento $w$ en $L^q$, siempre que $q$ se encuentre dentro de un rango específico relativo al exponente crítico.
• Explosión en Tiempo Finito (Régimen Subcrítico): Bajo el supuesto de que el término de forzamiento tiene masa total positiva ($\int_{\mathbb{R}^N} w(x) dx > 0$), los autores demuestran que no existen soluciones globales débiles si el exponente $p$ cae por debajo de un umbral crítico específico.
  • Para $-\alpha < \sigma \leq 0$ y $2s < N$, el exponente crítico se identifica como:
    $$p^* = \frac{N\alpha - 2s\sigma}{N\alpha - 2s(\alpha + \sigma)}.$$
    Si $1 < p < p^*$, las soluciones explotan en tiempo finito.
  • Si $\sigma > 0$, la explosión ocurre para todo $p > 1$.
• Existencia Global (Régimen Supercrítico):
  • Resultado de Datos Pequeños: Para $p \geq p^*$, se demuestra la existencia de soluciones globales si los datos iniciales $u_0$ y el término de forzamiento son suficientemente pequeños en sus respectivas normas críticas ($L^{p_c}$ y $L^{r_c}$).
  • Criterio de Existencia Robusto: Se establece un resultado de existencia global más general para el rango supercrítico. Este resultado relaja el requisito de pequeñez global, requiriendo en su lugar que los promedios locales de $u_0$ y $w$ sobre bolas grandes satisfagan condiciones específicas de decaimiento y pequeñez. Esto permite datos que pueden ser grandes en normas globales pero que decaen apropiadamente hacia el infinito.
• Identificación del Exponente Crítico: El artículo deriva explícitamente el umbral agudo de tipo Fujita $p^*$ para la ecuación de difusión totalmente fraccionaria con forzamiento creciente en el tiempo. Los autores señalan que este umbral difiere del caso no forzado ($p_F = 1 + \frac{2s\alpha}{N\alpha + 2s(1-\alpha)}$), demostrando que el forzamiento externo acelera el mecanismo de explosión al reducir efectivamente el exponente crítico.

Significancia y Reivindicaciones
El artículo afirma proporcionar el primer umbral agudo de tipo Fujita para ecuaciones de difusión fraccionaria espacio-temporales completas con forzamiento externo creciente en el tiempo. La significancia radica en:

1. Generalización: Extender el clásico fenómeno de Fujita (originalmente para la ecuación del calor) a modelos que combinan efectos de memoria (tiempo fraccionario) e interacciones espaciales de largo alcance (Laplaciano fraccionario) con fuentes externas.
2. Impacto del Forzamiento: Cuantificar cómo la ley de potencia temporal del término de forzamiento ($t^\sigma$) interactúa con los parámetros de difusión fraccionaria ($\alpha, s$) para alterar el exponente crítico. Los resultados indican que los términos de forzamiento pueden contrarrestar los efectos regularizadores de la difusión fraccionaria, expandiendo el rango de exponentes que conducen a la explosión.
3. Criterios de Existencia Refinados: Ir más allá de los supuestos estándar de pequeñez en normas globales para proporcionar criterios de existencia basados en la distribución espacial y los promedios locales de los datos, ofreciendo una comprensión más matizada de la solvabilidad global.

Los autores concluyen que sus hallazgos analíticos están respaldados por simulaciones numéricas en una dimensión, las cuales confirman visualmente la transición predicha entre la explosión y el comportamiento acotado en el exponente crítico derivado.