Nonlinear parabolic problem with time fractional derivative

Resumen Técnico: Problema Parabólico No Lineal con Derivada Fraccionaria y Potencial Singlar Doble

1. Planteamiento del Problema

El artículo estudia un problema de evolución parabólico no lineal que combina tres características complejas:

1. Derivada fraccionaria en el tiempo: Se utiliza la derivada de Liouville-Caputo de orden $\alpha \in (0, 1)$, denotada como $D^\alpha_t$.
2. Operador $p$-Laplaciano: La parte espacial está gobernada por el operador $\Delta_p u = \text{div}(|\nabla u|^{p-2}\nabla u)$, con $p > 2$.
3. Potencial Singlar Doble: El sistema incluye un potencial $W(x)$ que presenta singularidades tanto en el origen ($x=0$) como en la frontera del dominio $\partial\Omega$.

La ecuación principal (1.1) se formula en un dominio acotado estrellado $\Omega \subset \mathbb{R}^n$ ($n > p$) con condiciones de frontera de Dirichlet homogéneas:

$$\begin{cases} D^\alpha_t u - \Delta_p u - \mu W(x)|u|^{p-2}u = 0, & \text{en } Q = (0, T) \times \Omega, \\ u(t, x) = 0, & \text{en } [0, T] \times \partial\Omega, \\ u(0, x) = u_0(x), & \text{en } \Omega, \end{cases}$$

donde el potencial singular $W(x)$ se define como:
$$W(x) = |x|^{-p} \left[ 1 - \left( \frac{|x|}{\phi(x)} \right)^{\frac{n-p}{p}} \right]^{-p}$$
con $\phi(x)$ siendo una función homogénea de orden 0 que describe la geometría del dominio.

El objetivo central es determinar las condiciones bajo las cuales las soluciones débiles existen globalmente en el tiempo o experimentan una explosión en tiempo finito (blow-up), dependiendo de la constante de Hardy óptima $\Lambda(n, p) = \left(\frac{n-p}{p}\right)^p$.

2. Metodología

Los autores emplean una combinación de análisis funcional, teoría de ecuaciones diferenciales fraccionarias y desigualdades de Hardy. La estrategia se divide en cuatro etapas principales:

• Aproximación por Truncamiento: Dado que el potencial $W(x)$ es singular, se introduce un problema truncado (1.4) donde el potencial se reemplaza por $W_N(x) = \min\{N, W(x)\}$. Se demuestra la existencia de soluciones para este problema regularizado.
• Principio de Comparación para EDOs Fraccionarias: En la Sección 2, se establece un principio de comparación (Teorema 2.1) para ecuaciones diferenciales ordinarias fraccionarias (FODE). Esto es crucial para analizar el comportamiento asintótico de soluciones escalares que actúan como cotas.
• Principio de Comparación para EDPs: En la Sección 3.1, se extiende el principio de comparación al problema parabólico fraccionario (Teorema 3.1). La prueba utiliza una "identidad fundamental" (Lema 3.1, basada en Zacher) que relaciona la derivada fraccionaria de una función compuesta con la derivada de la función misma, permitiendo manejar el término no lineal $|u|^{p-2}u$.
• Estimaciones A Priori: Se derivan estimaciones de energía (Sección 3.2) para las soluciones truncadas $u_N$, utilizando la desigualdad de Hardy y propiedades de convexidad del operador $p$-Laplaciano.

3. Contribuciones Clave y Resultados Principales

El trabajo establece dos teoremas fundamentales que delinean el comportamiento de las soluciones débiles según la relación entre el parámetro $\mu$ y la constante de Hardy $\Lambda(n, p)$:

Teorema 1.1 (Existencia Global):

• Condición: Si $\mu < \Lambda(n, p)$ (el coeficiente del potencial es menor que la constante óptima de Hardy).
• Resultado: Para cualquier dato inicial $u_0 \in W^{1,p}_0(\Omega)$, la solución débil existe globalmente en el tiempo y no explota. Se prueban estimaciones de acotación uniformes en $L^2$ y $W^{1,p}$ para todo $t > 0$.
• Mecanismo: La disipación del operador $p$-Laplaciano domina el término de potencial singular, garantizando la estabilidad.

Teorema 1.2 (Explosión en Tiempo Finito):

• Condición: Si $\mu > \Lambda(n, p)$ y el dato inicial $u_0$ es positivo y acotado.
• Resultado: No existe solución global. La solución explota en un tiempo finito $T_m < \infty$.
• Mecanismo: El potencial singular es lo suficientemente fuerte para superar la disipación. Se construye una subsolución de la forma $v(t, x) = T(t)X(x)$, donde $X(x)$ es una autofunción positiva y $T(t)$ es una solución de una EDO fraccionaria que explota (Teorema 2.3). Por el principio de comparación, la solución original $u$ debe explotar al menos tan rápido como $T(t)$.

Análisis de EDOs Fraccionarias (Sección 2):
Un resultado independiente y significativo es la demostración rigurosa de que la ecuación $D^\alpha_t w = w^p$ con $p>1$ y $\alpha \in (0,1)$ tiene soluciones que explotan en tiempo finito, mejorando resultados previos de Shan y Lv [35] al proporcionar fórmulas explícitas para el tiempo de explosión $T_m$ en diferentes regímenes de parámetros.

4. Significado e Impacto

• Extensión de Teoría Clásica: El artículo generaliza resultados conocidos para ecuaciones parabólicas clásicas (con derivada temporal entera) al contexto de derivadas fraccionarias, un área de creciente importancia en física matemática (difusión anómala, materiales con memoria).
• Tratamiento de Singularidades Dobles: Aborda la complejidad de tener singularidades simultáneas en el origen y en la frontera, un escenario menos tratado que las singularidades puntuales internas.
• Umbral Óptimo: Identifica la constante de Hardy $\Lambda(n, p)$ como el umbral crítico exacto que separa la existencia global de la explosión, incluso en presencia de la memoria fraccionaria (derivada de orden $\alpha$).
• Aplicaciones: Los resultados son relevantes para modelar procesos físicos donde la difusión es anómala (subdifusión) y donde existen campos de fuerza intensos o defectos estructurales (singularidades), como en física molecular, cosmología cuántica y mecánica de fluidos en medios porosos.

En conclusión, el papel proporciona un marco analítico robusto para entender la dinámica de soluciones en sistemas no lineales con memoria y singularidades, estableciendo límites precisos para la estabilidad y la explosión de soluciones.