Flow of couple stress fluid through stenotic blood vessels

The effects of an axially symmetric mild stenosis on the flow of blood, when blood is represented by a couple stress fluid model, have been studied. It is found that, for a fixed stenosis size, the resistance to flow and wall shear stress increase as the couple stress parameter eta decreases from unity. A comparison of the results with those of the Newtonian case shows that the magnitude of resistance to flow and wall shear under a given set of conditions, is greater in the case of the couple stress fluid model. It is seen that even in the case of a mild stenosis (19% area reduction), resistance to flow and wall shear values are increased over those for no stenosis by 60% and 62%, respectively, when compared with the case of a Newtonian fluid.     

Blood Vessel Adaptation with Fluctuations in Capillary Flow Distribution

Throughout the life of animals and human beings, blood vessel systems are continuously adapting their structures – the diameter of vessel lumina, the thickness of vessel walls, and the number of micro-vessels – to meet the changing metabolic demand of the tissue. The competition between an ever decreasing tendency of luminal diameters and an increasing stimulus from the wall shear stress plays a key role in the adaptation of luminal diameters. However, it has been shown in previous studies that the adaptation dynamics based only on these two effects is unstable. In this work, we propose a minimal adaptation model of vessel luminal diameters, in which we take into account the effects of metabolic flow regulation in addition to wall shear stresses and the decreasing tendency of luminal diameters. In particular, we study the role, in the adaptation process, of fluctuations in capillary flow distribution which is an important means of metabolic flow regulation. The fluctuation in the flow of a capillary group is idealized as a switch between two states, i.e., an open-state and a close-state. Using this model, we show that the adaptation of blood vessel system driven by wall shear stress can be efficiently stabilized when the open time ratio responds sensitively to capillary flows. As micro-vessel rarefaction is observed in our simulations with a uniformly decreased open time ratio of capillary flows, our results point to a possible origin of micro-vessel rarefaction, which is believed to induce hypertension.     

Numerical Modeling of Interstitial Fluid Flow Coupled with Blood Flow through a Remodeled Solid Tumor Microvascular Network

Modeling of interstitial fluid flow involves processes such as fluid diffusion, convective transport in extracellular matrix, and extravasation from blood vessels. To date, majority of microvascular flow modeling has been done at different levels and scales mostly on simple tumor shapes with their capillaries. However, with our proposed numerical model, more complex and realistic tumor shapes and capillary networks can be studied. Both blood flow through a capillary network, which is induced by a solid tumor, and fluid flow in tumor’s surrounding tissue are formulated. First, governing equations of angiogenesis are implemented to specify the different domains for the network and interstitium. Then, governing equations for flow modeling are introduced for different domains. The conservation laws for mass and momentum (including continuity equation, Darcy’s law for tissue, and simplified Navier–Stokes equation for blood flow through capillaries) are used for simulating interstitial and intravascular flows and Starling’s law is used for closing this system of equations and coupling the intravascular and extravascular flows. This is the first study of flow modeling in solid tumors to naturalistically couple intravascular and extravascular flow through a network. This network is generated by sprouting angiogenesis and consisting of one parent vessel connected to the network while taking into account the non-continuous behavior of blood, adaptability of capillary diameter to hemodynamics and metabolic stimuli, non-Newtonian blood flow, and phase separation of blood flow in capillary bifurcation. The incorporation of the outlined components beyond the previous models provides a more realistic prediction of interstitial fluid flow pattern in solid tumors and surrounding tissues. Results predict higher interstitial pressure, almost two times, for realistic model compared to the simplified model.     

A model of oxygen exchange between an arteriole or venule and the surrounding tissue

A mathematical model is developed to study the effect of capillary convection on oxygen transport around segments of arterioles and venules that are surrounded by capillaries. These capillaries carry unidirectional flow perpendicular to the vessel. The discrete capillary structure is distributed in a manner determined by the capillary blood flow and capillary density. A nonlinear oxyhemoglobin dissociation curve described by the Hill equation is used in the analysis. Oxygen flux from the vessel is expressed as a relationship between Sherwood and Peclet numbers, as well as other dimensionless combinations involving parameters of the capillary bed. A numerical solution is obtained with a finite difference method. The numerical results obtained within the physiological range of parameters allow the prediction of longitudinal gradients of hemoglobin-oxygen saturation along the arterioles and venules.     

