Eric Forgoston

Robotic Tracking of Coherent Structures in Flows

Lagrangian coherent structures (LCSs) are separatrices that delineate dynamically distinct regions in general dynamical systems and can be viewed as the extensions of stable and unstable manifolds to general time-dependent systems. Identifying LCS in dynamical systems is useful for many applications, including oceanography and weather prediction. In this paper, we present a collaborative robotic control strategy that is designed to track stable and unstable manifolds in dynamical systems, including ocean flows. The technique does not require global information about the dynamics, and is based on local sensing, prediction, and correction. The collaborative control strategy is implemented with a team of three robots to track coherent structures and manifolds on static flows, a time-dependent model of a wind-driven double-gyre flow often seen in the ocean, experimental data that are generated by a flow tank, and actual ocean data. We present simulation results and discuss theoretical guarantees of the collaborative tracking strategy.

Converging towards the optimal path to extinction

Extinction appears ubiquitously in many fields, including chemical reactions, population biology, evolution and epidemiology. Even though extinction as a random process is a rare event, its occurrence is observed in large finite populations. Extinction occurs when fluctuations owing to random transitions act as an effective force that drives one or more components or species to vanish. Although there are many random paths to an extinct state, there is an optimal path that maximizes the probability to extinction. In this paper, we show that the optimal path is associated with the dynamical systems idea of having maximum sensitive dependence to initial conditions. Using the equivalence between the sensitive dependence and the path to extinction, we show that the dynamical systems picture of extinction evolves naturally towards the optimal path in several stochastic models of epidemics.

Initial-value problem for three-dimensional disturbances in a compressible boundary layer

An initial-value problem is formulated for a three-dimensional wave packet in a compressible boundary layer flow. The problem is solved using a Laplace transform with respect to time and Fourier transforms with respect to the streamwise and spanwise coordinates. The solution can be presented as a sum of modes consisting of continuous and discrete spectra of temporal stability theory. Two discrete modes, known as mode S and mode F, are of interest in high-speed flows since they may be involved in a laminar-turbulent transition scenario. The continuous and discrete spectrum are analyzed numerically for a hypersonic flow with Mach number M=5.6. The following features are revealed: (1) The synchronism of mode S with acoustic waves at a streamwise wave number α→0 is primarily two-dimensional; (2) at high angles of disturbance propagation, mode F is no longer synchronized with entropy and vorticity waves; (3) at high angles of disturbance propagation, the synchronism between mode S and mode F is not accompanied...

Robotic manifold tracking of coherent structures in flows

Tracking Lagrangian coherent structures in dynamical systems is important for many applications such as oceanography and weather prediction. In this paper, we present a collaborative robotic control strategy designed to track stable and unstable manifolds. The technique does not require global information about the fluid dynamics, and is based on local sensing, prediction, and correction. The collaborative control strategy is implemented on a team of three robots to track coherent structures and manifolds on static flows as well as a noisy time-dependent model of a wind-driven double-gyre often seen in the ocean. We present simulation and experimental results and discuss theoretical guarantees of the collaborative tracking strategy.

Set-based corral control in stochastic dynamical systems: Making almost invariant sets more invariant

We consider the problem of stochastic prediction and control in a time-dependent stochastic environment, such as the ocean, where escape from an almost invariant region occurs due to random fluctuations. We determine high-probability control-actuation sets by computing regions of uncertainty, almost invariant sets, and Lagrangian coherent structures. The combination of geometric and probabilistic methods allows us to design regions of control, which provide an increase in loitering time while minimizing the amount of control actuation. We show how the loitering time in almost invariant sets scales exponentially with respect to the control actuation, causing an exponential increase in loitering times with only small changes in actuation force. The result is that the control actuation makes almost invariant sets more invariant.

Accurate noise projection for reduced stochastic epidemic models.

