To summarize, we resolved the earlier observed mismatch between the theoretical predictions, based on the Euler equations, for the time evolution of a blast wave and results from hard sphere simulations in two and three dimensions. In this paper, we compare the results of direct numerical simulations of the Navier Stokes equation in two and three dimensions with the MD data, and show they match quantitatively. These results are consistent with those found recently for blast waves in one dimension [1]}, [2]}. We conclude that the order in which the scaling limit is taken is important in the blast problem. In the TvNS solution, the scaling limit is taken first – heat conduction and viscosity terms being irrelevant – before the exact solution is found. However, the correct limit to take is to find the solution with non-zero heat conduction and viscosity terms and then take the scaling limit.
