Navier Stokes tvd

shuo

I want to solve the Navier-Stokes equations using upwind and symmetric TVD schemes(Yee, Klopper.. 1990). The flux F have to be modified using the R whose columns are eigenvectors of dF/dU.

F_mod = 1/2*(F_i + F_i+1 + R*phi)

I have R for the inviscid F but not the viscous version. What I did instead was to only modify the inviscid fluxes and leave the viscous fluxes as is. The stresses themselves were calculated using central differencing.

I compared simulation of an supersonic impinging jet using EUler and N-S. The flow contours using the latter had 'strips' emanating from the jet where as the Eulers eq produced a smoother flowfield. Also the simulation using N-S blows up just before impingement. Any suggestions?

Cheers Shuo