Abstract Computational fluid dynamics (CFD) simulations are a vital part of the design process in the aerospace industry. Although reliable CFD results can be obtained with turbulence models, direct numerical simulation of complex bodies in three spatial dimensions (3D) is impracticable due to the massive amount of computational elements. For instance, a 3D direct numerical simulation of a turbulent boundary-layer over the wing of a commercial jetliner that resolves all relevant length scales using a serial CFD solver on a modern digital computer would take approximately 750 million years or roughly 20% of the earth’s age. Over the past 25 years, quantum computers have become the object of great interest worldwide as powerful quantum algorithms have been constructed for several important, computationally challenging problems that provide enormous speed-up over the best-known classical algorithms. In this paper, we adapt a recently introduced quantum algorithm for partial differential equations to Burgers’ equation and develop a quantum CFD solver that determines its solutions. We used our quantum CFD solver to verify the quantum Burgers’ equation algorithm to find the flow solution when a shockwave is and is not present. The quantum simulation results were compared to: (i) an exact analytical solution for a flow without a shockwave; and (ii) the results of a classical CFD solver for flows with and without a shockwave. Excellent agreement was found in both cases, and the error of the quantum CFD solver was comparable to that of the classical CFD solver.
Abstract Computational fluid dynamics (CFD) simulations are a vital part of the design process in the aerospace industry. Although reliable CFD results can be obtained with turbulence models, direct numerical simulation of complex bodies in three spatial dimensions (3D) is impracticable due to the massive amount of computational elements. For instance, a 3D direct numerical simulation of a turbulent boundary-layer over the wing of a commercial jetliner that resolves all relevant length scales using a serial CFD solver on a modern digital computer would take approximately 750 million years or roughly 20% of the earth’s age. Over the past 25 years, quantum computers have become the object of great interest worldwide as powerful quantum algorithms have been constructed for several important, computationally challenging problems that provide enormous speed-up over the best-known classical algorithms. In this paper, we adapt a recently introduced quantum algorithm for partial differential equations to Burgers’ equation and develop a quantum CFD solver that determines its solutions. We used our quantum CFD solver to verify the quantum Burgers’ equation algorithm to find the flow solution when a shockwave is and is not present. The quantum simulation results were compared to: (i) an exact analytical solution for a flow without a shockwave; and (ii) the results of a classical CFD solver for flows with and without a shockwave. Excellent agreement was found in both cases, and the error of the quantum CFD solver was comparable to that of the classical CFD solver.
As the energy sector transitions towards more sustainable forms of energy production, classical Computational Fluid Dynamics (CFD) increasingly struggles with the bottle- necks of simulating high complexity turbulent flows. While quantum computing (QC) theoretically offers an exponential speedup for large-scale linear systems, current Noisy Intermediate-Scale Quantum (NISQ) hardware lacks the error correction required for deep and complex circuits that mathematically ideal algorithms require. Furthermore, the inherent linearity of quantum mechanics seriously impedes the solving non-linear partial differential equations (PDEs), which are critical for fluid modeling. To bridge thisgap, thisthesisevaluatestheVariationalQuantumLinearSolver(VQLS),a shallow- circuit, hybrid algorithm, applied to the 1D viscous Burgers’ equation while utilizing the Cole-Hopf transformation to analytically linearize the equation. In both simulated and actual hardware, the VQLS algorithm is implemented across increasingly scaled topographies (2, 4, and 8 qubits). This work intends to identify the limitations of implementing quantum fluid dynamics in near-term quantum hardware.
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.
Abstract There is great interest in using quantum computers to efficiently simulate a quantum system’s dynamics as existing classical computers cannot do this. Little attention, however, has been given to quantum simulation of a classical nonlinear continuum system such as a viscous fluid even though this too is hard for classical computers. Such fluids obey the Navier–Stokes nonlinear partial differential equations, whose solution is essential to the aerospace industry, weather forecasting, plasma magneto-hydrodynamics, and astrophysics. Here we present a quantum algorithm for solving the Navier–Stokes equations. We test the algorithm by using it to find the steady-state inviscid, compressible flow through a convergent-divergent nozzle when a shockwave is (is not) present. We find excellent agreement between numerical simulation results and the exact solution, including shockwave capture when present. Finally, we compare the algorithm’s computational cost to deterministic and random classical algorithms and show that a significant speed-up is possible. Our work points to a large new application area for quantum computing with substantial economic impact, including the trillion-dollar aerospace industry, weather-forecasting, and engineered-plasma technologies.
This paper explores the potential of quantum computing for computational fluid dynamics (CFD) applications, with a focus on turbomachinery CFD. A harmonic balance solver is developed for Burgers’ equation, in which discrete Fourier transforms are approximated using a hybrid quantum-classical algorithm based on the quantum Fourier transform (QFT). Three novel algorithms are presented to estimate Fourier coefficients, providing complete knowledge of their amplitudes and phases, which is not possible with the standard QFT. Their behaviour is studied theoretically and numerically, under noiseless and noisy conditions. The algorithms, whose performance is limited by the extensive sampling required to achieve a small error on the Fourier coefficients, tolerate low levels of depolarising noise. The behaviour of the hybrid solver is investigated for a baseline case with a Reynolds number of 1000, 7 harmonics and 100 grid cells, using the best performing algorithm in a noiseless setting with up to 10 8 samples per QFT. Residuals decrease until errors introduced by statistical uncertainty dominate the total error on the solution. Nevertheless, the hybrid solutions match the classical one closely, with RMS residuals as low as 10 –4 . The impact of several solver parameters on convergence and solution quality is also assessed, including the effect of noise. Although the latter rapidly degrades the solution, the solver achieves satisfactory approximations to the classical solution with 0.001% of depolarising noise. Without seeking to demonstrate a quantum advantage, this work offers valuable insights into the opportunities and challenges of quantum computing, helping readers understand how to design, implement and study quantum algorithms, as well as evaluate their impact in CFD applications.
Fault-tolerant quantum computing is a promising technology to solve linear partial differential equations that are classically demanding to integrate. It is still challenging to solve non-linear equations in fluid dynamics, such as the Burgers equation, using quantum computers. We propose a novel quantum algorithm to solve the Burgers equation. With the Cole-Hopf transformation that maps the fluid velocity field u to a new field ψ, we apply a sequence of quantum gates to solve the resulting linear equation and obtain the quantum state ∣ψ⟩ that encodes the solution ψ. We also propose an efficient way to extract stochastic properties of u, namely the multi-point functions of u, from the quantum state of ∣ψ⟩. Our algorithm offers an exponential advantage over the classical finite difference method in terms of the number of spatial grids when a perturbativity condition in the information-extracting step is met.
Recent advances in quantum hardware and quantum computing algorithms promise significant breakthroughs in computational capabilities. Quantum computers can achieve exponential improvements in speed vs classical computers by employing principles of quantum mechanics like superposition and entanglement. However, designing quantum algorithms to solve the nonlinear partial differential equations governing fluid dynamics is challenging due to the inherent linearity of quantum mechanics, which requires unitary transformation. In this study, we first address in detail several challenges that arise when trying to deal with nonlinearity using quantum algorithms and then propose a novel pure quantum algorithm for solving a nonlinear Burgers' equation. We employed multiple copies of the state vector to calculate the nonlinear term, which is necessary due to the no-cloning theorem. By reusing qubits from the previous time steps, we significantly reduced the number of qubits required for multi-step simulations, from exponential/quadratic scaling in earlier studies to linear scaling in time in the current study. We also employed various advanced quantum techniques, including block-encoding, quantum Hadamard product, and the linear combination of unitaries, to design a quantum circuit for the proposed quantum algorithm. The quantum circuit was executed on quantum simulators, and the obtained results demonstrated excellent agreement with those from classical simulations.