The trouble with Hubble: Local versus global expansion rates in
  inhomogeneous cosmological simulations with numerical relativity
  In a fully inhomogeneous, anisotropic cosmological simulation performed by
solving Einstein's equations with numerical relativity, we find a local
measurement of the effective Hubble parameter differs by less than 1\% compared
to the global value. This variance is consistent with predictions from
Newtonian gravity. We analyse the averaged local expansion rate on scales
comparable to Type 1a supernova surveys, and find that local variance cannot
resolve the tension between the \citet{riess2018b} and \citet{planck2018a}
measurements.
