Boundary layer considerations in a multi-layer model for LDL accumulation

Abstract

Boundary layer effects for Low-Density Lipoprotein (LDL) concentration problems in a multi-layer artery model are analyzed in this work. Both a straight artery and aorta-iliac bifurcation are analyzed. Mass, momentum and species governing equations are based on the porous media theory and solved with the commercial finite-element based code COMSOL Multiphysics. For the straight artery, various inlet velocities, arterial sizes and intramural pressure values are investigated. Results are presented in terms of concentration profiles close to the lumen/endothelium interface and boundary layer thickness. It is shown that the boundary layer is affected by all of the three analyzed parameters. The results in this work will further clarify the concentration polarization effects imposed by the arterial wall.     

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.     

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.     

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.     

Internal elastic lamina affects the distribution of macromolecules in the arterial wall: a computational study

The internal elastic lamina (IEL), which separates the arterial intima from the media, affects macromolecular transport across the medial layer. In the present study, we have developed a two-dimensional numerical simulation model to resolve the influence of the IEL on convective-diffusive transport of macromolecules in the media. The model considers interstitial flow in the medial layer that has a complex entrance condition because of the presence of leaky fenestral pores in the IEL. The IEL was modeled as an impermeable barrier to both water and solute except for the fenestral pores that were assumed to be uniformly distributed over the IEL. The media were modeled as a heterogeneous medium composed of an array of smooth muscle cells (SMCs) embedded in a continuous porous medium representing the interstitial proteoglycan and collagen fiber matrix. Results for ATP and low-density lipoprotein (LDL) demonstrate a range of interesting features of molecular transport and uptake in the media that are determined by considering the balance among convection, diffusion, and SMC surface reaction. The ATP concentration distribution depends strongly on the IEL pore structure because ATP fluid-phase transport is dominated by diffusion emanating from the fenestral pores. On the other hand, LDL fluid-phase transport is only weakly dependent on the IEL pore structure because convection spreads LDL laterally over very short distances in the media. In addition, we observe that transport of LDL to SMC surfaces is likely to be limited by the fluid phase (surface concentration less than bulk concentration), whereas ATP transport is limited by reaction on the SMC surface (surface concentration equals bulk concentration).     

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.     

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.     

Concentration wave of a solute in an artery: the influence of curvature

Mass transport and diffusion phenomena in the arterial lumen are studied through a mathematical model. Blood flow is described by the unsteady Navier-Stokes equation and solute dynamics by an advection-diffusion equation, the convective field being provided by the fluid velocity. A linearization procedure over the steady state solution is carried out and an asymptotic analysis is used to study the effect of a small curvature with respect to the straight tube. Analytical and numerical solutions are found: the results show the characteristics of the long wave propagation and the role played by the geometry on the solute distribution and demonstrate the strong influence of curvature induced by the fluid dynamics.     

Interstitial flow through the internal elastic lamina affects shear stress on arterial smooth muscle cells

Interstitial flow through the tunica media of an artery wall in the presence of the internal elastic lamina (IEL), which separates it from the subendothelial intima, has been studied numerically. A two-dimensional analysis applying the Brinkman model as the governing equation for the porous media flow field was performed. In the numerical simulation, the IEL was modeled as an impermeable barrier to water flux, except for the fenestral pores, which were uniformly distributed over the IEL. The tunica media was modeled as a heterogeneous medium composed of a periodic array of cylindrical smooth muscle cells (SMCs) embedded in a fiber matrix simulating the interstitial proteoglycan and collagen fibers. A series of calculations was conducted by varying the physical parameters describing the problem: the area fraction of the fenestral pore (0. 001-0.036), the diameter of the fenestral pore (0.4-4.0 microm), and the distance between the IEL and the nearest SMC (0.2-0.8 microm). The results indicate that the value of the average shear stress around the circumference of the SMC in the immediate vicinity of the fenestral pore could be as much as 100 times greater than that around an SMC in the fully developed interstitial flow region away from the IEL. These high shear stresses can affect SMC physiological function.     

A steady-state filtration model for transluminal water movement in small and large blood vessels

It is now generally accepted that the intercellular cleft between adjacent endothelial cells is the primary pathway for the transluminal movement of water and small ions in the vasculature. A steady-state theoretical model has been developed to show quantitatively how the geometry of the intercellular cleft between adjacent endothelial cells is related to both the water movement and pressure distribution in the subendothelial space and to examine how the existence of a subendothelial interaction layer affects the hydraulic resistance of the media of vessels of varying wall thickness. The velocity and pressure fields in the media are described using porous matrix theory based on Darcy's law and a lubrication-type analysis is used to describe the flow in a variable geometry intercellular cleft. These two equations are solved simultaneously to determine the unknown pressure distribution beneath the endothelium and the flow in the arterial media. Application of this model shows that, when the tight junction in the cleft is 26 A or less, more than half of the total hydraulic resistance of the wall occurs across the endothelial cell monolayer, for a vessel whose wall thickness is less than 0.02 cm. This finding is in good agreement with the experimental findings of Vargas, et al. (1978) for rabbit aorta. Contrary to previous belief, the model shows that the filtration resistance of an arterial wall with intact endothelium does not scale linearly with wall thickness due to the highly nonlinear resistance of the endothelial interaction layer.     

