Documento completo qui
A pseudo-spectral code for the direct numerical simulation of incompressible Navier-Stokes equations for homogeneous isotropic turbulence is presented. The program is developed to run on graphical processing units (GPUs), exploiting the massive parallelism of modern GPU architectures. The pseudo- spectral method is based on Fourier transforms: the code works in wavenumbers domain, but non linearities are handled in physical domain instead of executing expensive complex convolution computations.
The GPU handles thousands of concurrent threads since it is composed of hundreds of cores, as opposed to a low number of threads concurrently running on a CPU. This characteristic makes operating on large data like matrices very efficient, providing a substantial speed-up compared to traditional CPU parallel code.