Equilibria of the Pressureless Euler System in Dimension 1

Resumo Técnico: Equilíbrios do Sistema de Euler Sem Pressão em Dimensão 1

Enunciado do Problema
Este artigo investiga o comportamento de longo prazo do sistema de Euler sem pressão unidimensional (PE), originalmente proposto por Zeldovich para modelar a formação de estruturas em grande escala. O sistema é dado pelas equações de continuidade e momento:
$$\partial_t \rho + \partial_y(\rho v) = 0, \quad \partial_t(\rho v) + \partial_y(\rho v^2) = 0$$
onde $\rho$ é uma distribuição de massa (medida de probabilidade com segundo momento finito) e $v$ é o campo de velocidade. Os autores focam especificamente em Soluções de Partículas Colantes (SPS), uma classe de soluções distributivas onde partículas de massa colam-se ao colidirem. Embora a existência de tais soluções tenha sido estabelecida via uma fórmula de projeção (Natile-Savaré), as condições sob as quais essas soluções convergem para um estado de equilíbrio, a natureza desse equilíbrio e o tempo necessário para o colapso permaneceram questões em aberto. Diferente do sistema Euler-Poisson atrativo, onde o equilíbrio é tipicamente uma única massa pontual, o sistema de Euler sem pressão permite uma rica variedade de distribuições de equilíbrio.

Metodologia
A análise baseia-se fortemente na fórmula de projeção derivada em [10], que traduz o problema do "espaço de posição" físico para o "espaço de massa". Ao utilizar o mapa de transporte ótimo $N_t: (0,1) \to \mathbb{R}$ (onde $N_t \# \mathcal{L} = \rho_t$), a evolução do sistema é caracterizada pela projeção da trajetória linear $N_0 + tV_0$ no cone convexo $K$ de funções não decrescentes em $L^2(0,1)$:
$$N_t = \text{proj}_K(N_0 + tV_0)$$
Aqui, $V_0 = v_0 \circ N_0$. Os autores utilizam ferramentas de análise convexa, especificamente propriedades do envelope convexo ($\text{conv}$) e da derivada à direita, para analisar o comportamento assintótico de $N_t$. Eles definem uma função primitiva $F(x) = \int_0^x V_0(w) dw$ e utilizam a invariância galileana do sistema para normalizar o momento total para zero.

Principais Contribuições e Resultados

1. Condição Necessária e Suficiente para o Equilíbrio:
   O artigo estabelece uma condição precisa, denotada como $(\star)$, baseada exclusivamente nos dados iniciais $(N_0, V_0)$, que determina se a SPS converge para um equilíbrio. A condição exige:

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

   Os autores provam que a SPS converge para um equilíbrio na métrica 1-Wasserstein se, e somente se, a condição $(\star)$ for satisfeita. Se $F(x) < 0$ para qualquer $x$, a solução diverge (a massa se espalha ao infinito).

2. Caracterização do Estado de Equilíbrio:
   Quando a condição $(\star)$ é satisfeita, o artigo fornece uma descrição completa da distribuição limite $\rho_\infty$. O equilíbrio é determinado pelos zeros da função $F$:

   • Regiões Estacionárias: Em qualquer intervalo onde $F$ é identicamente zero, a distribuição de massa permanece estacionária para todo o tempo ($N_t(x) = N_0(x)$).
   • Regiões de Colapso: Em qualquer intervalo $(a, b)$ onde $F(a)=F(b)=0$ e $F > 0$ estritamente no interior, a massa dentro desse intervalo colapsa para uma única medida de Dirac (massa pontual) localizada no centro de massa daquele "pedaço" inicial de massa.
     Isso implica que o equilíbrio é uma superposição de massa estacionária e massas pontuais colapsadas.

3. Tempo de Colapso:
   Os autores derivam uma fórmula explícita para o tempo de colapso ($T_c$), definido como o tempo mínimo $T$ tal que $\rho_t = \rho_\infty$ para todo $t \geq T$. Se a solução converge, $T_c$ é dado pelo supremo sobre todos os intervalos $(a,b)$ onde $F$ é positivo:
   $$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}$$
   Se este valor for infinito, a solução aproxima-se do equilíbrio assintoticamente, mas nunca o alcança em tempo finito.

4. Métricas de Convergência:
   Os resultados são formulados principalmente em termos da distância 1-Wasserstein ($W_1$). O artigo observa que, se o suporte da solução permanecer limitado (o que ocorre se o suporte inicial for limitado e $(\star)$ for satisfeito), a convergência em $W_p$ para qualquer $p$ segue-se. Os autores discutem a extensão para a convergência $W_2$, observando que, embora esta ocorra sob condições de monotonicidade específicas perto das fronteiras do intervalo de massa, uma condição geral dependente apenas dos dados iniciais para a convergência $W_2$ permanece uma dificuldade aberta devido à natureza global da operação de envelope convexo.

Significância e Alegações
O artigo afirma fornecer uma caracterização completa do comportamento de equilíbrio para soluções de partículas colantes ao sistema de Euler sem pressão em 1D usando apenas as condições iniciais. Ele resolve as três questões centrais relativas à existência de equilíbrio, a forma específica desse equilíbrio e o tempo de colapso. Os autores enfatizam que sua abordagem evita a necessidade de suposições de suporte compacto ao longo do tempo, que eram exigidas em literatura relacionada (ex: [8]). O trabalho distingue-se do sistema Euler-Poisson repulsivo, onde métodos preditivos como este falham devido ao comportamento patológico, destacando que a natureza impulsionada pelo momento do sistema de Euler sem pressão permite essa determinação precisa. Os resultados dependem da fórmula de projeção e da análise convexa, oferecendo um elo rigoroso entre a geometria da distribuição de velocidade inicial e o estado assintótico do sistema.