Comparison of LES of steady transitional flow in an idealized stenosed axisymmetric artery model with a RANS transitional model.

In this study, two different turbulence methodologies are investigated to predict transitional flow in a 75% stenosed axisymmetric experimental arterial model and in a slightly modified version of the model with an eccentric stenosis. Large eddy simulation (LES) and Reynolds-averaged Navier-Stokes (RANS) methods were applied; in the LES simulations eddy viscosity subgrid-scale models were employed (basic and dynamic Smagorinsky) while the RANS method involved the correlation-based transitional version of the hybrid k-ε/k-ω flow model. The RANS simulations used 410,000 and 820,000 element meshes for the axisymmetric and eccentric stenoses, respectively, with y(+) less than 2 viscous wall units for the boundary elements, while the LES used 1,200,000 elements with y(+) less than 1. Implicit filtering was used for LES, giving an overlap between the resolved and modeled eddies, ensuring accurate treatment of near wall turbulence structures. Flow analysis was carried out in terms of vorticity and eddy viscosity magnitudes, velocity, and turbulence intensity profiles and the results were compared both with established experimental data and with available direct numerical simulations (DNSs) from the literature. The simulation results demonstrated that the dynamic Smagorinsky LES and RANS transitional model predicted fairly comparable velocity and turbulence intensity profiles with the experimental data, although the dynamic Smagorinsky model gave the best overall agreement. The present study demonstrated the power of LES methods, although they were computationally more costly, and added further evidence of the promise of the RANS transition model used here, previously tested in pulsatile flow on a similar model. Both dynamic Smagorinsky LES and the RANS model captured the complex transition phenomena under physiological Reynolds numbers in steady flow, including separation and reattachment. In this respect, LES with dynamic Smagorinsky appeared more successful than DNS in replicating the axisymmetric experimental results, although inflow conditions, which are subject to caveats, may have differed. For the eccentric stenosis, LES with Smagorinsky coefficient of 0.13 gave the closest agreement with DNS despite the known shortcomings of fixed coefficients. The relaminarization as the flow escaped the influence of the stenosis was amply demonstrated in the simulations, graphically so in the case of LES.

Analysis of flow patterns in a patient-specific aortic dissection model.

Aortic dissection is the most common acute catastrophic event affecting the thoracic aorta. The majority of patients presenting with an uncomplicated type B dissection are treated medically, but 25% of these patients develop subsequent aneurysmal dilatation of the thoracic aorta. This study aimed at gaining more detailed knowledge of the flow phenomena associated with this condition. Morphological features and flow patterns in a dissected aortic segment of a presurgery type B dissection patient were analyzed based on computed tomography images acquired from the patient. Computational simulations of blood flow in the patient-specific model were performed by employing a correlation-based transitional version of Menter's hybrid k-epsilon/k-omega shear stress transport turbulence model implemented in ANSYS CFX 11. Our results show that the dissected aorta is dominated by locally highly disturbed, and possibly turbulent, flow with strong recirculation. A significant proportion (about 80%) of the aortic flow enters the false lumen, which may further increase the dilatation of the aorta. High values of wall shear stress have been found around the tear on the true lumen wall, perhaps increasing the likelihood of expanding the tear. Turbulence intensity in the tear region reaches a maximum of 70% at midsystolic deceleration phase. Incorporating the non-Newtonian behavior of blood into the same transitional flow model has yielded a slightly lower peak wall shear stress and higher maximum turbulence intensity without causing discernible changes to the distribution patterns. Comparisons between the laminar and turbulent flow simulations show a qualitatively similar distribution of wall shear stress but a significantly higher magnitude with the transitional turbulence model.

Analysis of flow disturbance in a stenosed carotid artery bifurcation using two-equation transitional and turbulence models.

In this study, newly developed two-equation turbulence models and transitional variants are employed for the prediction of blood flow patterns in a diseased carotid artery where the growth, progression, and structure of the plaque at rupture are closely linked to low and oscillating wall shear stresses. Moreover, the laminar-turbulent transition in the poststenotic zone can alter the separation zone length, wall shear stress, and pressure distribution over the plaque, with potential implications for stresses within the plaque. Following the validation with well established experimental measurements and numerical studies, a magnetic-resonance (MR) image-based model of the carotid bifurcation with 70% stenosis was reconstructed and simulated using realistic patient-specific conditions. Laminar flow, a correlation-based transitional version of Menter's hybrid k-epsilon/k-omega shear stress transport (SST) model and its "scale adaptive simulation" (SAS) variant were implemented in pulsatile simulations from which analyses of velocity profiles, wall shear stress, and turbulence intensity were conducted. In general, the transitional version of SST and its SAS variant are shown to give a better overall agreement than their standard counterparts with experimental data for pulsatile flow in an axisymmetric stenosed tube. For the patient-specific case reported, the wall shear stress analysis showed discernable differences between the laminar flow and SST transitional models but virtually no difference between the SST transitional model and its SAS variant.