Performance of a Two-Path Aliasing Free Calculation of a Spectral DNS Code

A direct numerical simulation (DNS) code was developed for solving incompressible homogeneous isotropic turbulence with high Reynolds numbers in a periodic box using the Fourier spectral method. The code was parallelized using the Message Passing Interface and OpenMP with two-directional domain decomposition and optimized on the K computer. High resolution DNSs with up to \(12288^3\) grid points were performed on the K computer using the code. Efficiencies of 3.84%, 3.14%, and 2.24% peak performance were obtained in double precision DNSs with \(6144^3\), \(8192^3\), and \(12288^3\) grid points, respectively. In addition, a two-path alias-free procedure is proposed and clarified its effectiveness for some number of parallel processes.