The non-linear evolution and spatial propagation of the thermal instability in accretion disks in galactic nuclei are investigated. For relevant black hole masses, M, and mass accretion rates, M0, parts of the disk are thermally unstable. The behavior of the unstable disks depends on the detailed prescription for the viscosity parameter a. For constant a, the thermal instability cannot propagate through a large radial extent, so the disk structure is only modified locally. We call this a "purr" type disk. On the other hand, if we prescribe α ∝ (h/r)n with n = 1 or 2, where h is the scale height of the disk, more extended regions in the disk suffer the thermal instability. This leads to dramatic outbursts of rapid accretion onto the hole, separated by longer intervals of quiescence. We call this a "roar" type disk. The rise time of a roar-type outburst is short, but the decay occurs more gradually on a viscous time scale. The repetition period is equation presented The duty cycles are short, τon/τoff ≤ 10-1. This allows high-luminosity quasars to be powered by a modest accretion source. The temperatures of the roar type disks are consistent with the blue bump observed in AGN spectra. Quiescent disks should be present in many galactic nuclei and should emit a characteristic thermal continuum at ∼ 103 K.