Boundary regularity for parabolic systems with nonstandard $(p,q)$-growth conditions in smooth convex domains

Resumen Técnico de "Regularidad en la Frontera para Sistemas Parabólicos con Condiciones de Crecimiento No Estándar (p, q) en Dominios Convexos Suaves"

Planteamiento del Problema
Este manuscrito investiga la regularidad en la frontera de las soluciones débiles locales de sistemas parabólicos no lineales de la forma:
$$\partial_t u_i - \text{div}(a(|Du|)Du_i) = f_i, \quad i = 1, \dots, N,$$
definidos en un cilindro espacio-temporales $\Omega_T = \Omega \times (0, T)$. El dominio $\Omega \subset \mathbb{R}^n$ ($n \ge 2$) es un conjunto $C^2$ acotado y convexo. Los coeficientes $a: \mathbb{R}_{>0} \to \mathbb{R}_{>0}$ son de tipo Uhlenbeck y satisfacen condiciones de crecimiento no estándar $(p, q)$ con $2 \le p \le q < p + \frac{4}{n+2}$. La inhomogeneidad $f$ pertenece a $L^{n+2+\sigma}(\Omega_T, \mathbb{R}^N)$ para algún $\sigma > 0$.

El objetivo principal es establecer una estimación de gradiente local $L^\infty$ (regularidad Lipschitz) para soluciones débiles $u$ que se anulan en la frontera lateral $(\partial \Omega)_T$, extendiendo los resultados previos de interior hacia la frontera en el contexto del crecimiento no estándar.

Metodología
La estrategia de la prueba sigue un enfoque de tres pasos, adaptando técnicas de la teoría de la regularidad elíptica y parabólica estándar al entorno de crecimiento no estándar:

1. Estimaciones A Priori para Soluciones Suaves:
   Los autores derivan primero una estimación cuantitativa del gradiente $L^\infty$ para soluciones clásicas suaves ($u \in C^3$) con coeficientes y datos suaves. Esto implica:

   • Establecer una desigualdad de tipo Caccioppoli que involucra segundas derivadas, utilizando la convexidad de $\Omega$ y las condiciones de Dirichlet homogéneas para descartar términos de frontera desfavorables mediante una identidad de geometría diferencial (Grisvard).
   • Derivar una desigualdad de tipo Hölder inversa. Este paso depende de la relación específica de los parámetros $q < p + \frac{4}{n+2}$, la cual asegura que la inmersión de Sobolev en el entorno parabólico (mejorando la integrabilidad de $p$ a $p + \frac{4}{n}$) sea suficiente para cerrar la estimación.
   • Ejecutar un procedimiento de iteración de Moser para obtener una cota $L^\infty$. Un componente crítico es el manejo del "déficit de escala" inherente a los sistemas de crecimiento $(p, q)$, lo que impide que la estimación sea homogénea.

2. Procedimiento de Aproximación:
   Dado que el resultado principal se aplica a soluciones débiles con datos y coeficientes no suaves, se construye un esquema de aproximación:

   • Se suavizan los coeficientes $a$ y el dato $f$ para obtener aproximaciones suaves $a_\varepsilon$ y $f_\varepsilon$.
   • Se añade un término de tipo $q$-Laplaciano no degenerado a los coeficientes para definir $b_\varepsilon$, asegurando condiciones de crecimiento $q$ estándar que garantizan la existencia de soluciones débiles únicas y regulares ($C^3$) $u_\varepsilon$.
   • Se demuestra una estimación de comparación entre las soluciones aproximantes $u_\varepsilon$ y la solución débil original $u$, demostando la convergencia fuerte de los gradientes $Du_\varepsilon \to Du$ en $L^p$.

3. Paso al Límite:
   La estimación cuantitativa del gradiente $L^\infty$ derivada para las aproximaciones suaves se aplica a $u_\varepsilon$. Al pasar al límite $\varepsilon \to 0$ y utilizar la semicontinuidad inferior de la norma $L^\infty$ y la convergencia de los datos, la estimación se transfiere a la solución débil original $u$.

Contribuciones Clave y Resultados
El resultado principal es el Teorema 1.3, el cual establece que para cualquier solución débil local $u$ que se anula en la frontera lateral, el gradiente $Du$ es esencialmente acotado en un entorno de la frontera. Específicamente, para cualquier cilindro parabólico $Q_\rho(z_0)$ centrado en un punto de la frontera $x_0 \in \partial \Omega$:
$$Du \in L^\infty(\Omega_T \cap Q_\rho(z_0), \mathbb{R}^{Nn}).$$
El artículo proporciona una estimación cuantitativa precisa (Ecuación 1.6):
$$\sup_{\Omega_T \cap Q_\rho(z_0)} |Du| \le C \left( \left( 1 + \|f\|_{L^{n+2+\sigma}}^{(n+2)(n+2+\sigma)/\sigma} \rho^{n+2} \right) \fint_{\Omega_T \cap Q_{2\rho}(z_0)} (1 + |Du|^p) \, dxdt \right)^{\frac{2p}{4+(p-q)(n+2)}}.$$

Las características técnicas clave de este resultado incluyen:

• Rango de Parámetros: El resultado se mantiene para el rango $2 \le p \le q < p + \frac{4}{n+2}$. Este rango coincide con los resultados recientes de interior para sistemas de crecimiento $(p, q)$ no estándar, pero está adaptado para el entorno de frontera.
• Déficit de Escala: El exponente $\frac{2p}{4+(p-q)(n+2)}$ representa el déficit de escala característico de los sistemas parabólicos con crecimiento no estándar. El artículo señala que este exponente tiende al infinito cuando $q$ se acerca al límite superior $p + \frac{4}{n+2}$, lo que resalta la nitidez de la condición.
• Regularidad del Dominio: La prueba requiere que $\Omega$ sea un dominio $C^2$. Esta regularidad es esencial para la aplicación de la identidad de geometría diferencial utilizada para manejar las integrales de frontera en la estimación de Caccioppoli.
• Regularidad de los Datos: Se asume que la inhomogeneidad $f$ pertenece a $L^{n+2+\sigma}$, lo cual es suficiente para la derivación de la estimación sin requerir la diferenciabilidad débil de $f$.

Significado y Reivindicaciones
El artículo posiciona su trabajo como la contraparte de frontera de los resultados existentes de regularidad de interior para sistemas parabólicos de crecimiento $(p, q)$ no estándar (citando [6, 42]) y extiende la teoría de regularidad de frontera conocida para sistemas parabólicos de crecimiento $p$ estándar (citando [5]) al caso no estándar.

Los autores afirman que este trabajo llena un vacío en la literatura respecto a la regularidad global (hasta la frontera) para sistemas parabólicos con crecimiento no estándar, un tema donde los resultados son significativamente más escasos que en el entorno elíptico. La estimación proporcionada es consistente con los resultados de frontera previos en el entorno elíptico $(p, q)$ [7] y el caso parabólico $(p, p)$ [5]. El artículo no pretende resolver el problema para el caso degenerado ($\mu=0$) o para estructuras no de tipo Uhlenbeck generales, señalando que estas últimas típicamente solo permiten la regularidad parcial. El enfoque permanece estrictamente en establecer la cota de Lipschitz local bajo las especificadas suposiciones estructurales y geométricas.