Self-similar solutions of the three-dimensional Muskat problem with surface tension

Resumen Técnico: Soluciones Autosimilares del Problema de Muskat Tridimensional con Tensión Superficial

Formulación del Problema
El artículo aborda el problema de Muskat de una fase en tres dimensiones para un fluido incompresible en un medio poroso con una frontera superior libre sujeta a tensión superficial. La interfaz $\Sigma_t$ se modela como la gráfica de una función $\eta(t, x)$ sobre $\mathbb{R}^2$. Las ecuaciones gobernantes combinan la ley de Darcy, la incompresibilidad y una condición de contorno dinámica donde el salto de presión se equilibra con la tensión superficial $\sigma \kappa(\eta)$. Este sistema se reduce a una ecuación parabólica cuasilineal de tercer orden para la altura de la interfaz:
$$\partial_t \eta = -\frac{\sigma}{\mu} G(\eta)\kappa(\eta),$$
donde $G(\eta)$ es el operador de Dirichlet–Neumann y $\kappa(\eta)$ es la curvatura media.

Los autores buscan soluciones autosimilares que emanan de datos iniciales cónicos no acotados, $\eta(0, x) = \varepsilon|x|$. Tales soluciones toman la forma $\eta_\varepsilon(t, x) = t^{1/3}U_\varepsilon(t^{-1/3}x)$, donde el perfil $U_\varepsilon$ satisface una ecuación estacionaria que involucra un término de dilatación no invariante por traslación $x \cdot \nabla$. El desafío central es la construcción de estas soluciones para $\varepsilon$ pequeño cuando los datos iniciales son no acotados y poseen una singularidad en el origen, situándose fuera del alcance de las teorías de Cauchy de datos pequeños estándar en espacios de Sobolev acotados o críticos.

Metodología y Estrategia
La prueba emplea un enfoque perturbativo, descomponiendo el perfil $U_\varepsilon$ en un término lineal principal y una corrección:
$$U_\varepsilon = \varepsilon u_* + v_\varepsilon,$$
donde $u_*$ es el "cono lineal regularizado" definido como $u_* = e^{-\frac{\sigma}{\mu}|D|^3}|x|$. La corrección $v_\varepsilon$ se busca en un espacio de funciones específico $X_r = \dot{H}^r(\mathbb{R}^2) \cap \dot{H}^{r+3}(\mathbb{R}^2)$ con $3/2 < r < 2$.

La metodología se apoya en tres pilares técnicos primordiales:

1. Análisis Lineal del Operador de Similitud:
   Se analiza el operador lineal $L = \frac{3\sigma}{\mu}|D|^3 - x \cdot \nabla + 1$ mediante su transformada de Fourier a lo largo de rayos. Los autores construyen un inverso por rayos $L^{-1}$ utilizando una estimación de Hardy ponderada. Este inverso gana tres derivadas en altas frecuencias mientras preserva la estructura cónica de bajas frecuencias, permitiendo la inversión de la ecuación del perfil módulo constantes.

2. Cancelación Cuadrática Favorable:
   Una dificultad crítica al manejar la no linealidad $G(\eta)\kappa(\eta)$ con datos no acotados es el comportamiento del término cuadrático. Los autores identifican una cancelación algebraica específica en el símbolo de la forma bilineal $B(f, g) = \frac{1}{2}(D_\eta G(0)[f]\Delta g + D_\eta G(0)[g]\Delta f)$. En la región de interacción de frecuencias "alta-alta hacia baja" (donde las frecuencias de entrada son grandes pero la frecuencia de salida es pequeña), el símbolo satisface una estimación $|b(\alpha, \beta)| \lesssim |\alpha+\beta||\beta|^3$. Esta ganancia adicional en la frecuencia de salida asegura la convergencia de la sumatoria homogénea de bajas frecuencias, la cual de otro modo divergiría para datos cónicos.

3. Estimaciones de Dirichlet–Neumann Adaptadas al Cono de Orden Finito:
   Las estimaciones estándar de Dirichlet–Neumann fallan para gráficas no acotadas porque las normas de Sobolev altas no controlan la energía de baja frecuencia. Los autores introducen normas de semi-espacio mixtas ($X_r$ e $Y_r$) que imponen pesos de Sobolev altos solo en las frecuencias horizontales positivas. Establecen estimaciones tímidas (tame estimates) para el operador de Dirichlet–Neumann y sus derivadas de forma, demostrando que el operador es $C^3$ en un entorno del perfil cónico. Esto permite el control de los términos residuales cúbicos en la expansión no lineal.

Resultados Clave
El artículo establece el siguiente teorema principal (Teorema 1.3):

• Existencia: Para $\varepsilon$ suficientemente pequeño, existe una solución autosimilar $\eta_\varepsilon(t, x) = t^{1/3}U_\varepsilon(t^{-1/3}x)$ para el problema de Muskat con tensión superficial.
• Estructura del Perfil: El perfil $U_\varepsilon$ es una perturbación del cono regularizado: $U_\varepsilon = \varepsilon u_* + v_\varepsilon$, donde la corrección $v_\varepsilon$ satisface $\|v_\varepsilon\|_{X_r} \leq C \varepsilon^2$.
• Expansión de Segundo Orden: Los autores identifican la corrección cuadrática principal $v^{(2)} = \frac{3\sigma}{\mu} L^{-1} B(u_*, u_*)$. La solución admite la expansión $U_\varepsilon = \varepsilon u_* + \varepsilon^2 v^{(2)} + O_{X_r}(|\varepsilon|^3)$ cuando $\varepsilon \to 0$.
• Traza Inicial y Regularización: La solución alcanza la traza inicial $\eta_\varepsilon(0, x) = \varepsilon|x|$ localmente uniforme. A pesar de la singularidad cónica en $t=0$, la solución es instantáneamente suavizada para cualquier $t > 0$; el perfil $U_\varepsilon$ es $C^1$, y la singularidad se "redondea" en una escala espacial de $O(t^{1/3})$.

Significancia y Reivindicaciones
El artículo afirma proporcionar la primera construcción de soluciones autosimilares para el problema de Muskat tridimensional con tensión superficial emanando de datos iniciales cónicos exactos. La significancia radica en superar las limitaciones de la teoría de Cauchy estándar, que típicamente requiere datos iniciales acotados o de Sobolev críticos. Al separar el componente cónico no acotado y utilizar un espacio de funciones especializado adaptado a la geometría del problema, los autores demuestran que el cono singular es instantáneamente regularizado por la tensión superficial.

El trabajo extiende resultados previos sobre problemas de Muskat impulsados por gravedad (donde el operador principal es de primer orden) al caso de la tensión superficial (operador de tercer orden). Los autores señalan que, si bien la cancelación algebraica utilizada es independiente de la dimensión, la elección específica de las normas de Sobolev y el rango del parámetro $r$ dependen de la dimensión espacial. El artículo concluye que la singularidad cónica no es una barrera para la buena adecuación (well-posedness) en el régimen autosimilar, siempre que el cono inicial sea suficientemente pequeño.