upwind for micronozzle flow
 

Good morning all, I'm running some simulations for micronozzle flow, and I can't understand which numerical method is best fitting: if I use a second order upwind scheme, I can see a weak shock in the divergent part of the nozzle, but if I use a first order upwind this shock becames weaker or disappear. Which solution should I consider correct? Thanks!