The shocking behaviour of moving fluids

Figure 1: The shock wave caused by a supersonic jet.

Originating authors are David Mumford and Christiane Rousseau.

Foreword: This vignette is more difficult than others. However, in a few pages, it tells you how to explain in simple terms one of the most difficult open problems at the beginning of the 21st century. The vignette contains enrichment material, that you can choose to read or skip. The editors of the Klein blog hesitated for a while posting this vignette. After testing it with teachers during two Klein workshops, who expressed that they enjoy being challenged by more difficult vignettes, they decided to test it on the blog. They are eager to hear your comments, and if some of you were motivated by this topic.

You have probably heard some planes breaking the sound barrier. What does that mean? It means that a shock wave in the atmosphere is created as in Figure 1. But what is a shock wave? Imagine heavy traffic on highways as a wave. A shock wave corresponds to collisions. To explain this, we develop our intuition with a 1D model: traffic on a one lane road at different speeds. You know that collisions could occur if drivers do not adjust their speed. The atmosphere is a fluid, and traffic is a rough model of 1D fluid which is convenient to develop our intuition. Under which conditions do shock waves or other singularities occur in fluids? A million dollar prize is offered for answering this question. This is what we are going to explain you.

1. The 1-D models

Burger’s equation

Let us construct a simple traffic model. We have cars at each position x along a line representing a road. At time t=0, the car at position x_0 has speed V_0(x_0). The car starting at x_0 will be at position x(t,x_0) at time t. Let V(x(t,x_0),t) be the speed of the car at position x(t,x_0) and time t. The hypothesis is that cars move at constant speed along a straight road. Hence, we have V(x(t,x_0),t)=V_0(x_0), which means that this function is constant in t. Hence, its derivative should be identically zero! By the chain rule for functions of several variables we have

    \begin{align*}\begin{split}\frac{d}{dt} V(x(t,x_0),t)&=\frac{\partial V}{\partial x} (x(t,x_0),t)\frac{\partial x}{\partial t}(t,x_0) +\frac{\partial V}{\partial t} (x(t,x_0),t)\\&= \frac{\partial V}{\partial x} (x(t,x_0),t)V_0(x_0) +\frac{\partial V}{\partial t} (x(t,x_0),t). \end{split}\end{align*}

Simply writing that the position on the road is given by x, and using that V(x(t,x_0),t)=V_0(x_0), the traffic equation becomes what is called Burger’s equation:

(1)   \begin{equation*} \frac{\partial V}{\partial t}+ \frac{\partial V}{\partial x} \cdot V=0.\end{equation*}

This equation says that each car maintains a constant speed as it moves long the road.

Another way of saying this is that if a car at position x_0 starts off with speed
V_0 = V(x_0,0), then at time t it will be at position x_0+tV_0 still with the same speed:
V(x_0+tV_0,t)=V_0 for all t.
Of course this is not a great idea! Collisions will occur. Figure 2 shows one such resulting chaos.

Figure 2: A shock wave in Burger’s equation for the initial speed given by V(x)=exp(-x²). Shown are the world lines of cars initially equally spaced on the road. The heavy black line is where the collisions happen.

Describing the shocks

What happens when we approach the shock? Suppose that at the beginning (t=0) two points x_0<y_0 had velocities V(x_0,0)>V(y_0,0). Let us denote the position of x_0 at time t by x_t. Similarly y_t denotes the position of y_0 at time t. Then x_t gets closer and closer to y_t and just before the shock x_t is very close to y_t. Yet, the difference of their velocities V(x_t,t)-V(y_t,t)=V(x_0,0)-V(y_0,0) has remained the same. This means that the mean slope of the function V(x,t) as a function of x for fixed t is becoming very large (see Figure 3 which shows the graph of V(x,t) as a function of x for several values of t). Hence, we see that \frac{\partial V} {\partial x} is becoming infinite when we approach the shock. This means that our differential equation (1) has a singularity at the shock.

Figure 3: In the same situation as figure 2, we show the graph of velocity as a function of space for six representative times. Note how the velocity becomes discontinuous when the shock wave develops: the slope of the velocity becomes infinite for fixed x, and also for fixed t.

2. Model real shocks or model how to avoid them?

