It is well known that rotating inviscid accretion flows with adequate injection parameters around black holes could form shock waves close to the black holes, after the flow passes through the outer sonic point and can be virtually stopped by the centrifugal force. We numerically examine such shock waves in 2D accretion flows with 10−5 to 106 Eddington critical accretion rates around a supermassive black hole with 106M⊙. As the results, the luminosities show QPO phenomena with modulations of a factor 2–3 and with quasi-periods of a few to several hours.