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

3D Flow Reconstruction Using Divergence Free Interpolation of Multiple 2D Contrast Enhanced Ultrasound Imaging Velocimetry Measurements. X. Zhou, V. Papadopoulou, C. H. Leow, P. E. Vincent, M. X. Tang. Accepted for publication in Ultrasound in Medicine and Biology.
Abstract: Quantification of 3D intravascular flow is valuable for studying arterial wall diseases but currently there is a lack of effective clinical tools for this purpose. Divergence Free Interpolation (DFI) using Radial Basis Function (RBF) is an emerging approach for full field flow reconstruction using experimental sparse flow field samples. Previous DFI reconstructs full field flow from scattered 3D velocity input obtained using Phase Contrast MRI with low temporal resolution. In this study, a new DFI algorithm is proposed to reconstruct full field flow from scattered 2D in-plane velocity vectors obtained using ultrafast contrast enhanced ultrasound (>1000 fps) and particle imaging velocimetry (Ultrasound PIV, or UIV). The full 3D flow field is represented by a sum of weighted divergence free RBFs in space. Due to the acquired velocity vectors being only in 2D and hence the problem being ill-conditioned, a regularized solution of the RBF weighting is achieved through Singular Value Decomposition (SVD) and the L-curve method. The effectiveness of the algorithm is demonstrated via numerical experiments for Poiseuille flow and helical flow with added noise, and it is shown that an accuracy of up to 95.6% can be achieved for Poiseuille flow (with 5% input noise). Experimental feasibility is also demonstrated by reconstructing full-field 3D flow from experimental 2D UIV measurements in a carotid bifurcation phantom. The method is typically faster for a range of problems compared to computational fluid dynamics, and has been demonstrated to be effective for the three flow cases.

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. Journal of Vascular Surgery, Volume 68, Issue 6, Pages 152S-163S, 2018.
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.



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.