Source author record

Mohammad Farazmand

Mohammad Farazmand appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
18topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

20 published item(s)

preprint2022arXiv

Model-assisted deep learning of rare extreme events from partial observations

To predict rare extreme events using deep neural networks, one encounters the so-called small data problem because even long-term observations often contain few extreme events. Here, we investigate a model-assisted framework where the training data is obtained from numerical simulations, as opposed to observations, with adequate samples from extreme events. However, to ensure the trained networks are applicable in practice, the training is not performed on the full simulation data; instead we only use a small subset of observable quantities which can be measured in practice. We investigate the feasibility of this model-assisted framework on three different dynamical systems (Rossler attractor, FitzHugh-Nagumo model, and a turbulent fluid flow) and three different deep neural network architectures (feedforward, long short-term memory, and reservoir computing). In each case, we study the prediction accuracy, robustness to noise, reproducibility under repeated training, and sensitivity to the type of input data. In particular, we find long short-term memory networks to be most robust to noise and to yield relatively accurate predictions, while requiring minimal fine-tuning of the hyperparameters.

preprint2022arXiv

Quantifying rare events in spotting: How far do wildfires spread?

Spotting refers to the transport of burning pieces of firebrand by wind which, at the time of landing, may ignite new fires beyond the direct ignition zone of the main fire. Spot fires that occur far from the original burn unit are rare but have consequential ramifications since their prediction and control remains challenging. To facilitate their prediction, we examine three methods for quantifying the landing distribution of firebrands: crude Monte Carlo simulations, importance sampling, and large deviation theory (LDT). In particular, we propose an LDT method that accurately and parsimoniously quantifies the low probability events at the tail of the landing distribution. In contrast, Monte Carlo and importance sampling methods are most efficient in quantifying the high probability landing distances near the mode of the distribution. However, they become computationally intractable for quantifying the tail of the distribution due to the large sample size required. We also show that the most probable landing distance grows linearly with the mean characteristic velocity of the wind field. Furthermore, defining the relative landed mass as the proportion of mass landed at a given distance from the main fire, we derive an explicit formula which allows computing this quantity as a function of the landing distribution at a negligible computational cost. We numerically demonstrate our findings on two prescribed wind fields.

preprint2022arXiv

Shape-morphing reduced-order models for nonlinear Schrodinger equations

We consider reduced-order modeling of nonlinear dispersive waves described by a class of nonlinear Schrodinger (NLS) equations. We compare two nonlinear reduced-order modeling methods: (i) The reduced Lagrangian approach which relies on the variational formulation of NLS and (ii) The recently developed method of reduced-order nonlinear solutions (RONS). First, we prove the surprising result that, although the two methods are seemingly quite different, they can be obtained from the real and imaginary parts of a single complex-valued master equation. Furthermore, for the NLS equation in a stationary frame, we show that the reduced Lagrangian method fails to predict the correct group velocity of the waves whereas RONS predicts the correct group velocity. Finally, for the modified NLS equation, where the reduced Lagrangian approach is inapplicable, the RONS reduced-order model accurately approximates the true solutions.

preprint2022arXiv

Stochastic compartmental models of COVID-19 pandemic must have temporally correlated uncertainties

Compartmental models are an important quantitative tool in epidemiology, enabling us to forecast the course of a communicable disease. However, the model parameters, such as the infectivity rate of the disease, are riddled with uncertainties, which has motivated the development and use of stochastic compartmental models. Here, we first show that a common stochastic model, which treats the uncertainties as white noise, is fundamentally flawed since it erroneously implies that greater parameter uncertainties will lead to the eradication of the disease. Then, we present a principled modeling of the uncertainties based on reasonable assumptions on the contacts of each individual. Using the central limit theorem and Doob's theorem on Gaussian Markov processes, we prove that the correlated Ornstein-Uhlenbeck process is the appropriate tool for modeling uncertainties in the infectivity rate. We demonstrate our results using a compartmental model of the COVID-19 pandemic and the available US data from the Johns Hopkins University COVID-19 database. In particular, we show that the white noise stochastic model systematically underestimates the severity of the Omicron variant of COVID-19, whereas the Ornstein-Uhlenbeck model correctly forecasts the course of this variant. Moreover, using an SIS model of sexually transmitted disease, we derive an exact closed-form solution for the asymptotic distribution of infected individuals. This analytic result shows that the white noise model underestimates the severity of the pandemic because of unrealistic noise-induced transitions. Our results strongly support the need for temporal correlations in modeling of uncertainties in compartmental models of infectious disease.

