[ayushmorx]

Hey Guys!

I have made a code to solve a very simple Burger's equation problem, using TVD 3 stage, upwind scheme. Wave velocity a=1. Since upwind is a perfect shift for this equation, I should be retaining the flow structure as it changes. Instead, it is dispersing, eventually becomes a sine wave. I have tried the single stage, and it works perfectly, so I am pretty sure there are some problems with my equation for TVD. I have referred quite a few books and have not been able to figure out the cause of the dispersion yet. I would really, really appreciate any help in this matter.

BTW ntime is the time in seconds and I am trying to solve the flow for 100 mesh size. Increasing the mesh size does retain the structure, but a change in temporal scheme should not significantly affect the structure of the wave according to me.

What do you think?