On the numerical modeling and optimization of a bottom-referenced heave-buoy array of wave energy converters