📣 Help Shape the Future of UKRI's Gateway to Research (GtR)

We're improving UKRI's Gateway to Research and are seeking your input! If you would be interested in being interviewed about the improvements we're making and to have your say about how we can make GtR more user-friendly, impactful, and effective for the Research and Innovation community, please email gateway@ukri.org.

Asymptotic approximation of the large-scale structure of turbulence in axisymmetric jets: a first principle jet noise prediction method

Lead Research Organisation: University of York
Department Name: Mathematics

Abstract

Ever since the jet age began in the 1950s, governments, scientists, and engineers have been acutely aware of the health effects created by aircraft noise--the prolonged exposure of which is highly damaging to human health. Increased noise pollution, for example, has been linked to cognitive impairment and behavioural issues in children, sleep disturbance (and consequent health issues therefrom) as well as the obvious hearing damage caused by the repeated intrusion of high levels of noise. The World Health Organization estimates that 1-million healthy life years are lost in Europe due to noise; this is mainly by cardiovascular disease via the persistent increase in stress level-with aviation noise being the largest contributor here. Moreover, the Aviation Environment Federation found that these issues place a £540M/year burden on UK government expenditure. While there has been tremendous progress in understanding aircraft noise, the doubling of flights in the past 20 years to a staggering 40 million (in the pre-Covid year 2019) has heightened the need for research into the physics of jet noise to uncover new reduced-order turbulence models. This proposal develops a novel mathematical model for jet flow turbulence using asymptotic analysis. The re-constructed turbulence structure will be used within a numerical code for fast noise prediction of a high-speed axisymmetric jet flow.

Fundamentally, a jet flow breaking down into turbulence creates pressure fluctuations that propagate away as sound. In 1952, Lighthill showed that the Navier-Stokes equations can be exactly re-arranged into a form where a wave operator acting on the pressure fluctuation, is equal to the double-divergence of the jet's Reynolds stress. When the auto-covariance of the Reynolds stress was assumed to be known for a fluid at rest, scaling properties of the acoustic spectrum were obtained such as the celebrated 8th power law. The generalized acoustic analogy formulated by Goldstein in 2003 advanced this idea by dividing the fluid mechanical variables into a steady base flow and its perturbation. The acoustic spectrum per unit volume is a tensor product of a propagator and the auto-covariance of the purely fluctuating Reynolds stress tensor. The propagator can be calculated by determining the Green's function of the Linearized Euler operator for an appropriate jet base flow however, as in Lighthill's theory, the auto-covariance tensor is assumed to be known, which invariably requires the use of Large-Eddy Simulation (LES) and experiments to obtain an approximate functional form for it. But LES data still uses immense computational resources and computing time when different nozzle operating points are needed for design optimization or when complex jets are considered.

What makes any alternative to modelling so complex is that the turbulence closure problem precludes a closed-form theory for the auto-covariance tensor. However, our recent work revealed that the peak noise can be accurately predicted when the propagator is determined at low frequencies that are of the same order as the jet spread rate (that is lesser than unity). This proposal, therefore, sets out an alternative, first-of-its-kind, analytical approach to determine the fluctuating Reynolds stress for a given mean flow solution. By solving the governing equations at this asymptotic scaling where the jet evolves temporally at the same rate it spreads in space, we determine the Large-Scale Turbulence (LST) structure in the jet. This approach is defined by a 2-dimensional system of equations for an axisymmetric jet and the computational time is expected to be an order-of-magnitude faster than LES. The LST-based solution of the Reynolds stress auto-covariance for peak jet noise will be compared to LES data provided by our project partners at several jet operating conditions. We aim to show that the LST model of turbulence provides accurate noise predictions and is a viable alternative to LES.

Publications

10 25 50
 
Description The aim of this project is to create a model for unsteady Reynolds stress source term by solving approximate forms of the nonlinear Euler equations.

