Equivalence of entropy solutions and gradient flows for pressureless 1D Euler systems

Resumo Técnico: Equivalência de Soluções de Entropia e Fluxos de Gradiente para Sistemas Euler 1D Sem Pressão

Enunciado do Problema
O artigo investiga soluções de distribuição para o problema de Cauchy para o sistema Euler unidimensional sem pressão:
$$\partial_t\rho + \partial_x(\rho v) = 0, \quad \partial_t(\rho v) + \partial_x(\rho v^2) = f[\rho],$$
onde $\rho$ é uma medida de probabilidade com momentos de segunda ordem limitados e $v$ é um campo de velocidade. O termo de força $f[\rho]$ é gerado pelo próprio fluido, com atenção específica dada ao sistema Euler–Poisson (atrativo ou repulsivo) e sistemas com confinamento quadrático e amortecimento linear. Um desafio central neste campo é a formação de singularidades (concentração de massa ou "blow-up") onde as soluções clássicas falham. O artigo aborda a necessidade de estender as soluções além dessas singularidades e de estabelecer uma continuação única e fisicamente relevante.

Dois frameworks de solução distintos foram desenvolvidos para tais sistemas:

1. Soluções Lagrangianas: Introduzidas por Brenier, Gangbo, Savaré e Westdickenberg [7], são formuladas como fluxos de gradiente $L^2$ no espaço de medidas de probabilidade, utilizando mapas de transporte ótimo e inclusões diferenciais envolvendo cones normais.
2. Soluções de Entropia: Estudadas por Nguyen e Tudorascu [36, 37] para o sistema Euler–Poisson, são definidas via leis de conservação escalar para a função de distribuição cumulativa, baseando-se na condição E de Oleĭnik para admissibilidade.

Embora ambos os frameworks produzam soluções únicas sob suposições específicas, sua equivalência não havia sido rigorosamente estabelecida para o sistema Euler de segunda ordem geral, particularmente na presença de concentração de massa e forças repulsivas onde as trajetórias das partículas podem se cruzar ou se separar.

Metodologia
Os autores estabelecem a equivalência entre esses dois conceitos construindo uma correspondência direta entre as estruturas geométricas que regem cada framework. O núcleo da metodologia envolve:

1. Transformação de Variáveis:

   • Para soluções Lagrangianas, a variável primária é o mapa de transporte ótimo $X(t, m)$, uma função não decrescente que leva a medida de Lebesgue para $\rho(t)$.
   • Para soluções de entropia, a variável primária é a função de distribuição $M(t, x)$, a inversa à direita de $X$.

2. Correspondência Geométrica:
   O artigo demonstra que os mecanismos que selecionam soluções únicas em ambos os frameworks são matematicamente idênticos:

   • Admissibilidade de Choque: A condição de que a velocidade do mapa de transporte $\dot{X}$ é a projeção $L^2$ da velocidade prescrita $U$ sobre o cone tangente do cone de funções não decrescentes ($\dot{X} = P_{TXK}U$) é mostrada como equivalente à condição de Rankine–Hugoniot para choques na função de distribuição $M$.
   • Princípio de Seleção: O requisito de que o resíduo $U - \dot{X}$ pertença ao cone normal $N_{X}K$ (garantindo que a solução permaneça dentro do cone admissível) é provado como equivalente à condição E de Oleĭnik para a lei de conservação escalar que rege $M$. Esta equivalência baseia-se em uma caracterização do cone normal por Natile e Savaré [35].

3. Derivação da Função de Fluxo:
   Os autores derivam uma lei de conservação escalar $\partial_t M + \partial_x U(t, M) = 0$, onde a função de fluxo $U(t, M)$ é a primitiva da velocidade prescrita da formulação Lagrangiana. Eles utilizam resultados de Golovaty e Nguyen [24] para lidar com funções de fluxo que não são continuamente Lipschitz, o que é necessário para o sistema Euler–Poisson.

Principais Contribuições e Resultados

• Teorema de Equivalência: O resultado principal é a prova de que as soluções Lagrangianas e as soluções de entropia são equivalentes para o sistema Euler 1D sem pressão. Especificamente, uma solução Lagrangiana $X$ gera uma solução de entropia $M$ (sua inversa generalizada), e vice-versa. Isso ocorre para termos de força $f[\rho]$ gerais que satisfazem condições específicas de limitação e continuidade, incluindo o sistema Euler–Poisson.
• Extensão Além do Blow-up: A equivalência permite que os autores definam soluções únicas para sistemas onde as soluções clássicas falham.
  • Sistema Euler–Poisson Repulsivo: O artigo estabelece que, para forças repulsivas ($\alpha < 0$), a massa pode se concentrar inicialmente, mas eventualmente se difundirá. O conceito de solução de entropia fornece um framework rigoroso para descrever esse comportamento, confirmando que partículas que colidem devido a velocidades iniciais "supercríticas" irão se separar, um comportamento ditado pela condição E de Oleĭnik.
  • Sistemas com Amortecimento e Confinamento: Para o sistema com amortecimento linear e confinamento quadrático (estudado em [14]), os autores provam a existência de soluções de entropia globalmente definidas. Eles mostram que, independentemente de ocorrer ou não a concentração de massa, o sistema converge assintoticamente para uma distribuição de massa uniforme com velocidade zero.
• Dissipação de Energia e Fisicalidade: Através de exemplos de dinâmica de partículas (Seções 6 e 7), o artigo ilustra que a solução de entropia/Lagrangiana é a continuação única "física". Em casos repulsivos, soluções alternativas "pegajosas" (onde as partículas permanecem presas após a colisão) violam a condição E de Oleĭnik e resultam em aumentos não físicos de energia. A solução de entropia garante que a energia seja não crescente.
• Estabilidade Assintótica: O artigo prova que, para sistemas amortecidos, a solução converge exponencialmente rápido para um estado de equilíbrio na métrica $T_2$. Para sistemas sem amortecimento, um resultado de estabilidade de Lyapunov é estabelecido, mostrando que as soluções permanecem próximas se os dados iniciais forem próximos.

Significância e Alegações
O artigo alega que a equivalência desses conceitos de solução é fundamental para entender o comportamento de longo prazo dos sistemas Euler sem pressão. Ao ligar a abordagem geométrica de fluxo de gradiente $L^2$ com a abordagem analítica de solução de entropia, os autores fornecem um framework unificado que:

1. Resolve ambiguidades na continuação da solução após a formação de singularidades, particularmente em regimes repulsivos onde as partículas podem colidir e se separar.
2. Generaliza resultados anteriores de equações de agregação de primeira ordem (ex: [4]) para sistemas de segunda ordem.
3. Fornece uma base rigorosa para o estudo de sistemas com amortecimento e confinamento além do tempo de vida das soluções clássicas.

Os autores enfatizam que o cenário unidimensional é essencial para esta equivalência devido à relação única entre mapas de transporte ótimo e funções de distribuição (inversas generalizadas) e o ordenamento das partículas. Eles observam que generalizar estes resultados para dimensões superiores é improvável de uma forma ingênua. O trabalho serve para unificar a intuição de "partículas pegajosas" com a teoria rigorosa de leis de conservação, oferecendo uma ferramenta robusta para analisar a complexa dinâmica de partículas em modelos de fluidos.