Particulate nature of blood determines macroscopic rheology: a 2-D lattice Boltzmann analysis

Historically, predicting macroscopic blood flow characteristics such as viscosity has been an empirical process due to the difficulty in rigorously including the particulate nature of blood in a mathematical representation of blood rheology. Using a two-dimensional lattice Boltzmann approach, we have simulated the flow of red blood cells in a blood vessel to estimate flow resistance at various hematocrits and vessel diameters. By including white blood cells (WBCs) in the flow, we also calculate the increase in resistance due to white cell rolling and adhesion. The model considers the blood as a suspension of particles in plasma, accounting for cell-cell and cell-wall interactions to predict macroscopic blood rheology. The model is able to reproduce the Fahraeus-Lindqvist effect, i.e., the increase in relative apparent viscosity as tube size increases, and the Fahraeus effect, i.e., tube hematocrit is lower than discharge hematocrit. In addition, the model allows direct assessment of the effect of WBCs on blood flow in the microvasculature, reproducing the dramatic increases in flow resistance as WBCs enter short capillary segments. This powerful and flexible model can be used to predict blood flow properties in any vessel geometry and with any blood composition.     

Atherosclerosis and flow in carotid arteries with authentic geometries

The influence of blood flow on the depositions and development of atherosclerotic lesions have been observed and described since the 19th century. Observations have shown that depositions correlate with regions of low wall shear stress. However, the exact correlations between depositions, vessel geometry and flow parameters are not yet known. The purpose of this study was the quantification of atherosclerosis risk factors in carotid bifurcation. This artery has attracted particular interest because lesions are often found in this bifurcation. Post mortem, the arteries are excised and vessel casts are produced. Afterwards, the arteries are analyzed morphometrically. The vessel casts are used for the assessment of some geometrical parameters. 31 carotid bifurcations were analyzed in this study. Eight vessel casts were digitized and rendered three-dimensional mathematical models of the arteries. These data were imported by the computational fluid dynamics program FLUENT. Further, the blood flow was reconstructed in a computer model based on the individual vessel geometry. The flow parameters, such as velocity, pressure and wall shear stress were computed. At the same time the geometrical parameters and wall alterations are known. This permits the comparison of the anatomical shape and its flow with the distribution and level of the wall alterations.     

A numerical model for studying the effect of plasma layer on the process of blood oxygenation in the pulmonary capillaries

A two layer model for the blood oxygenation in pulmonary capillaries is proposed. The model consists of a core of erythrocytes surrounded by a symmetrically placed plasma layer. The governing equations in the core describe the free molecular diffusion, convection, and facilitated diffusion due to the presence of haemoglobin. The corresponding equations in the plasma layer are based on the free molecular diffusion and the convective effect of the blood. According to the axial train model for the blood flow proposed by Whitmore (1967), the core will move with a uniform velocity whereas flow in the plasma layer will be fully developed. The resulting system of nonlinear partial differential equations is solved numerically. A fixed point iterative technique is used to deal with the nonlinearities. The distance traversed by the blood before getting fully oxygenated is computed. It is shown that the concentration of O2 increases continuously along the length of the capillary for a given ratio of core radius to capillary radius. It is found that the rate of oxygenation increases as the core to capillary ratio decreases. The equilibration length increases with a heterogeneous model in comparison to that in a homogeneous model. The effect of capillary diameters and core radii on the rate of oxygenation has also been examined.     

Radial artery wall shear stress evaluation in patients with arteriovenous fistula for hemodialysis access

