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

  • #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