Blood flow in the rabbit aortic arch and descending thoracic aorta.

The distribution of atherosclerotic lesions within the rabbit vasculature, particularly within the descending thoracic aorta, has been mapped in numerous studies. The patchy nature of such lesions has been attributed to local variation in the pattern of blood flow. However, there have been few attempts to model and characterize the flow. In this study, a high-order continuous Galerkin finite-element method was used to simulate blood flow within a realistic representation of the rabbit aortic arch and descending thoracic aorta. The geometry, which was obtained from computed tomography of a resin corrosion cast, included all vessels originating from the aortic arch (followed to at least their second generation) and five pairs of intercostal arteries originating from the proximal descending thoracic aorta. The simulations showed that small geometrical undulations associated with the ductus arteriosus scar cause significant deviations in wall shear stress (WSS). This finding highlights the importance of geometrical accuracy when analysing WSS or related metrics. It was also observed that two Dean-type vortices form in the aortic arch and propagate down the descending thoracic aorta (along with an associated skewed axial velocity profile). This leads to the occurrence of axial streaks in WSS, similar in nature to the axial streaks of lipid deposition found in the descending aorta of cholesterol-fed rabbits. Finally, it was observed that WSS patterns within the vicinity of intercostal branch ostia depend not only on local flow features caused by the branches themselves, but also on larger-scale flow features within the descending aorta, which vary between branches at different locations. This result implies that disease and WSS patterns in the vicinity of intercostal ostia are best compared on a branch-by-branch basis.

Effect of reverse flow on the pattern of wall shear stress near arterial branches.

Atherosclerotic lesions have a patchy distribution within arteries that suggests a controlling influence of haemodynamic stresses on their development. The distribution near aortic branches varies with age and species, perhaps reflecting differences in these stresses. Our previous work, which assumed steady flow, revealed a dependence of wall shear stress (WSS) patterns on Reynolds number and side-branch flow rate. Here, we examine effects of pulsatile flow. Flow and WSS patterns were computed by applying high-order unstructured spectral/hp element methods to the Newtonian incompressible Navier-Stokes equations in a geometrically simplified model of an aorto-intercostal junction. The effect of pulsatile but non-reversing side-branch flow was small; the aortic WSS pattern resembled that obtained under steady flow conditions, with high WSS upstream and downstream of the branch. When flow in the side branch or in the aortic near-wall region reversed during part of the cycle, significantly different instantaneous patterns were generated, with low WSS appearing upstream and downstream. Time-averaged WSS was similar to the steady flow case, reflecting the short duration of these events, but patterns of the oscillatory shear index for reversing aortic near-wall flow were profoundly altered. Effects of reverse flow may help explain the different distributions of lesions.

Endothelial nitric oxide production and transport in flow chambers: The importance of convection.

A computational model of Nitric Oxide (NO) production and transport within a parallel-plate flow chamber coated with endothelial cells is presented. The relationship between NO concentration and Wall Shear Stress (WSS) at the endothelium is investigated in detail. An increase in WSS is associated with two phenomena: enhanced NO production by the endothelial cells, and an increase in the velocity at which NO is convected out of the chamber. These two phenomena have opposite effects on endothelial NO concentration. In physiologically realistic cases, the balance between them is found to vary as WSS is raised, resulting in a complex non-monotonic dependence of endothelial NO concentration on WSS. Also, it is found that a NO concentration boundary layer develops within the chamber, leading to substantial spatial variations in NO concentration along the length of the device. Finally, the implications of a negative feedback mechanism (that affects NO production) are presented. The results emphasize the role of convection on NO transport within flow chambers, which has been overlooked or misinterpreted in most previous theoretical studies. It is hoped that the conclusions of this study can be used to aid accurate interpretation of related experimental data.