Upscaling the Navier-Stokes-Cahn-Hilliard model for incompressible multiphase flow in inhomogeneous porous media

Este trabalho apresenta a derivação rigorosa de um modelo macroscópico para o fluxo bifásico em meios porosos heterogêneos, obtido através da homogeneização das equações de Navier-Stokes e Cahn-Hilliard, incorporando formalmente o comportamento de molhabilidade e validando o modelo com simulações numéricas.

O essencial

Imagine que você está tentando entender como dois líquidos diferentes (como óleo e água) se misturam e se movem dentro de uma esponja muito complexa. Essa é a essência deste trabalho científico.

Aqui está uma explicação simples, usando analogias do dia a dia, do que os autores fizeram:

1. O Problema: A Esponja e o Labirinto

Pense em um meio poroso (como solo, rocha ou uma esponja) como um labirinto gigante feito de minúsculos túneis.

• Na escala microscópica (Pore-scale): Se você olhar de muito perto, através de um microscópio, você vê os líquidos desviando de cada pedrinha, colidindo e formando bolhas. É como tentar prever o movimento de cada gota de água em uma cachoeira. É extremamente difícil e lento de calcular para grandes áreas.
• Na escala macroscópica (Darcy-scale): Para engenheiros que querem saber quanto óleo sai de um poço ou como a água da chuva entra no solo, eles não querem saber o caminho de cada gota. Eles querem uma "média" geral, como se o líquido estivesse fluindo por um cano grosso e liso.

O problema é que as fórmulas tradicionais (chamadas de Lei de Darcy) são como regras de "chute": elas funcionam bem em casos simples, mas falham quando a coisa fica complexa, especialmente quando a "esponja" tem propriedades diferentes em lugares diferentes (heterogênea) ou quando a superfície da esponja "adora" um líquido e "odeia" o outro (molhabilidade).

2. A Solução: O Tradutor de Escalas

Os autores criaram um "tradutor matemático". Eles pegaram as leis físicas complexas que regem o movimento dos fluidos em cada micro-túnel (as equações de Navier-Stokes e Cahn-Hilliard) e usaram uma técnica chamada Média de Volume para "espremer" essa informação complexa e transformá-la em uma fórmula mais simples para a escala grande.

É como se você tivesse uma câmera de alta velocidade filmando uma multidão de pessoas correndo em um estádio cheio de obstáculos. Em vez de descrever o movimento de cada pessoa, você cria uma "média" que diz: "O fluxo geral de pessoas se move para o norte com velocidade X".

3. A Grande Inovação: O "Gosto" da Esponja

A parte mais brilhante deste trabalho é como eles lidaram com a molhabilidade (wetting).

• A Analogia: Imagine que a esponja é feita de um material que "ama" a água (hidrofílico) e "odeia" o óleo. Na física, isso significa que a água gruda nas paredes dos túneis e o óleo fica no meio.
• O Erro Antigo: Os modelos antigos tratavam essa interação de forma empírica (tentativa e erro), como se fosse um ajuste de botão aleatório.
• A Inovação: Os autores conseguiram incluir matematicamente esse "gosto" da esponja diretamente na equação principal. Eles mostraram que a força que puxa os líquidos (tensão superficial) não é apenas um número mágico, mas algo que pode ser calculado a partir das propriedades da superfície da esponja. É como se o modelo soubesse que a esponja é "pegajosa" para a água e ajustasse o fluxo automaticamente, sem precisar de palpites.

4. O Resultado: Um Modelo Unificado

Eles criaram uma única equação que funciona em dois níveis:

1. Se você olhar de perto, ela se comporta como a física real dos fluidos.
2. Se você olhar de longe, ela se transforma na Lei de Darcy (a regra usada na indústria).

É como ter um mapa que mostra tanto as ruas detalhadas de uma cidade quanto a visão de satélite de todo o país, e que consegue transitar entre os dois sem perder a precisão.

5. A Prova: Simulações Computacionais

