shock detection
 

hi!
i m making a code for higher order PVU scheme for 1-D Riemann problems, i m getting overshoots and undershoots in the test case of u_L=-19.59745, u_R=-19.59745, rho_L=rho_R=1.0, P_L=1000, P_R=0.01. interface=0.8, te=0.012, pliz let me know how to remove these overshoots and undershoots. i put condition to detect the shock as w.ge.0.5.AND.(u(i)-c(i).ge.u(i+1)-c(i+1).and.p(i).gt.p(i+1).and.rho(i).gt.rho(i+1))
where w is weight or gradient
thanks

Dear Khan,

You can achieve it using slope limiters, start with a max-min limiter for example. You can find lot of literature on the same.

Regards,

Ganesh