Monolithic kinetic algorithm for heterogeneous porous media systems using a continuous one-domain approach

Resumo Técnico: Algoritmo Cinético Monolítico para Sistemas de Meios Porosos Heterogêneos

Definição do Problema
A simulação do fluxo de fluidos em meios porosos heterogêneos, particularmente em sistemas que acoplam regiões de fluido livre com estruturas porosas, apresenta desafios significativos. As abordagens tradicionais frequentemente dependem de formulações de "dois domínios", que resolvem equações governantes distintas para regiões de fluido livre e poroso e requerem condições de contorno de interface explícitas (por exemplo, condições de Beavers-Joseph) para acoplá-las. Alternativamente, abordagens de "um domínio" utilizam um único conjunto de equações, mas historicamente enfrentaram dificuldades com a consistência em meios heterogêneos.

As formulações existentes do Método de Lattice Boltzmann (LBM) para meios porosos frequentemente sofrem de uma escala inconsistente da porosidade nas equações macroscópicas recuperadas. Especificamente, muitos modelos falham em escalar corretamente o gradiente de pressão e os termos de tensão viscosa com a porosidade quando a porosidade varia espacialmente. Além disso, muitas formulações expressam a tensão viscosa de Brinkman em termos de gradientes da velocidade intrínseca, em vez da velocidade superficial (Darcy), o que leva a imprecisões em domínios heterogêneos. Ademais, os modelos LBM padrão frequentemente carecem da flexibilidade para ajustar independentemente a viscosidade volumétrica efetiva ou a velocidade do som, limitando a eficiência computacional e a estabilidade.

Metodologia
Os autores propõem um modelo de Lattice Boltzmann monolítico de domínio único baseado em uma abordagem contínua de um domínio. A metodologia é construída sobre os seguintes pilares:

1. Modelo Macroscópico: A hidrodinâmica alvo é derivada utilizando o método de média de volume, resultando em um modelo generalizado de Darcy-Forchheimer-Brinkman. Este modelo é válido para todo o domínio computacional, tratando a transição entre regiões de fluido livre e poroso como uma variação contínua das propriedades físicas ($\epsilon$). As equações governantes incluem:

   • Uma equação de continuidade para a densidade média de volume ($\epsilon\rho$).
   • Uma equação de momento incorporando arrasto de Darcy, arrasto inercial de Forchheimer e um tensor de tensão viscosa efetiva ($\sigma_{eff}$).
   • Crucialmente, a tensão viscosa efetiva é expressa em termos dos gradientes da velocidade superficial ($v$), não da velocidade intrínseca ($u = v/\epsilon$), garantindo consistência com a formulação de Brinkman para porosidade variável.

2. Solver Cinético: O modelo utiliza uma única equação cinética generalizada em redes padrão (ex: D2Q9, D3Q19).

   • Distribuição de Equilíbrio: É empregada uma distribuição de equilíbrio discreta em forma de produto para garantir momentos corretos e invariância galileana. Isso inclui termos de correção específicos ($\Phi_{ext}$ e $\Phi_{\eta}$) para restaurar a isotropia no tensor de tensão viscosa e permitir o ajuste independente da viscosidade volumétrica.
   • Esquema de Forçamento: Um esquema de forçamento generalizado consistente (método Relax-to-Force) é usado para incorporar o arrasto de Darcy-Forchheimer, forças de corpo e as correções necessárias para gradientes de porosidade.
   • Escalonamento de Porosidade: Os campos localmente conservados são definidos como densidade e momento médios de volume, incorporando naturalmente o escalonamento de porosidade nos momentos cinéticos.
   • Avaliação de Velocidade: Um algoritmo iterativo é proposto para resolver a relação implícita entre a velocidade macroscópica e os termos de forçamento, lidando com o arrasto não linear de Darcy-Forchheimer e as correções viscosas.

3. Condições de Contorno: Os autores adaptam as condições de contorno de nós úmidos padrão (esquema de equilíbrio, extrapolação de não-equilíbrio e esquema de Tamm-Mott-Smith) para a formulação de meio poroso. Um procedimento específico é derivado para calcular a densidade do fluido intrínseco nas fronteiras usando continuidade de massa, levando em conta a presença do termo de forçamento e as variações de porosidade.

