Stochastic modeling and UQ for strongly nonlinear systems

A strongly nonlinear system is a chaotic dynamical system which exhibits energy transfer between various scales in time and space. The well-known examples of this behavior commonly appear in engineering and geophysical turbulence. In this project, we work on capturing relatively low-dimensional but statistically accurate models of these systems with the aid of data [1]. The end goal is to use these models, instead of very costly physics-based models, to characterize and predict extreme events (e.g. heat waves), study the statistical response to external output (e.g. climate trends), and synthesize long-time control strategies (e.g. mean drag reduction).

Characterizing the probability of extreme events for climate reanalysis data [1].
Related publications:
  1. H. Arbabi and T. Sapsis, Data-driven modeling of strongly nonlinear chaotic systems with non-Gaussian statistics, submitted, 2019.

Learning coarse-grained PDEs for particle-based or agent-based systems

Many physical, chemical and biological systems are described at a microscopic level (e.g. molecular reaction, response of neurons to stimula and vehicle platoons) but it is suspected that they may admit evolution equations at much larger length scales, e.g. in the form of PDEs. Discovering such coarse-grained evolution laws can lead to huge savings in computation-intensive tasks like prediction, design optimization and control. There are two challenges in machine learning these coarse-grained evolution laws from data: 1- the need for ample trainging data from microscopic simulations, and 2- ambiguity in the choice of coarse-grained dependent variable (e.g. density, pair probability etc.) in those equations. We combine modern tools of applied mathematics; such as equation-free numerics, optimal transport of probabilities, and manifold learning; to address these challenges.

We use unbalanced optimal transport + manifold learning on snapshots of particle distributions to discover the coarse-grained (PDE) variable and its relation to particle distribution moments [1].
Related publications:
  1. H. Arbabi and I. Kevrekidis, Particles to PDEs parsimoniously, submitted, 2020.
  2. H. Arbabi, J. Bunder, G. Samaey, A. Roberts and I. Kevrekidis, Linking Machine learning with multiscale numerics: data-driven discovery of homogenized equations, JOM, 2020

Data-driven analysis and control of high-dimensional nonlinear systems

Why do we hear so much about data? One reason is that we have a lot of it, but also in the context of dynamical systems, data is becoming instrumental in dealing with problems that have large dimensions and contain uncertainy. The traditional approach in dynamical systems theory has been centered around geometric analysis in the state space, but a lot of systems that we are interested in right now have either too many degrees of freedom to model directly (e.g. climate, turbulent flows and chemical kinetics) or lack reliable models for state-space analysis (e.g. biological and social networks). Data helps us detect and exploit low-dimensional structures in the former group, or find surrogate but amenable models for the latter. Two sets of techniques have contributed the most to data-driven analysis and control in the recent years: 1-operator-theoretic methods and 2-machine learning. Our goal is to push the boundaries of application for these methods, or invent others, that makes them applicable to problems ranging from turbulence and climate to cell motility and health monitoring.

Discovering the linearly evolving coordinates (i.e. Koopman eigenfunctions) on attractors of fluid flows and using them to extract quasi-periodic factor [2].
Related publications:
  1. H. Arbabi, M. Korda and I. Mezić, A data-driven Koopman model predictive control framework for nonlinear flows, submitted, conference version, preprint, 2018.
  2. H. Arbabi and I. Mezić, Study of dynamics in post-transient flows using Koopman mode decomposition, Physical Review Fluids, preprint , 2017.
  3. N. Govindarajan, H. Arbabi, L. van Blargian, T. Matchen, E. Tegling and I. Mezić, An operator-theoretic viewpoint to non-smooth dynamical systems: Koopman analysis of a hybrid pendulum, IEEE 55th Conference on Decision and Control (CDC), 2016, Preprint.

Numerical approximation of operator-theoretic objects for dynamical systems

