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.


'Towards Green Aviation with Python at Petascale' - Our simulations with PyFR on Piz Daint and Titan shortlisted for 2016 Gordon Bell Prize

'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

'Lifelines' - Our image of blood flow patterns in an arterio-venous fistulae wins prestigious BHF Reflections of Research award


Recent Papers

The Effect of Arterial Curvature on Blood Flow in Arterio-Venous Fistulae: Realistic Geometries and Pulsatile Flow. L. Grechy, F. Iori, R. Corbett, W. Gedroyc, N. Duncan, C. G. Caro, P. E. Vincent. Accepted for publication in Cardiovascular Engineering and Technology.
Abstract: Arterio-Venous Fistulae (AVF) are regarded as the 'gold standard' method of vascular access for patients with End-Stage Renal Disease (ESRD) who require haemodialysis. However, up to 60% of AVF do not mature, and hence fail, as a result of Intimal Hyperplasia (IH). Unphysiological flow and oxygen transport patterns, associated with the unnatural and often complex geometries of AVF, are believed to be implicated in the development of IH. Previous studies have investigated the effect of arterial curvature on blood flow in AVF using idealised planar AVF configurations and non-pulsatile inflow conditions. The present study takes an important step forwards by extending this work to more realistic non-planar brachiocephalic AVF configurations with pulsatile inflow conditions. Results show that forming an AVF by connecting a vein onto the outer curvature of an arterial bend does not, necessarily, suppress unsteady flow in the artery. This finding is converse to results from a previous more idealised study. However, results also show that forming an AVF by connecting a vein onto the inner curvature of an arterial bend can suppress exposure to regions of low wall shear stress and hypoxia in the artery. This finding is in agreement with results from a previous more idealised study. Finally, results show that forming an AVF by connecting a vein onto the inner curvature of an arterial bend can significantly reduce exposure to high WSS in the vein. The results are important, as they demonstrate that in realistic scenarios arterial curvature can be leveraged to reduce exposure to excessively low/high levels of WSS and regions of hypoxia in AVF. This may in turn reduce rates of IH and hence AVF failure.

High-Order Implicit Large-Eddy Simulations of Flow over a NACA0021 Aerofoil. J. S. Park, F. D. Witherden, P. E. Vincent. Accepted for publication in the AIAA Journal.
Abstract: In this study the GPU-accelerated solver PyFR is used to simulate flow over a NACA0021 aerofoil in deep stall at a Reynolds number of 270,000 using the high-order Flux Reconstruction (FR) approach. Wall-resolved Implicit Large Eddy Simulations (ILES) are undertaken on unstructured hexahedral meshes at fourth- and fifth-order accuracy in space. It was found that either modal filtering, or anti-aliasing via an approximate L2 projection, is required in order to stabilise simulations. Time-span averaged pressure coefficient distributions on the aerofoil, and associated lift and drag coefficients, are seen to converge towards experimental data as the simulation setup is made more realistic by increasing the aerofoil span. Indeed, the lift and drag coefficients obtained by fifth-order ILES with anti-aliasing via an approximate L2 projection agree better with experimental data than a wide range of previous studies. Stabilisation via modal filtering, however, is found to reduce solution accuracy. Finally, performance of various PyFR simulations is compared, and it is found that fifth-order simulations with anti-aliasing via an L2 projection are the most efficient. Results indicate that high-order FR schemes with anti-aliasing via an L2 projection are a good candidate for underpinning accurate wall-resolved ILES of separated, turbulent flows over complex engineering geometries.



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

Next Generation CFD: High-Order Accurate Simulations using Many-Core Platforms. P. E. Vincent. Swiss National Supercomputing Center, Lugano, Switzerland. August 2016.
PyFR: High-Order Accurate Cross-Platform Petascale Computational Fluid Dynamics with Python. F. D. Witherden, P. E. Vincent. NASA Ames, Moffett Field, CA, USA. May 2016.