# The Azimuth Project Blog - fluid flows and infinite dimensional manifolds (part 1) (Rev #12, changes)

Showing changes from revision #11 to #12: Added | Removed | Changed

## Or: waves that take the shortest path trough infinity

This page is a blog article in progress, written by Tim van Beek.

Water waves can do a lot of things that light waves, for example, cannot, like “breaking”: In mathematical models this difference shows up through the kind of partial differential equation (PDE) that models the waves:

• Light waves are modelled by linear equations while

• water waves are modelled by nonlinear equations.

Physicists like to point out that linear equations model things that do not interact, while nonlinear equations model things that interact with each other. In quantum field theory this is usually dubbed as “free fields” and “interacting fields”.

Some of these nonlinear PDE turn up as geodesic equations of infinite dimensional Riemannian manifolds that describe fluid flows. This is a fascinating observation, that is due to the Russian mathematician Vladimir Arnold. In this blog post I would like to talk a little bit about the involved concepts and show you a little toy example.

#### Fluid Flow modelled by Diffeomorphisms

The Euler viewpoint of fluids is that there are fluid packages or “particles”. The fluid flow is described by specifying where each package or particle is at a given time t. When we start with a time $t_0 = 0$ on a given manifold $M$, the flow of every fluid package is described by a path on $M$ parametrized by time, and for every time $t \gt t_0$ there is a diffeomorphism $g^t$ of $M$ defined by the requirement that it maps the initial position $y$ of every fluid package to the position $x = g^t(y)$ at time $t$. This picture is taken from the book

We will take as a model of the domain of the fluid flow a compact Riemannian manifold $M$ that is also a Lie group. A fluid flow, as pictured above, is then a path in the diffeomorphism group  DM D M. In order to apply geometric concepts in this situation, we will have to turn  DM D M or some closed subgroup of it into a manifold, which will be infinite dimensional.

The curvature of such a manifold can provide a great deal about the stability of fluid flows: On a manifold with negative curvature geodesics will diverge from each other. If we can model fluid flows as geodesics in a Riemannian manifold and calculate the curvature, we could try to infer a bound on weather forecasts. In fact, that is what Vladimir Arnold did!

If you never thought about manifolds in infinite dimensions, you may feel a little bit insecure as to how the concepts that you know from differential geometry can be generalized from finite dimensions. At least I felt this way when I first read about it. But it turns out that the part of the theory one needs to know in order to understand Arnold’s insight is not that scary, so I will talk a little bit about it next.

#### What You Should Know about Infinite Dimensional Smooth Manifolds

The basic strategy when handling finite dimensional, smooth, real manifolds is that you have a complicated manifold $M$, but also locally for every point $p \in M$ a neighborhood $U$ and an isomorphism (a “chart”) of $U$ to an open subset of the vastly simpler space $\mathbb{R}^n$, the “model space”. These isomorphisms can be used to transport concepts from $\mathbb{R}^n$ to $M$. In infinite dimensions it is however not that clear what kind of model space $E$ should be taken in place of $\mathbb{R}^n$. What structure should $E$ have?

Since we would like to differentiate, we should for example be able to define the differential of a curve in E:$E$:

$\gamma: \mathbb{R} \to E$

If we write down the differential quotient

 \gamma'(t_0) := \coloneqq \lim_{t \to 0} \frac{1}{t} (\gamma(t_0 +t) - \gamma(t_0))

we see that we need to be able to add elements, have a scalar multiplication, and a topology such that the algebraic operations are continuous, in order to make sense of this. Sets $E$ with this structure are called topological vector spaces.

A curve that has infinite many differentials in the sense above at all points is called a smooth curve, just as in the finite dimensional case.

So $E$ should at least be a topological vector space. We can, of course, put more structure on $E$ to make it “more similar” to $\mathbb{R}^n$, and choose as model space in ascending order of generality:

1) A Hilbert space which has a scalar product,

2) a Banach space that does not have a scalar product, but a norm,