We consider a stochastic susceptible-exposed-infected-recovered (SEIR) epidemiological model. Through the use of a normal form coordinate transform, we are able to analytically derive the stochastic center manifold along with the associated, reduced set of stochastic evolution equations. The transformation correctly projects both the dynamics and the noise onto the center manifold. Therefore, the solution of this reduced stochastic dynamical system yields excellent agreement, both in amplitude and phase, with the solution of the original stochastic system for a temporal scale that is orders of magnitude longer than the typical relaxation time. This new method allows for improved time series prediction of the number of infectious cases when modeling the spread of disease in a population. Numerical solutions of the fluctuations of the SEIR model are considered in the infinite population limit using a Langevin equation approach, as well as in a finite population simulated as a Markov process.

Maximal Sensitive Dependence and the Optimal Path to Epidemic Extinction

Extinction of an epidemic or a species is a rare event that occurs due to a large, rare stochastic fluctuation. Although the extinction process is dynamically unstable, it follows an optimal path that maximizes the probability of extinction. We show that the optimal path is also directly related to the finite-time Lyapunov exponents of the underlying dynamical system in that the optimal path displays maximum sensitivity to initial conditions. We consider several stochastic epidemic models, and examine the extinction process in a dynamical systems framework. Using the dynamics of the finite-time Lyapunov exponents as a constructive tool, we demonstrate that the dynamical systems viewpoint of extinction evolves naturally toward the optimal path.

Coherent Pattern Prediction in Swarms of Delay-Coupled Agents

We consider a general swarm model of self-propelling agents interacting through a pairwise potential in the presence of noise and communication time delay. Previous work has shown that a communication time delay in the swarm induces a pattern bifurcation that depends on the size of the coupling amplitude. We extend these results by completely unfolding the bifurcation structure of the mean field approximation. Our analysis reveals a direct correspondence between the different dynamical behaviors found in different regions of the coupling-time delay plane with the different classes of simulated coherent swarm patterns. We derive the spatiotemporal scales of the swarm structures, as well as demonstrate how the complicated interplay of coupling strength, time delay, noise intensity, and choice of initial conditions can affect the swarm. In particular, our studies show that for sufficiently large values of the coupling strength and/or the time delay, there is a noise intensity threshold that forces a transition of the swarm from a misaligned state into an aligned state. We show that this alignment transition exhibits hysteresis when the noise intensity is taken to be time dependent.

Noise, bifurcations, and modeling of interacting particle systems

We consider the stochastic patterns of a system of communicating, or coupled, self-propelled particles in the presence of noise and communication time delay. For sufficiently large environmental noise, there exists a transition between a translating state and a rotating state with stationary center of mass. Time delayed communication creates a bifurcation pattern dependent on the coupling amplitude between particles. Using a mean field model in the large number limit, we show how the complete bifurcation unfolds in the presence of communication delay and coupling amplitude. Relative to the center of mass, the patterns can then be described as transitions between translation, rotation about a stationary point, or a rotating swarm, where the center of mass undergoes a Hopf bifurcation from steady state to a limit cycle. Examples of some of the stochastic patterns will be given for large numbers of particles.

Escape Rates in a Stochastic Environment with Multiple Scales

We consider a stochastic environment with two time scales and outline a general theory that compares two methods to reduce the dimension of the original system. The first method involves the computation of the underlying deterministic center manifold followed by a naive replacement of the stochastic term. The second method allows one to more accurately describe the stochastic effects and involves the derivation of a normal form coordinate transform that is used to find the stochastic center manifold. The results of both methods are used along with the path integral formalism of large fluctuation theory to predict the escape rate from one basin of attraction to another. The general theory is applied to the example of a surface flow described by a generic, singularly perturbed, damped, nonlinear oscillator with additive, Gaussian noise. We show how both nonlinear reduction methods compare in escape rate scaling. Additionally, the center manifolds are shown to predict high pre-history probability regions of escape. The theoretical results are confirmed using numerical computation of the mean escape time and escape prehistory, and we briefly discuss the extension of the theory to stochastic control.