Lattice Boltzmann model for non-ideal compressible fluid dynamics

Este artigo apresenta um novo modelo de Boltzmann em rede que utiliza redes de primeiros vizinhos e termos de correção de quase-equilíbrio para simular escoamentos de fluidos compressíveis e não ideais com consistência termodinâmica e estabilidade numérica, validando com sucesso sua precisão por meio de simulações quantitativas de interações entre ondas de choque e gotas em números de Mach de até 1,47.

O essencial

Imagine que você está tentando simular como um fluido se comporta em um computador. Por muito tempo, os computadores foram excelentes em simular fluidos "ideais" — como a água fluindo suavemente em um rio ou o ar movendo-se lentamente ao redor de uma asa. Esses fluidos seguem regras simples e previsíveis.

Mas o que acontece quando o fluido está sob pressão e calor extremos, comportando-se como um gás denso que é quase um líquido, ou um líquido que é quase um gás? Este é o mundo dos fluidos compressíveis não ideais. Pense nisso como um fluido que está "estressado" e agindo de forma estranha, recusando-se a seguir as regras simples do mundo ideal. Isso ocorre em tecnologias avançadas como turbinas de CO2 supercrítico e ciclos energéticos orgânicos.

O problema é que as ferramentas computacionais existentes lutam com esses fluidos estressados. Elas ou travam, fornecem respostas erradas ou exigem tanta potência de computação que se tornam inúteis.

Este artigo apresenta uma nova e mais inteligente maneira de simular esses fluidos complicados usando um método chamado Método de Boltzmann em Rede (LBM). Aqui está como a nova abordagem dos autores funciona, explicada através de analogias simples:

1. O Sistema de "Duas Pistas"

A maioria dos métodos de simulação antigos tenta rastrear tudo (massa, velocidade, energia) com um único conjunto complicado de regras. Os autores perceberam que isso era como tentar dirigir um carro enquanto simultaneamente equilibra uma dúzia de bolas — fica bagunçado e instável.

Em vez disso, eles construíram um sistema de duas pistas:

• Pista A (A Multidão): Um conjunto de regras rastreia a densidade e a velocidade do fluido (quantas partículas existem e para onde estão indo).
• Pista B (A Energia): Um segundo conjunto separado de regras rastreia a energia total.
  Ao separar esses elementos, o computador não fica confuso. É como ter um controlador de tráfego dedicado para carros e outro separado para o combustível, garantindo que nenhum sistema derrube o outro.

2. O Atrator de "Quase-Equilíbrio"

Na física, os fluidos naturalmente desejam se estabelecer em um estado calmo chamado "equilíbrio". No entanto, nessas condições extremas, o fluido está constantemente sendo empurrado e puxado, então ele nunca se estabelece completamente.

Os autores inventaram um truque inteligente chamado atratore de "quase-equilíbrio".

• A Analogia: Imagine um cachorro perseguindo uma bola. A bola representa o "estado de calma perfeito". O cachorro representa o fluido. Em uma situação normal, o cachorro corre diretamente para a bola.
• O Problema: Neste fluido extremo, a bola continua se afastando ou mudando de forma. Se o cachorro apenas perseguir a bola cegamente, ele pode correr para fora de um penhasco (a simulação torna-se instável).
• A Solução: Os autores deram ao cachorro um "GPS" que prevê onde a bola estará uma fração de segundo depois, com base em como o vento (pressão) e o terreno (densidade) estão mudando. Esse alvo "deslocado" permite que o cachorro corra suavemente sem cair do penhasco. Isso garante que a simulação permaneça estável mesmo quando o fluido está se movendo muito rápido ou mudando de densidade rapidamente.

3. Corrigindo o Calor "Espúrio"

Quando os fluidos se movem rápido, eles geram calor. Em modelos computacionais padrão, o calor às vezes flui na direção errada ou cria calor "fantasma" falso que não existe na realidade.

• A Analogia: É como um termostato que acha que o quarto está congelando porque está medindo a corrente de ar de uma janela, e não a temperatura real do quarto.
• O Conserto: Os autores adicionaram um termo de "correção" específico às suas equações. Isso atua como um filtro que remove as correntes de ar falsas, garantindo que o calor flua exatamente como a física dita (Lei de Fourier), mesmo nessas condições extremas e não ideais.

4. Os Testes de "Tubo de Choque" e "Coluna Líquida"

Para provar que seu novo método funciona, eles não apenas fizeram matemática; realizaram testes extremos:

