Toward Quantum-Enabled Fluid Simulation: A Hybrid Surrogate for Lattice Boltzmann Collisions

This abstract has open access
Problem description and relevance

High-fidelity simulation of complex transport phenomena remains a central challenge in computational fluid dynamics. Direct Numerical Simulation (DNS) can resolve fluid behaviour with high accuracy, but its extremely large computational and memory requirements make it impractical for large-scale or industrial applications on classical hardware.

Quantum Computing is a potential way to reduce computational cost while keeping accuracy. However, quantum-assisted fluid dynamics still faces hardware limits, shallow circuit depths, and difficulty in representing nonlinear physical effects within quantum circuits. Consequently, many current methods must simplify physics, use fixed or constrained parameters, or remain limited to low-complexity flow regimes, which reduces their practical applicability.

The Lattice Boltzmann Method (LBM) is promising for quantum implementation due to its structured, locally defined update rules. Despite this, a major bottleneck arises with the collision operator. This component is inherently nonlinear and non-unitary, making it hard to implement directly on gate-based quantum computers.

The core challenge is both computational and physical. Efficiently and accurately capturing nonlinear fluid interactions at scale within a quantum framework is very difficult. Addressing this challenge has significant practical value: it would enable more efficient and accurate simulation of complex flows in engineering and science. For this reason, quantum-assisted formulations are motivated by the need to overcome the limits of classical computation while remaining compatible with near-term quantum hardware. This represents a key step toward practical quantum advantage in fluid dynamics simulation.

Submission ID :
27
Methodology :

The proposed approach adopts a hybrid quantum–classical machine learning framework in which a quantum machine learning (QML) surrogate model offloads the non-unitary collision dynamics to a parameterised quantum circuit acting as a function approximator.

The main inputs are encoded directly into the rotational parameters of single-qubit gates. Trainable scaling parameters are added to improve representational capacity. This extends the frequency spectrum that the circuit can learn in its implicit Fourier representation.

Each distribution function component is measured on a single qubit, mapping its expectation value to a collision outcome. This reduces hardware complexity, but requires multiple circuit runs to reconstruct a full distribution. Only three parameter sets are needed for two- and three-dimensional velocity sets. A classical post-processing stage restores physical consistency.

Design choices are evaluated with quantum hardware metrics such as expressibility, entanglement, and effective dimension. Next, these results are validated by evaluating different architectural configurations on a quantum simulator. Finally, training-dependent effects are systematically examined, including sensitivity to input range, interpolation, and extrapolation.

Practical demonstration :

The proposed quantum-assisted BGK collision model is validated on a gate-based quantum simulator. The demonstration shows that the quantum surrogate reproduces full BGK collision dynamics across benchmark cases. These include both two-dimensional (D2Q9) and three-dimensional (D3Q19) lattice configurations and a range of physical regimes.

Numerical results show the surrogate accurately reproduces global integral quantities and fine-scale flow structures. This is quantified using standard diagnostic measures: relative error, kinetic energy, and enstrophy. These measures capture both macroscopic conservation and small-scale dynamics.

The classical post-processing stage is systematically analysed. Its role is to restore physical consistency and to improve the stability of higher-order moments.

Application potential :

The proposed quantum-assisted BGK collision model is designed for scalability using a hybrid quantum–classical architecture. This framework partitions computational tasks according to their natural strengths. The quantum processor is solely responsible for evaluating the collision operator using a parameterised quantum circuit. Meanwhile, all streaming steps, macroscopic moment calculations, and global time evolution are handled efficiently on classical hardware. This separation ensures that only the non-linear tasks are managed by the quantum processor. The main objective is to develop an efficient formulation of computational fluid dynamics (CFD) methods for execution on quantum computing hardware. By doing so, the approach aims to exploit the computational capabilities of quantum systems to accelerate and enhance the simulation of complex fluid flows beyond classical limits.

Research Associate
,
Technical University of Munich
PhD Student
,
Technical University of Munich
,
Technical University of Munich
Siemens Digital Industries Software
Technical University of Munich
14 visits