Para provar que a teoria funcionava, eles usaram um computador para simular cenários reais:

• Fluxo em tubos: Verificaram se a velocidade do fluido batia com a teoria. (Bateu!)
• Deslocamento de fluidos: Simularam a injeção de um líquido para empurrar outro (como em recuperação de petróleo). O modelo conseguiu prever onde a frente de choque (a linha entre os dois líquidos) ia parar.
• Dedos Viscosos (Viscous Fingering): Quando um líquido menos viscoso (mais fino) empurra um mais viscoso (mais grosso), ele cria padrões que parecem dedos de uma mão. Eles mostraram que a "molhabilidade" da esponja pode impedir ou acelerar a formação desses "dedos". Se a esponja gosta do líquido que está entrando, os dedos se espalham mais; se não, eles ficam mais retos.

Resumo Final

Este trabalho é como construir uma ponte sólida entre a física microscópica complexa e a engenharia prática. Eles não apenas criaram uma nova fórmula, mas explicaram por que ela funciona, incorporando a "personalidade" das superfícies (molhabilidade) de forma rigorosa. Isso ajuda engenheiros a preverem com muito mais precisão como fluidos se movem em solos, rochas e materiais porosos, o que é crucial para coisas como limpeza de aquíferos, extração de petróleo e armazenamento de carbono.

Resumo técnico

Título: Escalonamento (Upscaling) do Modelo Navier-Stokes-Cahn-Hilliard para Escoamento Multifásico Incompressível em Meios Porosos Inhomogêneos

1. Problema e Contexto

O fluxo de fluidos multifásicos em meios porosos é fundamental em diversas aplicações de engenharia, como recuperação avançada de petróleo, remediação de águas subterrâneas e sequestro de CO2.

• Desafio: A modelagem direta na escala dos poros (DNS - Direct Numerical Simulation) é computacionalmente proibitiva para grandes escalas espaciais e temporais devido à complexidade geométrica dos meios porosos.
• Limitação Atual: As abordagens macroscópicas tradicionais baseadas na Lei de Darcy dependem fortemente de funções empíricas para permeabilidade relativa e pressão capilar. Essas funções frequentemente falham em capturar efeitos físicos complexos, como histerese, dependência da história de deslocamento e, crucialmente, os efeitos de molhabilidade (wetting) de forma rigorosa.
• Objetivo: Desenvolver um modelo macroscópico rigoroso, derivado diretamente das equações governantes na escala dos poros, que integre a dinâmica da interface e os efeitos de molhabilidade sem depender de ajustes empíricos arbitrários.

2. Metodologia

Os autores utilizam o método de Média Volumétrica (Volume Averaging Method) para escalar as equações do sistema Navier-Stokes-Cahn-Hilliard (NSCH) da escala dos poros para a escala de Darcy (REV - Representative Elementary Volume).

• Equações na Escala dos Poros:

  • O fluxo é governado pelas equações de Navier-Stokes para fluidos incompressíveis.
  • A evolução da interface entre as fases é descrita pela equação de Cahn-Hilliard, baseada em um funcional de energia livre que inclui energia de volume, energia interfacial e energia de parede (relacionada à molhabilidade).
  • A condição de contorno de molhabilidade é tratada via uma condição de Neumann não-linear na parede sólida.

• Processo de Escalonamento:

  1. Averaging (Média): As equações são integradas sobre um volume de controle (REV). As variáveis são decompostas em uma componente média intrínseca e uma componente de desvio espacial (decomposição de Gray).
  2. Teoremas de Média: Aplicação dos teoremas de média espacial e temporal para lidar com os termos de derivada e integração.
  3. Problemas de Fechamento (Closure): O processo gera termos não fechados (desvios espaciais). Para fechá-los:
     • Assumem-se restrições de escala de comprimento ($l \ll r_0 \ll L$).
     • Derivam-se equações para os desvios ($\hat{\phi}$, $\hat{u}$) dentro do REV.
     • Modelam-se os termos de desvio como funções das grandezas médias e coeficientes de transporte (tensores de dispersão, permeabilidade relativa, etc.).
  4. Incorporação da Molhabilidade: Um ponto central é a derivação formal de como a condição de contorno de molhabilidade na parede sólida se incorpora diretamente ao potencial químico médio na escala macroscópica.

