Self-similar solutions to the time-fractional Porous-Medium Equation

1. Planteamiento del Problema

El artículo estudia la existencia y propiedades de soluciones auto-similares no negativas para la Ecuación de Medio Poroso (PME) con derivada temporal fraccionaria (en el sentido de Caputo). La ecuación se define como:

$$\partial^\alpha_t u = \Delta(u^m)$$

donde:

• $t > 0, x \in \mathbb{R}^d$ ($d \ge 1$).
• $\alpha \in (0, 1)$ es el orden de la derivada fraccionaria de Caputo: $\partial^\alpha_t u(t, x) = \int_0^t \frac{(t-s)^{-\alpha}}{\Gamma(1-\alpha)} \partial_s u(s, x) ds$.
• $m > 0$ es el exponente de no linealidad.
• Las soluciones parten de una condición inicial de tipo delta de Dirac con masa finita $M > 0$: $u(0, \cdot) = M\delta_0$.
• Se busca soluciones que conserven la masa total $M$ para todo $t > 0$.

El objetivo principal es caracterizar los perfiles auto-similares $U(z)$, donde la solución tiene la forma $u(t, x) = t^{-a} U(t^{-b}|x|)$, para todo $d \ge 1$ y en el rango de exponentes $m > m_c = (d-2)_+/d$. Este rango es óptimo para la existencia de soluciones con masa finita.

2. Metodología

Los autores emplean una combinación de análisis asintótico, teoría de ecuaciones integrales y métodos de sub-soluciones y super-soluciones.

1. Reducción a Ecuación Integral:
   Mediante un análisis de escalado y sustitución de la forma auto-similar en la EDP, los autores derivan una ecuación integral para el perfil $U(z)$:
   $$U(z)^m = \int_z^{+\infty} K(z, \rho) U(\rho) d\rho$$
   Donde el núcleo $K(z, \rho)$ es no negativo y depende de los parámetros $\alpha, m, d$. Una contribución clave es la simplificación explícita de este núcleo, incluso en dimensiones $d=1$, expresándolo mediante funciones especiales (Beta incompleta) o integrales de potencias.

2. Método de Sub y Super-soluciones:
   Para demostrar la existencia, construyen pares de sub-soluciones y super-soluciones adecuadas que acotan el perfil buscado. Utilizan iteraciones de tipo Picard (o monotonía) para converger hacia la solución única del operador integral.

3. Análisis de Singularidades y Comportamiento Asintótico:
   Se estudia detalladamente el comportamiento del perfil cerca del origen ($z \to 0$) y en el infinito ($z \to \infty$). Se utilizan estimaciones precisas del núcleo $K$ y sus derivadas para determinar la regularidad y la tasa de decaimiento.

4. Límites Continuos:
   Se analizan rigurosamente los límites cuando los parámetros críticos varían:

   • $m \to 1$ (transición entre difusión lenta y rápida).
   • $\alpha \to 1$ (recuperación del caso clásico).
   • $m \to \infty$ (límite de "mesa").

3. Contribuciones Clave y Resultados Principales

El artículo establece resultados completos para tres regímenes distintos de difusión, generalizando resultados conocidos del caso clásico ($\alpha=1$) al caso fraccionario temporal.

A. Existencia y Unicidad

• Existencia: Se prueba la existencia de un perfil canónico $U_{\alpha, m}$ para todo $m > m_c$.
• Unicidad: Se demuestra la unicidad del perfil para $m \ge 1$. Para $m < 1$, se establece la unicidad de la familia de soluciones de masa finita que forman un "abanico asintótico" alrededor de la solución muy singular.
• Solución Débil: Se caracteriza la solución auto-similar como una "solución muy débil" de la EDP, satisfaciendo una formulación integral que incluye un término de memoria que recuerda la masa inicial (delta de Dirac) a través de la derivada fraccionaria.

B. Comportamiento según el Régimen de Difusión

