Ongoing issue with Pressure-Energy based schemes
This doesn't have very much information in it as it's just meant to track for me and @matthieu
We have just attempted to sort out some things in the drift, by drifting the smooth pressure correctly. Unfortunately this doesn't seem to have fixed things at all, but the code is still running.
We have identified the following three tests:
-
The Zeldovich with no viscosity, as we know the AV parameters are not the cause of the issue. Note we already 'pass' the comoving sodshock. This should help to see whether the cosmology factors in the hydro are correct. -
The Sedov blast with forced multi-dt; i.e. we force particles on the 'left' side of the box to have a time-step of some factor lower than those on the right. We can then look at the solutions in the four quadrants. I will probably do this in 2D. - Try fixing/re-implementing the Pressure-Entropy schemes to help see if it is a problem with the smooth pressure.