Wall shear stress (WSS) and pressure are important factors in the development of cerebral aneurysms. We aimed to develop a computational fluid dynamics simulator for flow in the complete circle of Willis to study the impact of variations in vessel radii and bifurcation angles on WSS and pressure on vessel walls. Methods-: Blood flow was modeled with Navier-Stokes equations as an incompressible newtonian fluid within rigid vessel walls. A model of the circle of Willis geometry was approximated as a network of tubes around cubic curves. Pulsatile inlet flow rates and constant outlet pressure were used as boundary conditions. Results-: The simulations confirmed that differences in vessel radii and asymmetric branch angles influence WSS magnitude and spatial distribution. High WSS occurred at locations where aneurysms are frequent and in anatomic variants known to be associated with an increased risk for aneurysm development. Conclusions-: Computational fluid dynamics analysis can be applied to the complete circle of Willis and should be used to study the pathophysiology of this complex vascular structure, including risk factors for aneurysm development. Further development of the method should include simulations with flexible vessel walls.