3) a Fréchet space that does not have a norm, but a metric,

4) a general topological vector space that need not be metrizable.

People talk accordingly of Hilbert, Banach and Fréchet manifolds. Since the space of smooth maps $C^{\infinity}(\mathbb{R}^n)$ of $\mathbb{R}^n \to \mathbb{R}$, for example, is not a Banach space but a Fréchet space, we should not expect that we can model diffeomorphism groups on Banach spaces, but on Fréchet spaces. So we will use the concept of Fréchet manifolds.

But if you are interested in a more general theory using locally convex topological vector spaces as model spaces, you can look it up here:

• Andreas Kriegl and Peter W. Michor: “The convenient setting of global analysis”, AMS 1999 (it is free for download here)

If you learn functional analysis, you will most likely start with operators on Hilbert spaces. One could say that the theory of topological vector spaces is about abstracting away as much structure from a Hilbert space and look what structure you need for important theorems to still hold true, like the open mapping/closed graph theorem. If you would like to learn more about this, my favorite book is this one:

• Francois Treves: “Topological Vector Spaces, Distributions and Kernels”, Dover Publications, 2006.

Since we replace the model space $\mathbb{R}^n$ with a Fréchet space $E$, there will be certain things that won’t work out as easily as for the finite dimensional $\mathbb{R}^n$, or not at all. The definition of a (topological) Fréchet manifold however is straight forward, you can keep it word for word. We will call a function between manifolds smooth if it maps smooth curves to smooth curves: In this way we get the definition of a smooth Fréchet manifold.

With tangent vectors, you may remember that there are several different ways to define them in the finite dimensional case, which turn out to be equivalent. Since there are situations in infinite dimensions where these definitions turn out to not be equivalent, I will be explicit and define tangent vectors in the “kinematic way”:

