The recent success of the fusion ignition experiment at Lawrence-Livermore National Laboratory has renewed excitement for inertial confinement fusion. Designing such experiments relies on computational modeling using radiation hydrodynamic codes run on massively parallel supercomputers which simulate hydrodynamic flows, radiation diffusion, and thermonuclear burn. Constructing a quantum algorithm for radiation hydrodynamics that provides a quantum speedup is of great interest. A recent quantum algorithm that solves the Navier-Stokes equations of hydrodynamics with a quadratic speedup marks a first step towards this goal. Here we take the next step and present a quantum algorithm for nonlinear radiation diffusion that also has a quadratic speedup. To verify the algorithm we consider radiation striking a cold, optically thick target that generates a Marshak wave in the target. Results of numerical simulation of the quantum algorithm are compared with those of a standard partial differential equation solver and excellent agreement is found.
The recent success of the fusion ignition experiment at Lawrence-Livermore National Laboratory has renewed excitement for inertial confinement fusion. Designing such experiments relies on computational modeling using radiation hydrodynamic codes run on massively parallel supercomputers which simulate hydrodynamic flows, radiation diffusion, and thermonuclear burn. Constructing a quantum algorithm for radiation hydrodynamics that provides a quantum speedup is of great interest. A recent quantum algorithm that solves the Navier-Stokes equations of hydrodynamics with a quadratic speedup marks a first step towards this goal. Here we take the next step and present a quantum algorithm for nonlinear radiation diffusion that also has a quadratic speedup. To verify the algorithm we consider radiation striking a cold, optically thick target that generates a Marshak wave in the target. Results of numerical simulation of the quantum algorithm are compared with those of a standard partial differential equation solver and excellent agreement is found.
In the radiation hydrodynamic simulations used to design inertial confinement fusion (ICF) and pulsed power experiments, nonlinear radiation diffusion tends to dominate CPU time. This raises the interesting question of whether a quantum algorithm can be found for nonlinear radiation diffusion which provides a quantum speedup. Recently, such a quantum algorithm was introduced based on a quantum algorithm for solving systems of nonlinear partial differential equations (PDEs) which provides a quadratic quantum speedup. Here, we apply this quantum PDE (QPDE) algorithm to the problem of a non-equilibrium Marshak wave propagating through a cold, semi-infinite, optically thick target, where the radiation and matter fields are not assumed to be in local thermodynamic equilibrium. The dynamics is governed by a coupled pair of nonlinear PDEs which are solved using the QPDE algorithm, as well as two standard PDE solvers: (i) Python's py-pde solver; and (ii) the KULL ICF simulation code developed at Lawrence-Livermore National Laboratory. We compare the simulation results obtained using the QPDE algorithm and the standard PDE solvers and find excellent agreement.
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.
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.
A practical and efficient way to use the power of quantum mechanics (quantum parallelism due to the superposition and entanglement of states) could enormously speed up numerical simulations of interest to computational physicists. The diffusion equation the nonlinear Burgers equations, the turbulent Navier-Stokes equations, and the many-body Schroedinger equation could all be solved much more rapidly if efficient quantum algorithms could be implemented for their solution. The presentations at the conference that are published here review the progress toward these goals.
Abstract Efficient simulation of a quantum system's dynamics is expected to be an important application area for quantum computers as existing classical computers cannot do this. However, quantum systems are not unique in being hard to simulate. For example, classical nonlinear continuum systems and fields are governed by nonlinear partial differential equations whose solution is also hard for classical computers. Solving such equations is essential for many economically important industries/applications such as the aerospace industry, weather‐forecasting, fiber‐optics communication, and plasma magneto‐hydrodynamics. This raises the question: can a quantum computer speed‐up solving these equations? In this Review, a new quantum algorithm is described for solving nonlinear partial differential equations for which the answer is yes. First, a new quantum algorithm is discussed for solving the Navier‐Stokes nonlinear partial differential equations which govern the flow of a viscous fluid. Its construction, verification, and computational cost are described, and it is shown that a significant quantum speed‐up is possible. Its generalization to a quantum algorithm for solving nonlinear partial differential equations is described. The Review closes with a discussion of next steps forward. These new quantum algorithms open up a large new application area for quantum computing with substantial economic impact, including the trillion‐dollar aerospace industry.
Simulating fluid dynamics on a quantum computer is intrinsically difficult due to the nonlinear and non-Hamiltonian nature of the Navier-Stokes equation (NSE). We propose a framework for quantum computing of fluid dynamics based on the hydrodynamic Schr\"odinger equation (HSE), which can be promising in simulating three-dimensional turbulent flows in various engineering applications. The HSE is derived by generalizing the Madelung transform to compressible or incompressible flows with finite vorticity and dissipation. Since the HSE is expressed as a unitary operator on a two-component wave function, it is more suitable than the NSE for quantum computing. The flow governed by the HSE can resemble a turbulent flow consisting of tangled vortex tubes with the five-thirds scaling of energy spectrum. We develop a prediction-correction quantum algorithm to solve the HSE. This algorithm is implemented for simple flows on the quantum simulator Qiskit with partial exponential speedup.