Resolvent Analysis Research Articles

Triglobal resolvent analysis of swept-wing wakes

Through triglobal resolvent analysis, we reveal the effects of wing tip and sweep angle on laminar separated wakes over swept wings. For the present study, we consider wings with semi-aspect ratios from $1$ to $4$ , sweep angles from $0^\circ$ to $45^\circ$ and angles of attack of $20^\circ$ and $30^\circ$ at a chord-based Reynolds number of $400$ and a Mach number of $0.1$ . Using direct numerical simulations, we observe that unswept wings develop vortex shedding near the wing root with a quasi-steady tip vortex. For swept wings, vortex shedding is seen near the wing tip for low sweep angles, while the wakes are steady for wings with high sweep angles. To gain further insights into the mechanisms of flow unsteadiness, triglobal resolvent analysis is used to identify the optimal spatial input–output mode pairs and the associated gains over a range of frequencies. The three-dimensional forcing and response modes reveal that harmonic fluctuations are directed towards the root for unswept wings and towards the wing tip for swept wings. The overlapping region of the forcing–response mode pairs uncovers triglobal resolvent wavemakers associated with self-sustained unsteady wakes of swept wings. Furthermore, we show that for low-aspect-ratio wings optimal perturbations develop globally over the entire wingspan. The present study uncovers physical insights on the effects of tip and sweep on the growth of optimal harmonic perturbations and the wake dynamics of separated flows over swept wings.

Reactive control of second Mack mode in a supersonic boundary layer with free-stream velocity/density variations

We consider closed-loop control of a two-dimensional supersonic boundary layer at $M=4.5$ that aims at reducing the linear growth of second Mack mode instabilities. These instabilities are first characterized with local spatial and global resolvent analyses, which allow us to refine the control strategy and to select appropriate actuators and sensors. After linear input–output reduced-order models have been identified, multi-criteria structured mixed $H_{2}$ / $H_{\infty }$ synthesis allows us to fix beforehand the controller structure and to minimize appropriate norms of various transfer functions: the $H_{2}$ norm to guarantee performance (reduction of perturbation amplification in nominal condition), and the $H_{\infty }$ norm to maintain performance robustness (with respect to sensor noise) and stability robustness (with respect to uncertain free-stream velocity/density variations). Both feedforward and feedback set-ups, i.e. with estimation sensor placed respectively upstream/downstream of the actuator, allow us to maintain the local perturbation energy below a given threshold over a significant distance downstream of the actuator, even in the case of noisy estimation sensors or free-stream density variations. However, the feedforward set-up becomes completely ineffective when convective time delays are altered by free-stream velocity variations of $\pm$ 5 %, which highlights the strong relevance of the feedback set-up for performance robustness in convectively unstable flows.

Resolvent analysis of a finite wing in transonic flow

Shock waves interacting with turbulent boundary layers on wings can result first in self-sustained flow unsteadiness and eventually in structural vibration. Due to its importance to modern wing design and aircraft certification, the transonic flow physics continue to be investigated intensively. Herein we focus the discussion on three main aspects. First, we assess a practical implementation of an iterative resolvent algorithm in the linear harmonic incarnation of an industrial computational fluid dynamics code for computing optimal forcing and response modes. This heavily relies on the efficient solution of large sparse linear systems of equations. Second, we showcase its application as a predictive tool to detect transonic buffet flow unsteadiness, well before a global stability analysis can first identify its dynamics through weakly damped eigenmodes, using the NASA common research model at wind-tunnel conditions. Third, we discuss its ability to uncover modal physics, not identifiable through global stability analysis, revealing higher-frequency wake and wingtip vortex modes while shedding some light on the elusive finite wing equivalent of the aerofoil buffet mode. We demonstrate that earlier computational limitations of resolvent analysis, when solving the truncated singular value decomposition using matrix-forming methods with direct matrix factorisation, have been overcome ready for industrial use.

Lee model and its resolvent analysis

We revisit the relativistic (2+1)-dimensional Lee model on flat space in light-front coordinates and on a space-time with a spatial section given by a compact manifold, in the usual canonical formalism. The simpler 2+1 dimension is chosen because renormalization is needed only for the mass difference but not required for the coupling constant and the wave function. The model is constructed non-perturbatively based on the resolvent formulation [B. T. Kaynak and O. T. Turgut, The relativistic Lee model on Riemannian manifolds, J. Phys. A: Math. Theor. 42(22) (2009) 225402]. The bound state spectrum is studied through its “principal operator” and bounds for the ground state energy are obtained. We show that the formal expression found indeed defines the resolvent of a self-adjoint operator–the Hamiltonian of the interacting system. Moreover, we prove an essential result that the principal operator corresponds to a self-adjoint holomorphic family of type-A, in the sense of Kato.

