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

Alberto Padovan

Alberto Padovan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynphysics.comp-phmath-phmath.DSmath.MPmath.OC

Trust snapshot

Quick read

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

6 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
preprint2025arXiv

Surface instabilities in laminar compressible boundary layers with sublimation

Surface patterns on ablating materials are observed in high-speed ground and flight tests, but the mechanisms behind their formation are not known. In this paper, the origin of surface patterns is investigated via a local linear stability analysis of compressible laminar boundary layers over a flat camphor plate. The effects of sublimation and conjugate heat transfer are included both on the baseflow and the linear fluctuations. This framework identifies one mode that fully characterizes the stability of the surface, which becomes unstable when the wall temperature exceeds that of an adiabatic wall, $T_{ad}$. These findings are consistent with experimental observations, where laminar flow conditions at adiabatic wall temperatures are found to be stable. The analysis also reveals that the nature of this surface mode varies as a function of the oblique angle $ψ= \tan^{-1}{β/α}$, where $α$ and $β$ are the streamwise and spanwise wavenumbers. Specifically, for baseflow temperatures below $\approx 1.15~T_{ad}$, the surface mode is most unstable at $ψ= 0$. Conversely, above $\approx 1.15~T_{ad}$ the surface is most unstable near the sonic angle $ψ_s = \cos^{-1}(1/M_e)$, which is the angle at which the normal Mach number equals one. Finally, a critical wavenumber is identified (i.e., one at which the temporal growth rate reaches a maximum) that is in good agreement with available experimental observations of turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Data-driven model reduction via non-intrusive optimization of projection operators and reduced-order dynamics

Computing reduced-order models using non-intrusive methods is particularly attractive for systems that are simulated using black-box solvers. However, obtaining accurate data-driven models can be challenging, especially if the underlying systems exhibit large-amplitude transient growth. Although these systems may evolve near a low-dimensional subspace that can be easily identified using standard techniques such as Proper Orthogonal Decomposition (POD), computing accurate models often requires projecting the state onto this subspace via a non-orthogonal projection. While appropriate oblique projection operators can be computed using intrusive techniques that leverage the form of the underlying governing equations, purely data-driven methods currently tend to achieve dimensionality reduction via orthogonal projections, and this can lead to models with poor predictive accuracy. In this paper, we address this issue by introducing a non-intrusive framework designed to simultaneously identify oblique projection operators and reduced-order dynamics. In particular, given training trajectories and assuming reduced-order dynamics of polynomial form, we fit a reduced-order model by solving an optimization problem over the product manifold of a Grassmann manifold, a Stiefel manifold, and several linear spaces (as many as the tensors that define the low-order dynamics). Furthermore, we show that the gradient of the cost function with respect to the optimization parameters can be conveniently written in closed-form, so that there is no need for automatic differentiation. We compare our formulation with state-of-the-art methods on three examples: a three-dimensional system of ordinary differential equations, the complex Ginzburg-Landau (CGL) equation, and a two-dimensional lid-driven cavity flow at Reynolds number Re = 8300.

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

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