Stochastic Compressible Euler Equations with Frictional Damping: Existence of $L^\infty$ Martingale Solutions and Asymptotic Porous Medium-Like Behavior

1. Definición del Problema

El estudio se centra en las ecuaciones de Euler compresibles isentrópicas unidimensionales con amortiguamiento friccional lineal, perturbadas por un ruido multiplicativo blanco en el tiempo. El sistema describe la evolución de la densidad $\rho$ y el momento $m = \rho u$ (donde $u$ es la velocidad del fluido) en un intervalo acotado $[0, 1]$.

Las ecuaciones gobernantes son:
$$\begin{cases} d\rho + \partial_x m \, dt = 0 \\ dm + \partial_x \left(\frac{m^2}{\rho} + p(\rho)\right) dt = -\alpha m \, dt + \sigma(x, \rho, m) \, dW(t) \end{cases}$$
donde:

• $\alpha > 0$ es el coeficiente de amortiguamiento friccional.
• $p(\rho) = \kappa \rho^\gamma$ es la ley de presión para un gas politrópico perfecto ($\gamma > 1$).
• $W(t)$ es un movimiento browniano unidimensional.
• $\sigma(x, \rho, m)$ es el coeficiente de ruido, que depende no linealmente del estado y satisface condiciones de soporte compacto y Lipschitz.

Condiciones de contorno e iniciales:

• Se imponen condiciones de contorno de Dirichlet para el momento: $m(t, 0) = m(t, 1) = 0$.
• Los datos iniciales $(\rho_0, m_0)$ pertenecen a $L^\infty(0, 1)$, con $\rho_0 \geq 0$ y masa total finita.

El objetivo principal es doble:

1. Establecer la existencia global en el tiempo de soluciones martingala en $L^\infty$.
2. Analizar el comportamiento asintótico (largo tiempo) de estas soluciones, demostrando su convergencia a un estado estacionario y su aproximación por la ecuación de medio poroso.

---

2. Metodología

Los autores emplean un enfoque analítico sofisticado que combina teoría de ecuaciones diferenciales estocásticas (SDEs), análisis funcional y teoría de compacidad compensada.

A. Esquema de Aproximación de Dos Niveles

Para probar la existencia, se utiliza un esquema de aproximación jerárquica:

1. Regularización Parabólica ($\epsilon$-nivel): Se introduce una viscosidad artificial $\epsilon \partial_{xx} U$ para transformar el sistema hiperbólico en uno parabólico estocástico. Esto permite obtener soluciones suaves y aplicar fórmulas de Itô.
2. Discretización Temporal ($\tau$-nivel): Se utiliza un esquema de división de operadores (splitting scheme) de tipo Lie-Trotter. En cada intervalo de tiempo $\tau$, el problema se descompone en:
   • Un subproblema determinístico (ecuaciones de Euler viscosas con amortiguamiento).
   • Un subproblema estocástico (donde la densidad se mantiene constante y solo el momento evoluciona bajo el ruido).
   • Se interpola linealmente entre las soluciones de estos subproblemas para construir una solución continua aproximada.

B. Técnicas de Compacidad y Límites

• Estimaciones Uniformes: Se derivan cotas $L^\infty$ uniformes en los parámetros de aproximación ($\epsilon, \tau$) utilizando invariantes de Riemann y regiones invariantes. Esto es crucial para manejar datos iniciales en $L^\infty$ y evitar la formación de vacíos (vacuum) no controlados.
• Teorema de Representación de Skorokhod: Se utiliza para pasar al límite $\tau \to 0$ y luego $\epsilon \to 0$, obteniendo convergencia casi segura en un espacio de probabilidad modificado.
• Compacidad Compensada y Medidas de Young: Para manejar la no linealidad en el límite $\epsilon \to 0$ (especialmente en los términos de flujo $\frac{m^2}{\rho}$), se emplea la teoría de medidas de Young y compacidad compensada estocástica (extendiendo trabajos de Feng-Nualart y Berthelin-Vovelle).

C. Análisis de Comportamiento Asintótico

Para estudiar el largo tiempo ($t \to \infty$), los autores:

1. Derivan estimaciones de momentos agudas para la entropía en el nivel de aproximación $\epsilon$.
2. Utilizan una transformación de variables escaladas exponencialmente para desacoplar la tasa de decaimiento.
3. Aplican el Lema de Gronwall en el contexto estocástico, combinado con resultados recientes de Yuskovych [58] sobre la convergencia casi segura de integrales estocásticas bajo ciertas condiciones de momentos.
4. Demuestran que el ruido, bajo condiciones de Lipschitz acotadas, no impide la disipación de energía causada por el amortiguamiento.

