Simulating the Hubbard Model with Equivariant Normalizing Flows