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:2603.00558

A Stable and General Quantum Fractional-Step Lattice Boltzmann Method for Incompressible Flows

Yang Xiao, Liming Yang, Chang Shu, Yinjie Du

quant-phphysics.flu-dynSubmitted Feb 28, 2026• Announced Mar 3, 2026• 0 citations
arXivHTMLPDF

Abstract

Quantum computing shows substantial potential in accelerating simulations and alleviating memory bottlenecks in computational fluid dynamics (CFD), owing to its inherent properties of superposition and entanglement. The lattice Boltzmann method (LBM), being largely algebraic in nature, has inspired the development of various quantum LBMs. However, most existing approaches fix the relaxation time at $τ$ = 1, thereby confining a given mesh resolution to simulations at a single Reynolds number. Although our earlier quantum lattice kinetic scheme (LKS) lifted this restriction, it suffers from instability at high Reynolds numbers. To address this challenge, we propose a quantum fractional-step LBM (FS-LBM). In this framework, the predictor step is implemented on a quantum circuit using the standard LBM formulation, while the corrector step is performed classically. The relaxation time is retained at $τ$ = 1 to ensure seamless compatibility with existing quantum LBMs. Benchmark simulations of representative two- and three-dimensional incompressible isothermal and thermal flows demonstrate that the quantum FS-LBM achieves accuracy and convergence orders consistent with its classical counterpart, while significantly outperforming the quantum LKS in both precision and stability. Notably, this work presents the first quantum LBM simulation of three-dimensional incompressible thermal flows.

Comments
51 pages, 18 figures

Related papers

Quantum Taylor series expansion and least squares-based fractional-step LBM for incompressible flows on non-uniform meshes

Liming Yang, Yang Xiao, Chang Shu +1 more

In this work, a quantum Taylor series expansion and Least squares-based Fractional-step Lattice Boltzmann Method (TLFS-LBM) is proposed to simulate incompressible flows on non-uniform meshes. Most existing quantum LBMs face two critical challenges: with the relaxation time τ fixed at 1, a fixed mesh resolution corresponds to only a single Reynolds number, and the streaming step, typically implemented through the quantum walk method, is limited to uniform meshes. To address these issues, we retain τ = 1 and decompose the quantum TLFS-LBM into predictor and corrector steps using a fractional-step method. The predictor step is implemented on a quantum circuit, but the resulting viscosity generally differs from the physical viscosity. A diffusion equation is therefore solved classically in the corrector step to restore the correct physical viscosity. Moreover, the streaming step is expressed in the form of A · f eq using Taylor series expansion and least-squares approximation, with singular value decomposition applied to factorize matrix A into unitary matrices for quantum implementation. Since A depends solely on mesh point coordinates, the proposed method is applicable to arbitrary mesh types. Validation through two-dimensional incompressible isothermal and thermal flow simulations confirms that the quantum TLFS-LBM achieves accuracy comparable to its classical counterpart, while extending the applicability of quantum LBM to non-uniform meshes and curved boundaries.

Chinese Journal of Aeronautics•Mar 26, 2026•3 citations

Simulating non-trivial incompressible flows with a quantum lattice Boltzmann algorithm

David Jennings, Kamil Korzekwa, Matteo Lostaglio +4 more

Quantum algorithms have been identified as a potential means to accelerate computational fluid dynamics (CFD) simulations, with the lattice Boltzmann method (LBM) being a promising candidate for realizing quantum speedups. Here, we extend the recent quantum algorithm for the incompressible LBM to account for realistic fluid dynamics setups by incorporating walls, inlets, outlets, and external forcing. We analyze the associated complexity cost and show that these modifications preserve the asymptotic scaling, and potential quantum advantage, of the original algorithm. Moreover, to support our theoretical analysis, we provide a classical numerical study illustrating the accuracy, complexity, and convergence of the algorithm for representative incompressible-flow cases, including the driven Taylor-Green vortex, the lid-driven cavity flow, and the flow past a cylinder. Our results provide a pathway to accurate quantum simulation of nonlinear fluid dynamics, and a framework for extending quantum LBM to more challenging flow configurations.

physics.flu-dynCross-listedarXiv:2512.05781•AIAA SciTech Forum•Dec 8, 2025•0 citations

Quantum parallelism meets lattice Boltzmann: A hybrid efficient framework for advection–diffusion simulation

Ruijie Xu, Gang Du

