Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2310.17973 · AVS Quantum Science

Lattice Boltzmann–Carleman quantum algorithm and circuit for fluid flows at moderate Reynolds number

Claudio Sanavio, Sauro Succi A

Published Apr 22, 2024• 40 citations
arXivHTMLPDF

Abstract

We present a quantum computing algorithm for fluid flows based on the Carleman-linearization of the Lattice Boltzmann (LB) method. First, we demonstrate the convergence of the classical Carleman procedure at moderate Reynolds numbers, namely, for Kolmogorov-like flows. Then we proceed to formulate the corresponding quantum algorithm, including the quantum circuit layout, and analyze its computational viability. We show that, at least for moderate Reynolds numbers between 10 and 100, the Carleman–LB procedure can be successfully truncated at second order, which is a very encouraging result. We also show that the quantum circuit implementing the single time-step collision operator has a fixed depth, regardless of the number of lattice sites. However, such depth is of the order of ten thousands quantum gates, meaning that quantum advantage over classical computing is not attainable today, but could be achieved in the near or mid-term future. The same goal for the multi-step version remains, however, an open topic for future research.

DOI
10.1116/5.0195549

Related papers

Carleman-lattice-Boltzmann quantum circuit with matrix access oracles

Claudio Sanavio, William Simon, Alexis Ralli +2 more

We apply Carleman linearization of the Lattice Boltzmann (CLB) representation of fluid flows to quantum emulate the dynamics of a 2D Kolmogorov-like flow. We assess the accuracy of the result and find a relative error of the order of $10^{-3}$ with just two Carleman iterates, for a range of the Reynolds number up to a few hundreds. We first define a gate-based quantum circuit for the implementation of the CLB method and then exploit the sparse nature of the CLB matrix to build a quantum circuit based on block-encoding techniques which makes use of matrix oracles. It is shown that the gate complexity of the algorithm is thereby dramatically reduced, from exponential to quadratic. However, due to the need of employing up to seven ancilla qubits, the probability of success of the corresponding circuit for a single time step is too low to enable multi-step time evolution. Several possible directions to circumvent this problem are briefly outlined.

arXiv:2501.02582•Jan 5, 2025•0 citations

Quantum computing for simulation of fluid dynamics

Claudio Sanavio, Sauro Succi A

The implementation of quantum algorithms for the simulation of classical fluid dynamics poses a fundamental challenge due to the nonlinearity of the fluid equations. In this work, we provide a pedagogical introduction to quantum computing algorithms for simulating classical fluids, with a special focus on the Carleman-Lattice Boltzmann algorithm, which has captured significant attention in the last couple of years. While this algorithm demonstrates satisfactory convergence to analytical solutions for systems at low-to-moderate Reynolds numbers, it also shows an exponential depth of the corresponding quantum circuit. As a result much further analysis is needed to assess the availability of the Carleman-Lattice Boltzmann method on a quantum computer.

arXiv:2404.01302•IntechOpen eBooks•May 8, 2024•6 citations

A Demonstration of Quantum Circuit Implementation for Obstacle Flow Using Carleman-Linearized Lattice Boltzmann Method

Kazumasa Ueno, Keita Kanno, Yasunori Lee

Fluid simulations, especially at high Reynolds numbers, are computationally expensive on classical computers, making them promising application targets for quantum computing. Recent studies have combined the lattice Boltzmann method (LBM) with Carleman linearization to design quantum algorithms for computational fluid dynamics (CFD). However, practical quantum-circuit implementations of these algorithms that incorporate non-periodic boundary conditions have not been fully explored. In this work, we implement a quantum algorithm for two-dimensional linearized fluid flow around an obstacle, using block-encoding of the linear-system matrix and quantum singular value transformation (QSVT) to solve it. Inflow, outflow, and no-slip boundary conditions are formulated as sparse matrix operations and efficiently embedded into quantum circuits using index-value encoding. We demonstrate logarithmic scaling of the required numbers of qubits and gates with respect to the number of lattice points, suggesting the potential feasibility of quantum-computational fluid dynamics simulations.

quant-pharXiv:2605.28135•May 28, 2026•0 citations

Carleman-lattice-Boltzmann quantum circuit with matrix access oracles

Claudio Sanavio, William Simon, Alexis Ralli +2 more

We apply Carleman linearization of the Lattice Boltzmann (CLB) representation of fluid flows to quantum emulate the dynamics of a two-dimensional Kolmogorov-like flow. We assess the accuracy of the result and find a relative error of the order of 10−3 with just two Carleman iterates for a range of the Reynolds number up to a few hundreds. We first define a gate-based quantum circuit for the implementation of the CLB method and then exploit the sparse nature of the CLB matrix to build a quantum circuit based on block-encoding techniques which makes use of matrix oracles. It is shown that the gate complexity of the algorithm is thereby dramatically reduced, from exponential to quadratic. However, due to the need of employing up to seven ancilla qubits, the probability of success of the corresponding circuit for a single time step is too low to enable multi-step time evolution. Several possible directions to circumvent this problem are briefly outlined.

Physics of Fluids•Mar 1, 2025•15 citations

