Numerical Study Of Hubbard Model For Iron-Based Superconductors