A parallelized eno procedure for direct numerical simulation of compressible turbulence

In this paper, a finite-difference based ENO (essentially nonoscillatory) procedure has been chosen for the direct numerical simulation (DNS) of compressible turbulence. The implementation of the ENO scheme follows the relatively efficient procedure in Shu et al. (1992), but the latter has been modified in the present paper to admit scalar conservation equations and to run on the iPSC/860 Paragon parallel supercomputer. DNS results with our procedure are in excellent agreement with pseudo-spectral and Padé approximation calculations in two and three dimensions. This is the case for a variety of initial conditions for compressible turbulence. The parallel algorithms presented are simple but quite efficient for DNS. with a speedup that approaches the theoretical value. Some of the attractive features include 1 ) minimum communication whereby a processor only communicates with two neighbors, 2) almost one hundred percent load balancing, 3) a checker-board approach to solve the Poisson equation reduces communication by a factor of approximately 2, and, 4) obtaining turbulence statistics is based on a 'global collect' approach, which is implemented to ensure that a single number, rather than a large matrix of numbers, is communicated between processors. The ENO code presented in this paper should be quite useful in its own right, while the parallel implementation should allow the simulation of fairly realistic problems.