• Validação Numérica:

  • Foi desenvolvido um método de Lattice Boltzmann (LBM) adaptado para resolver as equações médias derivadas.
  • O modelo foi testado em cenários de fluxo monofásico e bifásico.

3. Principais Contribuições

1. Modelo de Fluido Único Rigoroso: Diferente de modelos anteriores que tratam as fases separadamente (dois fluidos), este trabalho deriva um modelo de fluido único macroscópico. Isso elimina a necessidade de condições de contorno interfaciais empíricas explícitas.
2. Potencial Químico Médio Formal: A maior inovação é a incorporação formal do comportamento de molhabilidade no potencial químico médio ($\langle \mu \rangle_\alpha$). A pressão capilar e os efeitos de molhabilidade emergem naturalmente da derivação, em vez de serem inseridos via funções empíricas de pressão capilar.
3. Ponte entre Escalas: O modelo estabelece uma conexão teórica rigorosa entre a dinâmica microscópica (tensão superficial, ângulo de contato) e o comportamento macroscópico (Lei de Darcy generalizada, equação de Brinkman).
4. Generalização da Lei de Darcy: O modelo recupera a Lei de Darcy e a equação de Brinkman-Forchheimer como casos limites, mas com termos de tensão interfacial e viscosidade efetiva derivados termodinamicamente consistentemente.

4. Resultados

Os autores realizaram simulações numéricas para validar o modelo:

• Fluxo de Poiseuille em Meio Poroso Homogêneo: Os resultados numéricos concordaram excelentemente com a solução analítica, validando a parte hidrodinâmica do modelo (termos de arrasto e viscosidade efetiva).
• Fluxo em Canal Fluido-Poroso: O modelo conseguiu capturar corretamente a descontinuidade de tensão e a continuidade de velocidade na interface entre a região livre e a porosa, comparando-se bem com soluções analíticas para diferentes porosidades.
• Problema de Buckley-Leverett: A simulação da frente de deslocamento bifásico mostrou boa concordância com a solução analítica clássica, capturando a formação de choques e a saturação frontal, embora com pequenas oscilações numéricas típicas de esquemas de diferenças centrais.
• Instabilidade Viscosa (Viscous Fingering): Simulações qualitativas demonstraram a capacidade do modelo de capturar a formação de dedos viscosos. O estudo revelou que a molhabilidade (ângulo de contato) tem um efeito significativo na dispersão dos dedos; superfícies mais hidrofílicas tendem a dispersar os dedos incipientes e inibir sua formação, o que está em consonância com observações experimentais.

5. Significado e Conclusões

Este trabalho oferece uma estrutura teórica fundamental para a modelagem de fluxo multifásico em meios porosos. Ao derivar as equações macroscópicas a partir de princípios termodinâmicos e de conservação de massa/momento na escala dos poros, o modelo supera as limitações das abordagens puramente empíricas.

• Impacto: Permite prever o comportamento de fluxo em meios complexos onde a molhabilidade e a heterogeneidade são críticas, sem depender de calibração extensiva de parâmetros de pressão capilar.
• Futuro: Os autores apontam que alguns parâmetros de fechamento (coeficientes de dispersão e tensores de permeabilidade) ainda precisam ser melhor quantificados para geometrias gerais, e que a relação entre o novo termo de força interfacial e funções clássicas (como a função J de Leverett) merece análise futura.

Em resumo, o artigo fornece uma base teórica robusta e unificada que conecta a física da escala dos poros à escala de Darcy, oferecendo uma ferramenta poderosa para simulações mais precisas em engenharia de reservatórios e ciências ambientais.