• O Tubo de Choque: Eles simularam uma explosão súbita de pressão (uma onda de choque) movendo-se através de um gás denso. Em gases normais, essas ondas se comportam de uma maneira. Nesses gases "não ideais", as ondas podem fazer algo estranho: um "choque de rarefação" (uma onda que se espalha, mas ainda age como um choque agudo). Seu modelo previu com sucesso esse comportamento estranho, que os modelos mais antigos perderam.
• A Coluna Líquida: Eles simularam uma onda de choque de alta velocidade atingindo uma gota de líquido. Este é um teste muito difícil porque o choque salta, reflete e rasga a gota. Seu modelo lidou perfeitamente com a colisão, correspondendo aos experimentos do mundo real onde a gota de líquido se achata e se expande exatamente como deveria.

Por Que Isso Importa

Os autores afirmam que seu método é rápido, estável e preciso. Ele usa uma grade simples (como um tabuleiro de xadrez padrão) em vez de precisar de uma grade supercomplexa e esticada. Isso significa que os cientistas agora podem simular esses fluxos de fluidos extremos, de alta velocidade e não ideais em computadores padrão com alta precisão.

Em resumo: O artigo apresenta um novo "manual do motorista" para simulações computacionais que permite lidar com fluidos sob estresse extremo sem travar. Ao usar um sistema de duas pistas e um "GPS" inteligente para a energia do fluido, eles podem prever com precisão como esses fluidos complexos se comportam em cenários de alta velocidade, abrindo caminho para melhores projetos de sistemas energéticos avançados.

Resumo técnico

Resumo Técnico: Modelo de Boltzmann em Rede para Dinâmica de Fluidos Compressíveis Não Ideais

Enunciação do Problema
A dinâmica de fluidos compressíveis não ideais (NICFD) aborda o comportamento dos fluidos nos regimes próximo, trans e supercrítico, onde os fluidos desviam-se significativamente das leis dos gases ideais. Esses regimes são críticos para tecnologias energéticas modernas, como ciclos Rankine orgânicos e turbinas de CO₂ supercrítico. No entanto, dados experimentais nesses estados são escassos e difíceis de obter. Embora a Simulação Numérica Direta (DNS) seja necessária para compreender a física complexa (por exemplo, estruturas de choque, turbulência em vapores densos) e gerar bancos de dados de engenharia, as ferramentas numéricas existentes frequentemente lutam com consistência termodinâmica, estabilidade e a recuperação precisa dos limites hidrodinâmicos para equações de estado (EOS) não ideais. Especificamente, os métodos convencionais de Dinâmica dos Fluidos Computacional (CFD) enfrentam desafios na região espinodal termodinamicamente instável, onde o quadrado da velocidade do som adiabática torna-se negativo, exigindo tratamentos especiais como a regra da área igual de Maxwell. O Método de Boltzmann em Rede (LBM), embora bem-sucedido para fluxos incompressíveis e multifásicos, historicamente negligenciou a termodinâmica compressível não ideal em suas formulações cinéticas, frequentemente dependendo de regularização ad hoc ou estênceis estendidos que comprometem a eficiência ou a invariância de Galileu.

Metodologia
Os autores propõem um modelo cinético inovador e sua realização em Lattice Boltzmann, projetados para recuperar as equações compressíveis de Navier–Stokes–Korteweg (NSK) para fluidos não ideais. A metodologia é construída sobre os seguintes pilares:

1. Abordagem de Função de Distribuição Dupla (DDF): O modelo emprega duas funções de distribuição: $f$ para massa e momento, e $g$ para energia volumétrica ($\rho E$). Esse desacoplamento permite o uso de estênceis clássicos de primeiros vizinhos (especificamente D3Q27) sem exigir estênceis estendidos.
2. Atratores de Quase-Equilíbrio: O operador de colisão utiliza uma estrutura dupla de relaxação tipo BGK envolvendo um atrator de equilíbrio local ($f^{eq}, g^{eq}$) e um atrator de quase-equilíbrio deslocado ($f^\star_\lambda, g^\star_\lambda$).
   • O equilíbrio local é definido usando um parâmetro de temperatura de referência $\theta = P/\rho$ (trabalho de fluxo termodinâmico) em vez da temperatura termodinâmica $T$, garantindo a largura correta da distribuição de Maxwell independente da energia interna.
   • O equilíbrio deslocado introduz correções para a força de Korteweg (tensão superficial), viscosidade volumétrica e fluxo de calor. Ele emprega valores deslocados para velocidade de fluxo ($u^\star$), temperatura de referência ($\theta^\star$) e temperatura termodinâmica ($T^\star$).
3. Consistência Termodinâmica:
   • A temperatura de referência $\theta$ é definida como $P/\rho$, uma formulação que garante que o modelo permaneça válido para EOS não ideais (por exemplo, van der Waals).
   • Um termo de correção no fluxo de calor ($q^c_\lambda$) é introduzido para compensar componentes espúrios não-Fourier que surgem de gradientes de densidade em fluidos não ideais, garantindo a recuperação da lei de Fourier.
   • O deslocamento na temperatura de referência permite que a viscosidade volumétrica seja um parâmetro independente, ajustável e positivo-definido, contornando o problema da viscosidade volumétrica negativa que pode surgir em modelos BGK padrão para pressões não ideais.
