In this paper, we consider large eddy simulations (LES) for human stenotic carotids in presence of atheromasic plaque, a pathological condition where transitional effects to turbulence may occur, with relevant clinical implications such as plaque rupture. We provide a reference numerical solution obtained at high resolution without any subgrid scale model, to be used to assess the accuracy of LES simulations. In the context we are considering, ie, hemodynamics, we cannot refer to a statistically homogeneous, isotropic, and stationary turbulent regime; hence, the classical Kolmogorov theory cannot be used. For this reason, a mesh size and a time step are deemed fine enough if they allow to capture all the features of the velocity field in the shear layers developed after the bifurcation. To assess these requirements, we consider a simplified model of the evolution of a 2D shear layer, a relevant process in the formation of transitional effects in our case. Then, we compare the results of LES sigma model (both static and dynamic) and mixed LES models (where also a similarity contribution is considered). In particular, we consider a realistic scenario of a human carotid, and we use the reference solution as gold standard. The results highlight the accuracy of the LES sigma models, especially for the static model.