Source author record

Beverley J. McKeon

Beverley J. McKeon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2021arXiv

Kernel Learning for Robust Dynamic Mode Decomposition: Linear and Nonlinear Disambiguation Optimization (LANDO)

Research in modern data-driven dynamical systems is typically focused on the three key challenges of high dimensionality, unknown dynamics, and nonlinearity. The dynamic mode decomposition (DMD) has emerged as a cornerstone for modeling high-dimensional systems from data. However, the quality of the linear DMD model is known to be fragile with respect to strong nonlinearity, which contaminates the model estimate. In contrast, sparse identification of nonlinear dynamics (SINDy) learns fully nonlinear models, disambiguating the linear and nonlinear effects, but is restricted to low-dimensional systems. In this work, we present a kernel method that learns interpretable data-driven models for high-dimensional, nonlinear systems. Our method performs kernel regression on a sparse dictionary of samples that appreciably contribute to the underlying dynamics. We show that this kernel method efficiently handles high-dimensional data and is flexible enough to incorporate partial knowledge of system physics. It is possible to accurately recover the linear model contribution with this approach, disambiguating the effects of the implicitly defined nonlinear terms, resulting in a DMD-like model that is robust to strongly nonlinear dynamics. We demonstrate our approach on data from a wide range of nonlinear ordinary and partial differential equations that arise in the physical sciences. This framework can be used for many practical engineering tasks such as model order reduction, diagnostics, prediction, control, and discovery of governing laws.

preprint2020arXiv

From unsteady to quasi-steady dynamics in the streamwise-oscillating cylinder wake

The flow around a cylinder oscillating in the streamwise direction with a frequency, $f_f$, much lower than the shedding frequency, $f_s$, has been relatively less studied than the case when these frequencies have the same order of magnitude, or the transverse oscillation configuration. In this study, Particle Image Velocimetry and Koopman Mode Decomposition are used to investigate the streamwise-oscillating cylinder wake for forcing frequencies $f_f/f_s \sim 0.04-0.2$ and mean Reynolds number, $Re_0 = 900$. The amplitude of oscillation is such that the instantaneous Reynolds number remains above the critical value for vortex shedding at all times. Characterization of the wake reveals a range of phenomena associated with the interaction of the two frequencies, including modulation of both the amplitude and frequency of the wake structure by the forcing. Koopman analysis reveals a frequency spreading of Koopman modes. A scaling parameter and associated transformation are developed to relate the unsteady, or forced, dynamics of a system to that of a quasi-steady, or unforced, system. For the streamwise-oscillating cylinder, it is shown that this transformation leads to a Koopman Mode Decomposition similar to that of the unforced system.

preprint2020arXiv

Interaction of forced Orr-Sommerfeld and Squire modes in a low-order representation of turbulent channel flow

A resolvent-based reduced-order representation is used to capture time-averaged second-order statistics in turbulent channel flow. The recently-proposed decomposition of the resolvent operator into two distinct families related to the Orr-Sommerfeld and Squire operators [K. Rosenberg and B. J. McKeon, Efficient representation of exact coherent states of the Navier-Stokes equations using resolvent analysis, Fluid Dynamics Research 51, 011401 (2019)] results in dramatic improvement in the ability to match all components of the energy spectra and the $uv$ cospectrum. The success of the new representation relies on the ability of the Squire modes to compete with the vorticity generated by Orr-Sommerfeld modes, which is demonstrated by decomposing the statistics into contributions from each family. It is then shown that this competition can be used to infer a phase relationship between the two sets of modes. Additionally, the relative Reynolds number scalings for the two families of resolvent weights are derived for the universal classes of resolvent modes presented by Moarref et al. [R. Moarref, A. S. Sharma, J. A. Tropp, and B. J. McKeon, Model-based scaling of the streamwise energy density in high-Reynolds-number turbulent channels, Journal of Fluid Mechanics 734, 275 (2013)]. These developments can be viewed as a starting point for further modeling efforts to quantify nonlinear interactions in wall-bounded turbulence.

preprint2020arXiv

Self-sustained elastoinertial Tollmien-Schlichting waves

Direct simulations of two-dimensional plane channel flow of a viscoelastic fluid at Reynolds number Re = 3000 reveal the existence of a family of attractors whose structure closely resembles the linear Tollmien-Schlichting (TS) mode, and in particular exhibits strongly localized stress fluctuations at the critical layer position of the TS mode. At the parameter values chosen, this solution branch is not connected to the nonlinear TS solution branch found for Newtonian flow, and thus represents a solution family that is nonlinearly self-sustained by viscoelasticity. The ratio between stress and velocity fluctuations is in quantitative agreement for the attractor and the linear TS mode, and increases strongly with Weissenberg number, Wi. For the latter, there is a transition in the scaling of this ratio as Wi increases, and the Wi at which the nonlinear solution family comes into existence is just above this transition. Finally, evidence indicates that this branch is connected through an unstable solution branch to two-dimensional elastoinertial turbulence (EIT). These results suggest that, in the parameter range considered here, the bypass transition leading to EIT is mediated by nonlinear amplification and self-sustenance of perturbations that excite the Tollmien-Schlichting mode.