It has been extensively documented that changes in blood flow induce vascular remodeling and this phenomenon seems to be correlated to the shear forces imposed on the vessel wall by motion of blood. Wall shear stress, the tractive force that acts on the endothelium, has been shown to influence endothelial cell function. To study changes in wall shear stress that develop on the vessel wall upon changes of blood flow, we set up a technique that allows estimation of shear stress in the radial artery of patients on chronic hemodialysis therapy. The technique is based on color-flow Doppler examination of the radial artery before and after surgical creation of radiocephalic fistula for hemodialysis. Calculation of time function wall shear stress and blood flow rate in the radial artery is performed on the basis of arterial diameter, center-line velocity waveform and blood viscosity, using a numerical method developed according to Womersley's theory for pulsatile flow in tubes. The results presented confirm that the model developed is suitable for calculation of the wall shear stress that develops in the radial artery of patients before and after surgical creation of an arteriovenous fistula for hemodialysis. This methodology was developed for characterization of wall shear stress in the radial artery but may be well applied to other vessels that can be examined by echo-Doppler technique.     

Effects of couple stresses on the blood flow through an artery with mild stenosis

An analysis of the effects of couple stresses on the blood flow through thin artery in the presence of very mild stenosis has been carried out with the help of two nondimensional parameters, alpha (the length ratio parameter) and eta (the parameter characterizing the antisymmetric property of the couple stress tensor). It is shown that an increase in the couple stress (small value of alpha and eta), increases the resistance to the flow and the wall shear stress. These characteristics are further enhanced by the presence of the stenosis.     

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.     

Structural response of microcirculatory networks to changes in demand: information transfer by shear stress

Matching blood flow to metabolic demand in terminal vascular beds involves coordinated changes in diameters of vessels along flow pathways, requiring upstream and downstream transfer of information on local conditions. Here, the role of information transfer mechanisms in structural adaptation of microvascular networks after a small change in capillary oxygen demand was studied using a theoretical model. The model includes diameter adaptation and information transfer via vascular reactions to wall shear stress, transmural pressure, and oxygen levels. Information transfer is additionally effected by conduction along vessel walls and by convection of metabolites. The model permits selective blocking of information transfer mechanisms. Six networks, based on in vivo data, were considered. With information transfer, increases in network conductance and capillary oxygen supply were amplified by factors of 4.9 +/- 0.2 and 9.4 +/- 1.1 (means +/- SE), relative to increases when information transfer was blocked. Information transfer by flow coupling alone, in which increased shear stress triggers vascular enlargement, gave amplifications of 4.0 +/- 0.3 and 4.9 +/- 0.5. Other information transfer mechanisms acting alone gave amplifications below 1.6. Thus shear-stress-mediated flow coupling is the main mechanism for the structural adjustment of feeding and draining vessel diameters to small changes in capillary oxygen demand.     

Inert gas clearance from tissue by co-currently and counter-currently arranged microvessels

To elucidate the clearance of dissolved inert gas from tissues, we have developed numerical models of gas transport in a cylindrical block of tissue supplied by one or two capillaries. With two capillaries, attention is given to the effects of co-current and counter-current flow on tissue gas clearance. Clearance by counter-current flow is compared with clearance by a single capillary or by two co-currently arranged capillaries. Effects of the blood velocity, solubility, and diffusivity of the gas in the tissue are investigated using parameters with physiological values. It is found that under the conditions investigated, almost identical clearances are achieved by a single capillary as by a co-current pair when the total flow per tissue volume in each unit is the same (i.e., flow velocity in the single capillary is twice that in each co-current vessel). For both co-current and counter-current arrangements, approximate linear relations exist between the tissue gas clearance rate and tissue blood perfusion rate. However, the counter-current arrangement of capillaries results in less-efficient clearance of the inert gas from tissues. Furthermore, this difference in efficiency increases at higher blood flow rates. At a given blood flow, the simple conduction-capacitance model, which has been used to estimate tissue blood perfusion rate from inert gas clearance, underestimates gas clearance rates predicted by the numerical models for single vessel or for two vessels with co-current flow. This difference is accounted for in discussion, which also considers the choice of parameters and possible effects of microvascular architecture on the interpretation of tissue inert gas clearance.     