An end-to-end quantum algorithm for nonlinear fluid dynamics with bounded quantum advantage

David Jennings, Kamil Korzekwa, Matteo Lostaglio +3 more

Computational fluid dynamics (CFD) is a cornerstone of classical scientific computing, and there is growing interest in whether quantum computers can accelerate such simulations. The existing proposals for fault-tolerant quantum algorithms for CFD have almost exclusively been based on the Carleman embedding method, used to encode nonlinearities on a quantum computer. Here, we first identify several key challenges that must be addressed in order to establish quantum advantage using this method, including Carleman convergence, time-stepping requirements, condition-number scaling, and data extraction. Building on this analysis, we develop a novel algorithm for the incompressible lattice Boltzmann equation that addresses them, and provide its detailed analysis, including all potential sources of algorithmic complexity, as well as gate count estimates. We find that for an end-to-end problem, a modest quantum advantage may be preserved for selected observables in the high-error-tolerance regime. We establish a lower bound for the Reynolds number scaling of our quantum algorithm in dimension $D$ at Kolmogorov microscale resolution with $O(\mathrm{Re}^{\frac{3}{4}(1+\frac{D}{2})} \times q_M)$, where $q_M$ is a multiplicative overhead for data extraction with $q_M = O(\mathrm{Re}^{\frac{3}{8}})$ for the drag force. This upper bounds the scaling improvement over classical algorithms by $O(\mathrm{Re}^{\frac{3D}{8}})$. However, our numerical investigations suggest a lower speedup, with a scaling estimate of $O(\mathrm{Re}^{1.936} \times q_M)$ for $D=2$. Our results give robust evidence that small, but nontrivial, quantum advantages can be achieved in the context of CFD, and motivate the need for additional rigorous end-to-end quantum algorithm development.

quant-pharXiv:2512.03758•PRX Quantum•Dec 4, 2025•1 citation

Quantum algorithm for lattice Boltzmann (QALB) simulation of incompressible fluids with a nonlinear collision term

Wael Itani, Katepalli R. Sreenivasan, Sauro Succi A

We present a full quantum algorithm for the lattice Boltzmann method for simulating fluid flows, the only such algorithm to implement both the streaming and collision substeps as unitary operators using an efficient number of qubits. We use Hamiltonian simulation as the main route underlying the algorithm, and show that it uses qubits that scale logarithmically in the flow Reynolds number but gates that scale only polynomially. We motivate the discussion by a brief overview of existing attempts at simulating classical fluids on quantum computers and present a pedagogical discussion on assigning quantum operators to classical variables in both the streaming and the collision substeps, after highlighting the incompatibility of the latter with the implementations of the streaming step that exists in the literature. We use the Bhatnagar–Gross–Krook ansatz for the collision term, representing the relaxation toward an equilibrium distribution. For nonlinear collisions, we use Kowalski's framework that links the nonlinear dynamics of a system to the evolution of bosonic modes, assigning a Carleman linearization order to the truncation in the Fock space of the bosons. We present the qubit and gate complexities, in terms of the chosen accuracy and the Reynolds number. In the Appendix, we work out the details of implementing the operators of the truncated bosonic Fock space in terms of single-qubit gates as well as the error scaling for a general polynomial driving function.