Quantum computing offers a promising pathway to accelerate fluid dynamics simulations. In this work, we propose a hybrid quantum-classical lattice Boltzmann method (LBM) for solving the advection–diffusion equation, as a step toward an efficient and general quantum linear solver framework for fluid flows. Quantum circuits handle the computationally intensive collision and streaming steps, exploiting quantum parallelism and removing parameter constraints, while boundary corrections are managed classically to avoid unnecessary overhead. Various boundary conditions are incorporated within the quantum framework, enhancing the flexibility and applicability of the method. Numerical experiments with Gaussian hill initial conditions under periodic boundaries and zero initial conditions under mixed Dirichlet–Neumann boundaries, each in one and two dimensions, demonstrate the precision of the proposed approach. The computational complexity is analyzed in detail. The collision and streaming steps achieve a lower complexity than their classical LBM counterparts, although overall performance is still limited by the input and output costs. This work represents a critical step toward the practical quantum simulation of complex fluid flows.

Physics of Fluids•May 1, 2026•0 citations

Quantum algorithm for the lattice Boltzmann method with applications on real quantum devices

Antonio Bastida-Zamora, Ljubomir Budinski, Oskari Kerppo +8 more

We introduce a novel quantum algorithm for the lattice Boltzmann method (LBM) based on the one-step simplified LBM. The structure of the algorithm allows for more flexibility in modelling different physics in contrast to earlier quantum algorithms for the LBM, while retaining computational efficiency in terms of the gate and qubit complexity. The new algorithm has potential for full end-to-end quantum utility especially for linear problems. We discuss the implementation of examples in linear acoustics, as well as a nonlinear Navier-Stokes problem that was solved on an IBM QPU in a hybrid simulation loop.

quant-pharXiv:2603.02127•Mar 3, 2026•0 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

A quantum computing-based lattice Boltzmann method with a linearized non-equilibrium collision operator and modular circuit for practical flow simulation

Kang-Yang Zeng, Xiaodong Niu, Adnan Khan +2 more

Quantum computing (QC) has emerged as a potential computational tool in fluid dynamics. However, a practical QC protocol with end-to-end practical applications to fluid flow is still in nascent. We here propose a simple and practical hybrid QC algorithm based on lattice Boltzmann methods (LBM), which we term HQLBM, for fluid dynamics simulations. The existing QC-based LBMs are limited to simple low-Reynolds-number flow simulations due to their disability to handle the nonlinear collision term in LBM for QC. In this work, a linearized non-equilibrium collision, derived from kinetic theory, is presented for HQLBM simulating practical flows. The present approach ensures the unitary of quantum algorithms while keeping the collision relaxation parameter adjustable for simulating different Reynolds number flows. Most importantly, the linearized collision keeps the accuracy of the original nonlinear one in the lattice Boltzmann equation when it theoretically recovers its macroscopic counterpart—Navier–Stokes equations. The practicality of the proposed method is demonstrated by simulating two typical flows, including lid-driven and natural convection flows in a square cavity at different Reynolds and Rayleigh numbers, respectively. It is shown that the proposed HQLBM has good numerical accuracy compared to the conventional LBM, and lower complexity on comparisons to the QC-based LBM carried out on single-circuit QC even when it is performed on a Qiskit virtual quantum computer. This work offers a practical application of QC-based LBM for complex fluid dynamics problems.

