It's just hard to figure out what the functions are for a set of boundary conditions.
F = mx''(t) = mx''(-t) since d/dt x(-t) = -x'(-t), and d/dt (-x'(-t)) = x''(-t)
Navier-Stokes is only time-reversible if you ignore viscosity, because viscosity is velocity-dependent and you can already see signs of that being a problem in the derivation above (velocity pops out a minus sign under time reversal). From my reading the OP managed to derive viscous flow too, so there really is a break in time-symmetry happening somewhere.
I wonder if "t -> -t" is lost in the Boltzmann step or in the hydrodynamic step.
Besides that, I don't think anybody is really arguing that the correlations are actually lost after a collision, just that it's usually a good approximation to treat them as if they are.
Note that the solutions x(t) are not generally time symmetric. We aren't saying that x(t)=x(-t), we are saying that x(t) is a solution to the differential equation if and only if x(-t) is, which is a weaker statement.
It's hard to take full reversibility seriously given Newton's equations are not actually deterministic. If they're not deterministic, then they can't be fully reversible.
Of course maybe these non-deterministic regimes don't actually happen in realistic scenarios (like Norton's Dome), but maybe this is hinting at the fact that we need a better formalism for talking about these questions, and maybe that formalism will not be reversible in a specific, important way.