Global weak solutions to a compressible Navier--Stokes/Cahn--Hilliard system with singular entropy of mixing

Resumo Técnico: Soluções Globais Fracas para um Sistema Navier–Stokes Compressível/Cahn–Hilliard com Entropia de Mistura Singular

Enunciado do Problema
Este trabalho aborda a análise matemática de um sistema Navier–Stokes/Cahn–Hilliard compressível que modela a evolução de uma mistura binária de fluidos viscosos submetidos à separação de fases. O sistema é definido em um domínio limitado $\Omega \subset \mathbb{R}^3$ sobre um intervalo de tempo $(0, T)$. As equações governantes descrevem a conservação de massa, momento e a evolução da variável de fase (diferença de concentração de massa) acoplada a um potencial químico.

A principal novidade e o desafio deste estudo residem na escolha da densidade de energia livre. Diferente da literatura anterior sobre misturas compressíveis, que tipicamente emprega potenciais regulares com crescimento polinomial, este trabalho incorpora o potencial de Flory-Huggins (logarítmico), que é fisicamente relevante, e que é singular nos limites do intervalo físico $[-1, 1]$, impondo a restrição de que a variável de fase $c$ deve permanecer em $(-1, 1)$. O sistema está sujeito a condições de contorno de não-deslizamento para a velocidade e condições de Neumann homogêneas para a fase e o potencial químico.

Metodologia
Os autores utilizam um método de compacidade baseado em um esquema de regularização-aproximação para provar a existência de soluções fracas globais no tempo para dados iniciais arbitrariamente grandes. A metodologia central procede da seguinte forma:

1. Esquema de Aproximação: O potencial singular de Flory-Huggins $F(c)$ é aproximado por uma sequência de potenciais regulares $F_\varepsilon(c)$ definidos em toda a reta real, exibindo crescimento de no máximo ordem quadrática no infinito. Isso permite a aplicação de resultados de existência existentes para sistemas regularizados (especificamente aqueles estabelecidos por Abels e Feireisl para misturas compressíveis com potenciais regulares).
2. Estimativas Uniformes: Um componente crítico da prova é a derivação de limites uniformes para as soluções aproximadas $(\varrho_\varepsilon, u_\varepsilon, c_\varepsilon)$ que são independentes do parâmetro de regularização $\varepsilon$.
   • Estimativas de Energia: Desigualdades de energia padrão fornecem limites para a densidade, velocidade e o gradiente da variável de fase.
   • Novas Estimativas para Cahn-Hilliard Dependente da Densidade: Os autores derivam novas estimativas para o potencial químico $\mu_\varepsilon$ e o termo $\varrho_\varepsilon F'_\varepsilon(c_\varepsilon)$ sob suposições mínimas sobre a densidade (não-negativa, $\gamma$-integrável com $\gamma > 3/2$). Ao explorar o controle de $\nabla_x \mu_\varepsilon$ a partir do balanço de energia e utilizando uma função de teste específica envolvendo a diferença de massa total, eles estabelecem limites uniformes em $L^2$ para $\int_\Omega \varrho_\varepsilon |F'_\varepsilon(c_\varepsilon)| dx$.
   • Estimativas de Pressão: Para lidar com a compressibilidade, os autores adaptam a teoria de Lions-Feireisl. Eles derivam estimativas de integrabilidade melhoradas para a pressão $p(\varrho_\varepsilon)$ testando a equação do momento com uma função de teste específica baseada no operador de Bogovskii. Isso resulta em limites $L^q$ para a pressão, onde $q$ depende do expoente adiabático $\gamma$.
3. Compacidade e Passagem ao Limite: Usando os limites uniformes, os autores extraem subsequências convergentes. A convergência forte da densidade $\varrho_\varepsilon$ é estabelecida via equação de continuidade renormalizada e o argumento do fluxo viscoso efetivo. A convergência forte da variável de fase $c_\varepsilon$ é obtida adaptando argumentos de compacidade de estudos incompressíveis anteriores, aproveitando a convergência forte da densidade e os limites uniformes do potencial químico.

Principais Contribuições e Resultados
O resultado principal, enunciado no Teorema 1.1, estabelece a existência de soluções fracas globais no tempo para o sistema para $\gamma > 3/2$ em domínios limitados tridimensionais. A solução satisfaz as seguintes propriedades:

• Regularidade: A solução $(\varrho, u, c, \mu)$ possui classes de regularidade específicas, incluindo $\varrho \in C_{weak}([0, T]; L^\gamma(\Omega))$ e $c \in L^\infty(0, T; W^{1,2}(\Omega)) \cap L^2(0, T; W^{2, \frac{2\gamma}{\gamma+1}}(\Omega))$.
• Restrições Físicas: Uma conquista significativa é a prova de que a variável de fase $c$ satisfaz os limites físicos quase em toda parte. Especificamente:
  • $-1 \leq c \leq 1$ quase em toda parte em $(0, T) \times \Omega$.
  • $-1 < c < 1$ quase em toda parte no conjunto onde a densidade é positiva ($\{\varrho > 0\}$).
• Dados Iniciais Gerais: O resultado vale para dados iniciais com energia total finita, sem exigir que a densidade seja estritamente positiva ou limitada longe de zero, nem exigir que a concentração inicial esteja estritamente dentro de $(-1, 1)$ em todos os pontos (apenas que a massa relativa $M_c/M \in (-1, 1)$).

Significância e Alegações
O artigo afirma fornecer o primeiro resultado de existência para soluções fracas globais para o sistema Navier–Stokes/Cahn–Hilliard com o potencial de Flory-Huggins fisicamente relevante. Os autores destacam que trabalhos anteriores sobre misturas compressíveis foram limitados a potenciais regulares com crescimento polinomial, que falham em impor a restrição física $c \in [-1, 1]$.

A significância do trabalho reside em:

1. Preencher a Lacuna: Ele aborda a lacuna analítica entre modelos de campo de fase incompressíveis e compressíveis ao lidar com a entropia singular no regime compressível.
2. Suposições Mínimas: A derivação de estimativas para a equação de Cahn–Hilliard dependente da densidade não depende da positividade estrita ou do limite inferior da densidade, uma condição frequentemente exigida em estudos anteriores (por exemplo, no caso incompressível pelo segundo autor e Temam).
3. Consistência Física: A prova demonstra rigorosamente que a solução limite respeita os limites físicos da variável de fase, garantindo a consistência do modelo com a termodinâmica subjacente da separação de fases.

Os autores observam que sua abordagem depende de um método de compacidade e que a prova da convergência forte das densidades segue técnicas estabelecidas (teoria de Lions-Feireisl) adaptadas a este sistema acoplado específico. Eles também mencionam que os resultados podem ser estendidos para domínios bidimensionais com regularidade melhorada e que potenciais singulares mais gerais que satisfaçam condições de crescimento semelhantes podem ser tratados.