Consistent multiple-relaxation-time lattice Boltzmann method for the volume averaged Navier-Stokes equations

Este artigo propõe um novo método de Boltzmann em rede de múltiplos tempos de relaxamento consistente para as equações de Navier-Stokes médias por volume, que elimina velocidades espúrias e garante invariância galileana ao desacoplar a fração de vazio da densidade, permitindo a recuperação precisa das equações com erro de segunda ordem.

O essencial

Imagine que você está tentando prever como a água flui através de uma esponja gigante, ou como o sangue circula em um vaso cheio de células, ou ainda como a areia se move em um dreno. O problema é que, em escala microscópica, tudo é muito complicado: há milhões de partículas sólidas e líquidas interagindo de formas complexas.

Resolver essas equações uma por uma seria como tentar contar cada grão de areia de uma praia para prever a maré. Impossível e demorado demais.

É aqui que entra a Equação de Navier-Stokes Média por Volume (VANSE). Pense nela como uma "lente de zoom" que permite aos cientistas olhar para a esponja ou o sangue não como milhões de partículas individuais, mas como um fluido contínuo com "buracos" (os espaços vazios entre as partículas). É uma forma inteligente de simplificar a realidade para simulações computacionais.

O Problema: O "Fantasma" no Sistema

Os pesquisadores usavam um método chamado Método de Boltzmann de Rede (LBM) para simular isso. Imagine que o LBM é como um jogo de tabuleiro onde "partículas" de fluido pulam de casa em casa em uma grade.

O problema é que os métodos antigos (baseados em densidade) tinham um defeito grave: eles criavam "velocidades fantasmas".

• A Analogia: Imagine que você está tentando empurrar um carrinho de compras vazio em um supermercado. De repente, sem você fazer nada, o carrinho começa a acelerar sozinho para o lado, como se tivesse vida própria. Isso é o que acontecia nas simulações antigas: onde havia uma mudança brusca na quantidade de "buracos" (porosidade), o computador inventava um vento falso que não existia na realidade. Isso estragava todo o cálculo, tornando os resultados errados e instáveis.

Além disso, esses métodos antigos tinham dificuldade em lidar com fluidos muito viscosos (como mel ou óleo grosso), travando ou ficando instáveis.

A Solução: O "Mestre de Orquestra" (MRT)

Neste artigo, os autores (da Universidade Tsinghua e do Instituto de Tecnologia de Pequim) criaram uma nova versão desse método, chamada MRTLB-VANSE.

Vamos usar uma analogia para entender como eles consertaram o problema:

1. O Desacoplamento (A Quebra de Cadeia):
   Nos métodos antigos, a "densidade" do fluido e a "porosidade" (o quanto de espaço está ocupado por sólidos) estavam presas uma à outra de forma rígida. Se a porosidade mudava, a densidade era forçada a mudar de um jeito que criava erros.

   • A Solução: Eles criaram uma "equação de estado provisória". Pense nisso como colocar um amortecedor entre a densidade e a porosidade. Agora, quando a porosidade muda, a densidade não entra em pânico e não gera aqueles "ventos fantasmas". Eles separaram as variáveis para que uma não arrastasse a outra para o erro.

2. O Método de Relaxamento Múltiplo (MRT) - O Maestro:
   O método antigo usava um "relaxamento simples" (SRT), que é como ter um único maestro tentando controlar toda a orquestra ao mesmo tempo. Se um instrumento (um tipo de erro matemático) desafinava, o maestro não sabia como corrigir sem estragar o resto.

   • A Solução: Eles usaram o MRT (Relaxamento Múltiplo). Imagine agora uma orquestra onde cada seção (cordas, metais, percussão) tem seu próprio maestro sênior. O método MRT permite que o computador ajuste cada "tipo de erro" (como o estresse viscoso ou a pressão) de forma independente.
   • Eles adicionaram um "termo de penalidade" (uma correção matemática) que age como um "sistema de cancelamento de ruído". Se o computador começa a gerar um erro de velocidade fantasma, esse sistema detecta e cancela o erro instantaneamente, garantindo que a física do fluido permaneça realista.

