Black-hole universe: Time evolution

Time evolution of a black hole lattice toy model universe is simulated. The vacuum Einstein equations in a cubic box with a black hole at the origin are numerically solved with periodic boundary conditions on all pairs of faces opposite to each other. Defining effective scale factors by using the area of a surface and the length of an edge of the cubic box, we compare them with that in the Einstein-de Sitter universe. It is found that the behavior of the effective scale factors is well approximated by that in the Einstein-de Sitter universe. In our model, if the box size is sufficiently larger than the horizon radius, local inhomogeneities do not significantly affect the global expansion law of the Universe even though the inhomogeneity is extremely nonlinear.