Reduced-order Modeling on a Near-term Quantum Computer

TL;DR

This paper develops a quantum algorithm for reduced-order modeling of fluid flows using dynamic mode decomposition reformulated as an optimization problem, demonstrating quantum predictions for flow scenarios and analyzing computational complexity.

Contribution

It introduces a quantum approach to reduced-order modeling by reformulating DMD as an optimization problem solvable with quantum algorithms, advancing quantum fluid dynamics modeling.

Findings

Quantum-ROM predictions depend on bit precision and truncation level.

Quantum predictions are comparable to classical DMD results.

Analysis of quantum algorithm complexity and future prospects.

Abstract

Quantum computing is an advancing area of research in which computer hardware and algorithms are developed to take advantage of quantum mechanical phenomena. In recent studies, quantum algorithms have shown promise in solving linear systems of equations as well as systems of linear ordinary differential equations (ODEs) and partial differential equations (PDEs). Reduced-order modeling (ROM) algorithms for studying fluid dynamics have shown success in identifying linear operators that can describe flowfields, where dynamic mode decomposition (DMD) is a particularly useful method in which a linear operator is identified from data. In this work, DMD is reformulated as an optimization problem to propagate the state of the linearized dynamical system on a quantum computer. Quadratic unconstrained binary optimization (QUBO), a technique for optimizing quadratic polynomials in binary variables, allows for quantum annealing algorithms to be applied. A quantum circuit model (quantum approximation optimization algorithm, QAOA) is utilized to obtain predictions of the state trajectories. Results are shown for the quantum-ROM predictions for flow over a 2D cylinder at Re = 220 and flow over a NACA0009 airfoil at Re = 500 and $\alpha = 15^\circ{}$. The quantum-ROM predictions are found to depend on the number of bits utilized for a fixed point representation and the truncation level of the DMD model. Comparisons with DMD predictions from a classical computer algorithm are made, as well as an analysis of the computational complexity and prospects for future, more fault-tolerant quantum computers.