Efficient computation of global resolvent modes

From MaRDI portal
Publication:4993979

DOI10.1017/JFM.2021.364zbMATH Open1492.76062arXiv2008.10904OpenAlexW3165606180MaRDI QIDQ4993979FDOQ4993979


Authors: Eduardo Martini, Daniel Rodríguez, Aaron Towne, André V. G. Cavalieri Edit this on Wikidata


Publication date: 11 June 2021

Published in: Journal of Fluid Mechanics (Search for Journal in Brave)

Abstract: Resolvent analysis of the linearized Navier-Stokes equations provides useful insight into the dynamics of transitional and turbulent flows and can provide a model for the dominant coherent structures within the flow, particularly for flows with large gain separation. Direct computation of force and response modes using a singular value decomposition of the full resolvent matrix is feasible only for simple problems; despite recent progress, the cost of resolvent analysis for complex flows remains considerable. In this paper, we propose a new matrix-free method for computing resolvent modes based on integration of the linearized equations and the corresponding adjoint system in the time domain. Our approach achieves an order of magnitude speedup when compared to previous matrix-free time stepping methods by enabling all frequencies of interest to be computed simultaneously. Two different methods are presented: one based on analysis of the transient response, providing leading modes with fine frequency discretization; and another based on the steady-state response to a periodic forcing, providing optimal and suboptimal modes for a discrete set of frequencies. The methods are validated using a linearized Ginzburg-Landau equation and applied to the three dimensional flow around a parabolic body.


Full work available at URL: https://arxiv.org/abs/2008.10904




Recommendations




Cites Work


Cited In (21)

Uses Software





This page was built for publication: Efficient computation of global resolvent modes

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4993979)