Numerical-relativity simulations of long-lived remnants of binary neutron star mergers