For WP 1.1, the PI (Dr. Koshuriyan) derived the form of the non-linear Euler equations that we will use in WP2 & 3. An important decision we made was to keep the algebraic form of the non-linear terms on the right side of the Euler equations as simple as possible. That is, to retain only the gradient of the fluctuating Reynolds stress on the right hand side of the momentum equation. This approximation is accurate for jet noise prediction. See:

S. Stirrat, 2024, Unified approach to jet installation noise using mathematical modelling and spectral analysis of jet turbulence. PhD thesis supervised by Koshuriyan.

The fluctuating Reynolds stress is a stationary random function of space & time. Hence the "initial guess" for Reynolds stress cannot be deterministic and must also have this stationary random property. Therefore in WP1, we aimed to identify key spatial/temporal structures within the jet. This work was led by PDRA, Dr. Sadaf Arabi who spent the first year of the project learning different statistical methods and data-driven techniques.

One widely used data-driven technique that Dr. Arabi used is called Spectral Proper Orthogonal Decomposition (SPOD), which identifies the most energetic flow modes by constructing an optimal orthogonal basis that reduces the complexity of large temporally evolving datasets while preserving the essential features of the flow field.

For WP1.2, Dr. Sadaf Arabi used the SPOD algorithm to at Mach numbers (M = 0.5 & 0.9) using previously validated LES data (EP/R513349/136). The key findings are as follows:

1). All velocity components exhibit coherent structures in the unsteady flow field for both subsonic jets with similar energy content for all 3 components.

2). A small number of leading SPOD modes (less than 3) are sufficient to reconstruct the fluctuating Reynolds stress

3). Truncating the flow field to the potential core region (where the turbulence amplitude is greatest) allows a smaller number of modes to be retained in the reconstruction. These modes appear as coherent wavepackets.

See: Investigation of Orthogonal Modes in High-speed Axisymmetric Jet Turbulence: Instantaneous vs. Statistical Data (Arabi-Koshuriyan-Sescu). International Conference on Flow Dynamics, Nov-2024, Japan and journal paper (in prep) for Physics of Fluids.

Prior to the investigation in WP2, Dr Sadaf Arabi undertook a visit to the Polytechnic University of Madrid to use the linear Parabolized Stability code written by Prof. Rodriguez and to interact with his research team. Parabolized Stability approaches use linear stability theory to predict the growth rate of wave packet coherent structures present in the unsteady flow. Our preliminary results show that modal energy distribution obtained by linear PSE is similar to the fully non-linear SPOD calculation. This indicates that nonlinearity makes small contribution to the space/time coherent structure contained within the Reynolds stress.

PhD student Fay Oglethorpe has conducted preliminary work on WP2.1 by working out the coefficients of the nonlinear Euler equations using LES flow data.
Exploitation Route We hope that the numerical codes based on our low-order turbulence flow approximation can be used in both industry and academia in the following ways:

(1). Conduct fast jet noise prediction in isothermal (non-heated flows)

(2). Extend the analytical approach to more realistic heated flows.

(3). Within more complex aero-acoustic research problems - i.e. flows interacting with surfaces such as in the trailing edge noise problem.

(4). To conduct further research into the applicability of asymptotic approximations in turbulence.
Sectors Aerospace

Defence and Marine

URL https://eprints.whiterose.ac.uk/222207/1/ICFD_2024_S_POD_Arabi_Koshuriyan_Sescu_2024_.pdf
 
Title Orthogonal mode eduction from high speed jet flow turbulence data 
Description Spectral proper orthogonal decomposition is used to determine the energy content in unsteady flow data. It does this by expanding the flow field input into orthogonal modes. The spatial structure of these modes are in the form of coherent wave packets. This work fits into WP1 of project plan because the orthogonal modes can be used to determine a low-order model for the unsteady fluctuating Reynolds stress source term that is responsible for jet noise. The dataset is therefore the complete orthogonal mode spectral (i.e. the Fourier transform of the original input unsteady flow signal) data for the unsteady flow that generates the Reynolds stress and the orthogonal mode reconstruction for the Reynolds stress itself. 
Type Of Material Computer model/algorithm 
Year Produced 2024 
Provided To Others? Yes  
Impact The spectral orthogonal decomposition allows a low order reconstruction of the Fourier transform of the unsteady flow data to be obtained. This reconstruction is in terms of the spectral energy modes that appear are coherent wavepackets. Our work shows how many modes are needed to adequately reconstruct the spectral energy distribution of the fluctuating Reynolds stress source term. Since the latter appears on the right hand side of the non-linear Euler equations, a reduced order model can be used as a starting guess for an iterative numerical solution that determines it. 
URL https://eprints.whiterose.ac.uk/222207/1/ICFD_2024_S_POD_Arabi_Koshuriyan_Sescu_2024_.pdf
 
