N-body we can do, for example Illustris Project. To address the concern of this post, one would need non-linear general relativity simulation. That's more challenging.