Summary: | The Large Eddy Simulation (LES) technique is used to simulate the different types of Newtonian and non-Newtonian pulsatile blood flow in a constricted as well as in a dilated channel to gain insight of the transition-to-turbulent blood flow due to the arterial stenosis and aneurysm. In the stenosed model, a cosine shape stenosis is placed at the upper wall of a 3D channel which reduces the cross-sectional area, whereas the aneurysm which is also placed at the upper wall dilates the channel cross-sectional area. In LES, a top-hat spatial grid-filter is applied to the Navier-Stokes equations of motion to separate the large scale flows, which carry the majority of the energy, from the small scale known as sub-grid scale (SGS).The large scale flows are resolved fully while the unresolved SGS motions are modelled using two different dynamic models to determine the Smagorinsky constant at each time step. Initially, an additive sinusoidal pulsatile velocity profile is used at the inlet of the model stenosis to generate the unsteady oscillating flow and a comparison is made between the results obtained by the additive and non-additive pulsation. Secondly, the physiological pulsatile flow in the same model stenosis is investigated, where the physiological pulsation is generated at the inlet using the first four harmonics of the Fourier series of pressure pulse. A comparison between the LES and the coarse Direct Numerical Simulation (DNS) results is drawn and the effects of the various harmonics of pressure pulse, length and percentage of the stenosis on the flow field are examined. Transition-to-turbulent physiological flow through the model of a double stenosis and an aneurysm is also investigated. Finally, the physiological pulsatile flow in a model of single stenosis is investigated using the various non-Newtonian blood viscosity models and the results are compared with the Newtonian model. For the additive sinusoidal pulsation case the maximum ratio of the SGS to molecular viscosity is 0.709 and for the non-additive case is 0.78 while Re =2000. The shape of the post-stenotic re-circulation region is totally different between the additive and non-additive case. In the additive case the upper wall pressure drop is larger than the non-additive case. Due to the large amplitude of the oscillation, transition happens earlier and the peak turbulent kinetic energy occurs at the post-lip of the stenosis. The intensity of the turbulent kinetic energy is higher in the additive sinusoidal pulsation case than the physiological pulsation. The maximum contribution of the SGS motion to the large -scale motion is 37.4 percent for the first harmonic physiological pulsation while 97 percent contribution from the first four harmonics case for Re =2000. The centreline turbulent kinetic energy is slightly higher in the first harmonic case than the first four harmonics. For the higher area reduction of the stenosis, the stress drop at the upper wall, the maximum shear stress at the lower wall and the turbulent kinetic energy increased. The intensity of the shear stress and the turbulent kinetic energy decreased when the length of the stenosis is increased. The break frequency of the energy spectra found from -5/3 to -10/3 for the velocity fluctuations and from -5/3 to -7/3 for the pressure fluctuations. Due to the presence of the second stenosis, the stress drop, the adverse pressure gradient and the turbulent intensity of the flow enhance significantly. Inside the aneurysm a large re-circulation region exists and the flow is turbulent for a asymmetric aneurysm and maximum turbulent intensity occurs between the centre and the ending segment of the aneurysm. Owing to the effects of the non-Newtonian viscosity, the length of the post-stenotic re-circulation region increased as well as the streamwise velocity and the turbulent kinetic energy decreased.
|