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

Resumen Técnico: Simulación Cuántica de Flujo de Stokes mediante Schrödingerización y Compresibilidad Artificial

Planteamiento del Problema
La simulación de flujo de Stokes incompresible es crítica para aplicaciones en microfluídica, hidrodinámica de bajo número de Reynolds y transporte en medios porosos. Sin embargo, los métodos numéricos clásicos enfrentan un cuello de botella fundamental: la estructura de punto de silla del sistema de Stokes requiere la satisfacción exacta de la restricción de incompresibilidad ($\nabla \cdot u = 0$) a nivel discreto. Este acoplamiento conduce a sistemas lineales de gran tamaño y mal condicionados que son computacionalmente costosos de resolver, con requisitos de recursos que escalan de manera prohibitiva a medida que aumenta la dimensionalidad del problema. Si bien las técnicas de regularización clásica como la compresibilidad artificial existen, estas siguen siendo computacionalmente prohibitivas en entornos de alta dimensionalidad debido al escalamiento exponencial de los recursos requeridos.

Metodología
Los autores proponen un marco algorítmico cuántico unificado que aborda estos desafíos integrando la compresibilidad artificial con la técnica de Schrödingerización.

1. Regularización por Compresibilidad Artificial: Se relaja la estricta condición de incompresibilidad introduciendo una ecuación de evolución de pseudo-presión. Específicamente, la ecuación de continuidad $\nabla \cdot u = 0$ se reemplaza con $\nabla \cdot u = \varepsilon p$, donde $\varepsilon$ es un pequeño parámetro de compresibilidad artificial. Esto transforma el problema original de punto de silla en un problema de valor inicial acoplado que involucra la ecuación de cantidad de movimiento y una ecuación de evolución de presión parabólica. Los autores proporcionan estimaciones de error que demuestran una convergencia lineal ($O(\varepsilon)$) de la solución regularizada hacia la verdadera solución incompresible.

2. Discretización en Rejilla Escalonada (Staggered-Grid): Para asegurar la estabilidad numérica y un acoplamiento adecuado entre velocidad y presión, el sistema se discretiza utilizando un método de volúmenes finitos en rejilla escalonada. Este enfoque sitúa los componentes de velocidad y la presión en diferentes puntos nodales, evitando oscilaciones espurias comunes en las rejillas colocalizadas.

3. Mapeo de Schrödingerización: La innovación central reside en el mapeo de las ecuaciones diferenciales parciales (PDE) no hermíticas y disipativas resultantes en un sistema cuántico unitario. Mediante el empleo de una "transformación de fase deformada" ($w = e^{-q}u$) y una transformada de Fourier sobre la variable auxiliar $q$, los autores embeben la dinámica original en un espacio de Hilbert de mayor dimensión. Este embebido convierte el sistema disipativo en una ecuación de tipo Schrödinger ($\partial_t \hat{v} = i H_{Stokes} \hat{v}$) con un Hamiltoniano de Stokes hermítico $H_{Stokes}$.

4. Construcción del Circuito Cuántico: Los autores diseñan un circuito cuántico explícito para simular el operador de evolución temporal $U_{Stokes}(\tau) = \exp(i H_{Stokes} \tau)$.

   • El Hamiltoniano se descompone utilizando la fórmula de Lie-Trotter-Suzuki de primer orden en componentes elementales.
   • Los operadores de diferencia finita (desplazamiento y Laplaciano) se representan mediante codificaciones binarias y operadores de desplazamiento ($S_+, S_-$).
   • El circuito utiliza Transformadas de Fourier Cuánticas (QFT) para cambiar entre el espacio auxiliar $q$ y el espacio de momento $\eta$, permitiendo la simulación de la dinámica no unitaria original mediante evolución cuántica unitaria.
   • La solución se recupera mediante la QFT inversa y la medición proyectiva en el subespacio de la variable auxiliar donde $q > 0$.

Contribuciones Clave

• Marco Unificado: El artículo presenta el primer algoritmo cuántico explícito para las ecuaciones de Stokes que combina la compresibilidad artificial con la Schrödingerización, proporcionando una formulación coherente para el sistema acoplado de velocidad-presión.
• Diseño de Circuito Explícito: Los autores derivan una implementación detallada del circuito cuántico, incluyendo la descomposición del Hamiltoniano en puertas elementales (CNOT, Hadamard, Fase y puertas de rotación multicontroladas) y proporcionan un análisis riguroso del conteo de puertas.
• Análisis de Complejidad: Un análisis de complejidad riguroso establece una ventaja cuántica teórica. Los autores demuestran que, para problemas de alta dimensión, el método propuesto logra una aceleración exponencial en la dimensionalidad del problema en comparación con los métodos clásicos. Específicamente, mientras que la complejidad clásica escala como $O(\delta^{-(d+4)})$ bajo ciertas suposiciones, el enfoque cuántico evita esta dependencia exponencial de la dimensión $d$.
• Validación Numérica: La validez y escalabilidad del algoritmo se corroboran mediante simulaciones numéricas realizadas en la plataforma Qiskit, incluyendo pruebas de precisión en problemas de referencia y análisis de convergencia de parámetros.

Resultos
Las simulaciones numéricas en Qiskit validan la precisión del algoritmo. Los autores probaron el método en un sistema de Stokes acoplado con una solución analítica, demostrando que el error puntual disminuye a medida que aumentan los parámetros de discretización ($n_x, n_q$) y disminuyen el paso de tiempo ($dt$) y el parámetro de compresibilidad artificial ($\varepsilon$). Los resultados confirman la corrección de la discretización del operador diferencial y la estabilidad general del diseño del circuito cuántico.

Significancia y Reivindicaciones
El artículo afirma ofrecer un avance significativo en la dinámica de fluidos computacional al superar el cuello de botella de la dimensionalidad inherente a las simulaciones clásicas de flujo de Stokes. Al convertir un problema de punto de silla no unitario en una evolución cuántica unitaria, el método aprovecha el potencial de la computación cuántica para resolver problemas de dinámica de fluidos de alta dimensión que son actualmente intratables para las computadoras clásicas. Los autores enfatizan que su enfoque proporciona una "aceleración exponencial demostrable" en la dimensionalidad del problema, convirtiéndolo en un candidato viable para futuras aplicaciones en geometrías complejas e hidrodinámica de alta dimensión. El trabajo se posiciona como un paso fundacional hacia el traslado de tales algoritmos a hardware cuántico de corto plazo y, eventualmente, hacia la extensión del marco a las ecuaciones de Navier-Stokes no lineales.