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

Titolo: Algoritmo di Gas su Reticolo Quantistico Bidimensionale per Equazioni di Tipo Burgers Anisotrope

1. Problema e Contesto

L'applicazione del calcolo quantistico (QC) alla fluidodinamica computazionale (CFD) mira a superare i limiti dei metodi classici. Esistono diversi approcci per tradurre i modelli di gas su reticolo classici in computer quantistici:

• Codifica della base computazionale (CBE): Assegna un qubit per ogni bit classico, ma offre vantaggi quantistici limitati.
• Codifica delle ampiezze (Amplitude Encoding): Codifica le distribuzioni di probabilità nel vettore di stato, ma soffre della non unitarietà dell'operatore di collisione, richiedendo implementazioni multi-tempo e complicando la stima delle ampiezze.
• Modelli Quantistici Nativi (QLG): Il modello proposto da Yepez [25] rappresenta un'alternativa distinta e unitaria, dove l'operatore di collisione è intrinsecamente unitario. Questo permette di incorporare fenomeni dissipativi macroscopici (come la viscosità) all'interno di una dinamica microscopica reversibile.

Tuttavia, il modello di Yepez presenta sfide aperte: la sua simulabilità classica, la validità della sua formulazione in due dimensioni spaziali (2D) e la necessità di correggere la previsione analitica della viscosità. L'obiettivo di questo lavoro è rivedere, correggere ed estendere questo modello a 2D per simulare equazioni non lineari (tipo Burgers) mantenendo un numero minimo di risorse (qubit).

2. Metodologia

Gli autori adottano un approccio ibrido che combina derivazione analitica rigorosa e simulazione numerica:

• Modello Q-D1Q2 (1D): Partendo dal modello unidimensionale di Yepez, gli autori analizzano l'evoluzione del sistema composta da collisione (operatore unitario $\hat{C}$), misurazione degli osservabili (per aggiornare le popolazioni) e streaming classico.

  • Derivano una correzione analitica alla viscosità, utilizzando un'espansione di Taylor e un'espansione di Chapman-Enskog fino al secondo ordine nel passo spaziale ($\delta x \approx \epsilon$) e primo ordine nel tempo ($\delta t \approx \epsilon^2$), nel regime a basso numero di Mach.
  • Calcolano esplicitamente il termine $J_1 - J_0$ (derivata della collisione rispetto alle popolazioni), fondamentale per determinare il coefficiente di viscosità corretto.

• Generalizzazione Q-D2Q2 (2D): Viene introdotto un modello minimale a 2 dimensioni con solo due qubit per sito e due velocità discrete.

  • A differenza dei modelli classici (come FHP) che richiedono più velocità per la conservazione del momento, questo modello mantiene la struttura a due velocità per preservare la tracciabilità analitica.
  • Si assume un reticolo quadrato (o generico reticolo di Bravais) dove le popolazioni possono propagare sia in $x$ che in $y$.
  • Si deriva l'equazione differenziale alle derivate parziali (PDE) risultante, che descrive la dinamica della densità di massa $\rho$.

• Validazione Numerica:

  • Simulazioni 1D per confrontare la viscosità predetta analiticamente con quella sperimentale estratta dai dati.
  • Confronto con la soluzione analitica dell'equazione di Burgers (tramite trasformazione di Cole-Hopf) per verificare l'accuratezza della nuova formula di viscosità.
  • Simulazioni 2D confrontate con un metodo alle differenze finite classico (FDM) per validare la dinamica anisotropa.

3. Contributi Chiave

1. Correzione della Viscosità: Gli autori dimostrano che la formula della viscosità proposta in [25] è un'approssimazione valida solo a ordini inferiori. Derivano una nuova espressione esatta (Eq. 22 nel testo) che include termini di ordine superiore, migliorando significativamente l'accordo tra teoria e simulazione, specialmente nel regime di bassa viscosità.
   $$\nu = -\frac{\delta x^2}{\delta t} \frac{1}{2} \left( 1 - \frac{1}{\sin^2(\theta)\sqrt{\alpha^2 + 1}} \right)$$
   dove $\theta$ è un angolo di collisione e $\alpha$ dipende dagli angoli di Eulero.

