Spectral methods for the Navier-Stokes equations with one infinite and two periodic directions
Journal of Computational Physics
SIAM Journal on Scientific and Statistical Computing
Effects of the computational time step on numerical solutions of turbulent flow
Journal of Computational Physics
SIAM Journal on Numerical Analysis
Algorithm for solving tridiagonal matrix problems in parallel
Parallel Computing
Analysis of multigrid algorithms on massively parallel computers: architectural implications
Journal of Parallel and Distributed Computing
Journal of Computational Physics
Iterative methods for solving linear systems
Iterative methods for solving linear systems
Preconditioning for the Steady-State Navier--Stokes Equations with Low Viscosity
SIAM Journal on Scientific Computing
A comparison of optimal FFTs on torus and hypercube multicomputers
Parallel Computing
BoomerAMG: a parallel algebraic multigrid solver and preconditioner
Applied Numerical Mathematics - Developments and trends in iterative methods for large systems of equations—in memoriam Rüdiger Weiss
Parallel multigrid smoothing: polynomial versus Gauss--Seidel
Journal of Computational Physics
High order accurate solution of the incompressible Navier-Stokes equations
Journal of Computational Physics
Principles of Computational Fluid Dynamics
Principles of Computational Fluid Dynamics
Benchmarks on tera-scalable models for DNS of turbulent channel flow
Parallel Computing
Journal of Computational Physics
Hierarchical matrix preconditioners for the Oseen equations
Computing and Visualization in Science
A high-resolution code for turbulent boundary layers
Journal of Computational Physics
Breakdown-free ML(k)BiCGStab algorithm for non-Hermitian linear systems
ICCSA'05 Proceedings of the 2005 international conference on Computational Science and Its Applications - Volume Part IV
Journal of Computational Physics
Hi-index | 31.45 |
The emergence of ''petascale'' supercomputers requires us to develop today's simulation codes for (incompressible) flows by codes which are using numerical schemes and methods that are better able to exploit the offered computational power. In that spirit, we present a massively parallel high-order Navier-Stokes solver for large incompressible flow problems in three dimensions. The governing equations are discretized with finite differences in space and a semi-implicit time integration scheme. This discretization leads to a large linear system of equations which is solved with a cascade of iterative solvers. The iterative solver for the pressure uses a highly efficient commutation-based preconditioner which is robust with respect to grid stretching. The efficiency of the implementation is further enhanced by carefully setting the (adaptive) termination criteria for the different iterative solvers. The computational work is distributed to different processing units by a geometric data decomposition in all three dimensions. This decomposition scheme ensures a low communication overhead and excellent scaling capabilities. The discretization is thoroughly validated. First, we verify the convergence orders of the spatial and temporal discretizations for a forced channel flow. Second, we analyze the iterative solution technique by investigating the absolute accuracy of the implementation with respect to the different termination criteria. Third, Orr-Sommerfeld and Squire eigenmodes for plane Poiseuille flow are simulated and compared to analytical results. Fourth, the practical applicability of the implementation is tested for transitional and turbulent channel flow. The results are compared to solutions from a pseudospectral solver. Subsequently, the performance of the commutation-based preconditioner for the pressure iteration is demonstrated. Finally, the excellent parallel scalability of the proposed method is demonstrated with a weak and a strong scaling test on up to O(10^4) processing units and O(10^1^1) grid points.