Transonic Buffet Characteristics Under Conditions of Free and Forced Transition

Transonic buffet is commonly associated with self-sustained flow unsteadiness involving shock-wave/boundary-layer interaction over airfoils and wings. The phenomenon has been classified as either laminar or turbulent based on the state of the boundary layer immediately upstream of the shock foot, and distinct mechanisms for the two types have been suggested. The turbulent case is known to be associated with a global linear instability. Herein, large-eddy simulations are used for the first time to make direct comparisons of the two types by examining free- and forced-transition conditions. Corresponding simulations based on the Reynolds-averaged Navier–Stokes equations for the forced-transition case are also performed for comparison with the scale-resolving approach and for linking the findings with existing literature. Coherent flow features are scrutinized using both data-based spectral proper orthogonal decomposition of the time-marched results and operator-based global linear stability and resolvent analyses within the Reynolds-averaged Navier–Stokes framework. It is demonstrated that the essential dynamic features remain the same for the two buffet types (and for the two levels of the aerodynamic modeling hierarchy), suggesting that both types arise due to the same fundamental mechanism.

Parabolic resolvent modes for streaky structures in transitional and turbulent boundary layers

Resolvent analysis is a powerful tool in fluid mechanics, both for laminar and turbulent flows. Depending on application, computational cost is still a limiting factor. We develop a parabolic method for the calculation of low-frequency resolvent modes, associated with streaky structures, which can enable more than an order of magnitude speed up compared to usual global calculations. The method is applied to both laminar and turbulent boundary layers, unveiling physical mechanisms such as those associated with streak growth due to free-stream turbulence, in the laminar case, and the double peak in the frequency spectrum of velocity fluctuations, in a turbulent boundary layer.

Wing sweep effects on laminar separated flows

We reveal the effects of sweep on the wake dynamics around NACA 0015 wings at high angles of attack using direct numerical simulations and resolvent analysis. The influence of sweep on the wake dynamics is considered for sweep angles from $0^\circ$ to $45^\circ$ and angles of attack from $16^\circ$ to $30^\circ$ for a spanwise periodic wing at a chord-based Reynolds number of $400$ and a Mach number of $0.1$ . Wing sweep affects the wake dynamics, especially in terms of stability and spanwise fluctuations with implications on the development of three-dimensional (3-D) wakes. We observe that wing sweep attenuates spanwise fluctuations. Even as the sweep angle influences the wake, force and pressure coefficients can be collapsed for low angles of attack when examined in wall-normal and wingspan-normal independent flow components. Some small deviations at high sweep and incidence angles are attributed to vortical wake structures that impose secondary aerodynamic loads, revealed through the force element analysis. Furthermore, we conduct global resolvent analysis to uncover oblique modes with high disturbance amplification. The resolvent analysis also reveals the presence of wavemakers in the shear-dominated region associated with the emergence of 3-D wakes at high angles of attack. For flows at high sweep angles, the optimal convection speed of the response modes is shown to be faster than the optimal wavemakers speed suggesting a mechanism for the attenuation of perturbations. The present findings serve as a fundamental stepping stone to understanding separated flows at higher Reynolds numbers.

Improved convergence of the spectral proper orthogonal decomposition through time shifting

Spectral proper orthogonal decomposition (SPOD) is an increasingly popular modal analysis method in the field of fluid dynamics due to its specific properties: a linear system forced with white noise should have SPOD modes identical to response modes from resolvent analysis. The SPOD, coupled with the Welch method for spectral estimation, may require long time-resolved datasets. In this work, a linearised Ginzburg–Landau model is considered in order to study the method's convergence. Spectral proper orthogonal decomposition modes of the white-noise forced equation are computed and compared with corresponding response resolvent modes. The quantified error is shown to be related to the time length of Welch blocks (spectral window size) normalised by a convective time. Subsequently, an algorithm based on a temporal data shift is devised to further improve SPOD convergence and is applied to the Ginzburg–Landau system. Next, its efficacy is demonstrated in a numerical database of a boundary layer subject to bypass transition. The proposed approach achieves substantial improvement in mode convergence with smaller spectral window sizes with respect to the standard method. Furthermore, SPOD modes display growing wall-normal and spanwise velocity components along the streamwise direction, a feature which had not yet been observed and is also predicted by a global resolvent calculation. The shifting algorithm for the SPOD opens the possibility for using the method on datasets with time series of moderate duration, often produced by large simulations.

Oblique transition in hypersonic double-wedge flow

