Pulsatile flow in a 3D model of arterial double stenoses is investigated using a large eddy simulation (LES) technique. The computational domain that has been chosen is a simple channel with a biological-type stenosis formed eccentrically on the top wall. The pulsation was generated at the inlet using the first four harmonics of the Fourier series of the pressure pulse. The flow Reynolds numbers, which are typically suitable for a large human artery, are chosen in the present work. In LES, a top-hat spatial grid-filter is applied to the Navier-Stokes equations of motion to separate the large-scale flows from the sub-grid scale (SGS). The large-scale flows are then resolved fully while the unresolved SGS motions are modelled using a localized dynamic model. It is found that the narrowing of the channel causes the pulsatile flow to undergo a transition to a turbulent condition in the downstream region; as a consequence, a severe level of turbulent fluctuations is achieved in these zones. Transitions to turbulent of the pulsatile flow in the post stenosis are examined through the various numerical results, such as velocity, streamlines, wall pressure, shear stresses and root mean square turbulent fluctuations.