Recent advancements in quantum computing provide opportunities to address the challenges of computational resource limitations in computational fluid dynamics (CFD). This work presents a hybrid quantum-classical CFD framework that leverages quantum linear algorithms to practical flow simulations. By transforming nonlinear fluid dynamics problems into linear systems, the framework employs quantum linear algorithms to compute solutions, effectively integrating quantum and classical computing approaches. Furthermore, this framework employs subspace methods to map the original large-scale linear systems to small ones solved by quantum linear algorithms, thereby enabling the solution of large-scale problems using the currently limited quantum resources. We applied the famous Harrow–Hassidim–Lloyd (HHL) algorithm and variational quantum linear solver to simulate complex flows, including unsteady flows around a cylinder, simulation of an aircraft with turbulence model, and combustion flows. The performance and quantum resource consumption of these algorithms were evaluated in practical flow scenarios, demonstrating their effectiveness with an average relative error lower than 0.001%. We conducted quantum resource and convergence analyses to adapt this framework for use with near-term quantum computers. Our framework supports simulations of up to 9.4 × 106 grid cells, marking a significant advancement from previous quantum approaches limited to simple flows. This paper not only offers a way for employing quantum linear algorithms in solving complex flow problems but also provides insights into optimizing quantum algorithms for CFD applications, pushing forward the practical utilization of quantum CFD.
Recent advancements in quantum computing provide opportunities to address the challenges of computational resource limitations in computational fluid dynamics (CFD). This work presents a hybrid quantum-classical CFD framework that leverages quantum linear algorithms to practical flow simulations. By transforming nonlinear fluid dynamics problems into linear systems, the framework employs quantum linear algorithms to compute solutions, effectively integrating quantum and classical computing approaches. Furthermore, this framework employs subspace methods to map the original large-scale linear systems to small ones solved by quantum linear algorithms, thereby enabling the solution of large-scale problems using the currently limited quantum resources. We applied the famous Harrow–Hassidim–Lloyd (HHL) algorithm and variational quantum linear solver to simulate complex flows, including unsteady flows around a cylinder, simulation of an aircraft with turbulence model, and combustion flows. The performance and quantum resource consumption of these algorithms were evaluated in practical flow scenarios, demonstrating their effectiveness with an average relative error lower than 0.001%. We conducted quantum resource and convergence analyses to adapt this framework for use with near-term quantum computers. Our framework supports simulations of up to 9.4 × 106 grid cells, marking a significant advancement from previous quantum approaches limited to simple flows. This paper not only offers a way for employing quantum linear algorithms in solving complex flow problems but also provides insights into optimizing quantum algorithms for CFD applications, pushing forward the practical utilization of quantum CFD.
Computational Fluid Dynamics (CFD) plays a central role in mechanical system design by enabling the simulation of complex flow phenomena such as turbulence, phase transitions, and fluid–structure interactions. However, high-fidelity CFD simulations—particularly those involving turbulent and multiphase flows—are often limited by the computational cost of solving large-scale, nonlinear, and high-dimensional systems of equations. This study explores the potential of quantum-enhanced computational fluid dynamics, focusing on hybrid quantum–classical solvers as a next-generation approach to overcoming these limitations. By integrating quantum algorithms with classical numerical methods, hybrid frameworks aim to accelerate key CFD tasks such as linear system solving, optimization, and high-dimensional state exploration while preserving the stability and accuracy of established solvers. The paper examines how quantum principles such as superposition and entanglement can be leveraged to improve convergence rates and design-space exploration in turbulence modeling and multiphase flow simulations. Through a conceptual and methodological analysis, this work highlights the opportunities, current constraints, and future implications of quantum-enhanced CFD for mechanical system design, offering a forward-looking perspective on the convergence of quantum computing and fluid mechanics.
Quantum computing (QC) represents an emerging frontier with significant potential for transformative advancements in Computational Fluid Dynamics (CFD). This potential has been explicitly recognized by the CFD Vision 2030 report, emphasizing the strategic importance of exploring QC for future breakthroughs. Among quantum algorithms relevant to CFD, the Variational Quantum Linear Solver (VQLS) is particularly promising due to its suitability for solving linear and Hamiltonian-based problems, characteristic of fundamental CFD equations such as the pressure-Poisson equation. However, the practical implementation of VQLS for large-scale CFD problems faces critical scalability challenges, most notably the exponential increase in the number of quantum circuit evaluations and the barren plateau phenomenon, which severely hampers trainability and algorithm convergence. In this study, two novel strategies are introduced to address these challenges. First, we utilize a block encoding technique based on efficient oracles, thereby significantly reducing quantum circuit complexity and qubit overhead while accessing system matrix. Second, we incorporate a quantum-inspired evolutionary optimization (QIEO) approach, demonstrating enhanced training efficiency and effectively mitigating the barren plateau issue. These improvements are integrated within a Hybrid Quantum-Classical Method (HQCM) linear solver framework and demonstrated through solving a 2D steady-state heat equation, representative of the pressure-Poisson equations common in incompressible CFD simulations. The proposed methodologies yield promising results, showcasing improved circuit efficiency, enhanced stability during optimization, and accurate numerical solutions. This work lays a foundational step toward realizing practical, scalable quantum-assisted CFD solvers and sets clear pathways for future research aimed at further algorithmic and architectural enhancements.
We present a comprehensive review of quantum approaches for solving partial differential equations (PDEs) arising in computational fluid dynamics (CFD). We examine fully quantum approaches, including quantum linear system algorithms (QLSAs), ranging from the Harrow--Hassidim--Lloyd (HHL) algorithm to quantum singular value transformation (QSVT), Hamiltonian simulation, and quantum lattice Boltzmann methods (QLBMs), while emphasizing hybrid quantum--classical approaches, including quantum physics-informed neural networks (QPINNs) and amplitude-encoded variational PDE solvers. We focus on hardware-agnostic algorithms compatible with present noisy processors and emerging fault-tolerant architectures. For each framework, we analyze the mathematical formulation, algorithmic structure, and principal limitations. We also examine tensor-network (TN) representations, since CFD fields, differential operators, and geometrical information can often be encoded efficiently in low-rank form. The TN formalism bridges CFD discretizations and quantum states, operators, and circuits, enabling compact representations to be translated into tensor-programmable variational quantum algorithms (TP-VQAs). We further review benchmark problems, including Poisson, reaction, diffusion, and nonlinear model equations, and assess how well quantum algorithms capture key features of fluid dynamics. Our analysis highlights that potential quantum advantage is highly problem dependent and governed by condition number, representational complexity, state preparation, and measurement constraints. We outline capabilities, limitations, and challenges toward scalable quantum algorithms for CFD.
There has been significant progress in the development of quantum algorithms for solving linear systems of equations with a growing body of applications to Computational Fluid Dynamics (CFD) and CFD-like problems. This work extends previous work by developing a non-linear hybrid quantum-classical CFD solver and using it to generate fully converged solutions. The hybrid solver uses the SIMPLE CFD algorithm, which is common in many industrial CFD codes, and applies it to the 2-dimensional lid driven cavity test case. A theme of this work is the classical processing time needed to prepare the quantum circuit with a focus on the decomposition of the CFD matrix into a linear combination of unitaries (LCU). CFD meshes with up to 65x65 nodes are considered with the largest producing a LCU containing 32,767 Pauli strings. A new method for rapidly re-computing the coefficients in a LCU is proposed, although this reduces, rather than eliminates, the classical scaling issues. The quantum linear equation solver uses the Harrow, Hassidim, Lloyd (HHL) algorithm via a state-vector emulator. Test matrices are sampled from the classical CFD solver to investigate the solution accuracy that can be achieved with HHL. For the smallest 5x5 and 9x9 CFD meshes, full non-linear hybrid CFD calculations are performed. The impacts of approximating the LCU and the varying the number of ancilla rotations in the eigenvalue inversion circuit are studied. Preliminary timing results indicate that the classical computer preparation time needed for a hybrid solver is just as important to the achievement of quantum advantage in CFD as the time on the quantum computer. The reported HHL solutions and LCU decompositions provide a benchmark for future research. The CFD test matrices used in this study are available upon request.
The solution for non-linear, complex partial differential Equations (PDEs) is achieved through numerical approximations, which yield a linear system of equations. This approach is prevalent in Computational Fluid Dynamics (CFD), but it restricts the mesh size since the solution of the linear system becomes computationally intractable when the mesh resolution increases. The reliance on the ability of High-Performance Computers (HPC) to scale up and meet these requirements is myopic; such very high-fidelity simulations require a paradigm shift in computing. This paper presents an exploration of quantum methodologies aimed at achieving high accuracy in solving such a large system of equations. Leveraging recent works in Quantum Linear Solver Algorithms (QLSA) and variational algorithms suitable for Quantum Simulation in HPC, we aspire to push the boundaries of CFD-relevant problems that can be solved on hybrid quantum-classical framework. To this end, we consider the 2D, transient, incompressible, viscous, non-linear coupled Burgers equation as a test problem and investigate the accuracy of our approach by comparing results with a classical linear system of equation solvers, such as the Generalized Minimal RESidual method (GMRES). Through rigorous testing, our findings demonstrate that our quantum methods yield results comparable in accuracy to traditional approaches. Additionally, we demonstrate the accuracy, scalability, and consistency of our quantum method. Lastly, we present an insightful estimation of the resources our quantum algorithm needs to solve systems with nearly 2 billion mesh points.
Computational Fluid Dynamics simulations are crucial in industrial applications but require extensive computational resources, particularly for extreme turbulent regimes. While classical digital approaches remain the standard, quantum computing promises a breakthrough by enabling a more efficient encoding of large-scale simulations with a limited number of qubits. This work presents a practical numerical assessment of a hybrid quantum-classical approach to CFD based on the Lattice Boltzmann Method (LBM). The inherently non-linear LBM equations are linearized via a Carleman expansion and solved using the quantum Harrow Hassidim Lloyd algorithm (HHL). We evaluate this method on three benchmark cases featuring different boundary conditions, periodic, bounceback, and moving wall, using statevector emulation on high-performance computing resources. Our results confirm the validity of the approach, achieving median error fidelities on the order of 10−3 and success probabilities sufficient for practical quantum state sampling. Notably, the spectral properties of small lattice systems closely approximate those of larger ones, suggesting a pathway to mitigate one of HHL's bottlenecks: eigenvalue pre-evaluation.