O Resultado: Precisão e Estabilidade

Com essa nova abordagem, os pesquisadores provaram que:

• Sem Fantasmas: As "velocidades fantasmas" foram drasticamente reduzidas, mesmo em situações onde a porosidade muda bruscamente (como na fronteira entre areia e água).
• Precisão: O método é "segundo ordem", o que significa que, se você aumentar a resolução da simulação, o erro cai muito rápido, tornando os resultados extremamente precisos.
• Versatilidade: Funciona bem tanto para fluidos finos (como água) quanto para fluidos grossos (como mel), algo que os métodos antigos não conseguiam fazer bem.

Por que isso importa?

Essa melhoria é como passar de um mapa desenhado à mão, cheio de erros, para um GPS de alta precisão com satélites. Isso permite que engenheiros e cientistas simulem com muito mais confiança:

• Reatores químicos: Onde partículas sólidas e gases se misturam.
• Medicina: Entendendo melhor como o sangue flui em artérias com placas ou como células se movem.
• Meio Ambiente: Simulando como poluentes se movem no solo ou como a água flui em aquíferos.

Em resumo, os autores criaram um "super-herói" da simulação de fluidos que consegue lidar com cenários complexos e desordenados sem alucinar velocidades que não existem, tornando as previsões de engenharia muito mais seguras e confiáveis.

Resumo técnico

---

1. Problema e Motivação

As Equações de Navier-Stokes Médias por Volume (VANSE) são fundamentais para a simulação de sistemas multifásicos fluido-sólido em escala macroscópica (como leitos fluidizados, fluxo sanguíneo e meios porosos). Elas descrevem o comportamento do fluido considerando a fração de vazio (porosidade) sem resolver as interfaces detalhadas.

O método de Lattice Boltzmann (LB) é uma ferramenta popular para resolver essas equações devido à sua facilidade de implementação e paralelização. No entanto, os esquemas LB baseados em densidade existentes (especialmente os que utilizam o operador de colisão de tempo único, SRT) enfrentam dois problemas críticos:

1. Inconsistência Numérica: A recuperação das VANSE sofre de erros de alta ordem, levando a inconsistências no termo de gradiente de pressão e no tensor de tensão viscosa.
2. Velocidades Espúrias: Em campos de fração de vazio com grandes gradientes ou descontinuidades, os métodos anteriores geram "velocidades espúrias" (correntes não físicas) significativas, que podem exceder os limites de estabilidade do método e causar instabilidade numérica.

Embora métodos baseados em pressão tenham sido propostos recentemente para mitigar isso, eles ainda apresentam limitações em faixas amplas de viscosidade devido ao uso do operador SRT. Portanto, há uma necessidade urgente de desenvolver um esquema LB baseado em densidade que seja consistente, estável e livre de velocidades espúrias para VANSE.

2. Metodologia Proposta (MRTLB-VANSE)

Os autores propõem um novo método: o MRTLB-VANSE (Lattice Boltzmann de Múltiplos Tempos de Relaxamento para VANSE). A abordagem combina três inovações principais dentro de um esquema baseado em densidade:

A. Decuplagem da Fração de Vazio da Densidade

Para resolver a inconsistência no gradiente de pressão, os autores introduzem uma Equação de Estado Provisória ($p = \kappa \rho c_s^2$) na distribuição de equilíbrio da densidade.

• Isso permite modificar a distribuição de equilíbrio de densidade (usando uma expansão de Hermite de 3ª ordem) para decoplar a fração de vazio ($\phi$) da densidade ($\rho$) no operador de gradiente.
• A força de correção de pressão é redefinida para ser independente do gradiente de fração de vazio ($\nabla \phi$), eliminando oscilações não físicas nas descontinuidades.

B. Operador de Colisão MRT (Multiple-Relaxation-Time)

Em vez do operador BGK/SRT padrão, o método utiliza o operador MRT.

