  The spectrum of the Upsilon system is investigated using the Nonrelativistic
Lattice QCD approach to heavy quarks and ignoring light quark vacuum
polarization. We find good agreement with experiment for the Upsilon(1S),
Upsilon(2S), Upsilon(3S) and for the center of mass and fine structure of the
chi_b states. The lattice calculations predict b-bbar D-states with center of
mass at (10.20 +/- 0.07 +/- 0.03)GeV. Fitting procedures aimed at extracting
both ground and excited state energies are developed. We calculate a
nonperturbative dispersion mass for the Upsilon(1S) and compare with
tadpole-improved lattice perturbation theory.
