Skip to the Main Content

Note:These pages make extensive use of the latest XHTML and CSS Standards. They ought to look great in any standards-compliant modern browser. Unfortunately, they will probably look horrible in older browsers, like Netscape 4.x and IE 4.x. Moreover, many posts use MathML, which is, currently only supported in Mozilla. My best suggestion (and you will thank me when surfing an ever-increasing number of sites on the web which have been crafted to use the new standards) is to upgrade to the latest version of your browser. If that's not possible, consider moving to the Standards-compliant and open-source Mozilla browser.

October 25, 2021

The Kuramoto–Sivashinsky Equation (Part 2)

Posted by John Baez

I love the Kuramoto–Sivashinsky equation, beautifully depicted here by Thien An, because it’s one of the simplest partial differential equations that displays both chaos and a visible ‘arrow of time’. Now we’ll see that it’s also invariant under ‘Galilean transformations’ — that is, transformations into a moving frame of reference.

Last time I made some conjectures about the Kuramoto–Sivashinsky equation. Most notably, I conjecture that the ‘stripes’ you see above can appear or merge as time passes, but never disappear or split.

But I was quite confused for a while, because the Kuramoto–Sivashinsky equation comes in two forms. The integral form, which I spent most of my time discussing, is this:

h t+h xx+h xxxx+12(h x) 2=0 h_t + h_{x x} + h_{x x x x} + \frac{1}{2} (h_x)^2 = 0

where h(t,x)h(t,x) is a real function of two real variables, ‘time’ tt and ‘space’ xx. The derivative form is this:

u t+u xx+u xxxx+uu x=0u_t + u_{x x} + u_{x x x x} + u u_x = 0

where u(t,x)u(t,x) is again a real function of two real variables. If hh is a solution of the integral form,

u=h x u = h_x

is a solution of the derivative form. This is easy to see: just take the xx derivative of the equation with hh in it and see what you get.

By the way, beware of this: Thien An’s wonderful movie above shows the integral form of the equation, but she calls the function u.u.

Note that if hh has a local maximum as a function of x,x, then uu will go positive and then negative. This is important because it affects the definition of ‘stripe’. Here’s what stripes look in the derivative form of the Kuramoto–Sivashinsky equation, as computed by Steve Huntsman:

As you can see, the function uu goes negative and then positive as we increase xx moving through a stripe. This means that hh would have a local minimum.

Another thing puzzling me was that Steve Huntsman found solutions of the derivative form where the stripes in uu mostly move in the same direction as time passes:

This suggests that there’s some way to take a solution where the stripes aren’t moving much, and switch to a moving frame of reference to get a solution with moving stripes. And indeed I’ll show that’s true! But I’ll also show it’s not true for the integral form of the Kuramoto–Sivashinsky equation — at least not with periodic boundary conditions, which is what we’re using here.

Galilean transformations

Galileo was the guy who first emphasized that the laws of physics look the same in a moving frame of reference, but later Maxwell explained the idea very poetically:

Our whole progress up to this point may be described as a gradual development of the doctrine of relativity of all physical phenomena. Position we must evidently acknowledge to be relative, for we cannot describe the position of a body in any terms which do not express relations. There are no landmarks in space; one portion of space is exactly like every other portion, so that we cannot tell where we are. We are, as it were, on an unruffled sea, without stars, compass, sounding, wind or tide, and we cannot tell in which direction we are going. We have no log which we can case out to take a dead reckoning by; we may compute our rate of motion with respect to the neighboring bodies, but we do not know how these bodies may be moving in space.

All this was before Einstein. Before Einstein came along, a transformation into a moving frame of reference worked like this:

t t x x+vt\begin{array}{ccl} t &\mapsto& t \\ x &\mapsto& x + v t \end{array}

where vv is a constant, the velocity.

These are called ‘Galilean transformations’ And here’s something cool: the derivative form of the Kuramoto–Sivashinsky equation has Galilean transformations as symmetries! If u(t,x)u(t,x) is a solution, so is the boosted function

u(t,x)=u(t,x+vt)v u'(t,x) = u(t,x+v t) - v

The prime here does not mean derivative: uu' is just the name of our new boosted solution. To get this boosted solution, we do the obvious coordinate transformation into a moving frame of reference, but then a bit more: we must subtract the constant v.v.

Let’s see how this works! Suppose uu is a solution of the derivative form of the Kuramoto–Sivashinsky equation:

