Equilibria of the Pressureless Euler System in Dimension 1

Resumen Técnico: Equilibrios del Sistema de Euler sin Presión en Dimensión 1

Planteamiento del Problema
Este artículo investiga el comportamiento a largo plazo del sistema de Euler sin presión (PE) en una dimensión, propuesto originalmente por Zeldovich para modelar la formación de estructuras a gran escala. El sistema está dado por las ecuaciones de continuidad y de momento:
$$\partial_t \rho + \partial_y(\rho v) = 0, \quad \partial_t(\rho v) + \partial_y(\rho v^2) = 0$$
donde $\rho$ es una distribución de masa (medida de probabilidad con segundo momento finito) y $v$ es el campo de velocidad. Los autores se centran específicamente en las Soluciones de Partículas Pegajosas (SPS), una clase de soluciones distribucionales donde las partículas de masa se pegan al colisionar. Si bien la existencia de tales soluciones se establece mediante una fórmula de proyección (Natile-Savaré), las condiciones bajo las cuales estas soluciones convergen a un estado de equilibrio, la naturaleza de dicho equilibrio y el tiempo requerido para el colapso han permanecido como preguntas abiertas. A diferencia del sistema Euler-Poisson atractivo, donde el equilibrio es típicamente una única masa puntual, el sistema de Euler sin presión permite una rica variedad de distribuciones de equilibrio.

Metodología
El análisis se apoya fuertemente en la fórmula de proyección derivada en [10], que traslada el problema desde el "espacio de posición" físico al "espacio de masa". Al utilizar el mapa de transporte óptimo $N_t: (0,1) \to \mathbb{R}$ (donde $N_t \# \mathcal{L} = \rho_t$), la evolución del sistema se caracteriza mediante la proyección de la trayectoria lineal $N_0 + tV_0$ sobre el cono convexo $K$ de funciones no decrecientes en $L^2(0,1)$:
$$N_t = \text{proj}_K(N_0 + tV_0)$$
Aquí, $V_0 = v_0 \circ N_0$. Los autores emplean herramientas de análisis convexo, específicamente propiedades de la envolvente convexa ($\text{conv}$) y la derivada por la derecha, para analizar el comportamiento asintótico de $N_t$. Definen una función primitiva $F(x) = \int_0^x V_0(w) dw$ y utilizan la invariancia galileana del sistema para normalizar el momento total a cero.

Contribuciones Clave y Resultados

1. Condición Necesaria y Suficiente para el Equilibrio:
   El artículo establece una condición precisa, denotada como $(\star)$, basada únicamente en los datos iniciales $(N_0, V_0)$, que determina si la SPS converge a un equilibrio. La condición requiere:

   • Momento total cero: $\int_0^1 V_0(w) dw = 0$ (equivalente a $F(1)=0$).
   • Momento parcial no negativo: $\int_0^x V_0(w) dw \geq 0$ para todo $x \in (0,1)$ (equivalente a $F(x) \geq 0$).

   Los autores demuestran que la SPS converge a un equilibrio en la métrica 1-Wasserstein si y solo si se cumple la condición $(\star)$. Si $F(x) < 0$ para algún $x$, la solución diverge (la masa se dispersa hacia el infinito).

2. Caracterización del Estado de Equilibrio:
   Cuando la condición $(\star)$ se cumple, el artículo proporciona una descripción completa de la distribución límite $\rho_\infty$. El equilibrio está determinado por los ceros de la función $F$:

   • Regiones Estacionarias: En cualquier intervalo donde $F$ es identicamente cero, la distribución de masa permanece estacionaria para todo el tiempo ($N_t(x) = N_0(x)$).
   • Regiones de Colapso: En cualquier intervalo $(a, b)$ donde $F(a)=F(b)=0$ y $F > 0$ estrictamente en el interior, la masa dentro de ese intervalo colapsa hacia una medida de Dirac única (masa puntual) ubicada en el centro de masa de ese "trozo" inicial de masa.
     Esto implica que el equilibrio es una superposición de masa estacionaria y masas puntuales colapsadas.

3. Tiempo de Colapso:
   Los autores derivan una fórmula explícita para el tiempo de colapso ($T_c$), definido como el tiempo mínimo $T$ tal que $\rho_t = \rho_\infty$ para todo $t \geq T$. Si la solución converge, $T_c$ está dado por el supremo sobre todos los intervalos $(a,b)$ donde $F$ es positiva:
   $$T_c = \sup_{[a,b] \in I} \sup_{x \in (a,b)} \frac{\int_x^b N_0(w) dw + \frac{x-b}{b-a}\int_a^b N_0(w) dw}{\int_0^x V_0(w) dw}$$
   Si este valor es infinito, la solución se aproxima al equilibrio de forma asintótica pero nunca lo alcanza en un tiempo finito.

4. Métricas de Convergencia:
   Los resultados se enmarcan principalmente en términos de la distancia 1-Wasserstein ($W_1$). El artículo señala que si el soporte de la solución permanece acotado (lo cual ocurre si el soporte inicial es acotado y se cumple $(\star)$), la convergencia en $W_p$ para cualquier $p$ se sigue. Los autores discuten la extensión a la convergencia $W_2$, señalando que, si bien esta se cumple bajo ciertas condiciones de monotonicidad cerca de los límites del intervalo de masa, una condición general que dependa únicamente de los datos iniciales para la convergencia $W_2$ sigue siendo una dificultad abierta debido a la naturaleza global de la operación de la envolvente convexa.

Significancia y Reivindicaciones
El artículo afirma proporcionar una caracterización completa del comportamiento de equilibrio para las soluciones de partículas pegajosas al sistema de Euler sin presión en 1D utilizando únicamente las condiciones iniciales. Resuelve las tres cuestiones centrales relativas a la existencia de equilibrio, la forma específica de dicho equilibrio y el tiempo de colapso. Los autores enfatizan que su enfoque evita la necesidad de supuestos de soporte compacto en el tiempo, que eran requeridos en la literatura relacionada (por ejemplo, [8]). El trabajo se distingue del sistema Euler-Poisson repulsivo, donde los métodos predictivos de este tipo fallan debido a un comportamiento patológico, destacando que la naturaleza impulsada por el momento del sistema de Euler sin presión permite esta determinación precisa. Los resultados dependen de la fórmula de proyección y del análisis convexo, ofreciendo un vínculo riguroso entre la geometría de la distribución de velocidad inicial y el estado asintótico del sistema.