• O espaço de momentos permite o controle independente das taxas de relaxamento para diferentes modos (densidade, momento, tensão, etc.).
• Isso aumenta a estabilidade numérica, especialmente para viscosidades altas (onde o tempo de relaxamento seria baixo no esquema SRT).

C. Termo Fonte de Penalidade no Espaço de Momentos

A introdução da nova equação de estado provisória e da expansão de 3ª ordem gera erros no tensor de tensão viscosa recuperado, violando a invariância de Galileu.

• Para corrigir isso, os autores derivam e inserem um termo fonte de penalidade ($C$) no espaço de momentos durante o passo de colisão.
• Este termo elimina exatamente os erros numéricos indesejados no tensor de tensão viscosa, garantindo que as equações macroscópicas recuperadas sejam consistentes com as VANSE teóricas.

3. Validação Numérica e Resultados

O método foi validado através de uma análise de Chapman-Enskog (prova teórica de consistência) e diversas simulações numéricas:

• Fluxo em Meios Porosos Uniformes: O método recuperou com precisão as soluções analíticas para fluxos de Poiseuille e Couette em meios porosos, demonstrando capacidade de lidar com diferentes números de Darcy.
• Fluxo de Partículas Não Uniformes: Em cenários com gradientes espaciais e temporais de fração de vazio, o método mostrou-se consistente com soluções analíticas, corrigindo os erros de pressão observados em métodos anteriores.
• Teste de Campo Descontínuo (Sem Fluxo): Este é o teste crítico para velocidades espúrias. Em um campo com descontinuidade abrupta de fração de vazio e sem força externa:
  • Métodos anteriores (Zhang et al., Bukreev et al., Blais et al.) geraram velocidades espúrias significativas.
  • O MRTLB-VANSE reduziu drasticamente essas velocidades espúrias, mantendo o campo de fluxo quiescente, graças à preservação da isotropia da distribuição de densidade após o passo de streaming.
• Teste de Convergência (MMS - Método das Soluções Fabricadas):
  • Utilizando o MMS, os autores verificaram a ordem de convergência para campos de fração de vazio espaciais e espaço-temporais.
  • O método demonstrou convergência de segunda ordem precisa para a velocidade, tanto para fluxos com divergência zero quanto não zero.
  • Comparado ao método baseado em pressão de Fu et al. (que usa SRT), o MRTLB-VANSE manteve a precisão e estabilidade em viscosidades mais altas, onde o método SRT falhava devido à instabilidade.

4. Contribuições Principais

1. Consistência Teórica: Desenvolvimento de um esquema LB baseado em densidade que recupera as VANSE completas com consistência de segunda ordem, superando as limitações de esquemas anteriores.
2. Supressão de Velocidades Espúrias: A combinação da nova distribuição de equilíbrio e do termo de penalidade MRT elimina efetivamente as correntes espúrias em interfaces de fração de vazio descontínuas.
3. Robustez Numérica: O uso do operador MRT permite que o método funcione em uma faixa mais ampla de viscosidades e gradientes de porosidade, superando as limitações de estabilidade do operador SRT.
4. Eficiência Computacional: Ao evitar a necessidade de calcular gradientes de fração de vazio na força de correção de pressão, o método simplifica a implementação e reduz erros de discretização.

5. Significado e Perspectivas Futuras

Este trabalho representa um avanço significativo na modelagem de escoamentos multifásicos fluido-sólido. Ao fornecer um método LB baseado em densidade robusto e consistente, ele abre caminho para:

• Simulações mais confiáveis de processos industriais complexos (como combustão de coque, dissolução mineral e fraturamento capilar).
• A extensão do método para cenários tridimensionais (3D).
• O acoplamento com outros modelos mesoscópicos, como modelos de fases imiscíveis, transporte de espécies e modelos de partículas discretas (DEM), facilitando a simulação de processos líquido-gás-sólido complexos com alta eficiência computacional.

Em resumo, o MRTLB-VANSE resolve o dilema histórico entre a consistência teórica e a estabilidade numérica em simulações de meios porosos e multifásicos usando Lattice Boltzmann.