A multiphysics approach for modeling early atherosclerosis

This work is devoted to the development of a mathematical model of the early stages of atherosclerosis incorporating processes of all time scales of the disease and to show their interactions. The cardiovascular mechanics is modeled by a fluid–structure interaction approach coupling a non-Newtonian fluid to a hyperelastic solid undergoing anisotropic growth and a change of its constitutive equation. Additionally, the transport of low-density lipoproteins and its penetration through the endothelium is considered by a coupled set of advection–diffusion-reaction equations. Thereby, the permeability of the endothelium is wall-shear stress modulated resulting in a locally varying accumulation of foam cells triggering a novel growth and remodeling formulation. The model is calibrated and applied to an murine-specific case study, and a qualitative validation of the computational results is performed. The model is utilized to further investigate the influence of the pulsatile blood flow and the compliance of the artery wall to the atherosclerotic process. The computational results imply that the pulsatile blood flow is crucial, whereas the compliance of the aorta has only a minor influence on atherosclerosis. Further, it is shown that the novel model is capable to produce a narrowing of the vessel lumen inducing an adaption of the endothelial permeability pattern.     

Computational modeling of the arterial wall based on layer-specific histological data

Arterial walls typically have a heterogeneous structure with three different layers (intima, media, and adventitia). Each layer can be modeled as a fiber-reinforced material with two families of relatively stiff collagenous fibers symmetrically arranged within an isotropic soft ground matrix. In this paper, we present two different modeling approaches, the embedded fiber (EF) approach and the angular integration (AI) approach, to simulate the anisotropic behavior of individual arterial wall layers involving layer-specific data. The EF approach directly incorporates the microscopic arrangement of fibers that are synthetically generated from a random walk algorithm and captures material anisotropy at the element level of the finite element formulation. The AI approach smears fibers in the ground matrix and treats the material as homogeneous, with material anisotropy introduced at the constitutive level by enhancing the isotropic strain energy with two anisotropic terms. Both approaches include the influence of fiber dispersion introduced by fiber angular distribution (departure of individual fibers from the mean orientation) and take into consideration the dispersion caused by fiber waviness, which has not been previously considered. By comparing the numerical results with the published experimental data of different layers of a human aorta, we show that by using histological data both approaches can successfully capture the anisotropic behavior of individual arterial wall layers. Furthermore, through a comprehensive parametric study, we establish the connections between the AI phenomenological material parameters and the EF parameters having straightforward physical or geometrical interpretations. This study provides valuable insight for the calibration of phenomenological parameters used in the homogenized modeling based on the fiber microscopic arrangement. Moreover, it facilitates a better understanding of individual arterial wall layers, which will eventually advance the study of the structure–function relationship of arterial walls as a whole.     

Numerical analysis of the cooling effect of blood over inflamed atherosclerotic plaque

Atherosclerotic plaques with high likelihood of rupture often show local temperature increase with respect to the surrounding arterial wall temperature. In this work, atherosclerotic plaque temperature was numerically determined during the different levels of blood flow reduction produced by the introduction of catheters at the vessel lumen. The temperature was calculated by solving the energy equation and the Navier-Stokes equations in 2D idealized arterial models. Arterial wall temperature depends on three basic factors: metabolic activity of the inflammatory cells embedded in the plaque, heat convection due to luminal blood flow, and heat conduction through the arterial wall and plaque. The calculations performed serve to simulate transient blood flow reduction produced by the presence of thermography catheters used to measure arterial wall temperature. The calculations estimate the spatial and temporal alterations in the cooling effect of blood flow and plaque temperature during the measurement process. The mathematical model developed provides a tool for analyzing the contribution of factors known to affect heat transfer at the plaque surface. Blood flow reduction leads to a nonuniform temperature increase ranging from 0.1 to 0.25 degrees Celsius in the plaque/lumen interface of the arterial geometries considered in this study. The temperature variation as well as the Nusselt number calculated along the plaque surface strongly depended on the arterial geometry and distribution of inflammatory cells. The calculations indicate that the minimum required time to obtain a steady temperature profile after arterial occlusion is 6 s. It was seen that in arteries with geometries involving bends, the temperature profiles appear asymmetrical and lean toward the downstream edge of the plaque.     

Multiscale bio-chemo-mechanical model of intimal hyperplasia

We consider a computational multiscale framework of a bio-chemo-mechanical model for intimal hyperplasia. With respect to existing models, we investigate the interactions between hemodynamics, cellular dynamics and biochemistry on the development of the pathology. Within the arterial wall, we propose a mathematical model consisting of kinetic differential equations for key vascular cell types, collagen and growth factors. The luminal hemodynamics is modeled with the Navier–Stokes equations. Coupling hypothesis among time and space scales are proposed to build a tractable modeling of such a complex multifactorial and multiscale pathology. A one-dimensional numerical test-case is presented for validation by comparing the results of the framework with experiments at short and long timescales. Our model permits to capture many cellular phenomena which have a central role in the physiopathology of intimal hyperplasia. Results are quantitatively and qualitatively consistent with experimental findings at both short and long timescales.

     