preprint2019arXiv

A tale of two airfoils: resolvent-based modelling of an oscillator vs. an amplifier from an experimental mean

The flows around a NACA 0018 airfoil at a Reynolds number of 10250 and angles of attack of alpha = 0 (A0) and alpha = 10 (A10) are modelled using resolvent analysis and limited experimental measurements obtained from particle image velocimetry. The experimental mean velocity profiles are data-assimilated so that they are solutions of the incompressible Reynolds-averaged Navier-Stokes equations forced by Reynolds stress terms which are derived from experimental data. Spectral proper orthogonal decompositions (SPOD) of the velocity fluctuations and nonlinear forcing find low-rank behaviour at the shedding frequency and its higher harmonics for the A0 case. In the A10 case, low-rank behaviour is observed for the velocity fluctuations in two bands of frequencies. Resolvent analysis of the data-assimilated means identifies low-rank behaviour only in the vicinity of the shedding frequency for A0 and none of its harmonics. The resolvent operator for the A10 case, on the other hand, identifies two linear mechanisms whose frequencies are a close match with those identified by SPOD. It is also shown that the second linear mechanism, corresponding to the Kelvin-Helmholtz instability in the shear layer, cannot be identified just by considering the time-averaged experimental measurements as a mean flow due to the fact that experimental data are missing near the leading edge. The A0 case is classified as an oscillator where the flow is organised around an intrinsic instability while the A10 case behaves like an amplifier whose forcing is unstructured. For both cases, resolvent modes resemble those from SPOD when the operator is low-rank. To model the higher harmonics where this is not the case, we add parasitic resolvent modes, as opposed to classical resolvent modes which are the most amplified, by approximating the nonlinear forcing from limited triadic interactions of known resolvent modes.

preprint2016arXiv

Low-dimensional representations of exact coherent states of the Navier-Stokes equations from the resolvent model of wall turbulence

We report that many exact invariant solutions of the Navier-Stokes equations for both pipe and channel flows are well represented by just few modes of the model of McKeon & Sharma J. Fl. Mech. 658, 356 (2010). This model provides modes that act as a basis to decompose the velocity field, ordered by their amplitude of response to forcing arising from the interaction between scales. The model was originally derived from the Navier-Stokes equations to represent turbulent flows and has been used to explain coherent structure and to predict turbulent statistics. This establishes a surprising new link between the two distinct approaches to understanding turbulence.

preprint2016arXiv

On the correspondence between Koopman mode decomposition, resolvent mode decomposition, and invariant solutions of the Navier-Stokes equations

The relationship between Koopman mode decomposition, resolvent mode decomposition and exact invariant solutions of the Navier-Stokes equations is clarified. The correspondence rests upon the invariance of the system operators under symmetry operations such as spatial translation. The usual interpretation of the Koopman operator is generalised to permit combinations of such operations, in addition to translation in time. This invariance is related to the spectrum of a spatio-temporal Koopman operator, which has a travelling wave interpretation. The relationship leads to a generalisation of dynamic mode decomposition, in which symmetry operations are applied to restrict the dynamic modes to span a subspace subject to those symmetries. The resolvent is interpreted as the mapping between the Koopman modes of the Reynolds stress divergence and the velocity field. It is shown that the singular vectors of the resolvent (the resolvent modes) are the optimal basis in which to express the velocity field Koopman modes where the latter are not a priori known.

preprint2014arXiv

A foundation for analytical developments in the logarithmic region of turbulent channels

An analytical framework for studying the logarithmic region of turbulent channels is formulated. We build on recent findings (Moarref et al., J. Fluid Mech., 734, 2013) that the velocity fluctuations in the logarithmic region can be decomposed into a weighted sum of geometrically self-similar resolvent modes. The resolvent modes and the weights represent the linear amplification mechanisms and the scaling influence of the nonlinear interactions in the Navier-Stokes equations (NSE), respectively (McKeon & Sharma, J. Fluid Mech., 658, 2010). Originating from the NSE, this framework provides an analytical support for Townsend's attached-eddy model. Our main result is that self-similarity enables order reduction in modeling the logarithmic region by establishing a quantitative link between the self-similar structures and the velocity spectra. Specifically, the energy intensities, the Reynolds stresses, and the energy budget are expressed in terms of the resolvent modes with speeds corresponding to the top of the logarithmic region. The weights of the triad modes -the modes that directly interact via the quadratic nonlinearity in the NSE- are coupled via the interaction coefficients that depend solely on the resolvent modes (McKeon et al., Phys. Fluids, 25, 2013). We use the hierarchies of self-similar modes in the logarithmic region to extend the notion of triad modes to triad hierarchies. It is shown that the interaction coefficients for the triad modes that belong to a triad hierarchy follow an exponential function. The combination of these findings can be used to better understand the dynamics and interaction of flow structures in the logarithmic region. The compatibility of the proposed model with theoretical and experimental results is further discussed.

