On weak solutions for the stationary Cahn-Hillard-Navier-Stokes equations with singular potential

Resumo Técnico: Soluções Fracas para as Equações de Cahn-Hilliard-Navier-Stokes Estacionárias com Potencial Singular

Enunciado do Problema
Este artigo investiga a existência de soluções fracas para o sistema estacionário de Navier–Stokes–Cahn-Hilliard que descreve uma mistura fluida bifásica, compressível e tridimensional com interface difusa. O sistema é governado pela conservação de massa, momento e concentração, acoplados a uma equação de potencial químico. Uma característica distintiva do modelo é a inclusão de uma densidade de energia livre singular do tipo Flory–Huggins (logarítmica). Esta escolha garante que a fração mássica $c$ permaneça dentro do intervalo físico $(-1, 1)$ e permite a presença de estados de vácuo ($\rho = 0$).

O problema de valor de contorno específico envolve um domínio limitado $\Omega \subset \mathbb{R}^3$ com condições de Dirichlet para a velocidade e condições de Neumann para a concentração e o potencial químico. Os principais desafios matemáticos surgem de duas fontes:

1. O Cenário Estacionário: Diferentemente de problemas evolutivos, o sistema estacionário carece de uma desigualdade de energia a priori ($E(t) \leq E(0)$) impulsionada pela evolução temporal, o que tipicamente fornece o controle necessário sobre o potencial singular e o campo de concentração.
2. Potencial Singular e Vácuo: A singularidade logarítmica na energia livre impede regularizações padrão via expansão de Taylor, e a potencial degeneração da densidade perto de estados de vácuo complica os argumentos de compacidade necessários para a passagem ao limite.

Metodologia
Os autores estabelecem a existência através de um processo de limite de dois níveis aplicado a um sistema aproximado cuidadosamente construído.

1. Regularização Especializada do Potencial:
   Para abordar a falta de limites em $L^2$ sobre $c$ e o fracasso das expansões de Taylor padrão, os autores introduzem uma regularização especializada $f^\delta_2(c)$ para o termo logarítmico singular. Diferentemente de abordagens anteriores, esta aproximação é construída de modo que a diferença $f^\delta_2(c) - \frac{\theta_c}{2}c^2$ cresça apenas linearmente quando $|c| \to \infty$. Esta construção específica é projetada para cancelar o crescimento quadrático induzido pelo termo anti-difusivo, restaurando assim a compacidade necessária e permitindo estimativas uniformes.

2. Esquema de Aproximação de Dois Níveis:

   • Nível 1 (Limite $\epsilon$): Os autores introduzem primeiro uma viscosidade artificial e parâmetros de regularização ($\epsilon$) para construir um sistema regularizado. Usando um argumento de ponto fixo (Teorema do Ponto Fixo de Schaefer), eles provam a existência de soluções fortes para este sistema aproximado. Estimativas uniformes independentes de $\epsilon$ são derivadas, permitindo que o limite $\epsilon \to 0$ seja tomado, resultando em uma solução fraca para um sistema intermediário com energia livre regularizada.
   • Nível 2 (Limite $\delta$): O segundo passo envolve remover a regularização do potencial ($\delta \to 0$). Uma grande dificuldade aqui é controlar a integrabilidade em $L^{3/2}$ da densidade. Os autores empregam uma escolha específica de pressão artificial e um argumento de interpolação delicado. Este argumento explora a estrutura do potencial regularizado, especificamente a observação de que $\frac{\partial f^\delta_2}{\partial c} - \theta_c c = O(\ln \frac{1}{\delta})$ para $|c| > 1+\delta$. Isso permite fechar as estimativas para a norma da densidade e estabelecer limites uniformes independentes de $\delta$.

3. Convergência Forte e Limites Físicos:
   Passar ao limite $\delta \to 0$ requer provar a convergência forte da densidade para lidar com o termo de pressão não monotônico. Adaptando técnicas da teoria das equações de Navier–Stokes compressíveis estacionárias, os autores estabelecem a convergência forte da densidade. Além disso, provam que a concentração limite $c$ satisfaz os limites físicos estritos $-1 < c < 1$ quase em todo lugar no suporte da densidade ($\{\rho > 0\}$), garantindo que o termo singular esteja bem definido.

Resultados Principais
O teorema principal (Teorema 1.1) estabelece a existência de uma solução fraca $(\rho, u, \mu, c)$ para o problema estacionário sob as seguintes condições:

• Domínio: Um domínio limitado e simplesmente conexo com fronteira $C^2$.
• Forças: Forças externas $g_1, g_2$ são limitadas.
• Expoente Adiabático ($\gamma$):
  • $\gamma > 3/2$ se a força externa for irrotacional ($\nabla \times g_1 = 0$).
  • $\gamma > 5/3$ se a força externa for rotacional ($\nabla \times g_1 \neq 0$).
• Restrições de Massa: A massa total $m_1 > 0$ e a massa relativa $m_2$ satisfazem $m_2 \in (-m_1, m_1)$.

A solução satisfaz a equação de continuidade renormalizada, a formulação fraca das equações de momento e concentração, e a desigualdade de energia. Crucialmente, a concentração satisfaz $-1 < c(x) < 1$ quase em todo lugar onde $\rho(x) > 0$.

Significado e Afirmativas
Os autores afirmam que este trabalho fornece o primeiro resultado de existência para soluções fracas do sistema estacionário de Navier–Stokes–Cahn-Hilliard compressível que incorpora simultaneamente:

1. Um potencial de energia livre singular (tipo Flory–Huggins), que é fisicamente consistente para misturas binárias, mas matematicamente desafiador devido à falta de limites globais na concentração.
2. A presença de estados de vácuo ($\rho = 0$), o que introduz degeneração nas equações.

O artigo observa que, embora o intervalo de expoentes adiabáticos admissíveis $\gamma$ possa ser extensível com base em resultados para equações de Navier–Stokes compressíveis estacionárias de fluido único, o intervalo atual é limitado pela dificuldade de obter estimativas ponderadas para a concentração perto da fronteira devido à singularidade da energia livre. O trabalho representa um passo significativo na ponte entre modelos evolutivos (onde desigualdades de energia fornecem controle) e modelos estacionários com potenciais singulares.