Two-dimensional quantum lattice gas algorithm for anisotropic Burger-like equations

Resumo Técnico: Algoritmo de Gás em Rede Quântica Bidimensional para Equações do Tipo Burger Anisotrópicas

1. Problema e Contexto

A aplicação da Computação Quântica (QC) à Dinâmica dos Fluidos Computacional (CFD) busca superar as limitações dos métodos clássicos. Modelos de Gás em Rede Quântica (QLG) oferecem uma estrutura comum onde partículas ou funções de densidade de probabilidade evoluem através de etapas de colisão e streaming em uma rede discreta.
O modelo específico analisado neste trabalho, proposto originalmente por Yepez [25], é uma formulação unitária de QLG que reproduz comportamento hidrodinâmico sem codificar um modelo clássico, garantindo a unitariedade do operador de colisão. No entanto, existem questões em aberto:

• A precisão da viscosidade prevista na formulação original.
• A extensão do modelo para duas dimensões espaciais (2D) mantendo a simplicidade e a tratabilidade analítica.
• A capacidade de simular equações não lineares (como as de Burgers) com parâmetros hidrodinâmicos ajustáveis e viscosidades arbitrariamente baixas, algo difícil em métodos clássicos como o Lattice Boltzmann (LBM) devido a restrições de estabilidade numérica.

2. Metodologia

Os autores desenvolveram uma abordagem analítica e numérica baseada em três pilares principais:

• Derivação Analítica da Viscosidade (1D):

  • Revisaram o modelo 1D (Q-D1Q2) com dois qubits por sítio.
  • Realizaram uma expansão de Taylor e uma expansão de Chapman-Enskog até a segunda ordem no número de Knudsen ($\epsilon$).
  • Derivaram uma correção analítica para a viscosidade, refinando a expressão original de Yepez. A nova fórmula considera termos de ordem superior que foram negligenciados anteriormente, resultando em uma equação válida até $O(\epsilon^4)$ no regime de baixo número de Mach.
  • Demonstraram que o modelo simula a equação de Burgers: $\partial_t w + w\partial_x w = \nu \partial_{xx} w$, onde a velocidade macroscópica $w$ e a viscosidade $\nu$ dependem dos ângulos de colisão ($\theta, \zeta, \xi$).

• Generalização para 2D (Q-D2Q2):

  • Propuseram uma generalização mínima para 2D, mantendo apenas dois qubits por sítio e duas velocidades discretas.
  • Utilizaram uma rede quadrada (e discutiram outras redes de Bravais), permitindo que as populações fluam nas direções $x$ e $y$.
  • Derivaram a Equação Diferencial Parcial (PDE) geral para a densidade de massa $\rho$. O resultado é uma equação do tipo Burgers anisotrópica com um termo advectivo adicional e um tensor de difusão anisotrópico, dependendo da escolha dos vetores de velocidade.
  • A abordagem permite obter resultados analíticos fechados sem aproximações adicionais complexas, algo raro em esquemas de 2D com poucos qubits.

• Simulações Numéricas:

  • Implementaram o algoritmo para validar as previsões analíticas.
  • Calcularam a viscosidade experimental comparando com a solução analítica de Cole-Hopf para a equação de Burgers.
  • Testaram diferentes configurações de velocidade em 2D para observar a formação de choques e a anisotropia.

3. Contribuições Principais

1. Correção da Viscosidade: O trabalho fornece uma derivação rigorosa que corrige a expressão da viscosidade cinemática no modelo de Yepez. A nova fórmula (Eq. 22 no artigo) difere da original, especialmente em regimes de baixa viscosidade, e foi validada numericamente.
2. Modelo 2D Mínimo e Anisotrópico: Introduziram um modelo QLG 2D que utiliza apenas dois qubits por sítio. Embora isso leve a equações anisotrópicas (dependendo da escolha dos vetores de velocidade), o modelo é capaz de capturar a dinâmica não linear de Burgers em 2D, servindo como um passo intermediário crucial antes de modelos mais complexos (como FHP) que exigem mais recursos quânticos.
3. Viscosidade Arbitrariamente Baixa e Estabilidade: O modelo demonstra a capacidade de atingir viscosidades arbitrariamente pequenas sem perder a estabilidade numérica (graças à estrutura unitária e ao teorema H), ao contrário do LBM clássico de tempo único de relaxação, que se torna instável quando o tempo de relaxação se aproxima de 0.5. O custo é um tempo de relaxação mais longo (simulações mais longas) para atingir o equilíbrio.
4. Tratabilidade Analítica: A capacidade de obter soluções fechadas para as distribuições de equilíbrio e para a equação macroscópica derivada facilita a análise e a implementação em hardware quântico real.

4. Resultados

• Validação 1D: As simulações numéricas confirmaram a correção da viscosidade. A comparação entre a solução analítica com a viscosidade corrigida e os dados da simulação QLG mostrou um erro quadrático médio significativamente menor do que com a viscosidade original de Yepez (ver Fig. 7).
• Regime de Baixa Viscosidade: Foi demonstrado que, ao ajustar o ângulo de colisão $\theta$, é possível reduzir a viscosidade efetiva. No entanto, para $\theta$ muito próximo de 0 ou $\pi/2$, a termalização é lenta, exigindo mais passos de tempo para a formação do choque, o que explica discrepâncias em simulações de curta duração.
• Dinâmica 2D: As simulações em 2D mostraram a formação de choques em diferentes configurações de velocidade.
  • Configurações simétricas (ex: velocidades opostas em eixos ortogonais) resultaram em dinâmicas essencialmente 1D.
  • Configurações assimétricas (ex: uma velocidade em $x$, outra em $y$) geraram equações com termos de advecção cruzada e difusão anisotrópica, mas mantiveram o comportamento não linear de Burgers.
  • A comparação com o método de Diferenças Finitas (FDM) clássico mostrou boa concordância qualitativa e numérica até o momento da formação do choque, onde o FDM clássico tendeu a instabilizar, enquanto o QLG permaneceu estável.

5. Significado e Perspectivas Futuras

Este trabalho estabelece uma ponte importante entre a computação quântica nativa e a simulação de fluidos.

• Eficiência de Recursos: Ao demonstrar que é possível obter dinâmica não linear 2D com apenas dois qubits por sítio, o modelo oferece uma rota promissora para simulações em hardware quântico de escala intermediária (NISQ), onde o número de qubits é limitado.
• Paradigma Alternativo: O modelo foge da necessidade de codificar modelos clássicos (como FHP) em qubits, propondo um "gás em rede quântico nativo" com parâmetros hidrodinâmicos ajustáveis (viscosidade, velocidade do som) via ângulos de rotação quântica.
• Futuro: Os autores sugerem que o próximo passo lógico é introduzir um terceiro vetor de velocidade para incorporar a conservação de momento e, consequentemente, aproximar-se das equações de Navier-Stokes completas. Além disso, a exploração de esquemas de streaming puramente quânticos (sem medição clássica intermediária) é vista como uma direção crucial para realizar vantagens quânticas genuínas.

Em suma, o artigo refina a teoria de um modelo QLG existente, corrige seus parâmetros fundamentais e expande sua aplicabilidade para duas dimensões, provando ser uma ferramenta robusta e estável para simular dinâmica de fluidos não linear em futuros computadores quânticos.