
I built a GPU-resident CFD solver in CUDA and would love feedback on the architecture
I’ve been working on brae, a GPU-resident CFD solver for OpenFOAM-style cases.
The goal is to run existing finite-volume CFD cases on one GPU while avoiding the usual CPU-GPU transfer loop every iteration. The current path is simpleFoam-compatible: matrix assembly, pressure correction, turbulence updates, and sparse linear solves are designed to stay resident on the GPU.
Repo: https://github.com/simd-ai/brae
This is still early, but in the current benchmark it is around ~5× faster than a GPU-accelerated OpenFOAM setup on the same GPU, with validation error under 1%.
I’d really appreciate technical feedback on the CUDA architecture, benchmark methodology, or solver design.