BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Clarence W. Rowley

Clarence W. Rowley contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynmath.OCeess.SYphysics.comp-phSystems and ControlMachine Learningmath.DS

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
7topics
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

9 published item(s)

preprint2026arXiv

Resolvent4py: a parallel Python package for analysis, model reduction and control of large-scale linear systems

In this paper, we present resolvent4py, a parallel Python package for the analysis, model reduction and control of large-scale linear systems with millions or billions of degrees of freedom. This package provides the user with a friendly Python-like experience (akin to that of well-established libraries such as numpy and scipy), while enabling MPI-based parallelism through mpi4py, petsc4py and slepc4py. In turn, this allows for the development of streamlined and efficient Python code that can be used to solve several problems in fluid mechanics, solid mechanics, graph theory, molecular dynamics and several other fields.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Continuous-time balanced truncation for time-periodic fluid flows using frequential Gramians

Reduced-order models for flows that exhibit time-periodic behavior are critical for several tasks, including active control and optimization. One well-known procedure to obtain the desired reduced-order model in the proximity of a periodic solution of the governing equations is continuous-time balanced truncation. Within this framework, the periodic reachability and observability Gramians are usually estimated numerically via quadrature using the forward and adjoint post-transient response to impulses. However, this procedure can be computationally expensive, especially in the presence of slowly-decaying transients. Moreover, it can only be performed if the periodic orbit is stable in the sense of Floquet. In order to address these issues, we use the frequency-domain representation of the Gramians, which we henceforth refer to as frequential Gramians. First, these frequential Gramians are well-defined for both stable and unstable dynamics. In particular, we show that when the underlying system is unstable, these Gramians satisfy a pair of allied differential Lyapunov equations. Second, they can be estimated numerically by solving algebraic systems of equations that lend themselves to heavy computational parallelism and that deliver the desired post-transient response without having to follow physical transients. We demonstrate the method on a periodically-forced axisymmetric jet at Reynolds numbers Re=1250 and Re=1500. At the lower Reynolds number, the flow strongly amplifies subharmonic perturbations and exhibits vortex pairing about a Floquet-stable T-periodic solution. At the higher Reynolds number, the underlying T-periodic orbit is unstable and the flow naturally settles onto a 2T-periodic limit cycle characterized by pairing vortices. At both Reynolds numbers, we use a balanced reduced-order model to design a feedback controller and a state estimator to suppress vortex pairing.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optimal Control with Learning on the Fly: System with Unknown Drift

This paper derives an optimal control strategy for a simple stochastic dynamical system with constant drift and an additive control input. Motivated by the example of a physical system with an unexpected change in its dynamics, we take the drift parameter to be unknown, so that it must be learned while controlling the system. The state of the system is observed through a linear observation model with Gaussian noise. In contrast to most previous work, which focuses on a controller's asymptotic performance over an infinite time horizon, we minimize a quadratic cost function over a finite time horizon. The performance of our control strategy is quantified by comparing its cost with the cost incurred by an optimal controller that has full knowledge of the parameters. This approach gives rise to several notions of "regret." We derive a set of control strategies that provably minimize the worst-case regret; these arise from Bayesian strategies that assume a specific fixed prior on the drift parameter. This work suggests that examining Bayesian strategies may lead to optimal or near-optimal control strategies for a much larger class of realistic dynamical models with unknown parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optimizing Oblique Projections for Nonlinear Systems using Trajectories

Reduced-order modeling techniques, including balanced truncation and $\mathcal{H}_2$-optimal model reduction, exploit the structure of linear dynamical systems to produce models that accurately capture the dynamics. For nonlinear systems operating far away from equilibria, on the other hand, current approaches seek low-dimensional representations of the state that often neglect low-energy features that have high dynamical significance. For instance, low-energy features are known to play an important role in fluid dynamics where they can be a driving mechanism for shear-layer instabilities. Neglecting these features leads to models with poor predictive accuracy despite being able to accurately encode and decode states. In order to improve predictive accuracy, we propose to optimize the reduced-order model to fit a collection of coarsely sampled trajectories from the original system. In particular, we optimize over the product of two Grassmann manifolds defining Petrov-Galerkin projections of the full-order governing equations. We compare our approach with existing methods including proper orthogonal decomposition, balanced truncation-based Petrov-Galerkin projection, quadratic-bilinear balanced truncation, and the quadratic-bilinear iterative rational Krylov algorithm. Our approach demonstrates significantly improved accuracy both on a nonlinear toy model and on an incompressible (nonlinear) axisymmetric jet flow with $10^5$ states.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A computationally efficient approach for the removal of the phase shift singularity in harmonic resolvent analysis

