Researcher profile

Beverley J. McKeon

Beverley J. McKeon contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 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.