OATAO is an open access repository that collects the work of Toulouse researchers and makes it freely available over the web where possible. This is an author-deposited version published in: http://oatao. For this purpose, we perform the first simulations of blood flow in an anatomically accurate large human intra-cortical vascular network (~10000 segments), using a 1D non-linear model taking account of the complex rheological properties of blood flow in microcirculation. This model predicts blood pressure, blood flow and hematocrit distributions, as well as volumes of functional vascular territories, and regional flow at voxel and network scales. First, the influence of the prescribed boundary conditions (BCs) on the baseline flow structure is investigated, highlighting relevant lower-and upper-bound BCs. Independent of these BCs, large heterogeneities of baseline flow from vessel to vessel and from voxel to voxel, are demonstrated. These heterogeneities are controlled by the architecture of the intra-cortical vascular network. In particular, a correlation between the blood flow and the proportion of vascular volume occupied by arterioles or venules, at voxel scale, is highlighted. Then, the extent of venous contamination downstream to the sites of neuronal activation is investigated, demonstrating a linear relationship between the catchment surface of the activated area and the diameter of the intra-cortical draining vein.