Slow manifold

Definition

Consider the dynamical system

for an evolving state vector x → ( t ) and with equilibrium point x → ∗ . Then the linearization of the system at the equilibrium point is

The matrix A defines four invariant subspaces characterized by the eigenvalues λ of the matrix: as described in the entry for the center manifold three of the subspaces are the stable, unstable and center subspaces corresponding to the span of the eigenvectors with eigenvalues λ that have real part negative, positive, and zero, respectively; the fourth subspace is the slow subspace given by the span of the eigenvectors, and generalized eigenvectors, corresponding to the eigenvalue λ = 0 precisely. The slow subspace is a subspace of the center subspace, or identical to it, or possibly empty.

Correspondingly, the nonlinear system has invariant manifolds, made of trajectories of the nonlinear system, corresponding to each of these invariant subspaces. There is an invariant manifold tangent to the slow subspace and with the same dimension; this manifold is the slow manifold.

Stochastic slow manifolds also exist for noisy dynamical systems (stochastic differential equation), as do also stochastic center, stable and unstable manifolds. Such stochastic slow manifolds are similarly useful in modeling emergent stochastic dynamics, but there are many fascinating issues to resolve such as history and future dependent integrals of the noise.

Simple case with two variables

The coupled system in two variables x ( t ) and y ( t )

has the exact slow manifold y = x 2 on which the evolution is d x / d t = − x 3 . Apart from exponentially decaying transients, this slow manifold and its evolution captures all solutions that are in the neighborhood of the origin. The neighborhood of attraction is, roughly, at least the half-space y > − 1 / 2 .

Slow dynamics among fast waves

Edward Norton Lorenz introduced the following dynamical system of five equations in five variables to explore the notion of a slow manifold of quasi-geostrophic flow

Linearized about the origin the eigenvalue zero has multiplicity three, and there is a complex conjugate pair of eigenvalues, ± i . Hence there exists a three-dimensional slow manifold (surrounded by 'fast' waves in the X and Z variables). Lorenz later argued a slow manifold did not exist! But normal form arguments suggest that there is a dynamical system that is exponentially close to the Lorenz system for which there is a good slow manifold.

Eliminate an infinity of variables

In modeling we aim to simplify enormously. This example uses a slow manifold to simplify the 'infinite dimensional' dynamics of a partial differential equation to a model of one ordinary differential equation. Consider a field u ( x , t ) undergoing the nonlinear diffusion

with Robin boundary conditions

Parametrising the boundary conditions by b empowers us to cover the insulating Neumann boundary condition case b = 0 , the Dirichlet boundary condition case b = 1 , and all cases between.

Now for a marvelous trick, much used in exploring dynamics with bifurcation theory. Since parameter b is constant, adjoin the trivially true differential equation

Then in the extended state space of the evolving field and parameter, ( b , u ( x ) ) , there exists an infinity of equilibria, not just one equilibrium, with b = 0 (insulating) and u = constant, say u = a . Without going into details, about each and every equilibria the linearized diffusion has two zero eigenvalues and for a > 0 all the rest are negative (less than − π 2 a / 4 ). Thus the two-dimensional dynamics on the slow manifolds emerge (see emergence) from the nonlinear diffusion no matter how complicated the initial conditions.

Here one can straightforwardly verify the slow manifold to be precisely the field u ( x , t ) = a ( t ) ( 1 − b x 2 ) where amplitude a evolves according to

That is, after the initial transients that by diffusion smooth internal structures, the emergent behavior is one of relatively slow decay of the amplitude ( a ) at a rate controlled by the type of boundary condition (constant b ).

Notice that this slow manifold model is global in a as each equilibria is necessarily in the slow subspace of each other equilibria, but is only local in parameter b . We cannot yet be sure how large b may be taken, but the theory assures us the results do hold for some finite parameter b .

Perhaps the simplest nontrivial stochastic slow manifold

Stochastic modeling is much more complicated—this example illustrates just one such complication. Consider for small parameter ϵ the two variable dynamics of this linear system forced with noise from the random walk W ( t ) :

One could simply notice that the Ornstein–Uhlenbeck process y is formally the history integral

and then assert that x ( t ) is simply the integral of this history integral. However, this solution then inappropriately contains fast time integrals, due to the exp ⁡ ( s − t ) in the integrand, in a supposedly long time model.

Alternatively, a stochastic coordinate transform extracts a sound model for the long term dynamics. Change variables to ( X ( t ) , Y ( t ) ) where

then the new variables evolve according to the simple

In these new coordinates we readily deduce Y ( t ) → 0 exponentially quickly, leaving X ( t ) = ϵ W ( t ) undergoing a random walk to be the long term model of the stochastic dynamics on the stochastic slow manifold obtained by setting Y = 0 .

A web service constructs such slow manifolds in finite dimensions, both deterministic and stochastic.