Brinkman's term for the steady-state Navier-Stokes equations with prescribed flux rate or pressure drop in perforated pipes

Resumo Técnico: Homogeneização das Equações de Navier-Stokes em Regime Estacionário em Tubos Perfurados com Condições de Contorno Mistas

Enunciado do Problema
Este trabalho investiga o comportamento assintótico do movimento estacionário de um fluido viscoso e incompressível dentro de um tubo distorcido contendo um grande número de pequenas partículas sólidas. O domínio $\Omega_\varepsilon$ é um subconjunto perfurado de um domínio limitado e admissível $\Omega \subset \mathbb{R}^3$, obtido pela remoção de $N(\varepsilon)$ partículas de diâmetro $\varepsilon^3$ separadas por uma distância mútua $\varepsilon$. Esta configuração corresponde ao regime crítico ($\alpha = 3$) na teoria de homogeneização, onde o tamanho das partículas escala com o cubo da distância entre as partículas.

O escoamento do fluido é governado pelas equações de Navier-Stokes em regime estacionário com condições de contorno mistas. Especificamente:

1. Condições de não-deslizamento (no-slip) são impostas nas paredes físicas do tubo ($\Gamma_W$) e nas superfícies das partículas ($\partial K_n^\varepsilon$), onde as partículas podem se mover com velocidades constantes prescritas $\mu_n^\varepsilon$.
2. Condições de pressão de Bernoulli são aplicadas na entrada ($\Gamma_I$) e na saída ($\Gamma_O$). Estas envolvem a derivada normal da velocidade e a pressão de Bernoulli ($\Phi = p + \frac{1}{2}|u|^2$).
3. Dois cenários distintos são analisados:
   • Fluxo Líquido Prescrito: O fluxo total $F$ através da saída é fixo, enquanto a queda de pressão é uma constante desconhecida determinada pela solução.
   • Queda de Pressão Prescrita: A queda de pressão entre a entrada e a saída é fixa, enquanto o fluxo resultante é uma quantidade desconhecida.

O principal desafio matemático reside em derivar as equações efetivas (homogeneizadas) para o limite quando $\varepsilon \to 0$ sem impor suposições de pequenez sobre a força externa, a taxa de fluxo ou a queda de pressão. Isso é complicado pelo fato de que as estimativas de energia padrão para as equações de Navier-Stokes frequentemente falham em fornecer limites uniformes nestes problemas de valor de contorno misto sem restrições adicionais.

Metodologia
Os autores empregam o método da energia dentro do framework da teoria de homogeneização, seguindo a estratégia geral de Tartar. A prova procede em três etapas principais:

1. Limites Uniformes (estimativas independentes de $\varepsilon$):
   Estabelecer limites uniformes para as soluções $(u_\varepsilon, \Phi_\varepsilon)$ é o passo mais crítico e inovador. Diferente dos casos padrão onde estimativas de energia simples bastam, os autores não podem confiar apenas no teste da equação de momento com a própria solução devido às condições de contorno mistas e à queda de pressão desconhecida.

   • Em vez disso, eles utilizam um argumento de contradição baseado na lei de Bernoulli para soluções das equações de Euler estacionárias.
   • Eles assumem que as normas das soluções explodem conforme $\varepsilon \to 0$, normalizam a sequência e mostram que o limite satisfaz as equações de Euler.
   • Ao analisar o comportamento de contorno da pressão de Bernoulli limite nas paredes do tubo (usando resultados de regularidade local e a geometria específica da entrada/saída), eles demonstram que a pressão limite deve ser constante nas paredes. Isso leva a uma contradição com as condições de fluxo ou de queda de pressão, provando assim a limitação uniforme da sequência.
   • Para o problema de fluxo prescrito, uma suposição geométrica adicional é feita para que as partículas não se acumulem perto dos cilindros de entrada e saída, a fim de facilitar a análise de contorno da pressão limite.