We utilize resolvent and weakly nonlinear analyses in combination with direct numerical simulations (DNS) to identify mechanisms for oblique transition in a Mach $5$ hypersonic flow over an adiabatic slender double wedge. Even though the laminar separated flow is globally stable, resolvent analysis demonstrates significant amplification of unsteady external disturbances to the linearized flow equations. These disturbances are introduced upstream of the separation zone and they lead to the appearance of oblique waves further downstream. We demonstrate that the large amplification of oblique waves arises from the growth of fluctuation shear stress due to streamline curvature of the laminar base flow in the separated shear layer. This is in contrast to the attached boundary layers, where no such mechanism exists. We also use a weakly nonlinear analysis to show that the resolvent operator associated with linearization around the laminar base flow governs the evolution of steady reattachment streaks that arise from quadratic interactions of unsteady oblique waves. These quadratic interactions generate vortical excitations in the reattaching shear layer which lead to the formation of streaks in the recirculation zone and their subsequent amplification, breakdown and transition to turbulence downstream. Our analysis of the energy budget shows that deceleration of the base flow near reattachment is primarily responsible for amplification of steady streaks. Finally, we employ DNS to examine latter stages of transition to turbulence and demonstrate the predictive power of a weakly nonlinear input–output framework in uncovering triggering mechanisms for oblique transition in separated high-speed boundary layer flows.

Efficient global resolvent analysis via the one-way Navier–Stokes equations

Resolvent analysis is a powerful tool for modelling and analysing transitional and turbulent flows and, in particular, for approximating coherent flow structures. Despite recent algorithmic advances, computing resolvent modes for flows with more than one inhomogeneous spatial coordinate remains computationally expensive. In this paper we show how efficient and accurate approximations of resolvent modes can be obtained using a well-posed spatial marching method for flows that contain a slowly varying direction, i.e. one in which the mean flow changes gradually. First, we derive a well-posed and convergent one-way equation describing the downstream-travelling waves supported by the linearized Navier–Stokes equations. The method is based on a projection operator that isolates downstream-travelling waves. Integrating these one-way Navier–Stokes (OWNS) equations in the slowly varying direction, which requires significantly less CPU and memory resources than a direct solution of the linearized Navier–Stokes equations, approximates the action of the resolvent operator on a forcing vector. Second, this capability is leveraged to compute approximate resolvent modes using an adjoint-based optimization framework in which the forward and adjoint OWNS equations are marched in the downstream and upstream directions, respectively. This avoids the need to solve direct and adjoint globally discretized equations, therefore bypassing the main computational bottleneck of a typical global resolvent calculation. The method is demonstrated using the examples of a simple acoustics problem, a Mach 1.5 turbulent jet and a Mach 4.5 transitional zero-pressure-gradient flat-plate boundary layer. The optimal OWNS results are validated against corresponding global calculations, and the close agreement demonstrates the near-parabolic nature of these flows.

Sparsifying the resolvent forcing mode via gradient-based optimisation

We consider the use of sparsity-promoting norms in obtaining localised forcing structures from resolvent analysis. By formulating the optimal forcing problem as a Riemannian optimisation, we are able to maximise cost functionals whilst maintaining a unit-energy forcing. Taking the cost functional to be the energy norm of the driven response results in a traditional resolvent analysis and is solvable by a singular value decomposition (SVD). By modifying this cost functional with the $L_1$ -norm, we target spatially localised structures that provide an efficient amplification in the energy of the response. We showcase this optimisation procedure on two flows: plane Poiseuille flow at Reynolds number $Re=4000$ , and turbulent flow past a NACA 0012 aerofoil at $Re=23\,000$ . In both cases, the optimisation yields sparse forcing modes that maintain important features of the structures arising from an SVD in order to provide a gain in energy. These results showcase the benefits of utilising a sparsity-promoting resolvent formulation to uncover sparse forcings, specifically with a view to using them as actuation locations for flow control.

Low-frequency resolvent analysis of the laminar oblique shock wave/boundary layer interaction