Hemodynamics in the carotid artery bifurcation: a comparison between numerical simulations and in vitro MRI measurements

The presence of atherosclerotic plaques has been shown to be closely related to the vessel geometry. Studies on postmortem human arteries and on the experimental animal show positive correlation between the presence of plaque thickness and low shear stress, departure of unidirectional flow and regions of flow separation and recirculation. Numerical simulations of arterial blood flow and direct blood flow velocity measurements by magnetic resonance imaging (MRI) are two approaches for the assessment of arterial blood flow patterns. In order to verify that both approaches give equivalent results magnetic resonance velocity data measured in a compliant anatomical carotid bifurcation model were compared to the results of numerical simulations performed for a corresponding computational vessel model. Cross sectional axial velocity profiles were calculated and measured for the midsinus and endsinus internal carotid artery. At both locations a skewed velocity profile with slow velocities at the outer vessel wall, medium velocities at the side walls and high velocities at the flow divider (inner) wall were observed. Qualitative comparison of the axial velocity patterns revealed no significant differences between simulations and in vitro measurements. Even quantitative differences such as for axial peak flow velocities were less than 10%. Secondary flow patterns revealed some minor differences concerning the form of the vortices but maximum circumferential velocities were in the same range for both methods.     

Experimental assessment of wall shear flow in models

The blood flow immediately adjacent to the wall of a blood vessel or an artificial surface is of great interest. This flow defines the shear stress at the wall and is known to have a great physiological importance. The use of models is a viable method to investigate this flow. However, even in models the shear stress at the wall is difficult to assess. A new optical method is based on transparent models and uses particles in the model fluid, which are only visible near the wall. This is achieved with a model fluid having a defined opacity. This fluid obscures particles in the center of the models, but permits the observation and recording of particles close to the wall. The method has been applied for Hagen-Poiseuille flow and for the likewise well researched flow in a tube with a sudden expansion.     

Modeling of arterial stenosis and its applications to blood diseases

Blood flow in a stenosed tube has been modeled in the present studies. Blood flow is assumed to be represented by a couple stress fluid. Flow parameters such as velocity, resistance to flow, and shear stress distribution have been computed for different suspension concentrations (haematocrit), and for the blood diseases; polycythemia, plasma cell dyscrasias, and for Hb SS (sickle cell). The results have been compared with the case of normal blood and for other theoretical models. The importance of size effects in blood flow studies has been highlighted.     

Particle-fluid suspension model of blood flow through stenotic vessels with applications

The present study deals with the problem of blood flow through stenotic vessels when blood is represented by a particle-fluid suspension model, i.e. a suspension of red blood cells in plasma. The expressions for the dimensionless resistance to flow, the wall shear stress, and the shearing stress on the wall at the maximum height of the stenosis are derived. The results obtained in the analysis are discussed in brief, both qualitatively and quantitatively by comparison with other theories. It is observed that the magnitudes of the blood flow characteristics significantly increase with an increase in the red cell concentration. The importance of the decreasing vessel diameter is also pointed out. Finally, to observe the biological relevance of the analysis, the results obtained are used to compute the blood flow characteristics for normal and diseased blood using the experimental data from published literature and results are compared with those computed using the present theoretical approach.     

A viscoelastic model of blood capillary extension and regression: derivation, analysis, and simulation

This work studies a fundamental problem in blood capillary growth: how the cell proliferation or death induces the stress response and the capillary extension or regression. We develop a one-dimensional viscoelastic model of blood capillary extension/regression under nonlinear friction with surroundings, analyze its solution properties, and simulate various growth patterns in angiogenesis. The mathematical model treats the cell density as the growth pressure eliciting a viscoelastic response from the cells, which again induces extension or regression of the capillary. Nonlinear analysis captures two cases when the biologically meaningful solution exists: (1) the cell density decreases from root to tip, which may occur in vessel regression; (2) the cell density is time-independent and is of small variation along the capillary, which may occur in capillary extension without proliferation. The linear analysis with perturbation in cell density due to proliferation or death predicts the global biological solution exists provided the change in cell density is sufficiently slow in time. Examples with blow-ups are captured by numerical approximations and the global solutions are recovered by slow growth processes, which validate the linear analysis theory. Numerical simulations demonstrate this model can reproduce angiogenesis experiments under several biological conditions including blood vessel extension without proliferation and blood vessel regression.     