Physics of Fluids•Aug 1, 2025•6 citations
Quantum computing shows substantial potential in accelerating simulations and alleviating memory bottlenecks in computational fluid dynamics (CFD), owing to its inherent properties of superposition and entanglement. The lattice Boltzmann method (LBM), being largely algebraic in nature, has inspired the development of various quantum LBMs. However, most existing approaches fix the relaxation time at τττ = 1, thereby confining a given mesh resolution to simulations at a single Reynolds number. Although our earlier quantum lattice kinetic scheme (LKS) lifted this restriction, it suffers from instability at high Reynolds numbers. To address this challenge, we propose a quantum fractional-step LBM (FS-LBM). In this framework, the predictor step is implemented on a quantum circuit using the standard LBM formulation, while the corrector step is performed classically. The relaxation time is retained at τττ = 1 to ensure seamless compatibility with existing quantum LBMs. Benchmark simulations of representative two- and three-dimensional incompressible isothermal and thermal flows demonstrate that the quantum FS-LBM achieves accuracy and convergence orders consistent with its classical counterpart, while significantly outperforming the quantum LKS in both precision and stability. Notably, this work presents the first quantum LBM simulation of three-dimensional incompressible thermal flows.
In this work, a quantum Taylor series expansion and Least squares-based Fractional-step Lattice Boltzmann Method (TLFS-LBM) is proposed to simulate incompressible flows on non-uniform meshes. Most existing quantum LBMs face two critical challenges: with the relaxation time τ fixed at 1, a fixed mesh resolution corresponds to only a single Reynolds number, and the streaming step, typically implemented through the quantum walk method, is limited to uniform meshes. To address these issues, we retain τ = 1 and decompose the quantum TLFS-LBM into predictor and corrector steps using a fractional-step method. The predictor step is implemented on a quantum circuit, but the resulting viscosity generally differs from the physical viscosity. A diffusion equation is therefore solved classically in the corrector step to restore the correct physical viscosity. Moreover, the streaming step is expressed in the form of A · f eq using Taylor series expansion and least-squares approximation, with singular value decomposition applied to factorize matrix A into unitary matrices for quantum implementation. Since A depends solely on mesh point coordinates, the proposed method is applicable to arbitrary mesh types. Validation through two-dimensional incompressible isothermal and thermal flow simulations confirms that the quantum TLFS-LBM achieves accuracy comparable to its classical counterpart, while extending the applicability of quantum LBM to non-uniform meshes and curved boundaries.
Quantum algorithms have been identified as a potential means to accelerate computational fluid dynamics (CFD) simulations, with the lattice Boltzmann method (LBM) being a promising candidate for realizing quantum speedups. Here, we extend the recent quantum algorithm for the incompressible LBM to account for realistic fluid dynamics setups by incorporating walls, inlets, outlets, and external forcing. We analyze the associated complexity cost and show that these modifications preserve the asymptotic scaling, and potential quantum advantage, of the original algorithm. Moreover, to support our theoretical analysis, we provide a classical numerical study illustrating the accuracy, complexity, and convergence of the algorithm for representative incompressible-flow cases, including the driven Taylor-Green vortex, the lid-driven cavity flow, and the flow past a cylinder. Our results provide a pathway to accurate quantum simulation of nonlinear fluid dynamics, and a framework for extending quantum LBM to more challenging flow configurations.
Quantum computing offers a promising pathway to accelerate fluid dynamics simulations. In this work, we propose a hybrid quantum-classical lattice Boltzmann method (LBM) for solving the advection–diffusion equation, as a step toward an efficient and general quantum linear solver framework for fluid flows. Quantum circuits handle the computationally intensive collision and streaming steps, exploiting quantum parallelism and removing parameter constraints, while boundary corrections are managed classically to avoid unnecessary overhead. Various boundary conditions are incorporated within the quantum framework, enhancing the flexibility and applicability of the method. Numerical experiments with Gaussian hill initial conditions under periodic boundaries and zero initial conditions under mixed Dirichlet–Neumann boundaries, each in one and two dimensions, demonstrate the precision of the proposed approach. The computational complexity is analyzed in detail. The collision and streaming steps achieve a lower complexity than their classical LBM counterparts, although overall performance is still limited by the input and output costs. This work represents a critical step toward the practical quantum simulation of complex fluid flows.
We introduce a novel quantum algorithm for the lattice Boltzmann method (LBM) based on the one-step simplified LBM. The structure of the algorithm allows for more flexibility in modelling different physics in contrast to earlier quantum algorithms for the LBM, while retaining computational efficiency in terms of the gate and qubit complexity. The new algorithm has potential for full end-to-end quantum utility especially for linear problems. We discuss the implementation of examples in linear acoustics, as well as a nonlinear Navier-Stokes problem that was solved on an IBM QPU in a hybrid simulation loop.
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.
Quantum computing (QC) has emerged as a potential computational tool in fluid dynamics. However, a practical QC protocol with end-to-end practical applications to fluid flow is still in nascent. We here propose a simple and practical hybrid QC algorithm based on lattice Boltzmann methods (LBM), which we term HQLBM, for fluid dynamics simulations. The existing QC-based LBMs are limited to simple low-Reynolds-number flow simulations due to their disability to handle the nonlinear collision term in LBM for QC. In this work, a linearized non-equilibrium collision, derived from kinetic theory, is presented for HQLBM simulating practical flows. The present approach ensures the unitary of quantum algorithms while keeping the collision relaxation parameter adjustable for simulating different Reynolds number flows. Most importantly, the linearized collision keeps the accuracy of the original nonlinear one in the lattice Boltzmann equation when it theoretically recovers its macroscopic counterpart—Navier–Stokes equations. The practicality of the proposed method is demonstrated by simulating two typical flows, including lid-driven and natural convection flows in a square cavity at different Reynolds and Rayleigh numbers, respectively. It is shown that the proposed HQLBM has good numerical accuracy compared to the conventional LBM, and lower complexity on comparisons to the QC-based LBM carried out on single-circuit QC even when it is performed on a Qiskit virtual quantum computer. This work offers a practical application of QC-based LBM for complex fluid dynamics problems.