I think am already getting some results for inviscid burgers. I still need some insights.

I used Cflcof=dt(Smax)/dx, where 0

Hi Chinedu,

The optimal CFL should be close to 1 as, in linear problems, Godunov's method is the exact shift operator at CFL =1.

Now, in the case of Burgers' equation, the scheme will have 1st order level of dissipation so the numerical solution may only resemble the exact solution. Did you try refining your mesh to see if the numerical solution apprximates the exact solution?

Good luck,

Marco

Yes, I finally got a good approximation to the exact solution.