Quantum Simulation of Stokes Flow via Schrödingerisation and Artificial Compressibility

Resumo Técnico: Simulação Quântica de Fluxo de Stokes via Schrödingerização e Compressibilidade Artificial

Definição do Problema
A simulação de fluxo de Stokes incompressível é crítica para aplicações em microfluídica, hidrodinâmica de baixo número de Reynolds e transporte em meios porosos. No entanto, os métodos numéricos clássicos enfrentam um gargalo fundamental: a estrutura de ponto de sela do sistema de Stokes exige a satisfação exata da restrição de incompressibilidade ($\nabla \cdot u = 0$) no nível discreto. Esse acoplamento leva a sistemas lineares grandes e mal condicionados que são computacionalmente caros de resolver, com requisitos de recursos que escalam de forma proibitiva conforme a dimensionalidade do problema aumenta. Embora técnicas de regularização clássicas, como a compressibilidade artificial, existam, elas permanecem computacionalmente proibitivas em configurações de alta dimensão devido ao escalonamento exponencial dos recursos necessários.

Metodologia
Os autores propõem um framework algorítmico quântico unificado que aborda esses desafios integrando compressibilidade artificial com a técnica de Schrödingerização.

1. Regularização por Compressibilidade Artificial: A condição estrita de incompressibilidade é relaxada introduzindo uma equação de evolução de pseudo-pressão. Especificamente, a equação de continuidade $\nabla \cdot u = 0$ é substituída por $\nabla \cdot u = \varepsilon p$, onde $\varepsilon$ é um pequeno parâmetro de compressibilidade artificial. Isso transforma o problema original de ponto de sela em um problema de valor inicial acoplado, envolvendo a equação de momento e uma equação de evolução de pressão parabólica. Os autores fornecem estimativas de erro demonstrando convergência linear ($O(\varepsilon)$) da solução regularizada para a verdadeira solução incompressível.

2. Discretização em Grade Deslocada (Staggered-Grid): Para garantir a estabilidade numérica e o acoplamento adequado entre velocidade e pressão, o sistema é discretizado usando um método de volume finito em grade deslocada. Esta abordagem posiciona os componentes de velocidade e a pressão em diferentes pontos nodais, evitando oscilações espúrias comuns em grades colocalizadas.

3. Mapeamento por Schrödingerização: A inovação central reside no mapeamento das equações diferenciais parciais (EDPs) não-Hermitianas e dissipativas resultantes em um sistema quântico unitário. Ao empregar uma "transformação de fase deformada" ($w = e^{-q}u$) e uma transformada de Fourier sobre a variável auxiliar $q$, os autores incorporam a dinâmica original em um espaço de Hilbert de maior dimensão. Este embedding converte o sistema dissipativo em uma equação do tipo Schrödinger ($\partial_t \hat{v} = i H_{Stokes} \hat{v}$) com um Hamiltoniano Hermitiano $H_{Stokes}$.

4. Construção do Circuito Quântico: Os autores projetam um circuito quântico explícito para simular o operador de evolução temporal $U_{Stokes}(\tau) = \exp(i H_{Stokes} \tau)$.

   • O Hamiltoniano é decomposto usando a fórmula de Lie-Trotter-Suzuki de primeira ordem em componentes elementares.
   • Operadores de diferença finita (deslocamento e Laplaciano) são representados usando codificações binárias e operadores de deslocamento ($S_+, S_-$).
   • O circuito utiliza Transformadas de Fourier Quânticas (QFT) para alternar entre o espaço auxiliar $q$ e o espaço de momento $\eta$, permitindo a simulação da dinâmica original não-unitária via evolução quântica unitária.
   • A solução é recuperada via QFT inversa e medição projetiva no subespaço da variável auxiliar onde $q > 0$.

Principais Contribuições

• Framework Unificado: O artigo apresenta o primeiro algoritmo quântico explícito para as equações de Stokes que combina compressibilidade artificial com Schrödingerização, fornecendo uma formulação coerente para o sistema acoplado de velocidade-pressão.
• Design de Circuito Explícito: Os autores derivam uma implementação detalhada de circuito quântico, incluindo a decomposição do Hamiltoniano em portas elementares (CNOT, Hadamard, Fase e portas de rotação multi-controladas) e fornecem uma análise rigorosa da contagem de portas.
• Análise de Complexidade: Uma análise de complexidade rigorosa estabelece uma vantagem quântica teórica. Os autores demonstram que, para problemas de alta dimensão, o método proposto alcança um aceleração exponencial na dimensionalidade do problema em comparação com métodos clássicos. Especificamente, enquanto a complexidade clássica escala como $O(\delta^{-(d+4)})$ sob certas suposições, a abordagem quântica evita essa dependência exponencial da dimensão $d$.
• Validação Numérica: A validade e escalabilidade do algoritmo são corroboradas através de simulações numéricas realizadas na plataforma Qiskit, incluindo testes de precisão em problemas de referência e análises de convergência de parâmetros.

Resultados
Simulações numéricas no Qiskit validam a precisão do algoritmo. Os autores testaram o método em um sistema de Stokes acoplado com uma solução analítica, demonstrando que o erro pontual diminui conforme os parâmetros de discretização ($n_x, n_q$) aumentam e o passo de tempo ($dt$) e o parâmetro de compressibilidade artificial ($\varepsilon$) diminuem. Os resultados confirmam a correção da discretização do operador diferencial e a estabilidade geral do design do circuito quântico.

Significância e Alegações
O artigo afirma oferecer um avanço significativo na dinâmica de fluidos computacional ao superar o gargalo de dimensionalidade inerente às simulações clássicas de fluxo de Stokes. Ao converter um problema de ponto de sela não-unitário em uma evolução quântica unitária, o método aproveita o potencial da computação quântica para resolver problemas de dinâmica de fluidos de alta dimensão que são atualmente intratáveis para computadores clássicos. Os autores enfatizam que sua abordagem oferece uma "aceleração exponencial provável" na dimensionalidade do problema, tornando-a uma candidata viável para aplicações futuras em geometrias complexas e hidrodinâmica de alta dimensão. O trabalho posiciona-se como um passo fundamental para portar tais algoritmos para hardware quântico de curto prazo e, eventualmente, estender o framework para as equações de Navier-Stokes não-lineares.