We have modeled a situation in which all cars traveling along the road never change their speed. If the function V(x) is increasing with x, then the cars in the front are faster that the ones behind, and there is no problem. But if some slower cars are ahead of the faster ones, there will be a crash if no one modifies their speed. There are two ways to develop the model further. Coming back to our initial motivation, the supersonic shock wave, we can refine our equation to include collisions, a one-dimensional version of shocks. The underlying cause is the same for supersonic planes, namely the air in front of the plane cannot travel fast enough to escape the plane coming up from behind! Or we can refine our equation to allow for more intelligent driver behavior, having the faster car slow down and slower car speed up.

Smarter drivers

When a collision is imminent, the behavior of a driver with a brain would surely be influenced by his neighbors, tending to compromise his speed with that of his neighbors. The analogy is to think of a viscous fluid: if you have observed honey flowing, you may have remarked that it behaves as a whole, as if all parts were tied together: no drops are splashing around, they stay linked to the main stream. Water has viscosity about 10,000 times less and air 50 less than that — but not zero. So one possible correction to Burger’s equation which avoids collisions when the initial function speed V_0 is bounded will involve adding viscosity term. Such a term is just a multiple of the second derivative \frac{\partial^2 V}{\partial x^2}. Hence, a modified equation takes the form which is known as the viscous Burger’s equation:

(2)   \begin{equation*} \frac{\partial V}{\partial t}+ \frac{\partial V}{\partial x}\cdot V = \nu\frac{\partial^2 V}{\partial x^2},\end{equation*}

with \nu>0 being called the viscosity. The modified world lines can be seen in Figure 4 for the same function V(x) describing the initial speed as in Figure 2.

Figure 4: In the same situation as Figure 2, we add some viscosity. Note how the cars bunch up but do not collide any more. On the right, we see the graph of V(x,t).

Why a second derivative? Recall that the second derivative is the limit of

(3)   \begin{equation*} \frac{V(x+\delta x,t) - 2V(x,t) + V(x-\delta x,t)}{\delta x ^2}.\end{equation*}

You can see this because it is the first derivative of the first derivative so will be close to

    \begin{align*}\begin{split} & \frac{V'(x+\delta x/2, t)-V'(x-\delta x/2, t)}{\delta x} \\ & \approx \frac{1}{\delta x} \left(\frac{V(x+\delta x,t)-V(x,t)}{\delta x} - \frac{V(x,t)-V(x-\delta x,t)}{\delta x}\right). \end{split}\end{align*}

Let us rewrite this expression as

    \[\frac{2}{(\delta x)^2}\,\frac{\left(V(x+\delta x,t)-V(x,t)\right) +\left(V(x-\delta x,t)-V(x,t)\right)}{2}.\]

What means the second factor? It is a mean of two terms: the first is the difference of speed in the front and in the middle and the second is the difference of speed in the back and in the middle. The danger of a collision is when V(x-\delta x,t)>V(x,t)>V(x+\delta x,t). In that case, the first term is negative and the second is positive. If the first term wins, the mean is negative and expression (3) tells you that you should break, since the left hand side is the acceleration. If the second term wins, then the mean is positive and you should accelerate.
Just before a collision happens, there will be a discontinuity in velocity with a faster car about to hit a slower one.

Adding the viscosity removes the shocks in Burger’s equation! We can now move to tell you about the million dollar prizes.
Note that for x_0 fixed, the left hand side represents the acceleration, since it is the derivative of the function V(x(t,x_0), t) with respect to t. So the modified viscous equation says that if the middle car is slower than the average speed of its neighbors ahead and behind, then it accelerates. If it is faster, then it decelerates.

3. The challenges of modeling the real world flows

The real world flows being modeled
We have been cheating a bit to make nice clean math out of a definitely messier real world. In particular, we have neglected bringing in the density of the cars on the road. Better traffic models like the Lighthill-Whitham-Richards model deal with local average density D(x,t) and as well as local average speed V(x,t). But versions of Burger’s equation still dominate these refined models so our discussion was basically ok.

The motion of air around a plane, especially around its wings, and the motion of water around a boat are real world applications of much greater interest. Mathematicians have worked on modeling the latter ever since Euler in the 18th century. Air and water are fundamentally different because water is almost completely incompressible while air is highly compressible. The incompressibility of oil is used for example to transmit tremendous forces in machines like bulldozers and backhoes. The million dollar prize concerns the incompressible case in three dimensions.

The million dollar prizes

