The thermal suppression of heavy quark bound states represents an ideal observable for determining if one has produced a quark-gluon plasma in ultrarelativistic heavy-ion collisions. In recent years, a paradigm shift has taken place in the theory of quarkonium suppression due to new first principles calculations of the thermal widths of these states. These thermal widths are large, e.g. O(20-100 MeV) for the Upsilon(1S), and cause in-medium suppression of the states at temperatures below their traditionally defined disassociation temperatures. In order to apply the newly developed understanding to phenomenology, however, one must make detailed 3+1d dissipative hydrodynamical models of the plasma including the effects of finite shear viscosity. These effects include not only the modification of the time evolution of the temperature of the system, flow, etc., but also non-equilibrium modifications of the heavy quark potential itself. In this proceedings contribution, we briefly review the setup for these model calculations and present comparisons of theory with data from RHIC 200 GeV/nucleon Au-Au collisions, LHC 2.76 TeV/nucleon Pb-Pb, and LHC 5.02 TeV/nucleon Pb-Pb collisions as a function of number of participants, rapidity, and transverse momentum.

4 pages, 4 figures; Proceedings contribution to Quark Matter 2017, Feb 6-11 2017; v2 - some typos fixed and references updated