arXiv:2304.05915•Physics of Fluids•Jan 1, 2024•54 citations
We present a quantum computing algorithm for fluid flows based on the Carleman-linearization of the Lattice Boltzmann (LB) method. First, we demonstrate the convergence of the classical Carleman procedure at moderate Reynolds numbers, namely, for Kolmogorov-like flows. Then we proceed to formulate the corresponding quantum algorithm, including the quantum circuit layout, and analyze its computational viability. We show that, at least for moderate Reynolds numbers between 10 and 100, the Carleman–LB procedure can be successfully truncated at second order, which is a very encouraging result. We also show that the quantum circuit implementing the single time-step collision operator has a fixed depth, regardless of the number of lattice sites. However, such depth is of the order of ten thousands quantum gates, meaning that quantum advantage over classical computing is not attainable today, but could be achieved in the near or mid-term future. The same goal for the multi-step version remains, however, an open topic for future research.
We apply Carleman linearization of the Lattice Boltzmann (CLB) representation of fluid flows to quantum emulate the dynamics of a 2D Kolmogorov-like flow. We assess the accuracy of the result and find a relative error of the order of 10−310^{-3}10−3 with just two Carleman iterates, for a range of the Reynolds number up to a few hundreds. We first define a gate-based quantum circuit for the implementation of the CLB method and then exploit the sparse nature of the CLB matrix to build a quantum circuit based on block-encoding techniques which makes use of matrix oracles. It is shown that the gate complexity of the algorithm is thereby dramatically reduced, from exponential to quadratic. However, due to the need of employing up to seven ancilla qubits, the probability of success of the corresponding circuit for a single time step is too low to enable multi-step time evolution. Several possible directions to circumvent this problem are briefly outlined.
The implementation of quantum algorithms for the simulation of classical fluid dynamics poses a fundamental challenge due to the nonlinearity of the fluid equations. In this work, we provide a pedagogical introduction to quantum computing algorithms for simulating classical fluids, with a special focus on the Carleman-Lattice Boltzmann algorithm, which has captured significant attention in the last couple of years. While this algorithm demonstrates satisfactory convergence to analytical solutions for systems at low-to-moderate Reynolds numbers, it also shows an exponential depth of the corresponding quantum circuit. As a result much further analysis is needed to assess the availability of the Carleman-Lattice Boltzmann method on a quantum computer.
Fluid simulations, especially at high Reynolds numbers, are computationally expensive on classical computers, making them promising application targets for quantum computing. Recent studies have combined the lattice Boltzmann method (LBM) with Carleman linearization to design quantum algorithms for computational fluid dynamics (CFD). However, practical quantum-circuit implementations of these algorithms that incorporate non-periodic boundary conditions have not been fully explored. In this work, we implement a quantum algorithm for two-dimensional linearized fluid flow around an obstacle, using block-encoding of the linear-system matrix and quantum singular value transformation (QSVT) to solve it. Inflow, outflow, and no-slip boundary conditions are formulated as sparse matrix operations and efficiently embedded into quantum circuits using index-value encoding. We demonstrate logarithmic scaling of the required numbers of qubits and gates with respect to the number of lattice points, suggesting the potential feasibility of quantum-computational fluid dynamics simulations.
We apply Carleman linearization of the Lattice Boltzmann (CLB) representation of fluid flows to quantum emulate the dynamics of a two-dimensional Kolmogorov-like flow. We assess the accuracy of the result and find a relative error of the order of 10−3 with just two Carleman iterates for a range of the Reynolds number up to a few hundreds. We first define a gate-based quantum circuit for the implementation of the CLB method and then exploit the sparse nature of the CLB matrix to build a quantum circuit based on block-encoding techniques which makes use of matrix oracles. It is shown that the gate complexity of the algorithm is thereby dramatically reduced, from exponential to quadratic. However, due to the need of employing up to seven ancilla qubits, the probability of success of the corresponding circuit for a single time step is too low to enable multi-step time evolution. Several possible directions to circumvent this problem are briefly outlined.
Computational fluid dynamics (CFD) is a cornerstone of classical scientific computing, and there is growing interest in whether quantum computers can accelerate such simulations. The existing proposals for fault-tolerant quantum algorithms for CFD have almost exclusively been based on the Carleman embedding method, used to encode nonlinearities on a quantum computer. Here, we first identify several key challenges that must be addressed in order to establish quantum advantage using this method, including Carleman convergence, time-stepping requirements, condition-number scaling, and data extraction. Building on this analysis, we develop a novel algorithm for the incompressible lattice Boltzmann equation that addresses them, and provide its detailed analysis, including all potential sources of algorithmic complexity, as well as gate count estimates. We find that for an end-to-end problem, a modest quantum advantage may be preserved for selected observables in the high-error-tolerance regime. We establish a lower bound for the Reynolds number scaling of our quantum algorithm in dimension DDD at Kolmogorov microscale resolution with O(Re34(1+D2)×qM)O(\mathrm{Re}^{\frac{3}{4}(1+\frac{D}{2})} \times q_M)O(Re43​(1+2D​)×qM​), where qMq_MqM​ is a multiplicative overhead for data extraction with qM=O(Re38)q_M = O(\mathrm{Re}^{\frac{3}{8}})qM​=O(Re83​) for the drag force. This upper bounds the scaling improvement over classical algorithms by O(Re3D8)O(\mathrm{Re}^{\frac{3D}{8}})O(Re83D​). However, our numerical investigations suggest a lower speedup, with a scaling estimate of O(Re1.936×qM)O(\mathrm{Re}^{1.936} \times q_M)O(Re1.936×qM​) for D=2D=2D=2. Our results give robust evidence that small, but nontrivial, quantum advantages can be achieved in the context of CFD, and motivate the need for additional rigorous end-to-end quantum algorithm development.
We present a full quantum algorithm for the lattice Boltzmann method for simulating fluid flows, the only such algorithm to implement both the streaming and collision substeps as unitary operators using an efficient number of qubits. We use Hamiltonian simulation as the main route underlying the algorithm, and show that it uses qubits that scale logarithmically in the flow Reynolds number but gates that scale only polynomially. We motivate the discussion by a brief overview of existing attempts at simulating classical fluids on quantum computers and present a pedagogical discussion on assigning quantum operators to classical variables in both the streaming and the collision substeps, after highlighting the incompatibility of the latter with the implementations of the streaming step that exists in the literature. We use the Bhatnagar–Gross–Krook ansatz for the collision term, representing the relaxation toward an equilibrium distribution. For nonlinear collisions, we use Kowalski's framework that links the nonlinear dynamics of a system to the evolution of bosonic modes, assigning a Carleman linearization order to the truncation in the Fock space of the bosons. We present the qubit and gate complexities, in terms of the chosen accuracy and the Reynolds number. In the Appendix, we work out the details of implementing the operators of the truncated bosonic Fock space in terms of single-qubit gates as well as the error scaling for a general polynomial driving function.