Description Collaboration with Prof. Christophe Bogey (Ecolé Centralé de Lyon) 
Organisation Ecole Centrale de Lyon
Country France 
Sector Academic/University 
PI Contribution The collaboration with Prof. Bogey is an exciting opportunity bringing together applied mathematical modelling of jet noise that we have pioneered in my group with his expert Large-Eddy Simulation database. The Bogey datasets were obtained over the past 15 years and this represents an exciting opportunity to prove the robustness of our current approaches and also how the new project on turbulence reconstruction compares over a wide range of inflow conditions that the Bogey dataset has. Currently we have a plan in place on the order of calculations (i.e. which database to start from). My PhD student (Fay Oglethorpe) will do this as part of her doctoral work. She will begin her PhD in September 2024.
Collaborator Contribution Prof. Bogey has provided initial test data and code to begin the post-processing for the first jet in our plan for analysis of Bogey datasets. Fay Oglethorpe will use the initial data and code to familiarise herself with the Bogey database. This is scheduled to begin in September/October 2024.
Impact Currently underway. We are aiming for the November 2024 abstract submission for the American Institute of Aeronautics & Astronautics (2025) Conference.
Start Year 2024
 
Description Collaboration with Prof. Dani Rodriguez (Polytechnic University of Madrid) 
Organisation Technical University of Madrid
Country Spain 
Sector Academic/University 
PI Contribution The potential collaboration between PI (Koshuriyan) was established during the American Institute of Aeronautics & Astronautics (AIAA) Aeroacoustics meeting in June 2022. The low order asymptotic model of the turbulence velocity field (and subsequent jet noise source term reconstruction) that we will construct in this project will be compared to the so-called parabolized stability equations (PSE) that Prof. Rodriguez has focused a large part of his career on. Prof. Rodriguez has agreed to send us the codes he uses for the parabolized stability calculation. This is important because it will shows two things: (a). how our first-principles asymptotic theory compares against an ad-hoc engineering approach (b). where our approach fits within the realm of fully numerical turbulence simulations on the one hand and ad-hoc approaches like PSE. We have been conducting monthly Zoom meetings with Prof. Rodriguez and his Post-doc (Dr Ivan Padilla Montero), which I have organised. The project PDRA (Dr. Sadaf Arabi) visited the Polytechnic University of Madrid for 9 days in December 2024. I help organised this visit. The outcomes of the research undertaken are summarised below.
Collaborator Contribution Professor Rodriguez has guided Sadaf towards understanding his linear PSE code. This involved working through the derivation of the formulation and using his codes to work out the maximum energy modes of the flow turbulence. Dr. Sadaf Arabi began using his codes in months prior to her visit to Madrid. During the visit itself (10th-19th December). Dr. Montero guided her in using the linear PSE code. Since the code was written by Professor Rodriguez/Dr. Ivan Montero, they modified it appropriately in order to handle the existing LES mean flow fields we have. The latter were obtained in the EPSRC DTP Grant (ref: EP/R513349/136).
Impact This has proved to be a successful collaboration. We have used the PSE code to show how linear the maximum energy modes in the turbulence is. We are aiming to present a paper at AIAA Scitech 2026 based on this collaboration. Moreover, a journal paper at Physics of Fluids will be submitted by the end of the year. This paper is interdisciplinary in that it will involve theory developed at York showing how PSE is related to S-POD and the numerical computation that we performed using the Madrid PSE code.
Start Year 2024
 
