Search papers, labs, and topics across Lattice.
This paper details enhancements to the WaterLily.jl incompressible flow solver, focusing on the integration of MPI-based parallelism and an optimized geometric multigrid solver. The authors achieve near-ideal strong scalability and maintain over 85% weak-scaling efficiency, even with large grid sizes up to 1 billion cells. Key improvements include an adaptive under-relaxed red-black Gauss-Seidel smoother and anisotropic coarsening operators, significantly boosting the solver's performance on diverse hardware platforms.
Achieving over 96% inter-node weak scalability with a grid size of 1 billion cells transforms the capabilities of flow simulations in computational fluid dynamics.
We present recent performance-oriented developments in WaterLily, a scale-resolving incompressible flow solver written in pure Julia that runs seamlessly on CPUs and GPUs of any vendor. Supported by the newly added MPI-based parallelism, strong-scalability tests display a near-ideal linear trend, and weak-scaling efficiency is kept above 85% before node memory-concurrency contention dominates parallel performance. Inter-node weak scalability is sustained above 96% with grid size up to 1 billion cells. We further benchmark improvements to the geometric multigrid Poisson solver enabled by an adaptive under-relaxed red-black Gauss-Seidel smoother together with anisotropic coarsening operators.
