This has been up as a featured post for a while and nobody else has responded, so let me add something else. The following quote is from the paper I posted in Post #1,
"Thus the physical meaning of these calculations is simple: according to our approximation, it is not the average but the typical energy density that governs the expansion rate of the Universe. At high redshifts, where the distribution is fairly symmetric, the typical value of (mode of the PDF) is close to the average and the Universe evolves without backreaction. At late times skewness increases, the volume of the Universe is dominated by voids, and the typical value of is negative, thus effectively M < 1. High density regions, where metric perturbations are perhaps the largest, are inconsequential to this effect: what matters is the non-Gaussianity of the density distribution, in particular, the large volume fraction of low density regions, as advertised earlier."
This makes perfect sense to me. Early on, when the universe is of uniform density, the Friedmann equation works fine, but at late times the density becomes highly non-uniform. In today's universe, most of the volume is in regions of below average density. These regions expand faster because of their lower density, and so the apparent expansion of the universe accelerates. Can anyone offer a rebuttal and explain to me why this picture doesn't make sense?