Effect of the fluid-structure interactions on low-density lipoprotein transport within a multi-layered arterial wall

The effects of fluid-structure interactions (FSI) and pulsation on the transport of low-density lipoprotein (LDL) through an arterial wall are analyzed in this work. To this end, a comprehensive multi-layer model for both LDL transport as well as fluid-structure interaction (FSI) is introduced. The constructed model is analyzed and compared with the existing results in the limiting cases. Excellent agreement is found between the presented model and the existing results in the limiting cases. The presented model takes into account the complete multi-layered LDL transport while incorporating the FSI aspects to enable a comprehensive study of the deformation effect on the pertinent parameters of the transport processes within an artery. Since the flow inside an artery is time-dependent, the impact of pulsatile flow is also analyzed with and without FSI. A detailed analysis is presented to illustrate the consequence of different factors on the LDL transport in an artery.     

Computational analysis of coupled blood-wall arterial LDL transport

The transport of macromolecules, such as low density lipoproteins (LDLs), across the artery wall and their accumulation in the wall is a key step in atherogenesis. Our objective was to model fluid flow within both the lumen and wall of a constricted, axisymmetric tube simulating a stenosed artery, and to then use this flow pattern to study LDL mass transport from the blood to the artery wall. Coupled analysis of lumenal blood flow and transmural fluid flow was achieved through the solution of Brinkman's model, which is an extension of the Navier-Stokes equations for porous media. This coupled approach offers advantages over traditional analyses of this problem, which have used possibly unrealistic boundary conditions at the blood-wall interface; instead, we prescribe a more natural pressure boundary condition at the adventitial vasa vasorum, and allow variations in wall permeability due to the occurrence of plaque. Numerical complications due to the convection dominated mass transport process (low LDL diffusivity) are handled by the streamline upwind/Petrov-Galerkin (SUPG) finite element method. This new fluid-plus-porous-wall method was implemented for conditions typical of LDL transport in a stenosed artery with a 75 percent area reduction (Peclet number=2 x 10(8)). The results show an elevated LDL concentration at the downstream side of the stenosis. For the higher Darcian wall permeability thought to occur in regions containing atheromatous lesions, this leads to an increased transendothelial LDL flux downstream of the stenosis. Increased transmural filtration in such regions, when coupled with a concentration-dependent endothelial permeability to LDL, could be an important contributor to LDL infiltration into the arterial wall. Experimental work is needed to confirm these results.     

Effects of transmural pressure and wall shear stress on LDL accumulation in the arterial wall: a numerical study using a multilayered model

The accumulation of low-density lipoprotein (LDL) is recognized as one of the main contributors in atherogenesis. Mathematical models have been constructed to simulate mass transport in large arteries and the consequent lipid accumulation in the arterial wall. The objective of this study was to investigate the influences of wall shear stress and transmural pressure on LDL accumulation in the arterial wall by a multilayered, coupled lumen-wall model. The model employs the Navier-Stokes equations and Darcy's Law for fluid dynamics, convection-diffusion-reaction equations for mass balance, and Kedem-Katchalsky equations for interfacial coupling. To determine physiologically realistic model parameters, an optimization approach that searches optimal parameters based on experimental data was developed. Two sets of model parameters corresponding to different transmural pressures were found by the optimization approach using experimental data in the literature. Furthermore, a shear-dependent hydraulic conductivity relation reported previously was adopted. The integrated multilayered model was applied to an axisymmetric stenosis simulating an idealized, mildly stenosed coronary artery. The results show that low wall shear stress leads to focal LDL accumulation by weakening the convective clearance effect of transmural flow, whereas high transmural pressure, associated with hypertension, leads to global elevation of LDL concentration in the arterial wall by facilitating the passage of LDL through wall layers.     

Concentration polarization of atherogenic lipids in the arterial system

Nomenclature c, Normalized LDL concentration (C*/C0); C0, incoming (bulk) LDL concentration (gr/cm3); Cw, LDL concentration on the luminal surface (gr/cm3); ,wC time average value of LDL concentration on the luminal surface (gr/cm3); D, diffusion coef-ficient of LDL (cm2/s); Q, blood flow rate (mL/s); 0R, average internal radius of the artery (cm); Re, Reynolds number (002/Run); Sc, Schmidt number (/Dn); t, normalized time (00*/tuR); u, normalized axial velocity (0*/uu); 0u, time a…     

Large eddy simulation of LDL surface concentration in a subject specific human aorta