Description Research Collaboration with Tohoku University (Institute for Fluid Science) 
Organisation Tohoku University
Country Japan 
Sector Academic/University 
PI Contribution The collaboration with Prof. Yuji Hattori's group at the Institute for Fluid Science (IFS) at the University of Tohoku was initiated in 2017 by my collaborator Prof. Sescu. Each year since then we have visited or remotely presented. When there is mutual interest between IFS group headed by Prof. Hattori and our research papers, we have collaborate further or invistigated possible avenues for joint research proposals. In relation to my current EPSRC project, the research collaboration has resulted in my PDRA (Dr Sadaf Arabi) visiting Tohoku to present a paper at the 21st International Conference on Flow Dynamics. This was fully funded as described below.
Collaborator Contribution In relation to the current EPSRC project, the collaboration with IFS Tohoku has resulted in a Refereed Conference paper (see below). In context to that, the contributions made are as follows: Adrian Sescu -- modification of codes to post-process Large-Eddy Simulation (LES) data obtained in EPSRC DTP Grant (ref: EP/R513349/136). The code allowed us to preform a Spectral Proper Orthogonal Decomposition (S-POD) for use in WP1 of the present project. Yuji Hattori -- funding for presenting paper at ICFD Conference. Sadaf Arabi/Zamir Koshuriyan -- work on analysis and calculations for the S-POD of the LES data for WP1 of the project.
Impact Refereed Conference paper in 2024 Investigation of Orthogonal Modes in High-speed Axisymmetric Jet Turbulence: Instantaneous vs. Statistical Data. Authors: Sadaf Arabi, Zamir Koshuriyan*, Adrian Sescu Presented at the Twenty-first International Conference on Flow Dynamics, November 18 - 20, 2024, Sendai International Center, Sendai, Miyagi, Japan. - Japan, Sendai, Japan. Duration: 18 Nov 2024 ? 20 Nov 2024. https://www.ifs.tohoku.ac.jp/icfd/2024/ This research visit was fully funded by IFS at Tohoku University with total funding (Flight travel, hotel stay & conference registration fee) of approximately £2000. Previous sucessful research projects: Lagrange Multiplier-Based Optimal Control Technique for Streak Attenuation in High-Speed Boundary Layers. Omar Es Sahli, Adrian Sescu*, Mohammad Koshuriyan, Yuji Hattori, Makoto HirotaLagrange Multiplier-Based Optimal Control Technique for Streak Attenuation in High-Speed Boundary Layers. AIAA Journal, pp. 1-29, vol. 61, published online 2022. Es-Sahli, O., Sescu, A., Afsar, M. Z. and Hattori, Y. 2021. Investigation of Gortler vortices in high-speed boundary layers via an efficient numerical solution to the non-linear boundary region equations. Theoretical and Computational Fluid Dynamics (TCFD), Vol.36, pp. 237-249. https://doi.org/10.1007/s00162-021-00576-w
Start Year 2017
 
Description Industrial Engagement with BAE Systems PLC 
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 We have obtained £2K from EPSRC Impact Acceleration Account (IAA) at York. We have used the funding to bring BAE Engineers over here to York and for me to visit BAE.

In the first instance, 2 members of BAE Systems PLC came to visit my group at the University of York. The BAE Engineers were from various branches of the company in different parts of England. The aim of the visit was to develop a deeper of understanding of BAE's Engineering challenges in Acoustics as well as to present the mathematical modelling we have on in the areas of jet turbulence and aeroacoustics. The meeting was held on 27/11/23.

We also discussed the potential of the current project on turbulence reconstruction via first principle asymptotic analysis. We have planned future meetings to take place this year and early 2025.

The IAA funding end date is December 2025. The aim will be for the PI (Dr. Z. Koshuriyan) to visit BAE facilities in Barrow in Furness and Filton. During these visits, discussion about the outcomes from my NIA award will be made. We hope that BAE will provide support for EPSRC CASE PhD studentships to complement Koshuriyan's current EPSRC NIA award.
Year(s) Of Engagement Activity 2023