Finding the time for the first shock for a quasilinear first order PDE

In summary, the method of setting x_\xi = 0 can be used to determine the position of the first shock in PDEs of the form u_t + g(u)u_x = f(u), with initial condition u(x,0) = h(x). This method works by finding the intersections of neighbouring characteristics, which can be used to determine the first shock. However, in more complicated cases where x_\xi does not exist at a certain point, this method may not accurately determine the position of the first shock.
  • #1
BloonAinte
19
1
TL;DR Summary
To find a shock wave, do we always solve the equation ##x_{\xi}=0##? The PDEs I consider are of the form ##u_t + g(u) u_x = f(u)##, with initial condition ##u(x,0) = h(x)##.
To find a shock wave, do we always solve the equation ##x_{\xi}=0##? The PDEs I consider are of the form ##u_t + g(u) u_x = f(u)##, with initial condition ##u(x,0) = h(x)##. I have been looking at the solutions for problems in my homework sheet but this method was used with no explanation.

Why does this method work to find the formation of the first shock?

Thank you very much!
 
Physics news on Phys.org
  • #2
A shock forms when neighbouring charcteristics first intersect. Characteristics corresponding to neighbouring values of [itex]\xi[/itex] will intersect when [tex](x(\xi,t),t) = (x(\xi + \delta \xi, t), t)[/tex] which to first order in [itex]\delta \xi[/itex] requires [tex]\frac{\partial x}{\partial \xi} = 0.[/tex]
 
  • Like
Likes BloonAinte
  • #3
Thank you! Does this method always determine the first shock?

I have looked up an example to illustrate my doubt:

https://math.stackexchange.com/ques...-curve-and-sketch-characteristics-in-xt-plane

That question considers some PDE and obtains characteristics as follows:
IOfH2.png

In this question, the characteristics for ##x > 0.5## intersect those for ##x < -0.5##. These are not neighbouring characteristics, but they have intersections. However, ##t## is not minimal there, so it's not the position of the first shock. These also do not seem to solve ##x_{\xi} = 0##.

From this "investigation", I guess that this means that this method of setting ##x_{\xi}=0## will always give the position of the first shock, and this always occurs due to characteristics which are close together.

Does this sound accurate? Thank you!
 
  • #4
This example is complicated by the fact that [itex]x_\xi[/itex] does not exist at [itex]\xi = 0[/itex]. In this case, a characteristic in [itex]\xi > 0[/itex] intersects the characteristic [itex]x = 0[/itex] when [tex]
t = \frac{2\xi}{2 - \xi}.[/tex] This first happens at [tex]
\max \left\{ 0, \inf_{\xi \geq 0} \frac{2\xi}{2 - \xi} \right\} = 0.[/tex]
 
  • Like
Likes BloonAinte
  • #5
Thank you so much for all your help! :) I understand this more now ^^
 
  • Like
Likes jim mcnamara

Similar threads

Back
Top