1. Difusión Lenta ($m > 1$):

   • Soporte Compacto: Las soluciones tienen soporte compacto $[0, R]$.
   • Comportamiento en la Frontera Libre: Cerca del borde del soporte, el perfil se comporta como una potencia simple: $U(z) \sim (1-z)^{\frac{2-\alpha}{m-1}}$. Se proporciona la constante exacta de este comportamiento.
   • Límite de Mesa ($m \to \infty$): Las soluciones convergen a una función característica (forma de mesa) $\chi_{B_1}$, idéntica al caso clásico $\alpha=1$, perdiendo la singularidad en el origen.

2. Difusión Rápida Subcrítica ($m_c < m < 1$):

   • Soporte No Compacto: Las soluciones son positivas en todo $\mathbb{R}^d$.
   • Cola Pesada (Heavy Tails): A diferencia del caso lineal ($m=1$) que tiene colas exponenciales, aquí las soluciones tienen colas de tipo potencia: $U(z) \sim z^{-\frac{2}{1-m}}$.
   • Solución Muy Singular (VSS): Existe una solución auto-similar exacta de masa infinita $U^*(z) = c^* z^{-\frac{2}{1-m}}$. Las soluciones de masa finita forman un abanico asintótico alrededor de esta solución singular.
   • Se deriva una expansión asintótica de segundo orden para las colas.

3. Caso Lineal ($m = 1$):

   • Se revisa la solución conocida mediante la transformada de Fourier, que involucra funciones de Mittag-Leffler.
   • Se demuestra que el comportamiento de las soluciones fraccionarias converge al caso clásico cuando $\alpha \to 1$.

C. Comportamiento en el Origen ($z \to 0$)

Se establece un comportamiento preciso para $U^m$ cerca del origen, que depende de la dimensión $d$:

• Para $d \ge 2$, $U^m$ se comporta como el potencial newtoniano ($z^{2-d}$ o $\log z$), lo que implica que $\Delta(u^m)$ es una múltiplo de la delta de Dirac. Esto confirma que la derivada fraccionaria "recuerda" la condición inicial singular.
• Para $d = 1$, el comportamiento es acotado pero con una constante específica dependiente de $\alpha$ y $m$.

D. Continuidad en $m$

Aunque el perfil canónico $U_{\alpha, m}$ no es continuo en $m=1$ debido a la elección de normalización, se demuestra que las soluciones con masa fija $M$ son continuas respecto a $m$. Es decir, $\lim_{m \to 1} U_{\alpha, m, M}(z) = U_{\alpha, 1, M}(z)$.

4. Significado e Impacto

• Generalización Óptima: El trabajo extiende la teoría de soluciones auto-similares de la PME clásica al caso de tiempo fraccionario para todo el rango óptimo de exponentes $m > m_c$ y dimensiones $d \ge 1$.
• Nuevas Características Físicas: Destaca cómo la derivada fraccionaria temporal introduce efectos de memoria que modifican la regularidad en el origen y la estructura de las colas, manteniendo sin embargo la estructura de soporte compacto para $m>1$.
• Herramientas Analíticas: La derivación de la ecuación integral con un núcleo simplificado y su análisis mediante métodos de sub/super-soluciones proporciona un marco robusto para estudiar problemas no locales en tiempo.
• Validación Numérica: Se presenta un esquema numérico basado en iteraciones de Picard para el caso $m > 1$, validando teóricamente el comportamiento de las soluciones frente a variaciones de $\alpha$ y $m$.
• Problemas Abiertos: El artículo identifica claramente problemas futuros, como la unicidad para $m < 1$, el comportamiento de datos iniciales generales (atractores globales) y la extensión a la difusión muy singular ($m < m_c$) y problemas totalmente no locales.

En resumen, este artículo proporciona una teoría completa y rigurosa de las soluciones fundamentales (tipo Barenblatt) para la ecuación de medio poroso con tiempo fraccionario, cerrando brechas importantes en la comprensión de la difusión anómala en medios porosos.