We develop novel numerical methods and apply them to solve challenging fluid flow problems in various areas of science, engineering, and medicine. We are particularly interested in theoretical aspects of high-order numerical methods for unstructured grids, as well as their implementation for a range of modern hardware platforms.


'Step Inside a Jet Engine' - Results from our latest PyFR simulations of flow over low pressure turbine blades on show at the Imperial Fringe

'Implant may offer kidney patients easier dialysis' - Our latest work on suppressing unsteady flow in arterio-venous fistulae featured in the Times

'New Symmetric Quadrature Rules' - Checkout our latest paper on identification of symmetric quadrature rules for finite element methods

'Analysis of Tetrahedral Solution Points' - Checkout our latest paper on solution point placement for Flux Reconstrustion schemes on tetrahedra


Recent Papers

Heterogeneity in the Non-Planarity and Arterial Curvature of Arteriovenous Fistulae In Vivo. R. Corbett, L. Grechy, F. Iori, J. S. Crane, P. E. Herbert, P. Di Cocco, W. Gedroyc, P. E. Vincent, C. G. Caro, N. D. Duncan. Accepted for publication in the Journal of Vascular Surgery.
Abstract: Native arteriovenous fistulas (AVFs) for hemodialysis are susceptible to nonmaturation. Adverse features of local blood flow have been implicated in the formation of perianastomotic neointimal hyperplasia that may underpin nonmaturation. Whereas computational fluid dynamic simulations of idealized models highlight the importance of geometry on fluid and vessel wall interactions, little is known in vivo about AVF geometry and its role in adverse clinical outcomes. This study set out to examine the three-dimensional geometry of native AVFs and the geometric correlates of AVF failure.

A High-Order Cross-Platform Incompressible Navier-Stokes Solver via Artificial Compressibility with Application to a Turbulent Jet. N. A. Loppi, F. D. Witherden, A. Jameson, P. E. Vincent. Accepted for publication in Computer Physics Communications.
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 P-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 Re=10,000, 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 P-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.



PhD Studentship - Development of In-situ Visualisation and Analysis Technology for High-Fidelity Computational Fluid Dynamics
Summary: A PhD Studentship is currently available. The project, will involve addition of 'in-situ' visualisation, processing, and analysis technology to PyFR, an open-source high-order massively-parallel CFD platform, as well as its application to solve a range of challenging unsteady flow problems. Candidates should hold, or expect to obtain, an undergraduate degree in a numerate discipline. Previous programming experience is important (ideally Python, C++ and CUDA).


Recent Seminars

Towards Green Aviation with Python at Petascale. P. E. Vincent. Tokyo University of Science, Tokyo, Japan. December 2017.
Next Generation CFD: High-Order Accurate Simulations using Many-Core Platforms. P. E. Vincent. Swiss National Supercomputing Center, Lugano, Switzerland. August 2016.