Non-standard numerical methods applied to subsurface biobarrier formation models in porous media

Biofilm forming microbes have complex effects on the flow properties of natural porous media. Subsurface biofilms have the potential for the formation of biobarriers to inhibit contaminant migration in groundwater. Another example of beneficial microbial effects is the biotransformation of organic contaminants to less harmful forms, thereby providing an in situ method for treatment of contaminated groundwater supplies. Mathematical models that describe contaminant transport with biodegradation involve a set of coupled convection-dispersion equations with non-linear reactions. The reactive solute transport equation is one for which numerical solution procedures continue to exhibit significant limitations for certain problems of groundwater hydrology interest. Accurate numerical simulations are crucial to the development of contaminant remediation strategies. A new numerical method is developed for simulation of reactive bacterial transport in porous media. The non-standard numerical approach is based on the ideas of the ‘exact’ time-stepping scheme. It leads to solutions free from the numerical instabilities that arise from incorrect modeling of derivatives and reaction terms. Applications to different biofilm models are examined and numerical results are presented to demonstrate the performance of the proposed new method.     

Numerical modelling of mass transport in an arterial wall with anisotropic transport properties

Coronary artery disease results in blockages or narrowing of the artery lumen. Drug eluting stents (DES) were developed to replace bare metal stents in an effort to combat re-blocking of the diseased artery following treatment. The numerical models developed within this study focus on representing the changing trends of drug delivery from an idealised DES in an arterial wall with an anisotropic ultra-structure. To reduce the computational burden of solving coupled physics problems, a model reduction strategy was adopted. Particular focus has been placed upon adequately modelling the influence of strut compression as there is a paucity of numerical studies that account for changes in transport properties in compressed regions of the arterial wall due to stent deployment. This study developed an idealised numerical modelling framework to account for the changes in the directionally dependent porosity and tortuosities of the arterial wall as a result of radial strut compression. The results show that depending on the degree of strut compression, trends in therapeutic drug delivery within the arterial wall can be either increased or decreased. The study highlights the importance of incorporating compression into numerical models to better represent transport within the arterial wall and suggests an appropriate numerical modelling framework that could be utilised in more realistic patient specific arterial geometries.     

Dependence of leukocyte capture on instantaneous pulsatile flow

Atherosclerosis, an artery disease, is currently the leading cause of death in the United States in both men and women. The first step in the development of atherosclerosis involves leukocyte adhesion to the arterial endothelium. It is broadly accepted that blood flow, more specifically wall shear stress (WSS), plays an important role in leukocyte capture and subsequent development of an atherosclerotic plaque. What is less known is how instantaneous WSS, which can vary by up to 5 Pa over one cardiac cycle, influences leukocyte capture. In this paper we use direct numerical simulations (DNS), performed using an in-house code, to illustrate that leukocyte capture is different whether as a function of instantaneous or time-averaged blood flow. Specifically, a stenotic plaque is modeled using a computational fluid dynamics (CFD) solver through fully three-dimensional Navier-Stokes equations and the immersed boundary method. Pulsatile triphasic inflow is used to simulate the cardiac cycle. The CFD is coupled with an agent-based leukocyte capture model to assess the impact of instantaneous hemodynamics on stenosis growth. The computed wall shear stress agrees well with the results obtained with a commercial software, as well as with theoretical results in the healthy region of the artery. The analysis emphasizes the importance of the instantaneous flow conditions in evaluating the leukocyte rate of capture. That is, the capture rate computed from mean flow field is generally underpredicted compared to the actual rate of capture. Thus, in order to obtain a reliable estimate, the flow unsteadiness during a cardiac cycle should be taken into account.     

Transport of bacteria in porous media: II. A model for convective Transport and growth

A model is presented for the coupled processes of bacterial growth and convective transport of bacteria has been modeled using a fractional flow approach. The various mechanisms of bacteria retention can be incorporated into the model through selection of an appropriate shape of the fractional flow curve. Permeability reduction due to pore plugging by bacteria was simulated using the effective medium theory. In porous media, the rates of transport and growth of bacteria, the generation of metabolic products, and the consumption of nutrients are strongly coupled processes. Consequently, the set of governing conservation equations form a set of coupled, nonlinear partial differential equations that were solved numerically. Reasonably good agreement between the model and experimental data has been obtained indicating that the physical processes incorporated in the model are adequate. The model has been used to predict the in situ transport and growth of bacteria, nutrient consumption, and metabolite production. It can be particularly useful in simulating laboratory experiments and in scaling microbial-enhanced oil recovery or bioremediation processes to the field. (c) 1994 John Wiley & Sons, Inc.