Quantum Algorithm for Solving the Advection Equation using Hamiltonian Simulation

Peter Brearley,Sylvain Laizet
2024-06-19
Abstract:A quantum algorithm for solving the advection equation by embedding the discrete time-marching operator into Hamiltonian simulations is presented. One-dimensional advection can be simulated directly since the central finite difference operator for first-order derivatives is anti-Hermitian. Here, this is extended to industrially relevant, multi-dimensional flows with realistic boundary conditions and arbitrary finite difference stencils. A single copy of the initial quantum state is required and the circuit depth grows linearly with the required number of time steps, the sparsity of the time-marching operator and the inverse of the allowable error. Statevector simulations of a scalar transported in a two-dimensional channel flow and lid-driven cavity configuration are presented as a proof of concept of the proposed approach.
Quantum Physics
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to propose a quantum algorithm for solving the Advection Equation. Specifically, by embedding the discrete - time advancement operator into Hamiltonian simulation, the author designed a quantum algorithm that can effectively handle multi - dimensional industrially relevant flows. This algorithm is not only applicable to convection problems under one - dimensional periodic boundary conditions, but also extended to multi - dimensional flow problems with practical boundary conditions and arbitrary finite - difference templates. ### Background and Objectives of the Paper 1. **Importance of the Advection Equation** - The Advection Equation is a fundamental linear partial differential equation (PDE), widely used in multiple industries, such as oceanography, atmospheric science, geology, biomathematics, and heat exchanger design. - The form of the equation is: \[ \frac{\partial \phi}{\partial t}+u_j\frac{\partial \phi}{\partial x_j} = 0 \] where \(\phi\) represents the scalar field (e.g., temperature or concentration), and \(u_j\) is the \(j\) - th component of the advection velocity vector. 2. **Limitations of Existing Methods** - Classical computers face huge computational resource requirements when solving large - scale partial differential equations. - Quantum computing has become a new approach to solving this problem because of its potential computational advantages. ### Main Contributions of the Paper 1. **Design of the Quantum Algorithm** - A quantum algorithm based on Hamiltonian simulation was proposed, which can effectively solve the Advection Equation. - The core of the algorithm lies in embedding the discrete - time advancement operator into Hamiltonian simulation, thereby achieving efficient solution of multi - dimensional flow problems. 2. **Scope of Application of the Algorithm** - It is applicable to convection problems under one - dimensional periodic boundary conditions and also extended to multi - dimensional industrially relevant flow problems. - It can handle multi - dimensional flow problems with practical boundary conditions and arbitrary finite - difference templates. 3. **Performance Analysis of the Algorithm** - Through state - vector simulation, the effectiveness of the algorithm in two - dimensional channel flow and lid - driven cavity flow was verified. - The time complexity and error bounds of the algorithm were analyzed, demonstrating its efficiency under low numerical stability parameter \(r\). ### Experimental Results 1. **Two - Dimensional Laminar Channel Flow** - Spatial discretization was carried out using 64×64 grid points, with 12 qubits representing the solution and an additional auxiliary qubit, for a total of 13 qubits used. - The initial scalar field \(\phi\) was set as a sine wave in the horizontal direction. - Simulation results show that the algorithm can effectively capture the transport process of the scalar field. 2. **Lid - Driven Cavity Flow** - By simulating lid - driven cavity flow, the effectiveness of the algorithm in handling complex boundary conditions was further verified. ### Conclusion This paper successfully proposed an efficient quantum algorithm that can solve multi - dimensional Advection Equations and has demonstrated good performance in practical applications. This achievement lays the foundation for future applications of quantum computing in the field of fluid mechanics.