5

I have a question about the following Burgers' equation.

$u_t + (\frac12u^2)_x = 0 $ with $u(x,0) = \sin(x)$ on $[0,2\pi]$ and periodic boundary conditions.

When I studied this equation numerically, I notice that once the shock forms, the shock stays still in the same place, yet the magnitude of the function decreases as time passes. That is, the initial function $\sin(x)$ ranges from $-1$ to $1$, but as time passes, the magnitude of the maximum and minimum of the function decays.

Is this what is supposed to happen analytically?

EditPiAf
  • 20,898
mononono
  • 2,028

2 Answers2

5

The method of characteristics gives $u=f(x-ut)$ where $f(x) = \sin(x)$ represents the initial data $u(x,0)$. This implicit equation can be used to approach the solution, until a shock wave is met by the characteristics. Here is a sketch of the characteristics in the $x$-$t$ plane:

characteristic curves

As shown in this post, the shock is formed when characteristics intersect for the first time, i.e. at the breaking position and time \begin{aligned} x_B &= \arg \min f'(x) = (2n+1)\pi , \qquad (n\in \Bbb{Z})\\ t_B &= \frac{-1}{\min f'(x)} =1 \end{aligned} as shown in the figure. Setting $v=-u$, one notes the relationship $$ v(x_B - x, t) = \sin \left[x - v(x_B - x, t)\, t\right] $$ which proves the symmetry property $u(x_B - x, t) = -u(x, t)$.

By solving the implicit equation $u = f(x-ut)$, the method of characteristics yields the values $u_L$, $u_R$ of $u$ on the left and the right of the shock. The position $(x_s(t), t)$ of the shock wave is given by the Rankine-Hugoniot condition $$x'_s(t) = \tfrac{1}{2} (u_L+u_R)$$ with initial position $x_s(t_B) = x_B$. Therefore, using symmetry $u_L = -u_R$, the shock wave remains at the same place (static shock). Its magnitude $u_L$ increases and then vanishes in time, as shown in the figure below, where $u_L$ --obtained numerically-- is plotted with respect to $t$:

shock magnitude

EditPiAf
  • 20,898
  • Sorry if this is an old post, but could you describe how you calculated the Rankine-Hugoniot condition? I am confused on this as in the examples I have seen it is clear what $u_L$ and $u_R$ is, but less so here as the shock occurs within a single function. Here $f(u) = sin(x)$, but what do we use for $u_L$ and $u_R$? – CBBAM Jan 04 '22 at 22:13
  • 1
    @CBBAM Fact is that we can't find $u_L$ and $u_R$ analytically. Nevertheless, knowing the symmetry $u_L = -u_R$ is enough to deduce the shock trajectory. – EditPiAf Jan 05 '22 at 09:12
  • Thanks, that makes sense as to why I was beating my head over this for so long with no solution. – CBBAM Jan 05 '22 at 17:32
1

No, analytically that is not what happens. The equation in your question is discussed in Miller's Applied Asymptotic Analysis, and I quote from the middle of p. 78:

Once shocks form, there is no longer a global unique solution of the nonlinear initial-value problem.

It turns out that it's possible to modify the equation in such a way that the shocks are stable (though still move as time progresses, if I recall correctly) by adding a diffusion term. I definitely recommend reading section 3.6 of Miller's book.

  • 1
    @mononono Note that this chapter of Miller is about the method of characteristics (with implicit solution $u=\sin(x−ut)$), which is known to work only until a shock wave is encountered (see my answer for details). – EditPiAf Dec 05 '17 at 13:14
  • Is there an explicit formula of the analytic solution before shock is formed? I am trying to use this example to test the order of convergence of my numerical method and I need some reference. Thanks! – Ling Nov 19 '19 at 19:27