---

3. Contribuciones Clave

1. Primera prueba de convergencia casi segura (pathwise): A diferencia de la literatura previa que se limita a la existencia de medidas invariantes o convergencia estadística, este trabajo proporciona el primer resultado riguroso de convergencia casi segura de las soluciones a un estado estacionario constante para las ecuaciones de Euler estocásticas con amortiguamiento.
2. Soluciones en $L^\infty$ con datos iniciales $L^\infty$: Se establece la existencia de soluciones martingala en el espacio $L^\infty((0, \infty) \times (0, 1))$, lo cual es más fuerte que las soluciones en espacios de energía ($L^2$) o medidas de probabilidad consideradas en trabajos anteriores (como Berthelin y Vovelle).
3. Justificación Rigurosa de la Ley de Darcy Estocástica: Se demuestra que, asintóticamente, el momento del sistema estocástico converge a la ley de Darcy ($m \approx -\frac{1}{\alpha} \partial_x p(\rho)$) y la densidad converge al perfil de la ecuación de medio poroso determinista.
4. Nuevas Técnicas de Estimación de Momentos: Se desarrollan métodos novedosos para cuantificar y probar la decadencia de los efectos estocásticos mediante estimaciones de momentos superiores de la entropía, permitiendo controlar las integrales estocásticas en el límite.

---

4. Resultados Principales

Teorema 1.1: Existencia Global

Bajo suposiciones de datos iniciales acotados y un coeficiente de ruido con soporte compacto y Lipschitz continuo, existe una solución martingala $L^\infty$ débil de entropía global en el tiempo para el sistema (1.1). Esta solución satisface una desigualdad de entropía apropiada y las condiciones de contorno de Dirichlet en un sentido distribucional.

Teorema 1.2: Decaimiento Exponencial y Convergencia

Si la constante de Lipschitz del ruido ($A_0$) es suficientemente pequeña en relación con el coeficiente de amortiguamiento ($\alpha$), entonces:

• La solución $(\rho(t), m(t))$ converge casi seguramente y exponencialmente rápido a un estado estacionario constante $(\rho^*, 0)$, donde $\rho^*$ es la masa total inicial.
• La tasa de decaimiento es de la forma $C_0 e^{-rt}$ en la norma $L^2$.

Teorema 1.3: Aproximación por Ecuación de Medio Poroso

El comportamiento asintótico de la solución estocástica $(\rho_{IE}, m_{IE})$ está bien aproximado por la solución clásica de la ecuación de medio poroso determinista (PME):
$$\partial_t \rho = \partial_{xx} p(\rho), \quad m = -\frac{1}{\alpha} \partial_x p(\rho)$$
La diferencia entre la solución estocástica y la solución de la PME decae exponencialmente en el tiempo, validando que el ruido no altera la estructura asintótica fundamental del sistema, siempre que sea dominado por el amortiguamiento.

---

5. Significado e Impacto

Este trabajo es fundamental por varias razones:

• Puente entre Estocástico y Determinista: Proporciona una justificación rigurosa de que, a largo plazo, el comportamiento de un sistema de fluidos compresibles estocásticos con amortiguamiento está gobernado por la dinámica determinista de la ecuación de medio poroso. Esto es contraintuitivo, ya que el ruido blanco suele inducir comportamientos caóticos o difusivos diferentes; aquí, el amortiguamiento lineal domina la estadística a largo plazo.
• Avance en la Teoría de Ecuaciones Estocásticas: Supera las limitaciones de los resultados anteriores que solo garantizaban la existencia de medidas invariantes. La demostración de convergencia pathwise (casi segura) es un hito técnico significativo en el análisis de EDPs estocásticas hiperbólicas.
• Aplicabilidad Física: El modelo tiene relevancia directa en la dinámica de fluidos en medios porosos y, específicamente, en la modelización del flujo sanguíneo en arterias elásticas (donde las fluctuaciones aleatorias de los latidos cardíacos pueden modelarse como ruido). La demostración de que el sistema converge a un estado de equilibrio estable, incluso bajo perturbaciones aleatorias, es crucial para la comprensión de la hemodinámica a largo plazo.
• Metodología Robusta: El enfoque de dos niveles (viscosidad artificial + división de operadores) junto con las estimaciones de momentos de alta orden ofrece una plantilla metodológica que podría aplicarse a otros sistemas de conservación estocásticos con disipación.

En resumen, el artículo establece un marco teórico sólido para entender cómo el amortiguamiento y el ruido interactúan en sistemas de fluidos compresibles, demostrando que, bajo condiciones adecuadas, el sistema recupera un comportamiento determinista predecible en el límite de tiempo infinito.