In order to celebrate mathematics in the new millennium, Landon T. Clay established seven Prize Problems. The Prizes were conceived to record some of the most difficult problems with which mathematicians were grappling at the turn of the second millennium as well as to elevate in the consciousness of the general public the fact that in mathematics, the frontier is still open and abounds in important unsolved problems. A prize of 1 million dollars will be given for the solution of each problem — and one of them (the Poincaré conjecture) has been solved so far!

The 1 million dollar prize which interests us is concerned with the equations governing the flow of an incompressible fluid with viscosity in 3-dimensional space: these equations are called the Navier-Stokes equations, named after Claude-Louis Navier (1785–1836) and George Gabriel Stokes, (1819–1903) who worked on the motion of fluids. Except for the incompressibility which we will explain below, they are simply the three dimensional version of Burger’s equation with viscosity.

What’s the scoop? Most mathematicians conjecture that these Navier-Stokes equations never form shocks or other singularities like turbulence with eddies of all sizes, big to infinitesimal. Thus the equations should be solvable for time progressing as far as you want. If you prove that, the check is yours! But in science we are never sure of the conclusion before it is proved. Hence, if you can exhibit an example of initial conditions creating a shock or other singularity in finite time, the check is also yours. Why do mathematicians believe that no shock or other singularity will develop in Navier-Stokes equation? This was proven to be the case for 2-dimensional flows decades ago by a Russian woman, Olga Ladyzhenskaya. Since then, numerical simulations have not suggested that any singularities develop in the 3-dimensional case either. The simulation of solutions of partial differential equations like these is a huge computational challenge requiring sophisticated algorithms and programs running for long periods on very powerful computers.

Is that a contradiction with the shock wave of Figure 1? No! Air is compressible and the million dollar prize concerns incompressible flows.

The rest of this vignette is an enrichment for you, in two directions. In the first one, we explain how we allow singular solutions for Burger’s equation. (see section 5. READ MORE 1) If you have still remembered your multivariate calculus course, you can read how we modify Burger’s equation to model the movements of incompressible fluids in 3-dimensional space (see section 6. READ MORE 2).

4. Conclusion

What have we learnt? Many things. For instance, you could have imagined that modeling a fluid was so complex that you could not understand it. Now you have realized that it only involves the same basic principles of mechanics that you may have seen when studying forces, Newton laws, and conservation of energy. Also, the tools to express these physical laws simply involve the notions introduced in a multivariate calculus course. The form of the equations governing the motion of a fluid are very simple. Yet, there are many open questions about their solutions, even if they have been intensively studied, both theoretically and numerically, for more than 200 years.

If you have read the modeling of shocks, you have seen an important facet of the work of the mathematician. When solutions of a problem do not exist, the mathematician can create an object which is a solution. You are familiar with the creation of the complex numbers so as to add to \mathbb{R} the roots of any polynomial with real coefficients. Here, the mathematicians have created the distributions, thus allowing the partial differential equations (PDE) to have weak solutions.

But this has not sufficed to settle the million dollar prize described above. Let us quote Charles L. Fefferman describing the prize on the Clay website: “…, our understanding is at a very primitive level. Standard methods from PDE appear inadequate to settle the problem. Instead, we probably need some deep, new ideas.” Maybe one of your students will make the breakthrough some day…

5. READ MORE 1 – Going back to Burger’s equation

