This paper investigates intimal growth in arteries, induced by hemodynamical shear stress, through finite element simulation using the FEniCS computational environment. In our model, the growth of the intima depends on cross section geometry and shear stress. In this work, the arterial wall is modeled as three distinct layers: the intima, the media, and the adventitia, each with different mechanical properties. We assume that the cross section of the vessel does not change in the axial direction. We further assume that the blood flow is steady, non-turbulent, and unidirectional. Blood flow induces shear stress on the endothelium and stimulates the release of Platelet Derived Growth Factor (PDGF) which drives the growth. We simulate intimal growth for three distinct arterial cross section geometries. We show that the qualitative nature of intimal thickening varies depending on arterial geometry. For cross section geometries that are annular, the growth of the intima is uniform in the angular direction, and the endothelium stays circular as the intima grows. For non-annular cross section geometries, thicker intimas grow more compared to thinner ones, shear stress and intimal thickening are negatively correlated with the distance from the flow center, where the flow velocity is maximal. Over time, the maxima and minima of the curvature increase and decrease respectively, the PDGF concentration increases, and the lumen becomes more polygonal. The model provides a framework for coupling hemodynamics simulations to mathematical descriptions of atherosclerosis, both of which have been modeled separately in great detail.