preprint2013arXiv

Influence of a local change of depth on the behavior of bouncing oil drops

The work of Couder \textit{et al} (see also Bush \textit{et al}) inspired consideration of the impact of a submerged obstacle, providing a local change of depth, on the behavior of oil drops in the bouncing regime. In the linked videos, we recreate some of their results for a drop bouncing on a uniform depth bath of the same liquid undergoing vertical oscillations just below the conditions for Faraday instability, and show a range of new behaviors associated with change of depth. This article accompanies a fluid dynamics video entered into the Gallery of Fluid Motion of the 66th Annual Meeting of the APS Division of Fluid Dynamics.

preprint2013arXiv

Model-based scaling of the streamwise energy density in high-Reynolds number turbulent channels

We study the Reynolds number scaling and the geometric self-similarity of a gain-based, low-rank approximation to turbulent channel flows, determined by the resolvent formulation of McKeon & Sharma (2010), in order to obtain a description of the streamwise turbulence intensity from direct consideration of the Navier-Stokes equations. Under this formulation, the velocity field is decomposed into propagating waves (with single streamwise and spanwise wavelengths and wave speed) whose wall-normal shapes are determined from the principal singular function of the corresponding resolvent operator. Using the accepted scalings of the mean velocity in wall-bounded turbulent flows, we establish that the resolvent operator admits three classes of wave parameters that induce universal behavior with Reynolds number on the low-rank model, and which are consistent with scalings proposed throughout the wall turbulence literature. In addition, it was shown that a necessary condition for geometrically self-similar resolvent modes is the presence of a logarithmic turbulent mean velocity. We identify the scalings that constitute hierarchies of self-similar modes that are parameterized by the critical wall-normal location where the speed of the mode equals the local turbulent mean velocity. For the rank-1 model subject to broadband forcing, the integrated streamwise energy density takes a universal form which is consistent with the dominant near-wall turbulent motions. When the shape of the forcing is optimized to enforce matching with results from direct numerical simulations at low turbulent Reynolds numbers, further similarity appears. Representation of these weight functions using similarity laws enables prediction of the Reynolds number and wall-normal variations of the streamwise energy intensity at high Reynolds numbers (${Re}_τ\approx 10^3 - 10^{10}$).

preprint2011arXiv

A streamwise-constant model of turbulent pipe flow

A streamwise-constant model is presented to investigate the basic mechanisms responsible for the change in mean flow occuring during pipe flow transition. Using a single forced momentum balance equation, we show that the shape of the velocity profile is robust to changes in the forcing profile and that both linear non-normal and nonlinear effects are required to capture the change in mean flow associated with transition to turbulence. The particularly simple form of the model allows for the study of the momentum transfer directly by inspection of the equations. The distribution of the high- and low-speed streaks over the cross-section of the pipe produced by our model is remarkably similar to one observed in the velocity field near the trailing edge of the puff structures present in pipe flow transition. Under stochastic forcing, the model exhibits a quasi-periodic self-sustaining cycle characterized by the creation and subsequent decay of "streamwise-constant puffs", so-called due to the good agreement between the temporal evolution of their velocity field and the projection of the velocity field associated with three-dimensional puffs in a frame of reference moving at the bulk velocity. We establish that the flow dynamics are relatively insensitive to the regeneration mechanisms invoked to produce near-wall streamwise vortices and that using small, unstructured background disturbances to regenerate the streamwise vortices is sufficient to capture the formation of the high- and low-speed streaks and their segregation leading to the blunting of the velocity profile characteristic of turbulent pipe flow.

preprint2010arXiv

Amplification and Nonlinear Mechanisms in Plane Couette Flow

We study the input-output response of a streamwise constant projection of the Navier-Stokes equations for plane Couette flow, the so-called 2D/3C model. Study of a streamwise constant model is motivated by numerical and experimental observations that suggest the prevalence and importance of streamwise and quasi-streamwise elongated structures. Periodic spanwise/wall-normal (z-y) plane stream functions are used as input to develop a forced 2D/3C streamwise velocity field that is qualitatively similar to a fully turbulent spatial field of DNS data. The input-output response associated with the 2D/3C nonlinear coupling is used to estimate the energy optimal spanwise wavelength over a range of Reynolds numbers. The results of the input-output analysis agree with previous studies of the linearized Navier-Stokes equations. The optimal energy corresponds to minimal nonlinear coupling. On the other hand, the nature of the forced 2D/3C streamwise velocity field provides evidence that the nonlinear coupling in the 2D/3C model is responsible for creating the well known characteristic "S" shaped turbulent velocity profile. This indicates that there is an important tradeoff between energy amplification, which is primarily linear and the seemingly nonlinear momentum transfer mechanism that produces a turbulent-like mean profile.