2. Construção de Corretores e Operadores de Restrição:
   Para lidar com o domínio perfurado, os autores constroem corretores de Stokes $(\omega_{n,1}^\varepsilon, \omega_{n,2}^\varepsilon)$ no exterior das partículas e nas regiões anulares que as circundam. Estes corretores são usados para definir um operador de restrição $R_\varepsilon$.

   • Diferente de trabalhos anteriores no regime subcrítico, onde corretores capacitários escalares eram suficientes, o regime crítico exige estes corretores de Stokes vetoriais para dar conta da capacidade relativa não nula do enxame de partículas.
   • O operador de restrição $R_\varepsilon(\phi)$ mapeia funções de teste do domínio homogeneizado $\Omega$ para o domínio perfurado $\Omega_\varepsilon$, modificando-as próximo às partículas para satisfazer as condições de contorno nulas, enquanto preserva a propriedade de divergência nula.

3. Passagem ao Limite:
   Utilizando os limites uniformes e o operador de restrição, os autores passam ao limite na formulação fraca das equações de Navier-Stokes.

   • O termo convectivo não linear converge via convergência forte em $L^4$ (garantida pelos limites uniformes e imersões compactas).
   • A interação entre o fluido e as partículas gera um termo adicional na equação limite. Este termo surge do comportamento limite da soma das forças de arrasto exercidas pelas partículas, caracterizada pelas matrizes de resistência de Stokes $G_n$.

Contribuições Principais e Resultados
O artigo estabelece dois teoremas principais (Teoremas 2.4 e 2.5) descrevendo o limite homogeneizado do sistema:

• Emergência do Termo de Brinkman: As equações homogeneizadas no domínio limite $\Omega$ são as equações de Navier-Stokes aumentadas por um termo de Brinkman. A equação efetiva assume a forma:
  $$-\Delta u + (\nabla \times u) \times u + \nabla \Phi + Gu - J = f$$
  onde:

  • $Gu$ representa um termo de fricção proporcional à velocidade do fluido $u$, originado do arrasto coletivo das partículas (com $G$ sendo o limite da soma escalonada das matrizes de resistência de Stokes).
  • $J$ é um termo de fonte proveniente do fluxo limite das próprias partículas em movimento.
  • As condições de contorno na entrada e na saída permanecem do tipo misto envolvendo a pressão de Bernoulli.

• Extensão ao Regime Crítico com Dados Não Nulos: Os resultados estendem as descobertas de homogeneização anteriores (especificamente em [59] para o regime subcrítico) para o regime crítico ($\alpha=3$). Crucialmente, estes resultados sustentam-se sem qualquer suposição de pequenez sobre a força externa, taxa de fluxo ou queda de pressão. Esta é uma progressão significativa, visto que muitos resultados de homogeneização para Navier-Stokes dependem de suposições de dados pequenos para garantir unicidade ou existência de limites.

• Tratamento de Partículas em Movimento: O framework acomoda partículas com velocidades constantes não nulas $\mu_n^\varepsilon$, desde que sua energia cinética seja uniformemente limitada (Suposição 2.11). Isso permite uma modelagem mais realista das interações fluido-partícula em comparação com modelos de partículas fixas.

Significância
O artigo afirma que sua contribuição mais significativa é a derivação rigorosa da lei de Brinkman para as equações de Navier-Stokes estacionárias no regime crítico sob dados irrestritos. Os autores enfatizam que a derivação dos limites uniformes via argumento de contradição envolvendo a lei de Bernoulli é uma questão delicada, particularmente porque a estrutura complexa das perfurações deve ser considerada para garantir que o argumento se sustente ao nível de $\varepsilon$.

Ao construir com sucesso um verdadeiro operador de restrição que não modifica funções de teste que se anulam nas partículas e ao lidar com a convergência do termo de pressão em um cenário polidisperso com velocidades de partículas não nulas, o trabalho fornece uma base matemática robusta para modelar o escoamento de fluidos em tubos perfurados onde a densidade de partículas é alta o suficiente para induzir um arrasto do tipo Brinkman, mas não tão alta a ponto de resultar na lei de Darcy. Os resultados são diretamente aplicáveis a campos como hemodinâmica, engenharia de petróleo e aerodinâmica, onde o escoamento de fluidos através de arranjos de obstáculos sólidos é um fenômeno fundamental.