u t+u xx+u xxxx+uu x=0 u_t + u_{x x} + u_{x x x x} + u u_x = 0

Then defining

u(t,x)=u(t,x+vt)v u'(t,x) = u(t, x+v t) - v

we see the tt derivative of uu' has an extra vu xv u_x term:

u t(t,x)=u t(t,x+vt)+vu x(t,x+vt) u'_t(t,x) = u_t(t,x+v t) + v u_x(t,x+v t)

while its xx derivatives are simple:

u x(t,x)=u x(t,x+vt) u'_x(t,x) = u_x(t,x+v t)

and so on for the higher xx derivatives.

This lets us check that the boosted function uu' is also a solution:

u t+u xx+u xxxx+uu x = u t(t,x+vt)+vu x(t,x+vt) +u xx(t,x+vt)+u xxxx(t,x+vt) +(u(t,x+vt)v)u x(t,x+vt) = 0 \begin{array}{ccl} u'_t + u'_{x x} + u'_{x x x x} + u' u'_x &=& u_t(t,x+v t) + v u_x(t,x+v t) \\ \\ && + u_{x x}(t,x+v t) + u_{x x x x}(t,x+v t) \\ \\ && + (u(t,x+v t) - v) u_x(t,x+v t) \\ \\ &=& 0 \end{array}

Note how subtracting the constant vv when we transform uu exactly cancels the extra term vu xv u_x term we get when we transform the time derivative of u.u.

I got this idea from here:

thanks to someone on Twitter.

Does the integral form of the Kuramoto–Sivashinsky equation also have Galilean transformations as symmetries? Yes, but there’s a certain problem with them. Let’s see how it goes.

If h(t,x)h(t,x) is a solution, then so is this boosted function:

h(t,x)=h(t,x+vt)vx12tv 2 h'(t,x) = h(t,x+v t) - vx - \frac{1}{2}t v^2

But now the ‘fudge factors’ in the boosted function are more sneaky! We don’t just subtract a constant, as we did with the derivative form. We subtract a function that depends on both the space and time coordinates.

Let’s check. The first time derivative of hh' works like this:

h t(t,x)=h t(t,x+vt)+vh t(t,x+vt)12v 2 h'_t(t,x) = h_t(t,x+v t) + v h_t(t,x+v t) - \frac{1}{2} v^2

The first space derivative works like this:

h x(t,x)=h x(t,x+vt)v h'_x(t,x) = h_x(t,x+v t) - v

The second space derivative is simpler:

h xx(t,x)=h xx(t,x+vt) h'_{x x}(t,x) = h_{x x}(t,x+ v t)

and the higher space derivatives work the same way. Now suppose hh was a solution of the integral form of the Kuramoto–Sivashinsky equation:

h t+h xx+h xxxx+12(h x) 2=0 h_t + h_{x x} + h_{x x x x} + \frac{1}{2} (h_x)^2 = 0

We can use our formulas to check that hh' is a solution:

h t+h xx+h xxxx+12(h x) 2 = h t(t,x+vt)+vh t(t,x+vt)12v 2 +h xx(t,x+vt) +h xxxx(t,x+vt) +12(h x(t,x+vt)v) 2 = h x(t,x+vt)+h xx(t,x+vt) +h xxxx(t,x+vt)+12h x(t,x+vt) 2 = 0 \begin{array}{ccl} h'_t \! +\! h'_{x x}\! +\! h'_{x x x x} \!+ \! \frac{1}{2} (h'_x)^2 \! &=& h_t(t,x+v t) + v h_t(t,x+v t) \! - \!\frac{1}{2} v^2 \\ \\ && + h_{x x}(t,x+v t) \\ \\ && + h_{x x x x}(t,x+v t) \\ \\ && + \frac{1}{2}(h_x(t,x+v t) - v)^2 \\ \\ &=& h_x(t,x+v t) + h_{x x}(t,x+v t) \\ \\ && + h_{x x x x}(t,x+v t) + \frac{1}{2} h_x(t,x+v t)^2 \\ \\ &=& 0 \end{array}

The cancellations in the second step rely crucially on those sneaky fudge factors in the definition of the boosted solution h.h'. I chose those fudge factors precisely to make these cancellations happen.

Periodic boundary conditions

A lot of really interesting results about the Kuramoto–Sivashinsky equation involve looking at solutions that are periodic in space. For the derivative form this means