The influence of fluid shear stress on the remodeling of the embryonic primary capillary plexus

The primary capillary plexus in early yolk sacs is remodeled into matured vitelline vessels aligned in the direction of blood flow at the onset of cardiac contraction. We hypothesized that the influence of fluid shear stress on cellular behaviors may be an underlying mechanism by which some existing capillary channels remain open while others are closed during remodeling. Using a recently developed E-Tmod knock-out/lacZ knock-in mouse model, we showed that erythroblasts exhibited rheological properties similar to those of a viscous cell suspension. In contrast, the non-erythroblast (NE) cells, which attach among themselves within the yolk sac, are capable of lamellipodia extension and cell migration. Isolated NE cells in a parallel-plate flow chamber exposed to fluid shear stress, however, ceased lamellipodia extension. Such response may minimize NE cell migration into domains exposed to fluid shear stress. A two-dimensional mathematical model incorporating these cellular behaviors demonstrated that shear stress created by the blood flow initiated by the embryonic heart contraction might be needed for the remodeling of primary capillary plexus.     

Functional hierarchy of coronary circulation: direct evidence of a structure-function relation

The heart muscle is nourished by a complex system of blood vessels that make up the coronary circulation. Here we show that the design of the coronary circulation has a functional hierarchy. A full anatomic model of the coronary arterial tree, containing millions of blood vessels down to the capillary vessels, was simulated based on previously measured porcine morphometric data. A network analysis of blood flow through every vessel segment was carried out based on the laws of fluid mechanics and appropriate boundary conditions. Our results show an abrupt change in cross-sectional area that demarcates the transition from epicardial (EPCA) to intramyocardial (IMCA) coronary arteries. Furthermore, a similar pattern of blood flow was observed with a corresponding transition from EPCA to IMCA. These results suggest functional differences between the two types of vessels. An additional abrupt change occurs in the IMCA in relation to flow velocity. The velocity is fairly uniform proximal to these vessels but drops significantly distal to those vessels toward the capillary branches. This finding suggests functional differences between large and small IMCA. Collectively, these observations suggest a novel functional hierarchy of the coronary vascular tree and provide direct evidence of a structure-function relation.     

Numerical investigation of the non-Newtonian pulsatile blood flow in a bifurcation model with a non-planar branch

The pulsatile flow of non-Newtonian fluid in a bifurcation model with a non-planar daughter branch is investigated numerically by using the Carreau-Yasuda model to take into account the shear thinning behavior of the analog blood fluid. The objective of this study is to deal with the influence of the non-Newtonian property of fluid and of out-of-plane curvature in the non-planar daughter vessel on wall shear stress (WSS), oscillatory shear index (OSI), and flow phenomena during the pulse cycle. The non-Newtonian property in the daughter vessels induces a flattened axial velocity profile due to its shear thinning behavior. The non-planarity deflects flow from the inner wall of the vessel to the outer wall and changes the distribution of WSS along the vessel, in particular in systole phase. Downstream of the bifurcation, the velocity profiles are shifted toward the flow divider, and low WSS and high shear stress temporal oscillations characterized by OSI occur on the outer wall region of the daughter vessels close to the bifurcation. Secondary motions become stronger with the addition of the out-of-plane curvature induced by the bending of the vessel, and the secondary flow patterns swirl along the non-planar daughter vessel. A significant difference between the non-Newtonian and the Newtonian pulsatile flow is revealed during the pulse cycle; however, reasonable agreement between the non-Newtonian and the rescaled Newtonian flow is found. Calculated results for the pulsatile flow support the view that the non-planarity of blood vessels and the non-Newtonian properties of blood are an important factor in hemodynamics and may play a significant role in vascular biology and pathophysiology.