flux limiters
 

The TVD limiter definition psi(r+) given in Hirsch is

r+ = (u_i+2 - u_i)/(u_i+1 - u_i).

what if the denominator is zero as in the initial condition of flow througha shock tube?

Re: flux limiters
 

just set r+=0.0 at the initial state.