preprint2021arXiv

Investigating climate tipping points under various emission reduction and carbon capture scenarios with a stochastic climate model

We study the mitigation of climate tipping point transitions using an energy balance model. The evolution of the global mean surface temperature is coupled with the CO2 concentration through the green house effect. We model the CO2 concentration with a stochastic delay differential equation (SDDE), accounting for various carbon emission and capture scenarios. The resulting coupled system of SDDEs exhibits a tipping point phenomena: if CO2 concentration exceeds a critical threshold (around 478ppm), the temperature experiences an abrupt increase of about six degrees Celsius. We show that the CO2 concentration exhibits a transient growth which may cause a climate tipping point, even if the concentration decays asymptotically. We derive a rigorous upper bound for the CO2 evolution which quantifies its transient and asymptotic growths, and provides sufficient conditions for evading the climate tipping point. Combining this upper bound with Monte Carlo simulations of the stochastic climate model, we investigate the emission reduction and carbon capture scenarios that would avert the tipping point.

preprint2020arXiv

Mitigation of tipping point transitions by time-delay feedback control

In stochastic multistable systems driven by the gradient of a potential, transitions between equilibria is possible because of noise. We study the ability of linear delay feedback control to mitigate these transitions, ensuring that the system stays near a desirable equilibrium. For small delays, we show that the control term has two effects: i) a stabilizing effect by deepening the potential well around the desirable equilibrium, and ii) a destabilizing effect by intensifying the noise by a factor of $(1-τα)^{-1/2}$, where $τ$ and $α$ denote the delay and the control gain, respectively. As a result, successful mitigation depends on the competition between these two factors. We also derive analytical results that elucidate the choice of the appropriate control gain and delay that ensure successful mitigations. These results eliminate the need for any Monte Carlo simulations of the stochastic differential equations, and therefore significantly reduce the computational cost of determining the suitable control parameters. We demonstrate the application of our results on two examples.

preprint2020arXiv

Multiscale analysis of accelerated gradient methods

Accelerated gradient descent iterations are widely used in optimization. It is known that, in the continuous-time limit, these iterations converge to a second-order differential equation which we refer to as the accelerated gradient flow. Using geometric singular perturbation theory, we show that, under certain conditions, the accelerated gradient flow possesses an attracting invariant slow manifold to which the trajectories of the flow converge asymptotically. We obtain a general explicit expression in the form of functional series expansions that approximates the slow manifold to any arbitrary order of accuracy. To the leading order, the accelerated gradient flow reduced to this slow manifold coincides with the usual gradient descent. We illustrate the implications of our results on three examples.

preprint2019arXiv

Vortex boundaries as barriers to diffusive vorticity transport in two-dimensional flows

We put forward the idea of defining vortex boundaries in planar flows as closed material barriers to the diffusive transport of vorticity. Such diffusive vortex boundaries minimize the leakage of vorticity from the fluid mass they enclose when compared to other nearby material curves. Building on recent results on passive diffusion barriers, we develop an algorithm for the automated identification of such structures from general, two-dimensional unsteady flow data. As examples, we identify vortex boundaries as vorticity diffusion barriers in two flows: an explicitly known laminar flow and a numerically generated turbulent Navier--Stokes flow.

preprint2016arXiv

An adjoint-based approach for finding invariant solutions of Navier-Stokes equations

We consider the incompressible Navier--Stokes equations with periodic boundary conditions and time-independent forcing. For this type of flow, we derive adjoint equations whose trajectories converge asymptotically to the equilibrium and traveling wave solutions of the Navier--Stokes equations. Using the adjoint equations, arbitrary initial conditions evolve to the vicinity of a (relative) equilibrium at which point a few Newton-type iterations yield the desired (relative) equilibrium solution. We apply this adjoint-based method to a chaotic two-dimensional Kolmogorov flow. A convergence rate of 100% is observed, leading to the discovery of 21 new steady state and traveling wave solutions at Reynolds number Re=40. Some of the new invariant solutions have spatially localized structures that were previously believed to only exist on domains with large aspect ratios. We show that one of the newly found steady state solutions underpins the temporal intermittencies, i.e., high energy dissipation episodes of the flow. More precisely, it is shown that each intermittent episode of a generic turbulent trajectory corresponds to its close passage to this equilibrium solution.

