We investigate the flow dynamics and oscillatory behaviour of wall shear stress (WSS) vectors in intracranial aneurysms using high resolution numerical simulations. We analyse three representative patient-specific internal carotid arteries laden with aneurysms of different characteristics: (i) a wide-necked saccular aneurysm, (ii) a narrower-necked saccular aneurysm, and (iii) a case with two adjacent saccular aneurysms. Our simulations show that the pulsatile flow in aneurysms can be subject to a hydrodynamic instability during the decelerating systolic phase resulting in a high-frequency oscillation in the range of 20-50 Hz, even when the blood flow rate in the parent vessel is as low as 150 and 250 ml min 21 for cases (iii) and (i), respectively. The flow returns to its original laminar pulsatile state near the end of diastole. When the aneurysmal flow becomes unstable, both the magnitude and the directions of WSS vectors fluctuate at the aforementioned high frequencies. In particular, the WSS vectors around the flow impingement region exhibit significant spatio-temporal changes in direction as well as in magnitude.