The (kinematic) tangent vector space $T_p M$ of a Fréchet manifold $M$ at a point $p$ consists of all pairs $(p, c'(0))$ where $c$ is a smooth curve

$c: \mathbb{R} \to M \; \text{with} \; c(0) = p$

With this definition, the set of pairs $(p, c'(0)), p \in M$ forms a Fréchet manifold, the tangent bundle  TM T M, just as in finite dimensions.

The first serious (more or less) problem we encounter is the definition of the cotangent bundle: $\mathbb{R}^n$ is isomorphic to its dual vector space. This is still true for every Hilbert space (this is known as the Riesz representation theorem). It fails already for Banach spaces: The dual space will still be a Banach space, but a Banach space does not need to be isomorphic to its dual (those who are are called reflexive).

With Fréchet spaces things are even a little bit worse, because the dual of a Fréchet space (which is not a Banach space) is not even a Fréchet space! Since I did not know that and could not find a reference, I asked about this on mathoverflow here and promply got an answer. Mathoverflow is a really amazing platform for this kind of question!

So, if we naively define the cotangent space as in finite dimensions by taking the dual space of every tangent space, then the cotangent bundle won’t be a Fréchet manifold.

We will therefore have to be careful with the definition of differential forms for Fréchet manifolds and define it explicitly:

A differential form (a one form) $\alpha$ is a smooth map

$\alpha: T M \to \mathbb{R}$

where  TM T M is the tangent bundle, such that $\alpha$ restricts to a linear map on every tangent space $T_p M$.

Another pitfall is that theorems from multivariate calculus may fail in Fréchet spaces, like the existence and uniqueness theorem of Picard-Lindelöf for ordinary differential equations. Things are much easier in Banach spaces: If you take a closer look at multivariate calculus, you will notice that a lot of definitions and theorems actually make use of the Banach space structure of $\mathbb{R}^n$ only, so that a lot generalizes straight forward to infinite dimensional Banach spaces. But that is less so for Fréchet spaces.

It is nevertheless possible to define both integrals and differentials that behave much in the expected way. You can find a nice exposition of how this can be done in this paper:

• Richard S. Hamilton: The Inverse Function Theorem of Nash and Moser (Bulletin (New Series) of the American Mathematical Society Volume 7, Number 1, July 1982)

The story starts with the definition of the directional derivative that can be done just like in finite dimensions:

Let F and G be Fréchet spaces,$F$ and $G$ be Fréchet spaces, $U \subseteq F$ open and $P: U \to G$ a continuous map. The derivative of P at the point$P$ at the point $f \in U$ in the direction $h \in F$ is the map

 DP: D P: U \times F \to G
 DP(f) D P(f) h = \lim_{t \to 0} \frac{1}{t} ( P(f + t h) - P(f))

A simple, but nontrivial example is the operator

$P: C^{\infty}[a, b] \to C^{\infty}[a, b]$
 P(f) := \coloneqq ff' f f'

with the derivative

 DP(f) D P(f) h = f'h + fh' f h'

So for Fréchet manifolds we could define directly what it means to be smooth, instead of using the indirect way we did above, namely that a smooth map is one that maps smooth paths to smooth paths. But the latter works for more general topological vector spaces, too, so I thought I should mention it.

By now you should feel reasonably comfortable with the notion of a Fréchet manifold, so let us talk about the kind of gadget that Arnold used to describe fluid flows: Diffeomorphism groups that are both (infinite dimensional) Riemannian manifolds and Lie groups.

#### Reformulating the Geodesic Equation for an Invariant Metric

If $M$ is both a Riemannian manifold and a Lie group, it is possible to define the concept of left or right invariant metric . A left or right invariant metric d on$d$ on $M$ is one that does not change if we multiply the arguments with a group element:

A metric d is left invariant iff for all$d$ is left invariant iff for all $g, h_1, h_2 \in G$:

$d (h_1, h_2) = d(g h_1, g h_2)$

Similarly, d is right invariant iff:$d$ is right invariant iff:

$d (h_1, h_2) = d(h_1 g, h_2 g)$

How does one get a one-sided invariant metric?

Here is one possibility: If you take a Lie group $M$ off the shelf, you get two automorphisms for free, namely the left and right multiplication by a group element $g$:

$L_g, R_g: M \to M$
 L_g(h) := \coloneqq gh g h
 R_g(h) := \coloneqq hg h g

Pictorially speaking, you can use the differentials of these to transport vectors from the Lie algebra $\mathfrak{m}$ of $M$ - which is the tangential space at the group identity $T_{id}M$ - to any other tangent space $T_g M$. If you can define a scalar product on the Lie algebra, you can use this trick to transport the scalar product to all the other tangential spaces by left or right multiplication, which will get you a left or right invariant metric.

To be more precise, for every tangent vectors $U, V$ of a tangent space $T_{g} M$ there are unique vectors $X, Y$ that are mapped to $U, V$ by the differential of the right multiplication $R_g$, that is

 dR_g d R_g X = U \; \text{and} \; dR_g d R_g Y = V

and we can define the scalar product of $U$ and $V$ to have the value of that of $X$ and $Y$:

 \langle U, V \rangle := \coloneqq \langle X, Y \rangle

This works for the left multiplication $L_g$, too, of course.

For a one-sided invariant metric, the geodesic equation looks somewhat simpler than for general metrics. Let us take a look at that:

On a Riemannian manifold $M$ with tangential bundle  TM T M there is a unique connection, the Levi-Civita connection, with the following properties for vector fields  X, Y, Z \in TM T M:

$Z \langle X, Y \rangle = \langle \nabla_Z X, Y \rangle + \langle X, \nabla_Z Y \rangle \; \text{(compatibility with the metric)}$
$\nabla_X Y - \nabla_Y X = [X, Y] \; \text{(torsion freeness)}$

If we combine both formulas we get

$2 \langle \nabla_X Y, Z \rangle = X \langle Y, Z \rangle + Y \langle Z, X \rangle - Z \langle X, Y \rangle + \langle [X, Y], Z \rangle - \langle [Y, Z], X \rangle + \langle [Z, X], Y \rangle$

If the scalar products are constant along every flow, i.e. the metric is (left or right) invariant, then the first three terms on the right hand side vanish, so that we get

$2 \langle \nabla_X Y, Z \rangle = \langle [X, Y], Z \rangle - \langle [Y, Z], X \rangle + \langle [Z, X], Y \rangle$

This latter formula can be written in a more succinct way if we introduce the coadjoint operator. Remeber the adjoint operator defined to be

$ad_X Z = [X, Z]$

With the help of the scalar product we can define the adjoint of this operator:

 \langle ad^*_X Y, Z \rangle := \coloneqq \langle Y, ad_X Z \rangle = \langle Y, [X, Z] \rangle

Then the formula above for the covariant derivative can be written as

$2 \langle \nabla_X Y, Z \rangle = \langle ad_X Y, Z \rangle - \langle ad^*_Y X, Z \rangle - \langle ad^*_X Y, Z \rangle$

Since the inner product is nondegenerate, we can eliminate $Z$ and get

$2 \nabla_X Y = ad_X Y - ad^*_X Y - ad^*_Y X$

A geodesic curve is one whose tangent vector $X$ is transported parallel to itself, that is we have

$\nabla_X X = 0$

Using the formula for the covariant derivative for an invariant metric above we get

$\nabla_X X = - ad^*_X X = 0$

as a reformulation of the geodesic equation.

For time dependent dynamical systems, we have the time axis as an additional dimension and every vector field has $\partial_t$ as an additional summand. So, in this case we get as geodesic equation (again: for an invariant metric)

$\nabla_X X = \partial_t X - ad^*_X X = 0$

#### A Simple Example: the Circle

As a simple example we will look at the circle $S^1$ and its diffeomorphism group  DS^1 D S^1. The Lie algebra $Vec(S^1)$ of  DS^1 D S^1 can be identified with the space of all vector fields on $S^1$. If we sloppily identify $S^1$ with $\mathbb{R}/\mathbb{Z}$ with coordinate $x$, then we can write for vector fields $X = u(x) \partial_x$ and $Y = v(x) \partial_x$ the commutator

$[X, Y] = (u v_x - u_x v) \partial_x$

where $u_x$ is short for the derivative:

 u_x := \coloneqq \frac{d u}{d x}

And of course we have a scalar product via

$\langle X, Y \rangle = \int_{S^1} u(x) v(x) d x$

which we can use to define either a left or a right invariant metric on  DS^1 D S^1, by transporting it via left or right multiplication to every tangent space.

Let us evaluate the geodesic equation for this example. We have to calculate the effect of the coadjoint operator:

 \langle ad^*_X Y, Z \rangle := \coloneqq \langle Y, ad_X Z \rangle = \langle Y, [X, Z] \rangle

If we write for the vector fields $X = u(x) \partial_x$, $Y = v(x) \partial_x$ and $Z = w(x) \partial_x$, this results in

$\langle ad^*_X Y, Z \rangle = \int_{S^1} v (u w_x - u_x w) d x = - \int_{S^1} (u v_x + 2 u_x v) w d x$

where the last step employs integration by parts and uses the periodic boundary condition $f(x + 1) = f(x)$ for the involved functions.

So we get for the coadjoint operator

$ad^*_X Y = - (u v_x + 2 u_x v) \partial_x$

Finally, the geodesic equation

$\partial_t X + \nabla_X X = 0$

turns out to be

$u_t + 3 u u_x = 0$

A similar equation,

$u_t + u u_x = 0$

is known as Hopf or inviscid Burgers’ equation. It looks simple, but its solutions can produce behaviour that looks like turbulence, so it is interesting in its own right.

If we take a somewhat more sophisticated diffeomorphism group, we can get slightly more complicated and therefore more interesting partial differential equations like the Korteweg-de Vries equation. But since this post is quite long already, that topic will have to wait for another post.

category: blog, climate