Recent advancements of quantum technologies have triggered tremendous interest in exploring practical quantum advantage. The simulation of fluid dynamics, a highly challenging problem in classical physics but vital for practical applications, emerges as a potential direction. Here, we report an experiment on the digital simulation of unsteady flows with a superconducting quantum processor. The quantum algorithm is based on the Hamiltonian simulation using the hydrodynamic formulation of the Schrödinger equation. With the median fidelities of 99.97% and 99.67% for parallel single- and two-qubit gates respectively, we simulate the dynamics of a two-dimensional (2D) compressible diverging flow and a 2D decaying vortex with ten qubits. Note that the former case is an inviscid potential flow, and the latter one is an artificial vortical flow with an external body force. The experimental results well capture the temporal evolution of averaged density and momentum profiles, and qualitatively reproduce spatial flow fields with moderate noises. This work demonstrates the potential of quantum computing in simulating more complex flows, such as turbulence, for practical applications. Fluid dynamics simulation, a complex challenge in classical physics, is relevant for real-world applications and highlights the potential of quantum computing. The authors report an experiment for the digital simulation of unsteady flows on a superconducting quantum processor, and show that the results effectively capture the evolution of flow fields.
Recent advancements of quantum technologies have triggered tremendous interest in exploring practical quantum advantage. The simulation of fluid dynamics, a highly challenging problem in classical physics but vital for practical applications, emerges as a potential direction. Here, we report an experiment on the digital simulation of unsteady flows with a superconducting quantum processor. The quantum algorithm is based on the Hamiltonian simulation using the hydrodynamic formulation of the Schrödinger equation. With the median fidelities of 99.97% and 99.67% for parallel single- and two-qubit gates respectively, we simulate the dynamics of a two-dimensional (2D) compressible diverging flow and a 2D decaying vortex with ten qubits. Note that the former case is an inviscid potential flow, and the latter one is an artificial vortical flow with an external body force. The experimental results well capture the temporal evolution of averaged density and momentum profiles, and qualitatively reproduce spatial flow fields with moderate noises. This work demonstrates the potential of quantum computing in simulating more complex flows, such as turbulence, for practical applications. Fluid dynamics simulation, a complex challenge in classical physics, is relevant for real-world applications and highlights the potential of quantum computing. The authors report an experiment for the digital simulation of unsteady flows on a superconducting quantum processor, and show that the results effectively capture the evolution of flow fields.
Vortex interactions are commonly observed in atmospheric turbulence, plasma dynamics, and collective behaviors in biological systems. However, accurately simulating these complex interactions is highly challenging due to the need to capture fine-scale details over extended timescales, which places computational burdens on traditional methods. In this study, we introduce a quantum vortex method, reformulating the Navier–Stokes (NS) equations within a quantum mechanical framework to enable the simulation of multi-vortex interactions on a quantum computer. We construct the effective Hamiltonian for the vortex system and implement a spatiotemporal evolution circuit to simulate its dynamics over prolonged periods. By leveraging eight qubits on a superconducting quantum processor with gate fidelities of 99.97% for single-qubit gates and 99.76% for two-qubit gates, we successfully reproduce natural vortex interactions. Overall, we establish a framework that reformulates vortex dynamics into a normalized wavefunction representation compatible with quantum system unitary evolution, combined with the designed spatiotemporal encoding scheme, providing a concrete pathway toward leveraging quantum resources in fluid systems. Despite recent progress, simulations of complex fluid dynamics phenomena on quantum computers remain challenging. Here, the authors present a quantum friendly reformulation of vortex dynamics with a spatiotemporal encoding scheme and apply it to simulations on a superconducting quantum processor.
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.
Quantum simulation, as a practical application of noisy quantum computing, has aided the study of exotic quantum matters and the implementation of algorithms that outperform classical approaches. Superconducting qubits, one of the most promising candidates for realizing universal quantum computing, possess state-of-the-art features like easy integration of qubits, long coherence time, and high-fidelity single- and two-qubit gates. These characteristics have enabled applications of digital quantum simulation in the fields of physics, chemistry, and computer science. In this review, we first present the basic concepts of superconducting qubits, quantum gates, and digital quantum simulations. We also explore recent progress in digital quantum simulations using superconducting qubits, especially in relation to quantum chemistry, quantum matters, combinatorial optimization, and quantum machine learning. Finally, we address the current challenges of digital quantum simulation with superconducting qubits, and provide a perspective on the future of the field.
From fluid flow and transport to collective dynamics, numerical simulation of nonlinear partial differential equations underpins modern scientific computing. Extending this capability to quantum computers remains a longstanding challenge because nonlinear and non-Hermitian evolution is fundamentally incompatible with conventional Hamiltonian-based quantum simulation. Here we experimentally realize the time evolution of nonlinear fluid dynamics on a quantum processor using a hybrid variational framework for the viscous and inviscid Burgers equations. Our approach directly encodes the nonlinear dynamics into a variational optimization procedure, avoiding the enlarged linear embeddings and truncation overhead associated with Carleman linearization-based quantum algorithms. We further demonstrate convection-dominated dynamics corresponding to Reynolds numbers of order 102. We encode the governing evolution into parametrized quantum circuits and iteratively reconstruct the time-dependent field through quantum-classical optimization. By introducing a zero-noise extrapolation method without additional circuit-folding overhead, we accurately execute deep error-circuits with entangling-gate counts beyond those typical of Hadamard test circuits. We accurately reconstruct the time evolution across multiple timesteps despite hardware noise and finite device coherence. Our results constitute, to our knowledge, the first experimental realization of nonlinear time propagation on a quantum processor, extending quantum simulation beyond predominantly linear settings and establishing a route toward quantum computation for nonlinear continuum dynamics.
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 work aims to address the bottleneck issues of hardware resource limitation and decoherence error in the Hamiltonian simulation of quantum fluids, which are caused by the standard quantum Fourier transform and the evolution of momentum operators, resulting in excessively deep circuits and excessive two-qubit gates. We propose an approximate operator optimization scheme aimed at reducing the circuit depth in Hamiltonian evolution. The proposed scheme successfully reduces the depth of analog circuits from O(n2) to O(nlogn) or even O(n) by eliminating O(n2) redundant two-qubit entangling gates. In this work, the numerical experiments are implemented on a supercomputing-oriented quantum simulator, simulating two-dimensional unsteady divergent flow. Experimental results demonstrate that although the truncation of high-frequency qubit coupling terms introduces deterministic theoretical errors, scaling at O(n) for AQFT and O(n2) for momentum truncation, the optimized simulations successfully preserve the inherent macroscopic temporal evolution characteristics of the fluid in a 10-qubit simulation, achieving high correlation coefficients of r=0.933, r=0.941, and r=0.977 for density, X-momentum, and Y-momentum distributions respectively. Furthermore, we also analyzed the relationship between the algorithm truncation error and the hardware cumulative noise when the qubit number is extended to a higher level. This study proves that rationally adjusting truncation thresholds can establish an equilibrium point, preventing the hardware cumulative error from rapidly approaching 100% at the 20-30 qubit scale, providing a feasible engineering pathway for simulating complex fluid systems on real quantum devices in the future.