Drift of the 1D radiative shocks

I looked at the drift of the shock front for one computational time for all 18 models. One computational time is different for each model since the length scale is different for each model, so in order to compare apples to apples I decided that vdrift / vshock was the appropriate diagnostic.

The drift was not measurable for vs = 60 km/s. Even in the worst case, the drift velocity was only 0.065% of the shock velocity. Therefore, I think this effect is negligible. Obviously, if we ran the simulations for many computational times, the shock would eventually drift off the grid but that kind of simulation will take too long and will never be necessary. Still, there must be a very small error somewhere. I suspect it is somewhere in the chemistry or scaling, but again, the error is so small that we don't need to worry about it.

Attachments (1)

Download all attachments as: .zip

Comments

No comments.