Efficient Computation of Large-Scale Statistical Solutions to Incompressible Fluid Flows

Abstract

This work presents the development, performance analysis and subsequent optimization of a GPU-based spectral hyperviscosity solver for turbulent flows described by the three dimensional incompressible Navier-Stokes equations. The method solves for the fluid velocity fields directly in Fourier space, eliminating the need to solve a large-scale linear system of equations in order to find the pressure field. Special focus is put on the communication intensive transpose operation required by the fast Fourier transform when using distributed memory parallelism. After multiple iterations of benchmarking and improving the code, the simulation achieves close to optimal performance on the Piz Daint supercomputer cluster, even outperforming the Cray MPI implementation on Piz Daint in its communication routines. This optimal performance enables the computation of large-scale statistical solutions of incompressible fluid flows in three space dimensions.

References in corpus (6)

• HACC: Simulating Sky Surveys on State-of-the-Art Supercomputing Architectures
• Statistical solutions of hyperbolic conservation laws I: Foundations
• Statistical solutions of the incompressible Euler equations
• Computation of measure-valued solutions for the incompressible Euler equations
• On the convergence of the spectral viscosity method for the two-dimensional incompressible Euler equations with rough initial data
• Arbitrarily high-order (weighted) essentially non-oscillatory finite difference schemes for anelastic flows on staggered meshes