Resolvent analysis is used to study the low-frequency behaviour of the laminar oblique shock wave/boundary layer interaction (SWBLI). It is shown that the computed optimal gain, which can be seen as a transfer function of the system, follows a first-order low-pass filter equation, recovering the results of Touber & Sandham (J. Fluid Mech., vol. 671, 2011, pp. 417–465). This behaviour is understood as proceeding from the excitation of a single stable, steady global mode whose damping rate sets the time scale of the filter. Different Mach and Reynolds numbers are studied, covering different recirculation lengths $L$ . This damping rate is found to scale as $1/L$ , leading to a constant Strouhal number $St_{L}$ as observed in the literature. It is associated with a breathing motion of the recirculation bubble. This analysis furthermore supports the idea that the low-frequency dynamics of the SWBLI is a forced dynamics, in which background perturbations continuously excite the flow. The investigation is then carried out for three-dimensional perturbations for which two regimes are identified. At low wavenumbers of the order of $L$ , a modal mechanism similar to that of two-dimensional perturbations is found and exhibits larger values of the optimal gain. At larger wavenumbers, of the order of the boundary layer thickness, the growth of streaks, which results from a non-modal mechanism, is detected. No interaction with the recirculation region is observed. Based on these results, the potential prevalence of three-dimensional effects in the low-frequency dynamics of the SWBLI is discussed.

Identification and reconstruction of high-frequency fluctuations evolving on a low-frequency periodic limit cycle: application to turbulent cylinder flow

Oftentimes, turbulent flows exhibit a high-frequency turbulent component developing on a strong low-frequency periodic motion. In such cases, the low-frequency motion may strongly influence the spatio-temporal features of the high-frequency component. A typical example of such behaviour is the flow around bluff bodies, for which the high-frequency turbulent component, characterized by Kelvin–Helmholtz structures associated with thin shear layers, depends on the phase of the low-frequency vortex-shedding motion. In this paper, we propose extended versions of spectral proper orthogonal decomposition (SPOD) and of resolvent analysis that respectively extract and reconstruct the high-frequency turbulent fluctuation field as a function of the phase of the low-frequency periodic motion. These approaches are based on a quasi-steady (QS) assumption, which may be justified by the supposedly large separation between the frequencies of the periodic and turbulent components. After discussing their relationship to more classical Floquet-like analyses, the new tools are illustrated on a simple periodically varying linear Ginzburg–Landau model, mimicking the overall characteristics of a turbulent bluff-body flow. In this simple model, we in particular assess the validity of the QS approximation. Then, we consider the case of turbulent flow around a squared-section cylinder at a Reynolds number of $Re=22\,000$ , for which we show reasonable agreement between the extracted spatio-temporal fluctuation field and the prediction of QS resolvent analysis at the various phases of the periodic vortex-shedding motion.

Resolvent-based approach for pmb {H_2}-optimal estimation and control: an application to the cylinder flow

We consider estimation and control of the cylinder wake at low Reynolds numbers. A particular focus is on the development of efficient numerical algorithms to design optimal linear feedback controllers when there are many inputs (disturbances applied everywhere) and many outputs (perturbations measured everywhere). We propose a resolvent-based iterative algorithm to perform (i) optimal estimation of the flow using a limited number of sensors, and (ii) optimal control of the flow when the entire flow is known but only a limited number of actuators are available for control. The method takes advantage of the low-rank characteristics of the cylinder wake and provides full-dimensional solutions by implementing a terminal reduction technique based on resolvent analysis. Optimal feedback controllers are also obtained by combining the solutions of the estimation and control problems. We show that the performance of the estimators and controllers converges to the true global optima, indicating that the important physical mechanisms for estimation and control are of low rank.Graphic abstract

Self-similar mechanisms in wall turbulence studied using resolvent analysis

Self-similarity of wall-attached coherent structures in a turbulent channel at $Re_\tau =543$ is explored by means of resolvent analysis. In this modelling framework, coherent structures are understood to arise as a response of the linearised mean-flow operator to generalised frequency-dependent Reynolds stresses, considered to act as an endogenous forcing. We assess the self-similarity of both the wall-attached flow structures and the associated forcing. The former are educed from direct numerical simulation data by finding the flow field correlated with the wall shear, whereas the latter is identified using a frequency space version of extended proper orthogonal decomposition (Borée, Exp. Fluids, vol. 35, issue 2, 2003, pp. 188–192). The forcing structures identified are compared to those obtained using the resolvent-based estimation introduced by Towne et al. (J. Fluid Mech., vol. 883, 2020, A17). The analysis reveals self-similarity of both wall-attached structures – in quantitative agreement with Townsend's hypothesis of self-similar attached eddies – and the underlying forcing, at least in certain components.

Linear instability of a premixed slot flame: Flame transfer function and resolvent analysis

The response to forcing of a 2D laminar premixed slot flame is investigated by means of linear analysis, based on the compressible flow equations with a two-step reaction scheme for methane combustion. The flame transfer function (FTF) is computed from this linear model, in excellent agreement with reference nonlinear calculations. The input-output gain between externally applied forcing and the global heat release rate response is computed, and peaks in the gain are related to intrinsic thermoacoustic (ITA) modes. The receptivity of the flame to arbitrary flow forcing is characterised by the resulting amplitude of global heat release rate fluctuations. Linear resolvent analysis is used to identify optimal forcing structures and their associated flame response, leading to a discussion of the dominant mechanisms for the amplification of flow perturbations, which trigger flame oscillations. These seem to involve a resonance with ITA instability modes.

