Influence of Oscillating Flow on LDL Transport and Wall Shear Stress in the Normal Aortic Arch
J Soulis1 , *, G Giannoglou2, M Dimitrakopoulou 1, V Papaioannou1, S Logothetides3, D Mikhailidis4
Identifiers and Pagination:Year: 2009
First Page: 128
Last Page: 142
Publisher Id: TOCMJ-3-128
Article History:Received Date: 17/8/2009
Revision Received Date: 24/8/2009
Acceptance Date: 25/8/2009
Electronic publication date: 17/9/2009
Collection year: 2009
open-access license: This is an open access article licensed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted, non-commercial use, distribution and reproduction in any medium, provided the work is properly cited.
Lipid accumulation in the aortic wall is an important factor in the development of atherosclerosis. The Low Density Lipoprotein (LDL) at the surface of the endothelium in relation to Wall Shear Stress (WSS) in the normal human aortic arch under unsteady, normal flow and mass conditions was computationally analysed. Concave sides of the aortic arch exhibit, relatively to the convex ones, elevated LDL levels at the surface of the endothelium for all time steps. At the peak systolic velocity, the LDL level reaches a value 23.0% higher than that at entrance in the ascending-descending aorta region. The corresponding LDL levels at the surface of the endothelium for the near minimum entrance velocity instant reaches 26.0%. During the cardiac cycle, the highest area averaged normalized LDL taken up as compared to the lowest one is 0.69%. WSS plays an important role in the lipid accumulation. Low WSS regions are exposed to high LDL levels at the surface of the endothelium. Regions of elevated LDL levels do not necessarily co-locate to the sites of lowest WSS. The near wall paths of the velocities might be the most important factor for the elevated LDL levels at the surface of the endothelium.
Lipid accumulation in the aortic wall is an important factor in the development of atherosclerosis. Disturbed blood flow and occurrence of atherosclerotic plaques at certain sites of human arteries are strongly correlated . The flow of blood and the transport of macromolecules in the vascular system  are essential to understand atherogenesis . Biomechanical factors, such as Wall Shear Stress (WSS), blood viscosity and flow velocity, may be responsible for the localization and progression of atherosclerosis [4, 5]. Besides the spatial distribution of flow parameters, their temporal variation during the cardiac cycle has been proposed as an atherogenic factor . Complex configuration of the vessel promotes flow disturbances with low time-averaged WSS and high shear stress temporal oscillations during the cardiac cycle . The application of unsteady flow constitutes a challenging issue, since it demands computational power and analysis time. The WSS may affect endothelial permeability . Regional variations in the permeability of arterial endothelium may contribute to the localization of atheroscle- rosis [8, 9]. Luminal surface level of Low Density Lipoprotein (LDL) transport across the artery wall is considered to be important in atherogenesis [10, 11]. Elevated LDL co-localize with known sites of atherosclerotic plaque development, . Atherosclerotic regions are also low WSS regions and consequently it was difficult to determine whether it was low WSS, abnormal mass transfer or both that were contributing to atherosclerosis . Oscillating flow on LDL transport in the arterial wall showed that the steady flow assumption is inadequate and the instantaneous hemodynamic conditions have important influence on LDL transmural transport in arteries with disturbed and complicated flow patterns . The accumulation of LDL in vascular areas featuring a highly disturbed flow was examined, . Geometrical parameters such as curvature and variations of the luminal section strongly influence the LDL within the wall. Realizing that the research produced in macro-molecular transport in curved 3D blood vessels was rare, steady and unsteady flows and mass simulation was analyzed . The effect of near-wall blood flow velocity and plasma filtration velocity across the arterial wall on luminal surface concentration of LDL and the uptake of tritium-cholesterol were investigated . These results indicate that lipids accumulate at the luminal surface in areas where blood flow velocity and WSS are low and where the permeability of the endothelial layer is enhanced.
There have been very few studies investigating the fluid-wall coupled mass transport under oscillating flow conditions [12, 14, 16]. The current work analyzes the influence of oscillating flow and LDL mass transport patterns over the normal human aortic arch. The purpose is to elucidate the association of low WSS-elevated LDL levels at the surface of the endothelium in relation to atherosclerotic plaque localization over the physiological human aortic arch during the normal cardiac pulse. Emphasis is put into a) LDL, b) WSS, c) differentiation of the LDL and WSS between concave (outer) and convex (inner) aortic arch sides as well as into d) temporal differentiation of the LDL between end systolic and end diastolic periods. We demonstrate that due to the semi-permeable nature of the arterial walls, high LDL levels at the surface of the endothelium occur at certain regions of the normal aortic arch.
2. Materials and Methodology
The 3D geometry of the normal aortic arch computational model was generated using data compiled from several sources, Fig. (1). The computational model includes the ascending aorta, descending aorta, brachiocephalic artery, left common carotid artery and left subclavian artery.
Aortic arch geometry and non-structured computational grid including the: ascending aorta, descending aorta, brachiocephalic artery, left common carotid artery and left subclavian artery.
Average physiological human aortic arch velocity waveform. Computational analysis results refer to inlet flow velocities of a) 0.32 m/s (peak systole) and b) 0.005 m/s (near minimum entrance velocity).
Colour-coded contour plots of the Wall Shear Stress (N/m2) magnitude for the ascending-descending aorta at a) peak systole and b) near minimum entrance velocity, respectively.
Contour plots of the normalized LDL levels at the surface of the endothelium Cw/Co for the ascending-descending aorta of the aortic arch at a) peak systole and b) near minimum entrance velocity, respectively.
Contour plots of the normalized LDL levels at the surface of the endothelium Cw/Co for the convex side of the aortic arch at a) peak systole and b) near minimum entrance velocity, respectively.
Contour plots of the normalized LDL levels at the surface of the endothelium Cw/Co for the daughter vessels of the aortic arch at a) peak systole and b) near minimum entrance velocity, respectively.
Contour plots of the velocity magnitude (m/s) for an ascending-descending aorta section at a) peak systole and b) near minimum entrance velocity, respectively.
Contour plots of the strain rate (1/s) for an ascending-descending aorta section at a) peak systole and b) near minimum entrance velocity, respectively.
Strain rate (1/s) and molecular viscosity (kg/m-s) contour plots for ascending-descending aorta cross-sections at a) peak systole, and, b) near minimum entrance velocity, respectively.
Normalized LDL levels at the surface of the endothelium versus Wall Shear Stress (N/m2) for the aortic arch at a) peak systole, and, b) near minimum entrance velocity, respectively.
The assumptions made about the nature of the flow are that it is 3D, unsteady, laminar, isothermal, with no external forces applied on it while the aortic arch wall is comprised from non-elastic and semi-permeable material. The blood is considered to be non-Newtonian fluid obeying the power law. The flow equation is coupled with the convection-diffusion equation. The diffusion flux of LDL arises due to concentration gradients. In the absence of systematic and reliable experimental data for human arteries, we assume the molecular diffusivity D to be constant and equal to 15.0x10-12m2/s [9, 17]. The infiltration velocity Vw is set equal to 0.6x10-8 m/s . There is a limited amount of experimental data regarding the endothelial permeability K or else mass transfer coefficient. Specifications of K values are sometimes unreliable. For current analysis purposes this value is set equal to 2.0x10-10 m/s, .
A typical aortic arch blood flow waveform, shown in Fig. (2), is applied over the entire entrance cross-section. The waveform mimics typical aortic blood averaged flow velocity under resting condition. From the entire velocity waveform, lasting 800.0 msec, two characteristic time-points are selected for detail analysis. The first instant refers to peak systole occurring 123 msec after the systole onset with 0.32 m/s entrance velocity. The second instant refers to near minimum entrance velocity of 0.005 m/s corresponding to 522.75 msec. The inlet based Reynolds numbers range from 3.5 (minimum) to 475 (peak).
A uniform constant concentration Co of LDL is applied at the orifice of the ascending aorta. At the descending aorta, brachiocephalic artery, left common carotid artery and left subclavian artery outlets, the gradient of LDL concentration along the vessels is set equal to zero (zero flux, Newmann condition). Flow discharges are set analogous to the third power of the branching vessel inlet diameter according to Murray’s law .
The applied endothelium boundary condition at the semi-permeable aortic walls states that the net amount of LDL per unit area passing from it to the vessel wall is determined by the difference of the mass flow carried to the vessel wall by infiltration flow and the amount of flow which diffuses back to the main vessel flow. It is known that at arterial walls the LDL concentration values are higher than those of the vessel bulk flow. Furthermore, LDL concentration increases with increasing infiltration velocity.
The calculated LDL levels at the surface of the endothelium Cw are normalized with the inlet value of Co (=1.3 mg/ml). WSS contours show its magnitude and not its spatial direction.
3.1. WSS in the Ascending-Descending Aorta
WSS (N/m2) magnitude contour plots for the ascending-descending aorta of the aortic arch for the selected time instants are shown in Figs. (3a, b). The convex side of the aortic arch exhibits, relatively to the concave side, high WSS values. High WSS is encountered at the convex side of the end ascending and early descending aorta at all pulse wave times. Particular low WSS values appear along the concave parts of the ascending-descending aorta at regions located just downstream to the left subclavian artery.
3.2. LDL in the Ascending-Descending Aorta
Contour plots of the normalized luminal surface concentration Cw/Co for the ascending-descending aorta, at the selected characteristic time instants, are shown in Figs. (4a, b). At all time instants, the concave side of the aortic arch exhibits high LDL levels compared with the convex side. Low LDL levels at the surface of the endothelium are encountered at the end ascending and the early descending convex side of the aorta, while the main convex side encounters relatively high LDL levels Figs. (5aa, b). At peak systolic pressure, the lowest Cw/Co value is 1.13. High LDL values appear on the luminal surface of the wall along the concave parts of the ascending-descending aorta. Regions with elevated LDL levels at the surface of the endothelium occur at the concave ascending-descending aortic arch just downstream to the left subclavian artery.
3.3. LDL in the Brachiocephalic Artery, Left Common Carotid Artery and Left Subclavian Artery
Cw/Co contour plots for the brachiocephalic artery, left common carotid artery and left subclavian artery at the specific time instants are shown in Figs. (6a, b). Elevated LDL levels at the surface of the endothelium are observed in all vessels. This is due to low flow velocities which in turn gives rise to low WSS. The particular LDL levels are strongly affected by vessel curvature. As a rule, concave parts exhibit high luminal surface LDL. However, depending upon the particular bend side and the incoming flow velocity distribution, convex parts can exhibit high LDL levels.
Elevated LDL regions occur at the inlet surfaces of the brachiocephalic artery, left common carotid artery as well as at the downstream regions of the left subclavian artery.
3.4. Velocity and Strain Rate for an Ascending-Descending Aorta Transverse Section
Velocity magnitude and strain rate distributions for an ascending-descending aorta transverse section, at peak systolic pressure and near minimum entrance velocity instant, are shown in Figs. (7 and 8), respectively.
Peak systolic pressure velocities occur in regions located close to the concave part of the early descending aorta. This is due to the centrifugal forces acting on blood flow as it encounters the aorta bend. Low blood flow velocities are encountered in the concave ascending-descending aorta at regions located just downstream to left subclavian artery. Low strain rates regions cover large area, Fig. (8), shifting towards concave parts of the ascending-descending aorta.
3.5. Strain Rate and Molecular Viscosity at Ascending-Descending Aorta Cross-Sections
Strain rate and molecular viscosity contour plots at various aorta cross-sections for the extreme flow instant i.e. at peak systolic and at near minimum velocity entrance, are shown in Figs. (9a, b), respectively. At peak systolic pressure, in the upstream region of the descending aorta, increased molecular viscosity and low strain rates are shifted towards the concave side and cover a relatively wide area, Fig. (9a); at the same time instant, high strain rates appear at the convex side of the upstream descending aorta region. At near minimum entrance velocity instant, increased strain rates are also present in the above region. Throughout the aortic arch the molecular viscosity increases as the inlet flow velocity decreases, Fig. (9bb).
Typical relationships between LDL levels and WSS over the entire aortic arch surface at the extreme velocity instants, are shown in Figs. (10a, b). It is evident that the LDL levels at the surface of the endothelium increases with decreasing WSS. The concentration increases at a higher rate as the WSS values reduce to zero. WSS plays an important role in LDL levels at the surface of the endothelium. The dependence of LDL levels at the surface of the endothelium on WSS is higher at near minimum entrance velocity instant, Fig. 10b, than at peak velocity instant, Fig. (10a).
3.7. LDL Statistics
Table 1 shows the descriptive LDL statistics for the aortic arch. The total endothelial area of the computational analyzed aortic arch is 360.69 cm2. Results, at 5 specific time-points of the cardiac cycle, show the a) area averaged LDL levels at the surface of endothelium, b) normalized area averaged LDL levels at the surface of endothelium, c) amount of LDL taken up by all luminal aortic arch surfaces, and d) percentage LDL levels at the surface of the endothelium as compared to its lowest value taken up by all luminal aortic arch surfaces. It is worth noting that the maximum percentage increase of LDL levels at the surface of the endothelium taken up as compared with the lowest LDL at the surface of the endothelium within the cardiac cycle is 0.69 %.
Descriptive Statistics of a) Luminal Surface Concentration, b) Area Averaged Normalized Luminal Surface Concentration, c) Area Averaged Amount of LDL Taken up by all Luminal Aortic Arch Surfaces, and, d) Percentage Increase of LDL as Compared with the Lowest LDL Taken up by all Luminal Aortic Arch Surfaces at Specific Time-Points (Inlet Flow Velocities) of the Cardiac Cycle
|Time(msec)||Average Inlet Flow Velocity (m/s)||a) Area Averaged LDL Levels at the Surface of Endothelium (mg/ml)||b) Area Averaged Normalized LDL Levels at the Surface of the Endothelium, Cw/Co||c) LDL Taken up by all Luminal Surfaces (g/s)||d) LDL Taken up as Compared to the Lowest LDL %|
Spatial and temporal WSS variation is important in regulating the atheroprotective, normal physiology as well as the pathobiology and dysfunction of the vessel wall through complex molecular mechanisms that promote atherogenesis. The endothelial transport properties respond instantaneously to WSS. Wherever and whenever flow becomes abnormal, the balance between the mass flow of substances from main blood stream to the arterial wall and vice versa is disrupted. This unbalanced mass transfer probably initiates arterial wall disease.
In our analysis the hydraulic conductivity of the endothelium was considered to be invariant with WSS. We computationally analyzed the LDL distribution at the surface of the endothelium and tried to elucidate the elevated LDL regions using the healthy aortic arch. Elevated LDL may be responsible for the localization and progression of atherosclerosis. The LDL substance was considered as solute within plasma. This assumption is quite helpful for computational analysis and probably it is the first step for treating the LDL as discrete particle movement. It is known that High Density Lipoprotein (HDL) facilitates the removal of the LDL from the luminal wall. However, in our work we did not consider analyzing their interaction because of the computational complexity.
The problem of mass transfer quantification is crucial in the genesis and the progression of atherosclerosis . The preferred occurrence of atherosclerosis in low WSS regions is supported with extensive investigation of arterial fluid flow and mass transfer across the endothelium . The action of mass transport occurs within a very thin layer, which is located very close to the wall. Any perturbation of the flow quantities, as for example high or low strain rates or flow reversal, in this sensitive flow region initiates mass flow disturbance to and from the arterial wall. Elongation of the endothelial cells occurs in regions of high WSS having the longest cell axis oriented parallel to the flow direction . Conversely, in low WSS regions endothelial cells of polygonal shape occur and they are not oriented in any direction. It is also possible that these local alterations in endothelial cell morphology might yield different size permeability for the various blood constituents. The influences of WSS affecting the wall permeability on steady LDL transport in human right coronary arteries are starting to appear in the literature . They also recognized that it is vital to consider blood pulsatility when modelling the shear-dependent macromolecular transport in large arteries.
The relationship between WSS and LDL levels indicates that the LDL is elevated at locations where WSS is low. This is a rule for all time steps. At low WSS and as it approaches zero, the LDL levels at the surface of the endothelium are rapidly increasing, Figs. (10a, b). The same relationship shows that a lot of scattered points can be found at all applied waveform time instants. This is an indication that LDL levels at the surface of the endothelium are not only WSS dependent. The LDL-WSS distribution at near minimum entrance velocity, Fig. (10b), exhibits elevated LDL levels at the surface of the endothelium at low WSS values compared with peak velocity entrance time instant, Fig. (10a). This supports the hypothesis that low WSS results into high LDL levels at the surface of the endothelium.Under the action of hemodynamics the LDL particles move to reach a particular aortic arch site. Thereafter, it is the contact time and the interaction between LDL and aortic endothelial surface which really matters. The net amount of LDL passing into arterial wall depends upon the physical properties of the wall material. High luminal surface concentration does not necessarily denote that all the LDL molecules will be transported into the arterial wall. It is the permeability of the endothelium, which largely determines the final amount of mass passing through the wall.
Peak LDL regions of the ascending-descending aorta are located at the concave side of the descending aorta just downstream to left subclavian artery, while convex sides result in low LDL. All aortic arch walls are exposed to cholesterolemic environment although the applied mass and flow conditions refer to normal human geometry and normal mass-flow conditions.
The convex sides of the ascending–descending aortic arch exhibit, relatively to the concave side, small normalized LDL levels at all time steps. The low LDL region at peak systole, Fig. (6a), covers smaller endothelial area compared to minimum velocity entrance time instant, Fig. (6b). The mechanism of elevated LDL levels at the surface of the endothelium might also be attributed to the high curvature effects upon the flow pattern and the subsequent spatial accumulation of the circulating LDL as the mass transport equation specifies. It is the near wall paths of the velocities which might be the most important factor for elevated LDL levels at the surface of the endothelium for areas located either at the vicinity of bifurcations regions or at high curvature regions, Figs. (7, 8 and 9). These paths are strongly affected from the combined inlet velocity field (oscillation)-particular geometry. Subsequently, the developing shear stress flow field and particularly the molecular viscosities dictate the WSS which in turn alters the LDL at the surface according to the physical properties of the endothelium. Because of the flow complexity and the final flow pattern in the near to endothelium regions the concave side of the ascending-descending aortic arch exhibits, relatively to the convex one, elevated LDL levels at all time steps.
Regions of low WSS values do not necessary coincide with elevated areas of LDL levels at the surface of the endothelium. At all time instants the WSS pattern clearly exhibits its low regions. The pattern of LDL, shown in Fig. (4), does not exhibit its low regions as clear as the WSS does. At near minimum entrance velocity the LDL levels at the surface of the endothelium taken up as compared with the peak systolic LDL is 0.69 % higher, Table 1. It is known that atherosclerotic regions are localized into certain aortic arch areas. Atherosclerosis of the ascending aorta is a predictor of late stroke. High risk is predominantly linked to atheromas in its distal part and lesser curvature .Ascending aortic and arch atheroma detected by trans- esophageal echocardiography is an important new independent risk factor for cerebral ischemia . Furthermore, an atherosclerotic ascending aorta may be a marker of generalized atherosclerosis and thus of increased morbidity and mortality . Our results show that there is elevated LDL at the surface of the endothelium in regions of low WSS. At the concave side of the descending aorta just downstream to left subclavian artery the elevated LDL sites coincide with Farthing findings . How these findings influence local pathology (e.g. aortic dissection, aneurysm formation or aortic arch emboli) remains to be elucidated.
The assumption, that wall permeability to LDL is constant throughout needs to be altered. A novel aspect will be that the permeability is a function of either the WSS or spatial gradient of WSS. Arterial wall deformation needs to be taken into consideration, since the flow pattern-mass transport is seriously determined from the geometrical configuration of the artery. Arterial movement has to be taken into account since the applied acceleration force will yield a different flow field. Hypolipidaemic and antihypertensive agents may affect the localization of aortic arch atherosclerotic regions; this possibility should be investigated. Finally, the percentage of LDL deposition at the surface of the endothelium, in terms of the circulating LDL quantity, needs to be analyzed.
This computational study elucidates the LDL transport in relation to the WSS in the normal human aortic arch under oscillating flow conditions at rest. LDL levels at the surface of the arterial wall endothelium are flow-dependent, exhibiting increased permeation of LDL at specific regions. Concave sides of the aortic arch exhibit, relatively to the convex ones, elevated LDL levels at the surface of the endothelium at all times. Aortic arch walls are exposed to cholesterolemic environment although the applied mass and flow conditions refer to normal human geometry and normal mass-flow conditions. Low WSS regions are exposed to high LDL levels at the surface of the endothelium. Regions of high LDL levels at the surface of the endothelium do not necessarily co-locate to the sites of low WSS. For the peak entrance velocity, in the ascending-descending aorta region, the LDL level at the surface of the endothelium is 23.0% higher than that at entrance. The corresponding LDL levels at the surface of the endothelium for the near minimum entrance velocity reaches 26.0 %. The net amount of LDL mass per second taken up by all luminal aortic arch surfaces, also at peak velocity, is 1.17865 x10-11 g/s. The highest area averaged normalized LDL taken up as compared with the lowest one is 0.69 %. The near wall paths of the velocities might be the most important factor for elevated LDL levels at the surface of the endothelium.