A high-order cross-platform incompressible Navier-Stokes solver via artificial compressibility with application to a turbulent jet

File Description SizeFormat 
1-s2.0-S0010465518302248-main.pdfArticle In Press1.12 MBAdobe PDFView/Open
Title: A high-order cross-platform incompressible Navier-Stokes solver via artificial compressibility with application to a turbulent jet
Authors: Loppi, N
Witherden, F
Jameson, A
Vincent, PE
Item Type: Journal Article
Abstract: Modern hardware architectures such as GPUs and manycore processors are characterised by an abundance of compute capability relative to memory bandwidth. This makes them well-suited to solving temporally explicit and spatially compact discretisations of hyperbolic conservation laws. However, classical pressure-projection-based incompressible Navier–Stokes formulations do not fall into this category. One attractive formulation for solving incompressible problems on modern hardware is the method of artificial compressibility. When combined with explicit dual time stepping and a high-order Flux Reconstruction discretisation, the majority of operations can be cast as compute bound matrix–matrix multiplications that are well-suited for GPU acceleration and manycore processing. In this work, we develop a high-order cross-platform incompressible Navier–Stokes solver, via artificial compressibility and dual time stepping, in the PyFR framework. The solver runs on a range of computer architectures, from laptops to the largest supercomputers, via a platform-unified templating approach that can generate/compile CUDA, OpenCL and C/OpenMP code at runtime. The extensibility of the cross-platform templating framework defined within PyFR is clearly demonstrated, as is the utility of -multigrid for convergence acceleration. The platform independence of the solver is verified on Nvidia Tesla P100 GPUs and Intel Xeon Phi 7210 KNL manycore processors with a 3D Taylor–Green vortex test case. Additionally, the solver is applied to a 3D turbulent jet test case at , and strong scaling is reported up to 144 GPUs. The new software constitutes the first high-order accurate cross-platform implementation of an incompressible Navier–Stokes solver via artificial compressibility and -multigrid accelerated dual time stepping to be published in the literature. The technology has applications in a range of sectors, including the maritime and automotive industries. Moreover, due to its cross-platform nature, the technology is well placed to remain relevant in an era of rapidly evolving hardware architectures.
Issue Date: 25-Jun-2018
Date of Acceptance: 15-Jun-2018
URI: http://hdl.handle.net/10044/1/61494
DOI: https://dx.doi.org/10.1016/j.cpc.2018.06.016
ISSN: 0010-4655
Publisher: Elsevier
Journal / Book Title: Computer Physics Communications
Copyright Statement: © 2018 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
Sponsor/Funder: Engineering & Physical Science Research Council (EPSRC)
Engineering & Physical Science Research Council (EPSRC)
Funder's Grant Number: EP/K503381/1
Keywords: 01 Mathematical Sciences
02 Physical Sciences
08 Information And Computing Sciences
Nuclear & Particles Physics
Publication Status: Published online
Online Publication Date: 2018-06-25
Appears in Collections:Faculty of Engineering

Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.

Creative Commonsx