Transition, intermittency and phase interference effects in airfoil secondary tones and acoustic feedback loop

A large eddy simulation is performed to study secondary tones generated by a NACA0012 airfoil at angle of attack of $\alpha = 3^{\circ}$ with freestream Mach number of $M_{\infty} = 0.3$ and Reynolds number of $Re = 5 \times 10^4$ . Laminar separation bubbles are observed over the suction side and near the trailing edge, on the pressure side. Vortex shedding occurs aft of the suction side separation bubble, and vortex interaction results in merging or bursting such that coherent structures or turbulent packets are advected towards the trailing edge. This mechanism modulates the amplitude of the incident pressure signal, leading to different levels of noise emission. Despite the intermittent occurrence of laminar–turbulent transition, the noise spectrum depicts a main tone with multiple equidistant secondary tones. To understand the role of flow instabilities on the tones, the linearised Navier–Stokes equations are examined in their operator form through biglobal stability and resolvent analyses, and by time evolution of disturbances using a matrix-free method. These linear global analyses reveal amplification of disturbances over the suction side separation bubble. Spanwise-averaged pressure fluctuations elucidate aspects of the acoustic feedback loop mechanism in the nonlinear solutions. This feedback process is self-sustained by acoustic waves radiated from the trailing edge, which reach the most sensitive flow location near the leading edge, as identified by the resolvent analysis. Flow disturbances arising from secondary diffraction and phase interference among the most unstable frequencies computed in the eigenspectrum are shown to have an important role in the feedback loop.

Global stability and resolvent analyses of laminar boundary-layer flow interacting with viscoelastic patches

The attenuation of two-dimensional boundary-layer instabilities by a finite-length, viscoelastic patch is investigated by means of global linear stability theory. First, the modal stability properties of the coupled problem are assessed, revealing unstable fluid-elastic travelling-wave flutter modes. Second, the Tollmien–Schlichting instabilities over a rigid wall are characterised via the analysis of the fluid resolvent operator in order to determine a baseline for the fluid-structural analysis. To investigate the effect of the elastic patch on the growth of these flow instabilities, we first consider the linear frequency response of the coupled fluid-elastic system to the dominant rigid-wall forcing modes. In the frequency range of Tollmien–Schlichting waves, the energetic flow amplification is clearly reduced. However, an amplification is observed for higher frequencies, associated with travelling-wave flutter. This increased complexity requires the analysis of the coupled fluid-structural resolvent operator; the optimal, coupled, resolvent modes confirm the attenuation of the Tollmien–Schlichting instabilities, while also being able to capture the amplification at the higher frequencies. Finally, a decomposition of the fluid-structural response is proposed to reveal the wave cancellation mechanism responsible for the attenuation of the Tollmien–Schlichting waves. The viscoelastic patch, excited by the incoming rigid-wall wave, provokes a fluid-elastic wave that is out-of-phase with the former, thus reducing its amplitude.

Variational formulation of resolvent analysis

This work presents an alternative formulation of resolvent analysis based on the calculus of variations. The proposed formulation circumvents reliance on inversion of the linear operator and is inherently compatible with any arbitrary choice of norm. Several examples are presented which highlight the analytical and numerical advantages and limitations of the proposed formulation.

Hysteresis of two-dimensional flows around a NACA0012 airfoil at [formula omitted] and linear analyses of their mean flow.

Two-dimensional numerical simulations of the flow around a NACA0012 profile at Reynolds number Re=5000 show that unsteady periodic flows reach different saturated states when increasing or decreasing the angle of attack between 7∘ and 8∘. Within this range, the lift signal shows co-existing periodic states and period-doubling, as the wake undergoes a substantial change in character from the standard von-Kármán vortex street. Results of experiments in a water channel also indicate a change of the flow topology but at slightly lower angles of attack α=6°. A discussion of the discrepancy between numerical and experimental results is proposed in light of results about the three-dimensional transition of wake flows behind bluff bodies and airfoils. Finally, eigenvalue and resolvent analyses of time-averaged flows are used to investigate the two-dimensional transitions further. While a peak of energetic amplification is obtained at the frequency of a single periodic state, a double peak is observed for co-existing periodic states, the second one being at the frequency of the periodic state not used to compute the time-averaged flow. This behaviour also characterizes the resolvent analysis of the period-doubled states, although less pronounced.