The recently introduced harmonic resolvent framework is concerned with the study of the input-output dynamics of nonlinear flows in the proximity of a known time-periodic orbit. These dynamics are governed by the harmonic resolvent operator, which is a linear operator in the frequency domain whose singular value decomposition sheds light on the dominant input-output structures of the flow. Although the harmonic resolvent is a mathematically well-defined operator, the numerical computation of its singular value decomposition requires inverting a matrix that becomes exactly singular as the periodic orbit approaches an exact solution of the nonlinear governing equations. The very poor condition properties of this matrix hinder the convergence of classical Krylov solvers, even in the presence of preconditioners, thereby increasing the computational cost required to perform the harmonic resolvent analysis. In this paper we show that a suitable augmentation of the (nearly) singular matrix removes the singularity, and we provide a lower bound for the smallest singular value of the augmented matrix. We also show that the desired decomposition of the harmonic resolvent can be computed using the augmented matrix, whose improved condition properties lead to a significant speedup in the convergence of classical iterative solvers. We demonstrate this simple, yet effective, computational procedure on the Kuramoto-Sivashinsky equation in the proximity of an unstable time-periodic orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Inadequacy of Linear Methods for Minimal Sensor Placement and Feature Selection in Nonlinear Systems; a New Approach Using Secants

Sensor placement and feature selection are critical steps in engineering, modeling, and data science that share a common mathematical theme: the selected measurements should enable solution of an inverse problem. Most real-world systems of interest are nonlinear, yet the majority of available techniques for feature selection and sensor placement rely on assumptions of linearity or simple statistical models. We show that when these assumptions are violated, standard techniques can lead to costly over-sensing without guaranteeing that the desired information can be recovered from the measurements. In order to remedy these problems, we introduce a novel data-driven approach for sensor placement and feature selection for a general type of nonlinear inverse problem based on the information contained in secant vectors between data points. Using the secant-based approach, we develop three efficient greedy algorithms that each provide different types of robust, near-minimal reconstruction guarantees. We demonstrate them on two problems where linear techniques consistently fail: sensor placement to reconstruct a fluid flow formed by a complicated shock-mixing layer interaction and selecting fundamental manifold learning coordinates on a torus.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Distributed flexibility in inertial swimmers

We study a linear inviscid model of a passively flexible swimmer with distributed flexibility, calculating its propulsive performance and optimal distributions of flexibility. The frequencies of actuation and mean stiffness ratios we consider span a large range, while the mass ratio is fixed to a low value representative of swimmers. We present results showing how the trailing edge deflection, thrust coefficient, power coefficient, and efficiency vary with frequency, mean stiffness, and stiffness distribution. Swimmers with distributed flexibility have the same qualitative features as those with uniform flexibility. Significant gains in thrust can be made, however, by tuning the stiffness such that a resonant response is triggered, or by concentrating stiffness towards the leading edge if resonance cannot be triggered. To minimize power, the opposite is true. Meaningful gains in efficiency can be made at low frequencies by concentrating stiffness away from the leading edge, since doing so induces efficient travelling wave kinematics. We also speculate on the effects of a finite Reynolds number in the form of streamwise drag. The drag adds an offset to the net thrust produced by the swimmer, causing efficiency-maximizing distributions of flexibility to tend towards thrust-maximizing ones, representative of what is found in nature.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimal control with learning on the fly: a toy problem

We exhibit optimal control strategies for a simple toy problem in which the underlying dynamics depend on a parameter that is initially unknown and must be learned. We consider a cost function posed over a finite time interval, in contrast to much previous work that considers asymptotics as the time horizon tends to infinity. We study several different versions of the problem, including Bayesian control, in which we assume a prior distribution on the unknown parameter; and "agnostic" control, in which we assume nothing about the unknown parameter. For the agnostic problems, we compare our performance with that of an opponent who knows the value of the parameter. This comparison gives rise to several notions of "regret," and we obtain strategies that minimize the "worst-case regret" arising from the most unfavorable choice of the unknown parameter. In every case, the optimal strategy turns out to be a Bayesian strategy or a limit of Bayesian strategies.

0 citations0 reviewsNo code linkedOpen access