preprint2016arXiv

Defining Coherent Vortices Objectively from the Vorticity

Rotationally coherent Lagrangian vortices are formed by tubes of deforming fluid elements that complete equal bulk material rotation relative to the mean rotation of the deforming fluid volume. We show that initial positions of such tubes coincide with tubular level surfaces of the Lagrangian-Averaged Vorticity Deviation (LAVD), the trajectory integral of the normed difference of the vorticity from its spatial mean. LAVD-based vortices are objective, i.e., remain unchanged under time-dependent rotations and translations of the coordinate frame. In the limit of vanishing Rossby numbers in geostrophic flows, cyclonic LAVD vortex centers are precisely the observed attractors for light particles. A similar result holds for heavy particles in anticyclonic LAVD vortices. We also establish a relationship between rotationally coherent Lagrangian vortices and their instantaneous Eulerian counterparts. The latter are formed by tubular surfaces of equal material rotation rate, objectively measured by the Instantaneous Vorticity Deviation (IVD). We illustrate the use of the LAVD and the IVD to detect rotationally coherent Lagrangian and Eulerian vortices objectively in several two- and three-dimensional flows.

preprint2016arXiv

Dynamical indicators for the prediction of bursting phenomena in high-dimensional systems

Drawing upon the bursting mechanism in slow-fast systems, we propose indicators for the prediction of such rare extreme events which do not require a priori known slow and fast coordinates. The indicators are associated with functionals defined in terms of Optimally Time Dependent (OTD) modes. One such functional has the form of the largest eigenvalue of the symmetric part of the linearized dynamics reduced to these modes. In contrast to other choices of subspaces, the proposed modes are flow invariant and therefore a projection onto them is dynamically meaningful. We illustrate the application of these indicators on three examples: a prototype low-dimensional model, a body forced turbulent fluid flow, and a unidirectional model of nonlinear water waves. We use Bayesian statistics to quantify the predictive power of the proposed indicators.

preprint2016arXiv

Kinematics of fluid particles on the sea surface. Hamiltonian theory

We derive the John-Sclavounos equations describing the motion of a fluid particle on the sea surface from first principles using Lagrangian and Hamiltonian formalisms applied to the motion of a frictionless particle constrained on an unsteady surface. The main result is that vorticity generated on a stress-free surface vanishes at a wave crest when the horizontal particle velocity equals the crest propagation speed, which is the kinematic criterion for wave breaking. If this holds for the largest crest, then the symplectic two-form associated with the Hamiltonian dynamics reduces instantaneously to that associated with the motion of a particle in free flight, as if the surface did not exist. Further, exploiting the conservation of the Hamiltonian function for steady surfaces and traveling waves we show that particle velocities remain bounded at all times, ruling out the possibility of the finite-time blowup of solutions.

preprint2015arXiv

Polar rotation angle identifies elliptic islands in unsteady dynamical systems

We propose rotation inferred from the polar decomposition of the flow gradient as a diagnostic for elliptic (or vortex-type) invariant regions in non-autonomous dynamical systems. We consider here two- and three-dimensional systems, in which polar rotation can be characterized by a single angle. For this polar rotation angle (PRA), we derive explicit formulas using the singular values and vectors of the flow gradient. We find that closed level sets of the PRA reveal elliptic islands in great detail, and singular level sets of the PRA uncover centers of such islands. Both features turn out to be objective (frame-invariant) for two-dimensional systems. We illustrate the diagnostic power of PRA for elliptic structures on several examples.

preprint2014arXiv

Asymptotic dynamics of inertial particles with memory