The development of atherosclerosis is correlated to the accumulation of lipids in the arterial wall, which, in turn, may be caused by the build-up of low-density lipoproteins (LDL) on the arterial surface. The goal of this study was to model blood flow within a subject specific human aorta, and to study how the LDL surface concentration changed during a cardiac cycle. With measured velocity profiles as boundary conditions, a scale-resolving technique (large eddy simulation, LES) was used to compute the pulsatile blood flow that was in the transitional regime. The relationship between wall shear stress (WSS) and LDL surface concentration was investigated, and it was found that the accumulation of LDL correlated well with WSS. In general, regions of low WSS corresponded to regions of increased LDL concentration and vice versa. The instantaneous LDL values changed significantly during a cardiac cycle; during systole the surface concentration was low due to increased convective fluid transport, while in diastole there was an increased accumulation of LDL on the surface. Therefore, the near-wall velocity was investigated at four representative locations, and it was concluded that in regions with disturbed flow the LDL concentration had significant temporal changes, indicating that LDL accumulation is sensitive to not only the WSS but also near-wall flow.     

Concentration polarization of atherogenic lipids in the arterial system

The transport of atherogenic lipids (LDL) in a straight segment of an artery with a semi-permeable wall was simulated numerically. The numerical analysis predicted that a mass transport phenomenon called ’concentration polarization’ of LDL might occur in the arterial system. Under normal physiological flow conditions, the luminal surface LDL concentration was 5%–14% greater than the bulk concentration in a straight segment of an artery. The luminal surface LDL concentration at the arterial wall was flow-dependent, varying linearly with the filtration rate across the arterial wall and inversely with wall shear rate. At low wall shear rate, the luminal surface LDL concentration was very sensitive to changes in flow conditions, decreasing sharply as wall shear rate increased. In order to verify the numerical analysis, the luminal surface concentration of bovine serum albumin (as a tracer macromolecule) in the canine carotid artery was measured in vitro by directly taking liquid samples from the luminal surface of the artery. The experimental result was in very good agreement with the numerical analysis. The authors believe that the mass transport phenomenon of ‘concentration polarization’ may indeed exist in the human circulation and play an important role in the localization of atherosclerosis.     

Multiphysics simulation of blood flow and LDL transport in a porohyperelastic arterial wall model

Atherosclerosis localizes at a bend andor bifurcation of an artery, and low density lipoproteins (LDL) accumulate in the intima. Hemodynamic factors are known to affect this localization and LDL accumulation, but the details of the process remain unknown. It is thought that the LDL concentration will be affected by the filtration flow, and that the velocity of this flow will be affected by deformation of the arterial wall. Thus, a coupled model of a blood flow and a deformable arterial wall with filtration flow would be invaluable for simulation of the flow field and concentration field in sequence. However, this type of highly coupled interaction analysis has not yet been attempted. Therefore, we performed a coupled analysis of an artery with multiple bends in sequence. First, based on the theory of porous media, we modeled a deformable arterial wall using a porohyperelastic model (PHEM) that was able to express both the filtration flow and the viscoelastic behavior of the living tissue, and simulated a blood flow field in the arterial lumen, a filtration flow field and a displacement field in the arterial wall using a fluid-structure interaction (FSI) program code by the finite element method (FEM). Next, based on the obtained results, we further simulated LDL transport using a mass transfer analysis code by the FEM. We analyzed the PHEM in comparison with a rigid model. For the blood flow, stagnation was observed downward of the bends. The direction of the filtration flow was only from the lumen to the wall for the rigid model, while filtration flows from both the wall to the lumen and the lumen to the wall were observed for the PHEM. The LDL concentration was high at the lumenwall interface for both the PHEM and rigid model, and reached its maximum value at the stagnation area. For the PHEM, the maximum LDL concentration in the wall in the radial direction was observed at the position of 3% wall thickness from the lumenwall interface, while for the rigid model, it was observed just at the lumenwall interface. In addition, the peak LDL accumulation area of the PHEM moved about according to the pulsatile flow. These results demonstrate that the blood flow, arterial wall deformation, and filtration flow all affect the LDL concentration, and that LDL accumulation is due to stagnation and the presence of filtration flow. Thus, FSI analysis is indispensable.     

Structure of pulsatile flow in a model of elastic cerebral aneurysm