4. Discretização: O modelo utiliza um esquema explícito de segunda ordem de precisão derivado da integração trapezoidal ao longo das características. Emprega um estêncil padrão D3Q27 e trata termos fonte (força de Korteweg) via equilíbrio deslocado, evitando erros de discretização explícita de termos fonte.

Principais Contribuições

• Novo Framework Cinético: O artigo introduz um modelo cinético que recupera as equações hidrodinâmicas alvo (balanço de massa, momento e energia volumétrica) para fluidos compressíveis não ideais usando um estêncil cinético mínimo (primeiros vizinhos).
• Dissipação Positiva-Definida: Ao construir o modelo via atratores de quase-equilíbrio, os autores garantem que as taxas de dissipação de Navier–Stokes (viscosidade de cisalhamento e volumétrica) sejam positiva-definidas e invariantes de Galileu, um ponto de falha comum em tentativas anteriores de LBM para fluxos compressíveis não ideais.
• Robustez Termodinâmica: O modelo lida naturalmente com as regiões de instabilidade termodinâmica (espinodal) onde a velocidade do som torna-se imaginária no limite invíscido, pois o LBM herda a propagação ao longo de características discretas fixas, em vez de depender de solvers de EDP hiperbólicas que exigem tratamento especial nessas regiões.
• Validação de Fenômenos Não Clássicos: O estudo fornece a primeira validação quantitativa do LBM para interações choque-gota em números de Mach de até 1,47 em um contexto de fluido não ideal, demonstrando a capacidade do modelo de capturar estruturas de onda complexas.

Resultados
O modelo foi validado através de uma hierarquia de configurações cada vez mais complexas:

1. Estabilidade Linear e Dispersão: Simulações da velocidade do som, dissipação de onda de cisalhamento e condutividade térmica mostraram excelente concordância com soluções analíticas para um fluido de van der Waals ajustado às propriedades do nitrogênio em uma faixa de números de Mach e temperaturas reduzidas.
2. Consistência Multifásica: O modelo reproduziu com precisão as densidades de coexistência líquido-vapor e perfis de interface, demonstrando convergência de segunda ordem com o refinamento da malha e concordância com a regra da área igual de Maxwell.
3. Tubos de Choque Não Clássicos: Três casos de tubos de choque 1-D foram simulados para investigar a dinâmica de ondas não clássicas caracterizadas pelo derivado fundamental da dinâmica de gases ($\Gamma$).
   • Casos com $\Gamma < 0$ capturaram corretamente a formação de choques de rarefação (ondas de expansão afiadas) e ventiladores de compressão (ondas de compressão espalhadas), fenômenos distintos da dinâmica de gases clássica.
   • O modelo transitou com sucesso entre os regimes clássico ($\Gamma > 0$) e não clássico ($\Gamma < 0$) dentro de uma única simulação.
4. Interação Choque-Coluna Líquida: Foi realizada uma simulação 2D de uma onda de choque plana interagindo com uma coluna de líquido saturado próximo ao ponto crítico ($M_a = 1,47$).
   • O modelo capturou estruturas de onda complexas, incluindo choques transmitidos, choques refletidos, hastes de Mach e pontos triplos.
   • A comparação quantitativa da evolução da largura da coluna líquida contra dados experimentais e resultados numéricos anteriores mostrou boa concordância.
   • A análise termodinâmica revelou que a interação do choque leva o vapor circundante a um estado supercrítico, enquanto a coluna líquida permanece marginalmente subcrítica, capturando com precisão o forte acoplamento termodinâmico.

Significado e Alegações
Os autores afirmam que este trabalho estende a aplicabilidade dos métodos de Boltzmann em Rede para fluxos compressíveis não ideais de alta velocidade com um estêncil cinético mínimo. O significado primário reside em fornecer um framework termodinamicamente consistente e numericamente estável que não depende de estênceis estendidos ou regularização ad hoc. O modelo é apresentado como uma teoria cinética reduzida que visa o limite hidrodinâmico NSK, capaz de lidar com toda a gama de comportamentos não ideais, incluindo a região espinodal, sem as instabilidades numéricas frequentemente associadas ao CFD convencional nesses regimes. A simulação bem-sucedida de interações choque-gota em altos números de Mach serve como uma validação rigorosa, demonstrando o potencial do método para investigar fluxos complexos em tecnologias supercríticas e processos impulsionados por transições de fase.