Quantum Many-Body Simulations of the 2D Fermi-Hubbard Model in Ultracold Optical Lattices