This study examines the effect of aneurysmal wall elasticity on the structure of flow within an elastic aneurysm during pulsatile flow. We visualized flow structure in a model of an elastic saccular aneurysm located at the bifurcation of the anterior cerebral artery and extending to the anterior communicating artery, and measured changes in the diameter of the aneurysm wall during pulsatile flow using particle imaging velocimetry (PIV). We similarly measured these features during steady flow by PIV and found that dilation of the aneurysmal wall absorbed the dynamic energy within the aneurysm. Accordingly, aneurysm wall elasticity functions as a biocompatible reaction that relieves wall shear stress acting on the vascular wall during pulsatile flow, and should thus inhibit the development and rupture of an aneurysm.     

Effects of diffusion coefficients and struts apposition using numerical simulations for drug eluting coronary stents

In the context of drug eluting stent, we present two-dimensional numerical models of mass transport of the drug in the wall and in the lumen to study the effect of the drug diffusion coefficients in the three principal media (blood, vascular wall, and polymer coating treated as a three-compartment problem) and the impact of different strut apposition configurations (fully embedded, half embedded, and not embedded). The different conditions were analyzed in terms of their consequence on the drug concentration distribution in the arterial wall. We apply the concept of the therapeutic window to the targeted vascular wall region and derive simple metrics to assess the efficiency of the various stent configurations. Although most of the drug is dispersed in the lumen, variations in the blood flow rate within the physiological range of coronary blood flow and the diffusivity of the drug molecule in the blood were shown to have a negligible effect on the amount of drug in the wall. Our results reveal that the amount of drug cumulated in the wall depends essentially on the relative values of the diffusion coefficients in the polymer coating and in the wall. Concerning the strut apposition, it is shown that the fully embedded strut configuration would provide a better concentration distribution.     

Theoretical study on flow-dependent concentration polarization of low density lipoproteins at the luminal surface of a straight artery

It is suspected that physical and fluid mechanical factors play important roles in the localization of atherosclerotic lesions and intimal hyperplasia in man by affecting the transport of cholesterol in flowing blood to arterial walls. Hence, we have studied theoretically the effects of various physical and fluid mechanical factors such as wall shear rate, diffusivity of low density lipoproteins (LDL), and filtration velocity of water at the vessel wall on surface concentration of LDL at an arterial wall by means of a computer simulation of convective and diffusive transport of LDL in flowing blood to the wall of a straight artery under conditions of a steady flow. It was found that under normal physiologic conditions prevailing in the human arterial system, due to the presence of a filtration flow of water at the vessel wall, flow-dependent concentration polarization (accumulation or depletion) of LDL occurs at a blood/endothelium boundary. The surface concentration of LDL at an arterial wall takes higher values than that in the bulk flow in that vessel, and it is affected by three major factors, that is, wall shear rate, gamma w, filtration velocity of water at the vessel wall, Vw, and the distance from the entrance of the artery, L. It increases with increasing Vw and L, and decreasing gamma w hence the flow rate. Thus, under certain circumstances, the surface concentration of LDL could rise locally to a value which is several times higher than that in the bulk flow, or drop locally to a value even lower than a critical concentration for the maintenance of normal functions and survival of cells forming the vessel wall. These results suggest the possibility that all the vascular phenomena such as the localization of atherosclerotic lesions and intimal hyperplasia, formation of cerebral aneurysms, and adaptive changes of lumen diameter and wall structure of arteries and veins to certain changes in hemodynamic conditions in the circulation are governed by this flow-dependent concentration polarization of LDL which carry cholesterol.     

Coupled computational analysis of arterial LDL transport -- effects of hypertension

Hypertension, a risk factor for atherosclerosis, increases the uptake of low density lipoproteins (LDL) by the arterial wall. Our objective in this work was to use computational modeling to identify physical factors that could be partially responsible for this effect. Fluid flow and mass transfer patterns in the lumen and wall of an arterial model were computed in a coupled manner, replicating as closely as possible previous experimental studies in which LDL uptake into the artery wall was measured in straight, excised arterial segments. Under conditions of both flow and no-flow, simulations predicted an increase in concentration polarization of LDL at the artery wall when arterial pressure was increased from 120 to 160 mmHg. However, this led to only a slight increase in mean LDL concentration within the arterial wall. However, if the permeability of the endothelium to LDL was allowed to vary with intra-arterial pressure, then the simulations predicted that the uptake of LDL would be enhanced 1.9-2.6 fold at higher pressure. The magnitude of this increase was consistent with experimental data. We conclude that the concentration polarization effects, enhanced by elevated intra-arterial pressure, cannot explain the increase in LDL uptake seen under hypertensive conditions. Instead, the data are most consistent with a pressure-linked increase in endothelial permeability to LDL.     