2. Generalizzazione Minimale 2D: Viene proposto un modello 2D che utilizza solo due velocità (due qubit). Questo porta a equazioni di tipo Burgers anisotrope. La forma generale dell'equazione derivata è:
   $$\partial_t \rho + \mathbf{a} \cdot \nabla \rho + \mathbf{b} \cdot [(1-\rho)\nabla \rho] - \nabla \cdot (\mathbf{D} \nabla \rho) = 0$$
   dove $\mathbf{D}$ è un tensore di diffusione anisotropa. A seconda della scelta dei vettori di velocità, si possono ottenere termini di advezione aggiuntivi o termini incrociati ($\partial_{xy}$).

3. Stabilità e Viscosità Arbitrariamente Piccola: Il modello QLG nativo soddisfa un teorema H, garantendo stabilità incondizionata. A differenza dei modelli LBM classici (Single-Relaxation-Time) che diventano instabili per viscosità troppo basse (quando il tempo di rilassamento $\tau \to 0.5$), il modello QLG permette di ottenere viscosità arbitrariamente piccole controllando l'angolo di collisione $\theta$, a costo di tempi di rilassamento più lunghi.

4. Risultati

• Conferma della Correzione Analitica: Le simulazioni numeriche mostrano che i dati sperimentali seguono la nuova formula di viscosità (Eq. 22) con un errore molto inferiore rispetto alla formula originale di Yepez (Eq. 19), specialmente per valori di $\theta$ lontani dai limiti estremi.
• Accordo con la Soluzione di Burgers: Il confronto con la soluzione analitica dell'equazione di Burgers (tramite la trasformazione di Cole-Hopf) conferma che il modello QLG riproduce correttamente la fisica del sistema, con un errore quadratico medio (L2) inferiore rispetto alla previsione precedente.
• Dinamica 2D: Le simulazioni 2D mostrano la formazione di shock e la propagazione di onde non lineari. L'analisi conferma che, con due sole velocità, la dinamica è intrinsecamente anisotropa. Tuttavia, scegliendo opportunamente i vettori di velocità (es. uno orizzontale e uno verticale), è possibile eliminare i termini di derivata mista, ottenendo un'equazione di Burgers con un termine di advezione aggiuntivo.
• Stabilità vs. Risoluzione: È stato osservato che viscosità molto piccole richiedono tempi di simulazione più lunghi per raggiungere l'equilibrio termico (thermalization), ma il modello rimane numericamente stabile anche quando i metodi classici (FDM) divergono alla formazione dello shock.

5. Significato e Prospettive

Questo lavoro stabilisce un ponte fondamentale tra la meccanica quantistica e la fluidodinamica classica:

• Efficienza delle Risorse: Dimostra che è possibile catturare dinamiche non lineari complesse (Burgers) e potenzialmente Navier-Stokes con un numero minimo di qubit (2 per sito), rendendo il modello promettente per l'esecuzione su hardware quantistico attuale (NISQ).
• Flessibilità dei Parametri: Offre un controllo diretto e continuo sulla viscosità tramite parametri di collisione, una caratteristica difficile da ottenere in modo stabile nei modelli LBM classici.
• Percorso verso Navier-Stokes: Sebbene il modello attuale sia anisotropo e non conservi il momento in tutte le direzioni con solo due velocità, la struttura analitica aperta suggerisce che l'aggiunta di una terza velocità (un terzo qubit) potrebbe permettere di recuperare la conservazione del momento e avvicinarsi alla dinamica di Navier-Stokes completa.
• Simulazione Nativa: Il modello rappresenta un metodo numerico "quantistico nativo", che potrebbe offrire vantaggi computazionali reali rispetto alle simulazioni classiche, specialmente per problemi di fluidodinamica turbolenta o ad alta viscosità controllata.

In conclusione, l'articolo fornisce una base teorica solida e correzioni analitiche essenziali per l'uso pratico degli algoritmi di gas su reticolo quantistico, aprendo la strada a simulazioni fluidodinamiche più efficienti e scalabili su computer quantistici.