u(t,x+L)=u(t,x) u(t,x+L) = u(t,x)

It’s easy to see that if uu obeys this equation, so does its boosted version:

u(t,x)=u(t,x+vt)v u'(t,x) = u(t,x+v t) - v

After all, a quick calculation shows

u(t,x+L)=u(t,x+L+vt)v=u(t,x+vt)v=u(t,x+L) u'(t,x+L) = u(t,x+L+v t) - v = u(t,x+v t) - v = u'(t,x+L)

So for any spatially periodic solution with some stripes that are basically standing still, there’s a spatially periodic boosted version where they’re moving at velocity v.v. You can think of a spatially periodic solution as a function on the cylinder. In the boosted version, the stripes spiral around this cylinder like the stripes on a barber pole!

But for the integral form this doesn’t work! Suppose hh is a solution of the integral form that is periodic in space:

h(t,x+L)=h(t,x) h(t,x+L) = h(t,x)

Now its boosted version is defined by

h(t,x)=h(t,x+vt)vx12tv 2 h'(t,x) = h(t,x+ v t) - vx - \frac{1}{2}t v^2

and this is not periodic in space:

h(t,x+L) = h(t,x+L+vt)v(x+L)12tv 2 = h(t,x+vt)v(x+L)12tv 2 h(t,x+vt)=h(t,x) \begin{array}{ccl} h'(t,x+L) &=& h(t,x+L+v t) - v(x+L) - \frac{1}{2}t v^2 \\ \\ &=& h(t,x+v t) - v(x+L) - \frac{1}{2}t v^2 \\ \\ &\ne& h(t,x+v t) = h'(t,x) \end{array}

A conserved quantity

The derivative form of the Kuramoto–Shivashinsky equation has an unexpected symmetry under Galilean transformations, or boosts, even for periodic solutions. It also has an interesting conserved quantity, namely

0 Lu(t,x)dx \displaystyle{ \int_0^L u(t,x) \, d x }

To see this, note that

ddt 0 Lu(t,x)dx = 0 Lu t(t,x)dx = 0 L(u xx+u xxxx+uu x)dx = (u x+u xxx+12u 2)| 0 L = 0 \begin{array}{ccl} \displaystyle{ \frac{d}{dt} \int_0^L u(t,x) \, d x } &=& \displaystyle{ \int_0^L u_t(t,x) \, d x } \\ \\ &=& \displaystyle{ - \int_0^L \left( u_{x x} + u_{x x x x} + u u_x \right) \, d x } \\ \\ &=& \displaystyle{ \left. - \left( u_x + u_{x x x} + \frac{1}{2} u^2 \right) \right|_0^L } \\ \\ &=& 0 \end{array}

Note that if we boost a solution, replacing u(t,x)u(t,x) by u(t,x+vt)v,u(t,x+v t) - v, we also subtract a constant from this conserved quantity, namely vL.v L. So, it seems that the conserved quantity

0 Lu(t,x)dx \displaystyle{ \int_0^L u(t,x) \, d x }

is a way of measuring ‘how much the stripes are moving, on average’.

I’ll venture a guess that when

0 Lu(t,x)dx=0 \displaystyle{ \int_0^L u(t,x) \, d x = 0 }

the stripes are stationary, on average. One reason I’m willing to guess this is that this equation has another meaning, too.

I mentioned that we can get solutions uu of the derivative form from solutions hh of the integral form by differentiating them, like this:

u=h x u = h_x

Does every solution of the derivative form arise this way from a solution of the integral form? Given uu we can define hh by

h(t,x)= 0 xu(t,y)dy \displaystyle{ h(t,x) = \int_0^x u(t,y) \, d y }

and we can check it obeys the integral form of the Kuramoto–Sivashinsky equation. I’ve been showing you lots of calculations and we’re probably both getting tired of them, so I won’t actually go through this check. But here’s the interesting part. Suppose we restrict attention to spatially periodic solutions! If uu is spatially periodic, hh will be so iff

0 Lu(t,x)dx=0 \displaystyle{ \int_0^L u(t,x) \, d x = 0}

since we have

h(t,x+L) = 0 x+Lu(t,y)dy = 0 xu(t,y)dy+ x x+Lu(t,y)dy = h(t,x) \begin{array}{ccl} h(t,x+L) &=& \displaystyle{ \int_0^{x+L} u(t,y) \, d y } \\ \\ &=& \displaystyle{ \int_0^x u(t,y) \, d y + \int_x^{x+L} u(t,y) \, d y } \\ \\ &=& h(t,x) \end{array}