Effect of the endothelial glycocalyx layer on arterial LDL transport under normal and high pressure

To quantitatively investigate the role of the endothelial glycocalyx layer (EGL) in protecting the artery from excessive infiltration of atherogenic lipids such as low density lipoproteins (LDLs), a multilayer model with the EGL of an arterial segment was developed to numerically simulate the flow and the transport of LDLs under normal and high pressure. The transport parameters of the layers of the model were obtained from the hydrodynamic theory, the stochastic theory, and from the literature. The results showed that the increase in the thickness of the EGL could lead to a sharp drop in LDL accumulation in the intima. A partial damage to the EGL could compromise its barrier function, hence leading to enhanced infiltration/accumulation of LDLs within the wall of the arterial model. Without the EGL, hypertension could lead to a significantly enhanced LDL transport into the wall of the model. However, the intact EGL could protect the arterial wall from hypertension so that the LDL concentration in the intima layer was almost the same as that under normal pressure conditions. The results also showed that LDL concentration within the arterial wall increased with Φ (the fraction of leaky junctions) on the intima layer. The increase in LDL concentration with Φ was much more dramatic for the model without the EGL. For instance, without the EGL, a Φ of 0.0005 could lead LDL concentration within the arterial wall to be even higher than that predicted for the EGL intact model with a Φ of 0.002. In conclusion, an intact EGL with a sufficient thickness may act as a barrier to LDL infiltration into the arterial wall and has the potential to suppress the hypertension-driven hike of LDL infiltration/accumulation in the arterial wall.     

Computational study of LDL mass transport in the artery wall

A two-dimensional (2D) numerical simulation of convective–diffusive transport of LDL in the artery wall, coupled with the wall shear stress gradient (WSSG)-dependent LDL consumption of smooth muscle cells (SMCs) is presented. SMCs are modeled as an array of solid cylindrical pillars embedded in a continuous porous media which represents the interstitial proteoglycan and collagen fiber matrix. The internal elastic lamina (IEL), which separates the artery media from the intima, is modeled as an impermeable barrier to both water and LDL except for the fenestral pores that are assumed to be uniformly distributed over the IEL. The predictions demonstrate a range of interesting features of LDL transport and uptake in the media. For cells immediately below the fenestral pores, LDL uptake of SMCs is highly dependent on WSSG. Moreover, the rate of LDL consumption by SMCs is also affected by the diameter of the fenestral pore. This will be helpful in understanding the involvement of transmural transport processes in the initiation and development of atherosclerosis.     

The effect of the endothelial glycocalyx layer on concentration polarisation of low density lipoprotein in arteries

It has been postulated that a flow-dependent (and hence spatially varying) low density lipoprotein (LDL) concentration polarisation layer forms on the luminal surface of the vascular endothelium. Such a layer has the potential to cause heterogeneity in the distribution of atherosclerotic lesions by spatially modulating the rate of LDL transport into the arterial wall. Theoretical analysis suggests that a transmural water flux which is spatially heterogeneous at the cellular scale can act to enhance LDL concentration polarisation in a shear dependent fashion. However, such an effect is only observed if a relevant Peclet number (i.e. the ratio of LDL convection to LDL diffusion) is of order unity or greater. Based on the diffusivity of LDL in blood plasma, such a Peclet number is found to be far less than unity, implying that the aforementioned enhancement and shear dependence will not occur. However, this conclusion ignores the existence of the endothelial glycocalyx layer (EGL), which may inhibit the diffusion of LDL near the luminal surface of the endothelium, and hence raise any Peclet number associated with the transport of LDL. The present study numerically investigates the effect of the EGL, as well as a heterogeneous transmural water flux, on arterial LDL concentration polarisation. Particular attention is paid to measures of LDL concentration polarisation thought relevant to the rate of transendothelial LDL transport. It is demonstrated that an EGL is unlikely to cause any additional shear dependence of such measures directly, irrespective of whether or not LDL can penetrate into the EGL. However, it is found that such measures depend significantly on the nature of the interaction between LDL and the EGL (parameterised by the height of the EGL, the depth to which LDL penetrates into the EGL, and the diffusivity of LDL in the EGL). Various processes may regulate the interaction of LDL with the EGL, possibly in a flow dependent and hence spatially non-uniform fashion. It is concluded that any such processes may be as important as vascular scale flow features in terms of spatially modulating transendothelial LDL transport via an LDL concentration polarisation mechanism.     

Numerical simulation of the flow field and mass transport pattern within the coronary artery

