Simulations of the equatorial thermosphere anomaly: Field-aligned ion drag effect