2

I have been struggling on this problem for a while now so here it is:

I am looking to work out when the shockwaves occur in the traffic flow model given below:

A traffic flow governed by the equations below is observed at $t = 0$ to have a dip in the density distribution at $x = 0$ given by:

$$ \rho (x,t=0) = \left\{\begin{aligned} &a &&: x \lt 0\\ &a- \frac {a}{2}(1+x)&&: -1<x< 0\\ &a- \frac {a}{2}(1-x)&&: 0<x<1\\ &a &&: 1<x\\ \end{aligned} \right.$$

where $0<a<1$

Where the equation of characteristics is given by $x-(1-2 \rho)t=x_0$

I believe I have the correct sketch of this function however I cannot figure out how to find the equations and then the solution of the shocks.

Any explanation and help would be great.

juper
  • 267
  • 1
  • 10

1 Answers1

2

Here, we are dealing with the macroscopic traffic-flow model by Lighthill-Whitham-Richards (LWR), which consists in a scalar hyperbolic conservation law $$ \frac{\partial}{\partial t} \rho + \frac{\partial}{\partial x} Q(\rho) = 0 \, , $$ where the flux $Q(\rho)=\rho\, v(\rho)$ depends only on the density of cars $\rho$. The simplest expression for the car velocity of the form $v(\rho) = 1 - \rho$ was introduced by Greenshields. The car density is constant along the characteristics, which are given by $dx/dt = Q'(\rho) = 1 - 2\rho$. The characteristic curve issued from $x_0$ is a straight line in $x$-$t$ space, $x(t) = (1-2\rho_0(x_0))t + x_0$, on which the car density equals its initial value $$ \rho_0(x_0) = \left\lbrace \begin{aligned} &a &&\text{if}\quad x_0<-1\\ &a(1-x_0)/2 &&\text{if}\quad {-1}<x_0<0\\ &a(1+x_0)/2 &&\text{if}\quad 0<x_0<1\\ &a &&\text{if}\quad 1<x_0\, . \end{aligned}\right. $$ Let us assume that $0< a < 1/2$. The slope of the characteristics is larger in the slit* than outside. Since the slope is positive, the characteristics cross on the right of the slit: a shock-wave arises. On the left of the shock trajectory $(x(t),t)$, we have the data coming from the slit. From the equation of characteristics, we deduce $$ \rho_0(x_0) = -\frac{1}{2}\left( \frac{x(t) - x_0}{t} -1 \right) = \frac{a}{2} \left(1 \mp x_0\right) . $$ Solving with respect to $x_0$, one obtains $$ x_0 = \frac{x(t) - t + at}{ 1 \pm at} \quad\text{and}\quad \rho_0(x_0) = \frac{a}{2} \frac{1 \mp (x(t) - t)}{1 \pm at} \, , $$ depending whether $(x(t),t)$ has $x_0$ in $[-1,0]$ or in $[0,1]$. On the right, we have the data from outside the slit: $\rho_0(x_0) = a$. Therefore, the Rankine-Hugoniot condition writes $$ x'(t) = \frac{Q\!\left(\frac{a}{2} \frac{1 \mp (x(t) - t)}{1 \pm at}\right) - Q(a) }{\frac{a}{2} \frac{1 \mp (x(t) - t)}{1 \pm at} - a}\, , $$ which is a differential equation with initial condition $x(0) = 1$. Solving this equation gives the shock trajectory. The case $1/2<a<1$ is tackled in a similar manner.


*slit : area where ${-1}<x_0<1$.

EditPiAf
  • 20,898
  • this is an excellent explanation, thank you! Is there a way to calculate the shocks for a arbitrary value of $a$? Also ,what do you mean by the slit? – juper May 01 '17 at 22:52
  • I'm struggling to work out where the $\frac{a}{2} \frac{1 \mp (x(t) - t)}{1 \pm at}$ has come from. I understand that the shock occurs where the characteristics in the slit intersect other characteristics but where has this equation come from? Thank you for helping me by the way :) – juper May 01 '17 at 23:03
  • Thank you. Just a clarification, in my question I had the density in the slit to be $a- \frac {a}{2}(1+x)$ but in your $\rho$ you don't have the $a-$ part. Why is this? – juper May 02 '17 at 11:43
  • okay thank you, this has helped a lot! – juper May 02 '17 at 21:51