In this study, the steady and pulsatile flow field with mass transport analysis in an anatomically correct model of coronary artery is simulated numerically using a specific patient data from a 64-multislice computed tomography scanner. It is assumed that the blood flow is laminar and that the Navier-Stokes equations of motion are applied. Downstream of the bifurcation, a strong skewing occurs towards the flow divider walls as a result of branching. For the low-density lipoprotein (LDL) transport analysis where a specific boundary condition at the arterial walls is applied, LDL is generally elevated at locations where shear stress distribution is low, but it does not co-locate at whole domain. This numerical simulation gives an insight, as well as detailed quantitative data, of haemodynamic conditions in the left coronary artery as well as mass transfer patterns for a specific patient.     

Analysis of hemodynamic fluid phase mass transport in a separated flow region

The mass transfer behavior in the recirculation region downstream of an axisymmetric sudden expansion was examined. The Reynolds number, 500, and Schmidt number, 3200, were selected to model the mass transfer of molecules, such as ADP, in the arterial system. In a first step the transient mass transport applying zero diffusive flux at the wall was analyzed using experiments and two computational codes. The two codes were FLUENT, a commercially available finite volume method, and FTSP, a finite element code developed at Graz University of Technology. The comparison of the transient wall concentration values determined by the three methods was excellent and provides a measure of confidence for computational mass transfer calculations in convection dominated, separated flows. In a second step the effect of the flow separation on the stationary mass transport applying a permeability boundary condition at the water-permeable wall was analyzed using the finite element code FTSP. The results show an increase of luminal ADP surface concentration in the upstream and in the downstream tube of the sudden expansion geometry in the range of six and twelve percent of the bulk flow concentration. The effect of flow separation in the downstream tube on the wall concentration is a decrease of about ten percent of the difference between wall concentration and bulk concentration occurring at nearly fully developed flow at the downstream region at a distance of 66 downstream tube diameters from the expansion. The decrease of ADP flux into the wall is in the range of three percent of the flux at the downstream region.     

Flow-tissue interaction with compliance mismatch in a model stented artery

The insertion of an endovascular prosthesis is known to have a thrombogenic effect that is also a consequence of the interaction between the flowing blood and the stented arterial segment; in fact the prosthesis induces a compliance mismatch and a possible small expansion along the vessel that eventually gives rise to an anomalous distribution of wall shear stresses. The fluid dynamics inside a rectilinear elastic vessel with compliance and section variation is studied here numerically. A recently introduced perturbative approach is employed to model the interaction between the fluid and the elastic tissue; this approximate technique is first validated by comparison with a complete solution within a simple one-dimensional model of the same system. Then it is applied to an axisymmetric model in order to evaluate the flow dynamics and the distribution of wall shear stress in the stented vessel. Compliance mismatch is shown to produce more intense negative wall shear stresses in the stented segment while rapid variations of wall shear stress are found at the stent ends. These effects are enhanced when the prosthesis is accompanied by a small increase of the vessel lumen.     

Effect of non-Newtonian and pulsatile blood flow on mass transport in the human aorta

To investigate the effects of both non-Newtonian behavior and the pulsation of blood flow on the distributions of luminal surface LDL concentration and oxygen flux along the wall of the human aorta, we numerically compared a non-Newtonian model with the Newtonian one under both steady flow and in vivo pulsatile flow conditions using a human aorta model constructed from MRI images. The results showed that under steady flow conditions, although the shear thinning non-Newtonian nature of blood could elevate wall shear stress (WSS) in most regions of the aorta, especially areas with low WSS, it had little effect on luminal surface LDL concentration (c(w)) in most regions of the aorta. Nevertheless, it could significantly enhance c(w) in areas with high luminal surface LDL concentration through the shear dependent diffusivity of LDLs. For oxygen transport, the shear thinning non-Newtonian nature of blood could slightly reduce oxygen flux in most regions of the aorta, but this effect became much more apparent in areas with already low oxygen flux. The pulsation of blood flow could significantly reduce c(w) and enhance oxygen flux in these disturbed places. In most other regions of the aorta, the oxygen flux was also significantly higher than that for the steady flow simulation. In conclusion, the shear shining non-Newtonian nature of blood has little effect on LDL and oxygen transport in most regions of the aorta, but in the atherogenic-prone areas where luminal surface LDL concentration is high and oxygen flux is low, its effect is apparent. Similar is for the effect of pulsatile flow on the transport of LDLs. But, the pulsation of blood flow can apparently affect oxygen flux in the aorta, especially in areas with low oxygen flux.