Recent experimental and numerical observations have shown the significance of the Basset--Boussinesq memory term on the dynamics of small spherical rigid particles (or inertial particles) suspended in an ambient fluid flow. These observations suggest an algebraic decay to an asymptotic state, as opposed to the exponential convergence in the absence of the memory term. Here, we prove that the observed algebraic decay is a universal property of the Maxey--Riley equation. Specifically, the particle velocity decays algebraically in time to a limit that is $\mathcal O(ε)$-close to the fluid velocity, where $0<ε\ll 1$ is proportional to the square of the ratio of the particle radius to the fluid characteristic length-scale. These results follows from a sharp analytic upper bound that we derive for the particle velocity. For completeness, we also present a first proof of existence and uniqueness of global solutions to the Maxey--Riley equation, a nonlinear system of fractional-order differential equations.

preprint2014arXiv

How coherent are the vortices of two-dimensional turbulence?

We use recent developments in the theory of finite-time dynamical systems to objectively locate the material boundaries of coherent vortices in two-dimensional Navier--Stokes turbulence. We show that these boundaries are optimal in the sense that any closed curve in their exterior will lose coherence under material advection. Through a detailed comparison, we find that other available Eulerian and Lagrangian techniques significantly underestimate the size of each coherent vortex.

preprint2014arXiv

The Maxey-Riley Equation: Existence, Uniqueness and Regularity of Solutions

The Maxey--Riley equation describes the motion of an inertial (i.e., finite-size) spherical particle in an ambient fluid flow. The equation is a second-order, implicit integro-differential equation with a singular kernel, and with a forcing term that blows up at the initial time. Despite the widespread use of the equation in applications, the basic properties of its solutions have remained unexplored. Here we fill this gap by proving local existence and uniqueness of weak solutions. For certain initial velocities between the particle and the fluid, the results extend to strong solutions. We also prove continuous differentiability of the weak and strong solutions with respect to their initial conditions. This justifies the search for coherent structures in inertial flows using the Cauchy--Green strain tensor.

preprint2013arXiv

Attracting and repelling Lagrangian coherent structures from a single computation

Hyperbolic Lagrangian Coherent Structures (LCSs) are locally most repelling or most attracting material surfaces in a finite-time dynamical system. To identify both types of hyperbolic LCSs at the same time instance, the standard practice has been to compute repelling LCSs from future data and attracting LCSs from past data. This approach tacitly assumes that coherent structures in the flow are fundamentally recurrent, and hence gives inconsistent results for temporally aperiodic systems. Here we resolve this inconsistency by showing how both repelling and attracting LCSs are computable at the same time instance from a single forward or a single backward run. These LCSs are obtained as surfaces normal to the weakest and strongest eigenvectors of the Cauchy-Green strain tensor.

preprint2013arXiv

Detecting invariant manifolds, attractors, and generalized KAM tori in aperiodically forced mechanical systems

We show how the recently developed theory of geodesic transport barriers for fluid flows can be used to uncover key invariant manifolds in externally forced, one-degree-of-freedom mechanical systems. Specifically, invariant sets in such systems turn out to be shadowed by least-stretching geodesics of the Cauchy-Green strain tensor computed from the flow map of the forced mechanical system. This approach enables the finite-time visualization of generalized stable and unstable manifolds, attractors and generalized KAM curves under arbitrary forcing, when Poincare maps are not available. We illustrate these results by detailed visualizations of the key finite-time invariant sets of conservatively and dissipatively forced Duffing oscillators.

preprint2013arXiv

Shearless transport barriers in unsteady two-dimensional flows and maps

We develop a variational principle that extends the notion of a shearless transport barrier from steady to general unsteady two-dimensional flows and maps defined over a finite time interval. This principle reveals that hyperbolic Lagrangian Coherent Structures (LCSs) and parabolic LCSs (or jet cores) are the two main types of shearless barriers in unsteady flows. Based on the boundary conditions they satisfy, parabolic barriers are found to be more observable and robust than hyperbolic barriers, confirming widespread numerical observations. Both types of barriers are special null-geodesics of an appropriate Lorentzian metric derived from the Cauchy--Green strain tensor. Using this fact, we devise an algorithm for the automated computation of parabolic barriers. We illustrate our detection method on steady and unsteady non-twist maps and on the aperiodically forced Bickley jet.