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

Identifying Eigenmodes of Averaged Small-Amplitude Perturbations to Turbulent Channel Flow. A. S. Iyer, F. D. Witherden, S. I. Chernyshenko, P. E. Vincent. Accepted for publication in Journal of Fluid Mechanics.
Abstract: Eigenmodes of averaged small-amplitude perturbations to a turbulent channel flow - which is one of the most fundamental canonical flows - are identified for the first time via an extensive set of high-fidelity GPU-accelerated direct numerical simulations. While the system governing averaged small-amplitude perturbations to turbulent channel flow remains unknown, the fact such eigenmodes can be identified constitutes direct evidence that it is linear. Moreover, while the eigenvalue associated with the slowest-decaying anti-symmetric eigenmode mode is found to be real, the eigenvalue associated with the slowest-decaying symmetric eigenmode mode is found to be complex. This indicates that the unknown linear system governing the evolution of averaged small-amplitude perturbations cannot be self-adjoint, even for the case of a uni-directional flow. In addition to elucidating aspects of the flow physics, the findings provide guidance for development of new unsteady Reynolds-averaged Navier-Stokes turbulence models, and constitute a new and accessible benchmark problem for assessing the performance of existing models, which are used widely throughout industry.

Optimal Runge-Kutta Schemes for Pseudo Time-Stepping with High-Order Unstructured Methods. B. C. Vermeire, N. A. Loppi, P. E. Vincent. Journal of Computational Physics, Volume 383, Pages 55-71, 2019.
Abstract: In this study we generate optimal Runge-Kutta (RK) schemes for converging the Artificial Compressibility Method (ACM) using dual time-stepping with high-order unstructured spatial discretizations. We present optimal RK schemes with between s=2 and s=7 stages for Spectral Difference (SD) and Discontinuous Galerkin (DG) discretizations obtained using the Flux Reconstruction (FR) approach with solution polynomial degrees of k=1 to k=8. These schemes are optimal in the context of linear advection with predicted speedup factors in excess of 1.80x relative to a classical RK44 scheme. Speedup factors of between 1.89x and 2.11x are then observed for incompressible Implicit Large Eddy Simulation (ILES) of turbulent flow over an SD7003 airfoil. Finally, we demonstrate the utility of the schemes for incompressible ILES of a turbulent jet, achieving good agreement with experimental data. The results demonstrate that the optimized RK schemes are suitable for simulating turbulent flows and can achieve significant speedup factors when converging the ACM using dual time-stepping with high-order unstructured spatial discretizations.



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.