Imperial College London

Professor Peter Vincent

Faculty of EngineeringDepartment of Aeronautics

Professor of Computational Fluid Dynamics



+44 (0)20 7594 1975p.vincent




211City and Guilds BuildingSouth Kensington Campus






BibTex format

author = {Loppi, N and Witherden, F and Jameson, A and Vincent, PE},
doi = {10.1016/j.cpc.2018.06.016},
journal = {Computer Physics Communications},
pages = {193--205},
title = {A high-order cross-platform incompressible Navier-Stokes solver via artificial compressibility with application to a turbulent jet},
url = {},
volume = {233},
year = {2018}

RIS format (EndNote, RefMan)

AB - 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. Moreove
AU - Loppi,N
AU - Witherden,F
AU - Jameson,A
AU - Vincent,PE
DO - 10.1016/j.cpc.2018.06.016
EP - 205
PY - 2018///
SN - 0010-4655
SP - 193
TI - A high-order cross-platform incompressible Navier-Stokes solver via artificial compressibility with application to a turbulent jet
T2 - Computer Physics Communications
UR -
UR -
VL - 233
ER -