Home

Welcome

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.

News

'Turbulent Channel Flow' - Checkout our latest paper on identifying eigenmodes of averaged small-amplitude perturbations to turbulent channel flow

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

Link:
'Implant may Offer Kidney Patients Easier Dialysis' - Our latest work on suppressing unsteady flow in arterio-venous fistulae featured in the Times

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

Link:

Recent Papers

PyFR v2.0.3: Towards Industrial Adoption of Scale-Resolving Simulations. F. D. Witherden, P. E. Vincent, W. Trojak, Y. Abe, A. Akbarzadeh, S. Akkurt, M. Alhawwary, L. Caros, T. Dzanic, G. Giangaspero, A. S. Iyer, A. Jameson, M. Koch, N. Loppi, S. Mishra, R. Modi, G. Saez-Mischlich, J. S. Park, B. C. Vermeire, L. Wang. Accepted for publication in Computer Physics Communications.
Abstract: PyFR is an open-source cross-platform computational fluid dynamics framework based on the high-order Flux Reconstruction approach, specifically designed for undertaking high-accuracy scale-resolving simulations in the vicinity of complex engineering geometries. Since the initial release of PyFR v0.1.0 in 2013, a range of new capabilities have been added to the framework, with a view to enabling industrial adoption. In this work, we provide details of these enhancements as released in PyFR v2.0.3, including improvements to cross-platform performance (new backends, extensions of the DSL, new matrix multiplication providers, improvements to the data layout, use of task graphs) and improvements to numerical stability (modal filtering, anti-aliasing, artificial viscosity, entropy filtering), as well as the addition of prismatic, tetrahedral and pyramid shaped elements, improved domain decomposition support for mixed element grids, improved handling of curved element meshes, the addition of an adaptive time-stepping capability, the addition of incompressible Euler and Navier-Stokes solvers, improvements to file formats and the development of a plugin architecture. We also explain efforts to grow an engaged developer and user community and provided a range of examples that show how our user base is applying PyFR to solve a wide range of fundamental, applied and industrial flow problems. Finally, we demonstrate the accuracy of PyFR v2.0.3 for a supersonic Taylor-Green vortex case, with shocks and turbulence, and provided latest performance and scaling results on up to 1024 AMD Instinct MI250X accelerators of Frontier at ORNL (each with two GCDs) and up to 2048 Nvidia GH200 GPUs of Alps at CSCS. We note that absolute performance of PyFR accounting for the totality of both hardware and software improvements has, conservatively, increased by almost 50x over the last decade.

Link:
Effects of Free-Stream Eddies on Triangular Airfoils for Martian Rotorcraft. L. Caros, O.R.H. Buxton, P. E. Vincent. Accepted for publication in AIAA Journal.
Abstract: This study uses Direct Numerical Simulations (DNS) with PyFR (www.pyfr.org) to investigate the impact of free-stream disturbances on the performance of non-conventional triangular airfoils for Martian rotorcraft. First, triangular airfoils optimized for steady uniform free- stream conditions at a 12 degree angle of attack are subjected to unsteady upstream conditions containing free-stream eddies with various target turbulence intensities and characteristic length scales. Results show that the free-stream eddies lead to a breakdown of coherent suction-side vortices and a reduction of aerodynamic forces, with the magnitude of this effect increasing with higher turbulence intensities and larger length scales. Building on this, DNS are then used to optimize triangular airfoils at a 12 degree angle of attack for unsteady upstream conditions containing free-stream eddies. The new designs still leverage an unsteady vortex roll-up mechanism to generate lift. However, their shapes are different to those obtained when optimizing for steady uniform inflow conditions, with their apexes closer to the trailing edge. This downstream apex position allows for the smaller and weaker vortices that result from the interaction of the free-stream eddies with the separated shear layer to cover a larger suction surface area. Taken together, the results demonstrate the importance of considering free-stream disturbances when designing non-conventional airfoils for use on a Martian rotorcraft.

Link:

Recent Seminars

PyFR: Taking Scale-Resolving Simulations from Academia to Industry. P. E. Vincent NASA Ames, Moffett Field, CA, USA, December 2024.
Aerodynamic Optimisation of Aerofoils for Martian Rotorcraft Using Direct Numerical Simulations. P. E. Vincent AIAA Journal Seminar Series, October 2024.

Openings

PhD Studentship in Aeronautics - High-Fidelity Simulation of Titan/Mars Entry Vehicles with PyFR
Summary: Next-generation Entry, Descent, and Landing (EDL) systems for Titan and Mars must safely slow down increasingly large payloads. One particular challenge occurs during the transonic phase of descent, where the spacecraft is subject to aerodynamic instabilities that can cause uncontrolled oscillations, posing a significant risk of mission failure. This project will further develop the GPU-accelerated computational fluid dynamics flow solver PyFR - implementing improved shock capturing approaches and a full 6-DOF free-flight capability - and use it to study dynamic stability in the transonic phase of descent. The work will be undertaken in collaboration with Texas A&M University and NASA Ames.

Advert: