Extracting coherent structures and low-order modelling using Optimal Mode Decomposition

Lead Research Organisation: Imperial College London
Department Name: Aeronautics

Abstract

Since nearly 25% of UK emissions are generated by transport, improving the energy efficiency of aviation and road vehicles is vital if the UK is to meet legally binding targets to significantly reduce greenhouse gas emissions by 2050. One promising approach to achieving this aim is to improve the aerodynamic performance of ground and air vehicles using active flow control. In active flow control, measurements are taken from the flow in real time and, in response to those measurements, action is taken to positively modify the flow. For example, movable surfaces on an aircraft's wing or on the rear of a road vehicle may be dynamically adjusted to reduce aerodynamic drag and hence decrease fuel consumption.

To implement an active control strategy typically requires an approximating flow model, which allows a beneficial control action to be computed from a given flow measurement. However, a significant challenge is that flows of practical interest have highly complicated and intricate dynamics, while computational and theoretical restrictions mean that the approximating models must be simple, or low-dimensional. One method of bridging this gap is to identify coherent structures in a flow, that is, spatial features which are of dynamical importance to its evolution and performance. If a small set of representative coherent structures can be identified, they can be used as the building blocks for a flow model and therefore facilitate a successful flow control strategy.

The aim of this project is to develop the Optimal Mode Decomposition (OMD) algorithm, a method which systematically extracts dynamically important coherent structures from ensembles of experimental or numerical fluid flow data. The systematic nature of the data extraction is important since it provides a tractable method of analysing the increasing large-scale sets of fluid flow data that are now available. The proposed research will thoroughly benchmark the performance of the OMD algorithm across a range of fundamental and challenging fluid flows, assess the suitability of the extracted structures to be used for flow modelling and estimation and, finally, undertake a theoretical study to explain any observed behaviour. This research will therefore develop a methodology to underpin the application of flow control for improved energy efficiency and consequently addresses a current and fundamentally important environmental and economic issue.

Planned Impact

The transportation sector generates roughly 20% of EU countries' CO2 emissions. Although emissions from other industrial sectors, such as energy generation, are expected to decrease over the coming decades, a predicted increase in demand means this is not the case for transportation. Consequently, it is of great importance to improve the fuel efficiency of, for example, road vehicles and aviation in order to address key climate change targets. Since aerodynamic performance is fundamental to the fuel and energy efficiency of road and air vehicles, understanding and controlling the behaviour of fluid flows presents an important, but technically challenging, pathway towards achieving this aim.

This research seeks to address this challenge by developing a methodology - the Optimal Mode Decomposition algorithm - that will help researchers transition from the ability to accurately measure or observe a fluid flow, to being able to model, predict and improve its behaviour. Therefore, the impact of this research will be to produce a tool that enables the application of flow control technologies that can be used to improve fuel efficiency in transportation. This promises economic benefits through reduced fuel consumption and associated cost savings. More importantly, improvements in fuel efficiency will have environmental and societal benefits ranging from direct improvements to air quality to the vital medium-term impact of helping to mitigate the effects of climate change. Outcomes of the proposed research will be presented at an industry-sponsored workshop towards to highlight the developed methodologies outside of the academic community, and also at the Imperial College Festival, a major public outreach event, in order to motivate flow modelling and control as a fascinating and practically important challenge to the next generation of researchers.
 
Description The fundamental purpose of this grant is to develop a computational method, the Optimal Mode Decomposition (OMD) algorithm, which facilitates the analysis of complex, high-dimensional data arising from fluid mechanics experiments and simulations. It extracts coherent structures which are the building blocks for low-dimensional models of fluid flows that can be used to design control strategies to improve flow efficiency.

The main findings of the research project to date are:

Benchmarking (Work Package 1): The OMD algorithm has been shown to be able to extract known coherent structures from fluid flows with a range of dynamic complexity and has been shown to outperform (in terms of correctly identifying the underlying system's frequency and growth rate) the widely employed Dynamic Mode Decomposition Algorithm when non-linearity is present in the underlying dynamics (employing a standard Stuart-Landau model) and for the canonical flow of the cylinder wake. This evidence suggests that the OMD algorithm has the potential to generate more accurate flow models than possible using existing methodologies.

Application to Noisy Data (Work Package 2): The situation was investigated in which the underlying data ensemble was corrupted by assuming either spatial or temporal down-sampling. This corresponds to analyzing flow data in which information may be unavailable at certain spatial locations in the flow or that the flow could not be sampled at as high a temporal rate as would have been desired. Both situations are experienced in experimental practice. The first key finding is that the OMD algorithm is effective when applied to spatially down-sampled data. Conversely, too aggressive temporal down-sampling (past the Nyquist criterion) was found to produce systematic errors in the extracted coherent structures and their associated dynamics. Subsequently, a key finding was the OMD algorithm could be modified to mitigate the effects of temporal down-sampling. This resulted in the development of a new algorithm, irregularly-sampled OMD (is-OMD), which has been demonstrated to perform well on both carefully controlled flows and highly-complex experimental flow data. Furthermore, new theoretical results were developed, proving that the is-OMD algorithm has rigorous convergence properties. This work has recently been published in the Journal of Fluid Mechanics (Krol J, Wynn A, 2017, Dynamic reconstruction and data reconstruction for subsampled or irregularly sampled data, Journal of Fluid Mechanics, 825, pp 133-166).

Use of OMD modes in flow estimation (Work Package 3): Being able to infer the velocity field of a flow using only information on its boundary is of fundamental practical importance: if one knows what a flow is doing, this opens the door to improving its behaviour and measurements can only typically be taken on the boundary of a flow due to physical limitations on the placement of sensors. Turbulent flow past an axisymmetric bluff body was analysed using the OMD algorithm and it was found that certain coherent structures identified by OMD could be determined from pressure sensors on the base of the body using stochastic estimation. It has been discovered that OMD modes perform better at extracting wake structures whose temporal dynamics have more prominent dominant frequencies than modes extracted using standard methodologies (e.g. Proper Orthogonal Decomposition). This temporal prominence is beneficial when constructing low-order flow models from the modes for use in flow estimation. This work is ongoing and is being taken forward by the PhD student funded by Imperial College in support of this grant.

A further development currently being pursued by the PhD student is to use OMD to provide spare approximate models for turbulent flows. Sparse identification is useful and necessary, since modal flow decomposition methods typically produce a large number of flow structures, only a few of which may be of interest to the user (or mathematical modeller). We have developed a new addition to the OMD algorithm which uses the graph-theoretical notion of Maximal Cliques to collect together similar modes, and subsequently represent each collection by a single flow structure. By applying this method iteratively, we are able to vastly reduce the number of degrees of freedom necessary to represent important features in a turbulent flow. Preliminary results was published and presented at the 11th International Symposium on Turbulence and Shear Flow Phenomena (TSFP11) in Southampton, UK, in July 2019 (http://www.tsfp-conference.org/proceedings/2019/103.pdf). Subsequently, this work was invited to be submitted for publication as a journal article as part of a special issue of the International Journal of Heat and Fluid Flow. This work has now been published as a journal article (Beit Sadi M, Krol J, Wynn A, Data-driven feature identification and sparse representation of turbulent flows, International Journal of Heat and Fluid Flow, vol 88, 2021, 108766).
Exploitation Route A new algorithm, is-OMD, has been developed which is an extension of the OMD algorithm allowing it to be applied to temporally downsampled data. This can be implemented by other researchers and practitioners to analyse fluid flow data. The academic work initiated by this grant is being taken forward in the work of a PhD student hired with supporting funding for the grant.

Recently, the research supported by this grant underpinned a successful funding application to support a forthcoming project whose aim is to enable active drag reduction for turbulent boundary layers ([EnAble]: Developing and Exploiting Intelligent Approaches for Turbulent Drag Reduction, EPSRC, EP/T021144/1, £1.35m FEC). The purpose of the grant will be to use computational optimization to discover novel forcing mechanisms for power-efficient drag reduction for boundary layer flows. Optimization of the flow control strategy will be achieved using a Bayesian optimization approach, with Optimal Mode Decomposition used to analyse the modifications made to flow structures by the control: this will allow new physical mechanisms for drag reduction to be understood and discovered. This project is in collaboration with Airbus who have expressed an interest in using any successful strategies developed in the proposal for to enable drag reduction (and, hence, increase fuel-efficiency) on aircraft wings and fuselages.
Sectors Aerospace, Defence and Marine,Transport

 
Description The Optimal Model Decomposition (OMD) algorithm is currently being used in work supported by the leading aircraft manufacturer Airbus to help understand a new method for aircraft drag reduction. Compelling evidence suggests that injecting low-velocity fluid into a turbulent boundary layer can reduce its friction drag and enable significant fuel savings and emissions reductions. Designing the optimal version of this flow control strategy requires an understanding of the key physical changes that boundary layer forcing has on a fluid, and the OMD algorithm is allowing us to understand these changes. While this work is currently at a low Technology Readiness Level, there is a concrete pathway to application in the aerospace sector through Airbus's involvement. Wider impact of this work in the aerospace sector has also been achieved via the PI's invitation to lecture engineers from the Aerospace manufacturer COMAC on flow control and modal decomposition methods in fluid mechanics. Popularisation of the OMD algorithm achieved though this project has led to the OMD algorithm being widely used and well-cited (now 139 times) in the fluid mechanics community. This success has allowed the central concept of the algorithm, that a flow's evolution over a short time can be approximated by a subspace projection followed by a low-order dynamical mapping, to be generalised by the fluids community to fit within the framework of machine learning (where the projection steps have been generalised to auto-encoder and decoder neural networks). This means that OMD and related algorithms can take advantage of the latest advances in machine intelligence and may therefore expect to have continuing and lasting impact in both fluid mechanics and data-science.
First Year Of Impact 2021
Sector Aerospace, Defence and Marine
Impact Types Economic

 
Description PhD Studentship, Imperial College London
Amount £63,000 (GBP)
Organisation Imperial College London 
Department Department of Aeronautics
Sector Academic/University
Country United Kingdom
Start 10/2017 
End 03/2021
 
Description [EnAble]: Developing and Exploiting Intelligent Approaches for Turbulent Drag Reduction
Amount £498,547 (GBP)
Funding ID EP/T021144/1 
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 10/2020 
End 09/2023
 
Description Presentation to Siemens-Gamesa Renewable Energy 
Form Of Engagement Activity A formal working group, expert panel or dialogue
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Industry/Business
Results and Impact During 2018, I gave presentations to the wind turbine manufacturer Siemens-Gamesa Renewable Energy concerning the potential for flow control and modelling techniques to improve efficiency for offshore wind turbine arrays. The Optimal Mode Decomposition algorithm formed a key part of the approach, facilitating the modelling procedure upon which the proposed techniques were based.

The meeting led to Siemens-Gamesa supporting a potential future collaboration via an EPSRC-funded research project ((EP/S036873/1, in review, Siemens provided a letter of support) co-written with a colleague at Imperial College, Dr. O. Buxton. If successful, we will work with field data from Siemens to improve models which predict wind-turbine wakes (and therefore provide the tools to enable greater efficiency of wind farms).
Year(s) Of Engagement Activity 2018
 
Description Technical Meeting with Airbus to discuss flow control strategies. 
Form Of Engagement Activity A formal working group, expert panel or dialogue
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Industry/Business
Results and Impact A representative of Airbus, Germany, visited Imperial College in January 2020 to discuss Airbus' involvement in a forthcoming EPSRC-funded grant to research new methods of drag reduction in turbulent boundary layers. (EP/T021144/1). The research in the optimal mode decomposition algorithm supported by this project will feed into the new grant as a tool to understand the physics underpinning any new flow control mechanisms discovered as part of the project. Airbus's involvement provides a possible pathway to implementation of this research.
Year(s) Of Engagement Activity 2020
 
Description Training sessions for engineers from Chinese national aerospace manufacturer COMAC 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Industry/Business
Results and Impact A training session on modern methods of Flow Control, detailing the data-driven modelling approaches developed in this grant, was delivered to senior engineers from the aircraft manufacturer COMAC.
Year(s) Of Engagement Activity 2020
 
Description Visit to Vestas Technology, Isle of Wight, December 2017 
Form Of Engagement Activity A formal working group, expert panel or dialogue
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Industry/Business
Results and Impact A visit was made to Vestas Technology's manufacturing and research plant on the Isle of Wight in December 2017. Vestas Technology is a leading wind-turbine manufacturer and I, along with another colleague from Imperial College London Dr Oliver Buxton, informed them of potential applications of flow control technology that could be used to make their turbines operate more efficiently. A component of this proposal involved flow modelling and analysis using the Optimal Mode Decomposition algorithm developed as part of this EPSRC proposal. The visit has led to ongoing discussions regarding future work between ourselves and Vestas.
Year(s) Of Engagement Activity 2017