Allowing singular solutions
But shocks do exist in practice! As a scientist, we cannot just stop our analysis here. A “singularity” is a very rich phenomenon and the mathematician has no choice but to put on his or her mathematical glasses and think more deeply about this behavior in the context of the equation. This approach is called introducing generalized or `weak’ solutions that are not smooth. In our case, cars with distinct speeds are colliding. Think of this as a huge pile-up of ever more cars speeding in from the left with the pile-up moving fast enough to collide with ever more cars on the right. Thus we want to allow V(x,t) be discontinuous in x, having one limit on the left and another from the right.

But the equation uses derivatives of V. To deal with these, we need an idea which goes back to the 19th century and introduce a science-fiction kind of function called a delta function.

Figure 5: Heaviside and delta functions.

A short digression on \delta functions.
Consider the Heaviside function

    \[H(x)= \begin{cases} 0,&x<a,\\ b, &x=a,\\ 1, &x>a, \end{cases}\]

with graph in Figure 5.
Its derivative is the delta function at the point a. Note that this derivative is given by

    \[H'(x)= \delta(x)= \begin{cases} 0,&x<a,\\ +\infty, &x=a,\\ 0, &x>a, \end{cases}\]

so it is not a function in the general sense, but only a “generalized function”.
Yet, if the fundamental theorem of calculus were to remain true, we should have

    \[\int_{-\infty}^{+\infty} \delta(x)\;dx= \int_{-\infty}^{+\infty} H'(x)\;dx = H(x)\Bigr|_{-\infty}^{+\infty} = 1-0 =1,\]

i.e. the delta function is infinite exactly enough so that the integral \int_{\mathbb R} \delta(x) dx is equal to 1, that is the area under the “graph” of \delta is 1! This trick was used widely by the physicists in the 19th century and provided correct answers, thus challenging the mathematicians in search of rigor. It could finally be made rigorous in the 1940s by Laurent Schwartz, who called these new beasts, distributions.
To make the thing work, he said distributions aren’t functions f(x) with values at points
x, but they are “things” determined instead by their average values, which you think of as the integral of the distribution times a weighting function. A real life example is given by images made by cameras: you can measure light hitting a sensor but the sensor has a size; you can’t measure the light hitting an infinitesimal point.

Back to describing our shock wave, referring to Figure 2, we must let V be discontinuous across the bold slanting line at the top right. Burger’s equation will still hold if you use delta functions and make the shock itself move at the average of the speeds of the cars colliding from the left and the right. In fact, look near a point (x,t) on the locus of shocks. Let s be the speed of the shock wave itself so that the shock wave is tangent to a line x=st+a. The equation of this line can also be written as x-st=a. Hence parallel lines are given by x-st= a' with a'\in \mathbb R. On a small scale we can suppose that the speed is (approximately) constant on each side of the shock, and that it takes the value V_L on the left of the shock, and V_R on the right. So we need to show that

(4)   \begin{equation*} s=\frac{V_L+V_R}{2}.\end{equation*}

Since V is constant on each side of the shock, then V is constant along the lines x-st= a' with a'\in \mathbb R, and hence, that V(x,t) depends only on the value of x-st. If we let y=x-st, this allows to write V(x,t)=\widetilde{V}(y) for some (one variable!) function
\widetilde{V}. (The graph of this function could look like one of the back curves in Figure 3.) Hence, \widetilde{V} is of the form cH(y)+d, where H is a Heaviside function at a, and c=V_R-V_L.

Note that V\cdot \frac{\partial V}{\partial x} =\frac{\partial(\tfrac12 V^2)}{\partial x}, which allows to write Burger’s equation as

(5)   \begin{equation*} \frac{\partial V}{\partial t} + \frac{\partial(\tfrac12 V^2)}{\partial x} = 0. \end{equation*}

Since \tfrac12 V^2 has a jump of size \tfrac12 (V_R^2-V_L^2) = \tfrac12 (V_R-V_L)(V_R+V_L), the function \tfrac12 \widetilde{V} is of the form \tfrac12 (V_R^2-V_L^2)H(y)+ e. Hence,

    \[\frac{d\left(\tfrac12 \widetilde{V}^2\right)}{dy}=\tfrac12 (V_R^2-V_L^2)\delta(y).\]

Since \tfrac12 V^2(x,t) = \tfrac12 \widetilde{V}^2(y) then, by the chain rule,

(6)   \begin{equation*}\frac{\partial \left( \tfrac12V^2\right)}{\partial x} = \frac{d \left( \tfrac12\widetilde{V}^2\right)}{dy} \frac{\partial y}{\partial x}= \tfrac12 (V_R^2-V_L^2)\delta(y).\end{equation*}


(7)   \begin{equation*}\frac{\partial V}{\partial t}=\frac{d \widetilde{V}}{dy} \frac{\partial y}{\partial t}=(-s)\frac{d \widetilde{V}}{dy}= -s(V_R-V_L)\delta(y).\end{equation*}

Replacing (6) and (7) in (5) yields

    \[\left[\tfrac12 (V_R-V_L)(V_R+V_L) -s(V_R-V_L)\right]\delta(y)=0,\]

which can only be true if (4) is satisfied.

6. READ MORE 2 – Modeling the flow of a fluid in 3 dimensions
We want to model the flow of water, or more generally, the flow of any incompressible liquid. We already said that the water at each point x\in \mathbb R^3 is moving with speed given by a vector V(x,t). Already the three components of V are three unknown functions that we need to find. But actually we will need a fourth function, the pressure p(x,t). The pressure is that force that enables earth-moving equipment to do their work: a pump pressurizes a reservoir of oil which is linked by hoses to the pistons which move the buckets, shovels, etc. The pressure in a fluid comes from the force exerted on an element of a fluid from its neighbors. In general liquids like water and oil are incompressible, i.e. their volume does not change when compressed while gases are compressible. Now, to model the movement of a fluid, start with the same idea as traffic: assume the water tries to retain its velocity as time progresses but viscosity is causing the velocity at each point to average out with that of its neighbors. This can be expressed by equation (2) if we interpret its terms vectorially.

For the viscosity term we need the three dimensional way of adding a second derivative term. Now this term is supposed to say some function changes in the direction of being an average of its neighbors. On a road, you had two neighbors, one behind you and one in front. In three dimensions, thinking of this small cube around a point in the water, you have six neighbors, one to the left, one to the right, one in front, one behind, one above and one below. The effect of this is that instead of adding one term \tfrac{\partial^2 V}{\partial x^2}, you need to add three terms, namely:

    \[\triangle V := \frac{\partial^2 V}{\partial x_1^2} + \frac{\partial^2 V}{\partial x_2^2} + \frac{\partial^2 V}{\partial x_3^2}\]

an expression known as the laplacian of V.

Another problem is that the flow given by this equation may increase or decrease the density of the water, contradicting the fact that water is incompressible. What we need is that at any point, if we imagine a small cube around the point, the amount of water coming into the cube equals the amount of water going out. There’s a very simple way of ensuring this: the divergence of V must be zero at all times. If x=(x_1,x_2,x_3) and V=(V_1,V_2,V_3), then

    \[\text{div}(V) := \frac{\partial V_1}{\partial x_1} + \frac{\partial V_2}{\partial x_2} + \frac{\partial V_3}{\partial x_3} = 0.\]

If we think of the derivative as a limit, the vanishing of the divergence says:

    \begin{align*} \begin{split}& V_1(x_1+\delta,x_2,x_3)-V_1(x_1-\delta,x_2,x_3) \\ + &V_2(x_1,x_2+\delta,x_3)-V_2(x_1,x_2-\delta,x_3) \\ +& V_3(x_1,x_2,x_3+\delta)- V_3(x_1,x_2,x_3-\delta)\to 0,\end{split}\end{align*}

which amounts to adding up the water outflow and minus the water inflow at each of the six faces of the small cube and making it zero.

Finally to ensure that, under the flow, the divergence of V stays zero, we need to add an additional force, which is caused by the water molecules pushing against each other, refusing to compress. The resistance to compression is measured by the pressure and the gradient of the pressure is a force which accelerates each bit of the water so as to maintain constant density.

Euler’s and Navier-Stokes’ equations
We are now ready to put all this together in the equations of an incompressible flow. The equation has four unknown functions which must be solved for, the three components V and the pressure p. The pressure must be there to ensure incompressibility. The Navier-Stokes equation is the system of four equations:

    \[\boxed{\begin{aligned} \frac{\partial V}{\partial t} + V \cdot \frac{\partial V}{\partial x}&= -\nabla p + \nu\triangle V \\ \text{div}(V) &=0 \end{aligned}}\]

(Here \nabla p is just the gradient of p, namely the vector (\frac{\partial p}{\partial x_1}, \frac{\partial p}{\partial x_2}, \frac{\partial p}{\partial x_3}). The minus sign is justified by the fact that the left hand side represents the acceleration. If the pressure increases then x decelerates and we get a negative acceleration.) The special case where the viscosity \nu=0 is called Euler’s equation. In dimension 1, incompressibility implies that the speed is constant, and hence there are no shocks. It has been proved that there are also no shocks in dimension 2 and, given any initial speed function, there is a unique solution that exists for infinite time. In dimension 3, without the taming effect of viscosity, the likelihood that singularities develop appears much greater — but nobody knows for sure.

7. Bibliography

1. A Mathematical Introduction to Fluid Mechanics, Alexandre Chorin and Jerrold Marsden, Springer 1993.

2. Wikipedia article’_equation.

3. Nonlinear Partial Differential Equations for Scientists and Engineers, Lokenath Debnath, Birjauser-Boston, 2004.

This post is also available in: German, Arabic

This entry was posted in Mathematics Within the Last 100 Years. Bookmark the permalink.

Leave a Reply

Your email address will not be published. Required fields are marked *