2

I am trying to plot the characteristic lines for the inviscid Burgers equation which is $$u_t +uu_x=0.$$

From what I understand, with the initial condition $u(x,0)=f(x)$ and using the method of characteristics, the solution can be obtained implicitly as

$$u(x,t)=f(\xi)$$

where $\xi=x-ut$ and $\frac{dx}{dt}=f(\xi)$ describes the characteristics curve.

Then to plot the characteristics onto the $(x,t)$ plane, one could rearrange to get

$$t=\frac{1}{u}(x-\xi)=\frac{1}{f(\xi)}(x-\xi) $$

Assuming the above is correct, then consider the following initial condition:

\begin{equation} u(x,0)=f(x)= \left \{ \begin{aligned} &a^2-x^2, && \text{if}\ |x|)\leq a \\ &0, && \text{if}\ |x|)\geq a \\ \end{aligned} \right. \end{equation}

which gives the solution \begin{equation} u(x,t)=f(\xi)= \left \{ \begin{aligned} &a^2-\xi^2, && \text{if}\ |x|)\leq a \\ &0, && \text{if}\ |x|)\geq a \\ \end{aligned} \right. \end{equation}

where $\xi=x-ut=x-(a^2-\xi^2)t.$

The characteristics equation would then be $$t=\frac{1}{a^2-\xi^2}(x-\xi)$$

and substituting $\xi$

$$t=\frac{1}{a^2-(x-ut)^2}(ut) \tag{1}$$

which shouldn't be correct since we are expecting the characteristics to be straight lines.

Furthermore, I think I have started to enter into a "loop" with the substitutions. May I ask what have I done wrong?

An example of the characteristic lines onto the $x-t$ plane can be seen in Inviscid Burgers equation


This example was obtained from Debnath, Example 5.2.1.

EditPiAf
  • 20,898
mathnoob
  • 115

1 Answers1

0

Indeed, we have $u=f(x-ut)$. It follows from the method of characteristics which gives the set of lines $$ x(t) = x_0 + f(x_0) t , \quad\quad x_0\in \Bbb R , $$ where the function $f = u(\cdot ,0)$ represents the initial data $x \mapsto (a^2-{x}^2)\,\Bbb I_{|x|\leq a}$ and $\Bbb I$ is the indicator function (see e.g. this post for the case $a=1$). Along these curves, $u(t) = f(x_0)$ is constant and equal to its initial value. The equation for $x(t)$ above can be rewritten as $x_0 = x - u t$, and thus, we get $u=f(x-ut)$. However, starting from this equation to recover the equation of the characteristic lines may not be a good idea.

To solve for $u(x,t)$, we must distinguish two cases here:

  • The first case is $|x-ut|\geq a$ for which $u=0$. Therefore, if $|x|\geq a$, we have $u=0$.

  • The second case is $|x-ut|\leq a$ for which $u=a^2 - (x-ut)^2$. Solving the quadratic equation, we find the expression of $ u $ for $|x|\leq a$.

Finally, $$ u(x,t) = \frac{2 t x - 1 + \sqrt{1 - 4 t x + 4 a^2 t^2} }{2 t^2} \Bbb I_{|x|\leq a} . $$ Note that characteristics intersect at the breaking time $t_b = 1/(2a)$, where the classical solution becomes multi-valued: a shock wave is generated.

EditPiAf
  • 20,898
  • May I know what is the exact function used to plot the characteristic curves? That is my question which I did not make clear. So I am interested in $x=x_0+f(x_0)t$; what is $f(x_0)$? Is it $a^2-x_{0}^{2}$ ? Also, do you know why my attempt in finding the equation responsible for the characteristic lines was wrong? – mathnoob Aug 26 '19 at 20:35
  • @mathnoob answer edited. Yes this is what you want to draw. The problem with your approach is that it was similar to "watching a movie from the end to the beginning"! – EditPiAf Aug 27 '19 at 06:58