Broadly, I am interested in reduced order modeling, where one uses information about a large-scale dynamical system to find a much smaller system that well approximates the input-output mapping of the original system. This smaller system can then be simulated for a much smaller computational cost. The "information" used to construct such models can be an internal description of the system, or various types of input-output data.
Of particular interest to me is the case where one cannot assume to have an internal description of the system, but must construct reduced order models directly from data. Constructing such models poses intriguing numerical challenges, and much of my work thus far has been in tackling these challenges.
When constructing reduced order models from input-output data, one can choose to work with information in the frequency domain or time domain. Frequency domain techniques have many attractive properties, including bounds on the mismatch between the outputs of the reduced and true system for any input. Unfortunately, frequency information can at times be costly or impossible to obtain.
My work provides a numerically robust method to learn frequency information from time-domain data. This work enables the powerful tools of frequency-based reduced order modeling to be used in the case where only time-domain data is available.
The article is published in SISC, andMatLab code is available.
Finding reduced order models that are optimal in an appropriate sense has classically required access to a realization of the full order model. The Transfer Function Iterative Rational Krylov Algorithm (TF-IRKA) constructs optimal ROMs from frequency data, but in some practical scenarios, this information may be difficult or impossible to obtain.
In this work, we provide the necessary analysis to produce optimal reduced order models directly from a single time-domain simulation of the full-order model.
The article is published in SISC, and MatLab code is available.
Data-driven reduced order modeling is a tool to construct models of physical phenomena when access to governing equations or a state-space model is not possible, but input-output data is abundant. While there are many available methods capable of producing models that fit the given data to high accuracy, such methods typically do not preserve underlying differential structures which are important for interpretation of the models. A common structure to appear in mechanical and acoustic systems is an explicit dependence on the second time derivative. Inspired by the popular Adaptive Anderson-Antoulas (AAA) algorithm, we develop a structure-preserving second-order AAA algorithm which is capable of learning highly accurate second-order surrogate models from frequency response measurements.
An arXiv preprint is available, as well as the MatLab code.