where in the last step, the first integral is h(t,x)h(t,x) by definition and the second is zero because

0 Lu(t,x)dx=0 \displaystyle{ \int_0^L u(t,x) \, d x = 0 }

and uu is periodic with period LL.

Summary

Limiting ourselves to spatially periodic solutions, we see:

  • The derivative form of the Kuramoto–Sivashinsky equation has symmetry under boosts.

  • The integral form does not.

  • Solutions of the integral form correspond precisely to solutions of the derivative form with

0 Lu(t,x)dx=0 \displaystyle{ \int_0^L u(t,x) \, d x = 0 }

  • Boosting a solution of the derivative form with

0 Lu(t,x)dx=0 \displaystyle{ \int_0^L u(t,x) \, d x = 0 }

gives a solution where this quantity is not zero.

Given all this, I think the good theorems about spatially periodic solutions of the Kuramoto–Sivashinsky equation will focus on solutions of the integral form, or equivalently solutions of the differential form with

0 Lu(t,x)dx=0 \displaystyle{ \int_0^L u(t,x) \, d x = 0 }

Epilogue: symmetries of the heat equation

At first I was shocked that the Kuramoto–Sivashinsky equation had Galilean transformations as symmetries, because it’s a relative of the heat equation, and the heat equation obviously doesn’t have Galilean transformations as symmetries: you can’t get a solution of the heat equation where a wave of heat moves along at constant velocity as it spreads out.

But then I remembered that the heat equation

u t=u xx u_t = u_{x x}

is very similar to Schrödinger’s equation:

u t=iu xx u_t = i u_{x x}

And Schrödinger’s equation does have Galilean transformations as symmetries, since it describes a free particle! Physicists know how they work. You don’t just replace u(t,x)u(t,x) with u(t,x+vt),u(t,x+v t), you also multiply it by a suitable complex phase that depends on tt and xx. That is, you multiply u(t,x+vt)u(t,x+v t) by the exponential of some imaginary-valued function of tt and x.x.

So then I realized that the heat equation does have Galilean symmetries! They work just as for Schrödinger’s equation, but now you have to multiply u(t,x+vt)u(t,x+v t) by the exponential of some real-valued function of tt and x.x.

This grows exponentially as you move in one direction in space, so people thinking about heat often don’t think about such solutions. Yes, you can get a ‘wave of heat that moves along at constant velocity as it spreads out’, but it’s like a huge tsunami wave rolling in from infinity!

I was very proud of myself for having discovered this weird Galilean invariance of the heat equation, and I posed it as a puzzle on Twitter. But then Tom Price pointed out a paper about it:

It turns out the Galilean invariance of the heat equation, and much more, was discovered by Sophus Lie in 1882!

The so-called ‘Schrödinger group’ also includes certain fractional linear transformations of the plane. So this raises the question of whether the Kuramoto–Sivashinsky equation has even more symmetries than those in the Galilei group (translations, spatial reflections and Galilean transformations).

But it seems this too has been studied:

This discusses the analogue of the Kuramoto–Sivashinsky equation in two dimensions of space and one of time. It uses some techniques to compute its symmetry group, but I don’t see any sign of invariance under a big group like the ‘Schrödinger group’.

Posted at October 25, 2021 3:07 AM UTC

TrackBack URL for this Entry:   https://golem.ph.utexas.edu/cgi-bin/MT-3.0/dxy-tb.fcgi/3362

2 Comments & 0 Trackbacks

Jared

The heat equation is equivalent to Burgers equation u t+uu x=u xxu_t + u u_x =u_{xx} under the Cole Hopf transformation and Burgers enjoys the same Galilean invariance as Kuramoto-Sivashinsky, so that is one way to see it.

Posted by: Jared on October 29, 2021 10:10 PM | Permalink | Reply to this

Jared

Tried to post this once but it failed for some reason.

The heat equation can be mapped the Burger’s equation u t+uu x=u xxu_t + u u_x = u_{xx}. Burger’s equation enjoys the same Galilean invariance as (derivative) Kuramoto-Sivashinsky, so from that point of view it is not surprising that the heat equation has a Galilean invariance.

Posted by: Jared on October 29, 2021 10:27 PM | Permalink | Reply to this

Post a New Comment