The goal of this project is to devise and analyze numerical algorithms that enable data-driven analysis and control of dynamical systems. For example, Dynamic Mode Decomposition (DMD) is a linear data decomposition algorithm that has become the primary tool for computing the Koopman eigenfunctions (i.e. linearly evolving coordinates) in simple nonlinear systems. In previous works, we have established rigorous results on the connection between this algorithm [1,2], and some other signal processing methods [3], with operator-theoretic viewpoint for nonlinear systems. We are interested in improving the performance of numerical toolbox for more efficient and robust data-driven modeling, for example by including optimal sampling of state and inputs.

Combination of linear algebra and ergodic theory makes a strong tool for operator-theoretic analysis of dynamical systems [1]: data-driven computation of an orthogonal function basis for observables on Lorenz attractor via Hankel-SVD (top) and asymptotic phase (i.e. what initial conditions synchronize) in limit-cycling systems via Hankel-DMD (bottom).
Related publications:
  1. H. Arbabi and I. Mezić, Ergodic theory, Dynamic Mode Decomposition and computation of spectral properties of the Koopman operator, SIAM Journal on Applied Dynamical Systems, 2017.
  2. N. I. Mezić and H. Arbabi, On the Computation of Isostables, Isochrons and Other Spectral Objects of the Koopman Operator Using the Dynamic Mode Decomposition, NOLTA, 2017.
  3. H. Arbabi and I. Mezić, H. Arbabi and I. Mezić, Study of dynamics in post-transient flows using Koopman mode decomposition, Physical Review Fluids, preprint , 2017.

Mixing and transport in fluid flows

Mixing is a crucial part of many processes related to fluids like movement of oil slicks in the ocean, processing of polymers and metals, and biodmedical testing. Despite its prevalence, the process of mixing in flows, even in the absence of diffusion and particle-size effects, is not fully understood. Most of rigorous analysis in the field of advective mixing comes from the theory of chaotic advection which treats the fluid flow as a dynamical system with trajectories that describe the paths of tracers. A strong theme of our work is to understand mixing by combining the fluid mechanical apsect (i.e. Navier-Stokes dynamics) with the methods of chaotic advection. For example, this research has led to extension of classical Prandtl-Batchelor theory to unsteady 2D flows [1] which predicts that at high Reynolds the rotational structures in the flow resemble rigid body rotation. Interstingly, resonance arguments from chaotic advection show that this constant-vorticity condition implies that mixing in the core of those flows would be weaker than wall-adjacent areas, which is radically different from mixing in low-Reynolds flows [2].

On a more applied side, we collaborated with a group in Croatia to design an efficient algorithm for search on the ocean surface. The new algorithm [3], which combines ideas from statistical mechanics and ocean mixing, has shown great promise in a computational replication of search for the missing Malaysian Airliner MH370. Another exciting collaboration was formed with scientists in Marine Science Institute to devise a framework for prediction of oil spill movement on the ocean [4]. This problem was motivated by the 2015 Refugio oil spill in northern Santa Barbara county.

2D rotational flows at high Reynolds develop constant-vorticity cores [2] (top row ) which leads to rigid-body kinematics. This rigid-body type of motion leads to weaker mixing in the core (hypergraphs in the middle row) [1]. Bottom row shows quantitaive characterization of mixing in cavity flow which exhibits this behavior.
Related publications:
  1. H. Arbabi and I. Mezić, Prandtl-Batchelor theorem for flows with quasi-periodic time dependence, Journal of Fluid Mechanics, preprint, 2019.
  2. H. Arbabi and I. Mezić, Spectral analysis of mixing in 2D high-Reynolds flows, submitted, preprint, 2019.
  3. S. Ivić, B. Crnković, H. Arbabi, S. Loire, P. Clary and I. Mezić, Search strategy in a complex and dynamic environment: the MH370 case, submitted, preprint, 2020
  4. B. Emery, L. Washburn, I. Mezić, S. Loire, H. Arbabi, C. Ohlmann, J. Harlan, Oil Spill Trajectories from HF Radars: Applied Dynamical Systems Methods vs. a Lagrangian Stochastic Model, AGU Fall Meeting, Poster Session, 2016.