Principais Contribuições

• Formulação de Um Domínio Consistente: O artigo apresenta a primeira formulação LBM que recupera um tensor de tensão viscosa de Brinkman com escalonamento de porosidade totalmente consistente, expresso em termos da velocidade superficial. Isso permite que o modelo lide com porosidade espacialmente variável sem condições de interface artificiais.
• Correção de Termos de Pressão e Viscosidade: Os autores derivam e implementam correções de forçamento específicas para retificar inconsistências nos termos de gradiente de pressão e tensão viscosa que surgem em meios heterogêneos, garantindo que o modelo recupere o limite macroscópico correto.
• Ajuste Independente de Parâmetros: O modelo permite que a velocidade do som e a viscosidade volumétrica efetiva sejam especificadas independentemente da viscosidade de cisalhamento. Isso aumenta a eficiência computacional (permitindo velocidades de som mais altas) e a estabilidade numérica (ajustando a viscosidade volumétrica).
• Algoritmo Monolítico: A abordagem utiliza uma única equação e solver cinético para todo o domínio, simplificando a implementação para geometrias complexas envolvendo tanto fluido livre quanto regiões porosas.

Resultados e Validação
O modelo foi validado através de uma série de testes de referência (benchmarks) e estudos numéricos:

• Análise de Consistência: A análise espectral temporal linear confirmou que o modelo recupera a velocidade do som efetiva correta e as taxas de dissipação para modos de cisalhamento e acústicos. A análise demonstrou que o equilíbrio em forma de produto com as correções propostas preserva a invariância galileana e a isotropia, ao contrário de formulações de equilíbrio de segunda ordem anteriores.
• Fluxos de Referência:
  • Fluxo de Couette: Simulações em um canal poroso homogêneo sob gravidade mostraram excelente concordância com soluções numéricas de referência para ambos os perfis de velocidade e densidade em vários números de Reynolds e de Darcy.
  • Fluxo de Poiseuille: O modelo capturou com precisão o fluxo em canais homogêneos e, crucialmente, em canais com perfis de porosidade heterogêneos. Ele previu corretamente efeitos de canalização de fluxo, onde o fluido flui preferencialmente através de regiões de alta porosidade.
• Demonstrações Numéricas:
  • Fluxo ao Redor de um Cilindro Permeável: O modelo simulou com sucesso o fluxo confinado ao redor de um cilindro poroso, identificando cinco regimes distintos baseados no número de Reynolds e na porosidade do núcleo. Capturou fenômenos como estabilização de esteira, desprendimento de recirculação e instabilidades de esteira distante, mostrando boa concordância qualitativa com a literatura.
  • Cavidade Impulsionada por Tampa com Obstáculo Poroso: O modelo simulou uma cavidade com um obstáculo poroso interno, demonstrando como a permeabilidade afeta a formação de vórtices, a penetração do fluxo e os padrões de recirculação. Os resultados destacaram a transição do comportamento sólido para o comportamento permeável conforme a porosidade aumenta.

Significância e Alegações
Os autores afirmam que este trabalho fornece um framework robusto e unificado para simular o transporte em meios porosos heterogêneos. Ao eliminar a necessidade de condições de contorno de interface explícitas e garantir o escalonamento consistente de todos os termos hidrodinâmicos, o modelo é particularmente adequado para:

• Estudos de design em estágio inicial, em escala de dispositivo, de estruturas porosas projetadas com porosidade espacialmente variável.
• Simulação de interações complexas entre fluidos livres e corpos permeáveis (ex: células de combustível, sistemas de filtragem, combustão porosa).
• Captura de fenômenos de transporte em sistemas onde a transição entre fluxo livre e meio poroso é contínua ou onde interfaces nítidas são difíceis de resolver.

O artigo conclui que, embora o estudo atual se concentre na validação 2D, a formulação é inerentemente multidimensional e pronta para extensão para geometrias 3D. Identifica-se que trabalhos futuros são necessários para incluir o balanço de energia (transferência de calor), misturas multicomponentes e fluxos multifásicos para abordar plenamente o transporte acoplado em sistemas de engenharia realistas.