A Quantum Spectral Solver for Periodic Incompressible Stokes Flow

Abstract

We present a quantum spectral solver for the steady incompressible Stokes equations on a two-dimensional periodic domain. The method uses the Quantum Fourier Transform as a coherent change of basis and exploits the resulting spectral structure of the Stokes operator: the Laplacian becomes diagonal, while incompressibility is enforced mode by mode through a Helmholtz projection. In two dimensions, this projection is realized by a mode-dependent rotation from Cartesian velocity components to longitudinal--transverse coordinates, followed by component-conditioned inverse-Laplacian scaling. The velocity and pressure fields are encoded as quantum states over Fourier modes and physical components, and the corresponding spectral factors are implemented through polynomially encoded amplitude blocks. The construction extends recent quantum spectral methods in computational mechanics to an incompressible flow operator with explicit pressure--velocity splitting and divergence-free projection. The approach is also compatible with multiscale finite-element architectures in which quantum parallelism can simultaneously update all representative volume element (RVE) states. Numerical verification includes a steady vortex, a regularized periodic force-dipole benchmark, and an RVE-inspired Kolmogorov-like fluctuation benchmark. The latter illustrates how the circuit can recover a homogenized kinetic-energy observable without reconstructing the full velocity field, consistent with the role of averaged quantities in multiscale flow calculations. Under the standard assumptions of efficient state preparation and observable estimation, the circuit has polylogarithmic dependence on the grid resolution, with the polynomial degree and tile count appearing as explicit approximation and implementation parameters.