A problem in the mathematical biophysics of blood circulation: II. Relation between pressure and flow of a viscous fluid in an elastic distensible tube

In an elastic distensible tube, like a blood vessel, the radius is determined by the equality of the hydrostatic pressure and the elastic forces. If a viscous fluid flows through such a tube, there is a pressure drop along the line of flow. This results in a variation of the radius of the tube along the axis. An approximate expression, valid within a limited range of values, is derived for the radius of the tube as a function of the distance along the axis. Another approximate expression is derived for the relation between pressure drop and total flow in such a case. For sufficiently high rates of flow the pressure drop does not vary linearly with the flow, as in the usual poiseuille's law, but more rapidly.     

Heterogeneous mechanics of the mouse pulmonary arterial network

Individualized modeling and simulation of blood flow mechanics find applications in both animal research and patient care. Individual animal or patient models for blood vessel mechanics are based on combining measured vascular geometry with a fluid structure model coupling formulations describing dynamics of the fluid and mechanics of the wall. For example, one-dimensional fluid flow modeling requires a constitutive law relating vessel cross-sectional deformation to pressure in the lumen. To investigate means of identifying appropriate constitutive relationships, an automated segmentation algorithm was applied to micro-computerized tomography images from a mouse lung obtained at four different static pressures to identify the static pressure–radius relationship for four generations of vessels in the pulmonary arterial network. A shape-fitting function was parameterized for each vessel in the network to characterize the nonlinear and heterogeneous nature of vessel distensibility in the pulmonary arteries. These data on morphometric and mechanical properties were used to simulate pressure and flow velocity propagation in the network using one-dimensional representations of fluid and vessel wall mechanics. Moreover, wave intensity analysis was used to study effects of wall mechanics on generation and propagation of pressure wave reflections. Simulations were conducted to investigate the role of linear versus nonlinear formulations of wall elasticity and homogeneous versus heterogeneous treatments of vessel wall properties. Accounting for heterogeneity, by parameterizing the pressure/distention equation of state individually for each vessel segment, was found to have little effect on the predicted pressure profiles and wave propagation compared to a homogeneous parameterization based on average behavior. However, substantially different results were obtained using a linear elastic thin-shell model than were obtained using a nonlinear model that has a more physiologically realistic pressure versus radius relationship.     

A framework for incorporating 3D hyperelastic vascular wall models in 1D blood flow simulations

We present a novel framework for investigating the role of vascular structure on arterial haemodynamics in large vessels, with a special focus on the human common carotid artery (CCA). The analysis is carried out by adopting a three-dimensional (3D) derived, fibre-reinforced, hyperelastic structural model, which is coupled with an axisymmetric, reduced order model describing blood flow. The vessel transmural pressure and lumen area are related via a Holzapfel–Ogden type of law, and the residual stresses along the thickness and length of the vessel are also accounted for. After a structural characterization of the adopted hyperelastic model, we investigate the link underlying the vascular wall response and blood-flow dynamics by comparing the proposed framework results against a popular tube law. The comparison shows that the behaviour of the model can be captured by the simpler linear surrogate only if a representative value of compliance is applied. Sobol’s multi-variable sensitivity analysis is then carried out in order to identify the extent to which the structural parameters have an impact on the CCA haemodynamics. In this case, the local pulse wave velocity (PWV) is used as index for representing the arterial transmission capacity of blood pressure waveforms. The sensitivity analysis suggests that some geometrical factors, such as the stress-free inner radius and opening angle, play a major role on the system’s haemodynamics. Subsequently, we quantified the differences in haemodynamic variables obtained from different virtual CCAs, tube laws and flow conditions. Although each artery presents a distinct vascular response, the differences obtained across different flow regimes are not significant. As expected, the linear tube law is unable to accurately capture all the haemodynamic features characterizing the current model. The findings from the sensitivity analysis are further confirmed by investigating the axial stretching effect on the CCA fluid dynamics. This factor does not seem to alter the pressure and flow waveforms. On the contrary, it is shown that, for an axially stretched vessel, the vascular wall exhibits an attenuation in absolute distension and an increase in circumferential stress, corroborating the findings of previous studies. This analysis shows that the new model offers a good balance between computational complexity and physics captured, making it an ideal framework for studies aiming to investigate the profound link between vascular mechanobiology and blood flow.

     

Pulsatile blood flow, shear force, energy dissipation and Murray's Law

Background  

Murray's Law states that, when a parent blood vessel branches into daughter vessels, the cube of the radius of the parent vessel is equal to the sum of the cubes of the radii of daughter blood vessels. Murray derived this law by defining a cost function that is the sum of the energy cost of the blood in a vessel and the energy cost of pumping blood through the vessel. The cost is minimized when vessel radii are consistent with Murray's Law. This law has also been derived from the hypothesis that the shear force of moving blood on the inner walls of vessels is constant throughout the vascular system. However, this derivation, like Murray's earlier derivation, is based on the assumption of constant blood flow.     

Cell-free plasma layer in cerebral microvessels.

Two diameters of vessel and red cell column in cerebral microvessels (> 29.8 microns in diameter) of cat were measured together with red cell velocity, using a two fluorescent tracer method. A fluorescein isothiocyanate (FITC)-labeled red cell was adopted as a flow tracer to measure the cell velocity with a dual window technique. Based on the fluorescence image, the red cell column diameter was measured. Plasma was stained with rhodamine-B isothiocyanate (RITC)-labeled dextran to measure the vessel diameter. The thickness of the cell-free plasma layer could be determined from the difference of the two diameters. The obtained thickness of the cell-free layer was not described by a simple function of vessel diameter or red cell velocity; it was dependent on the pseudo shear rate defined by the ratio of cell velocity to vessel radius. The layer thickness increased with a decrease in the pseudo shear rate.     

Experimental determination of wall shear rate in canine carotid arteries perfused in vitro

The mathematical model of Hung (Tsai and Hung, 1984) is employed to determine the wall shear rate acting on canine carotid arteries perfused in vitro. Model equations for pulsatile flow in a deformable vessel are coupled with experimental data of dynamic pressure drop, flow rate, vessel radius and radial wall motion. Derived quantities, e.g. velocity profiles and wall shear, are obtained for vessels exposed to 'normotensive' hemodynamics, 'hypertension' simulations and perfusions in which the compliance of the vessel wall is deliberately altered. Our results indicate that wall shear varies markedly as a function of the hemodynamic environment. The effects of vessel radius vs flow rate on the development of wall shear are also demonstrated. It is found that convective processes correlate with the magnitude of wall shear in the 'hypertension' simulations. The present findings and complementary published data may explain, at least in part, the variations in vessel wall transport and endothelial cell biology we observe as a function of the hemodynamic environment. For example we have documented that the exposure of canine carotids to 'hypertensive' (vs 'normotensive') hemodynamics is associated with an increased flux of lipoproteins (LDL) into the intima and luminal media. Alternations in wall compliance, on the other hand, profoundly influence endothelial shape, orientation and cytoskeletal array.     

A model of pulsatile flow in a uniform deformable vessel.

Simulations of blood flow in natural and artificial conduits usually require large computers for numerical solution of the Navier-Stokes equations. Often, physical insight into the fluid dynamics is lost when the solution is purely numerical. An alternative to solving the most general form of the Navier-Stokes equations is described here, wherein a functional form of the solution is assumed in order to simplify the required computations. The assumed forms for the axial pressure gradient and velocity profile are chosen such that conservation of mass is satisfied for fully established pulsatile flow in a straight, deformable vessel. The resulting equations are cast in finite-difference form and solved explicitly. Results for the limiting cases of rigid wall and zero applied pressure are found to be in good agreement with analytical solutions. Comparison with the experimental results of Klanchar et al. [Circ. Res. 66, 1624-1635 (1990]) also shows good agreement. Application of the model to realistic physiological parameter values provides insight as to the influence of the pulsatile nature of the flow field on wall shear development in the presence of a moving wall boundary. Specifically, the model illustrates the dependence of flow rate and shear rate on the amplitude of the vessel wall motion and the phase difference between the applied pressure difference and the oscillations of the vessel radius. The present model can serve as a useful tool for experimentalists interested in quantifying the magnitude and character of velocity profiles and shearing forces in natural and artificial biologic conduits.     

Factors influencing distribution of cerebral gas embolism

The results of a series of experiments in baboons centrally cooled to a brain temperature of 18 °C have been described. These studies were characterized by 50% reduction in blood flow to cerebellum after cooling, relative increase in distribution to cerebellar cortex as opposed to white matter, and reduced cerebral oxygen extraction.The physical factors of importance are the viscosity of the blood, chiefly the effect to hematocrit, surface tension, the surface tension coefficient between gas and blood, the radius of the vessel, aud the velocity of the bubble.     

Experimental determination of wall shear rate in canine carotid arteries perfused in vitro

The mathematical model of Hung (Tsai and Hung, 1984) is empolyed to determine the wall shear rate acting on canine carotid arteries perfused in vitro. Model equations for pulsatile flow in a deformable vessel are coupled with experimental data of dynamic pressure drop, flow rate, vessel radius and radial wall motion. Derived quantities, e.g. velocity profiles and wall shear, are obtained for vessels exposed to ‘normotensive’ hemodynamics, ‘hypertension’ simulations and perfusions in which the compliance of the vessel wall is deliberately altered. Our results indicate that wall shear varies markedly as a function of the hemodynamic environment. The effects of vessel radius vs flow rate on the development of wall shear are also demonstrated. It is found that convective processes correlate with the magnitude of wall shear in the ‘hypertension’ simulations.The present findings and complementary published data may explain, at least in part, the variations in vessel wall transport and endothelial cell biology we observe as a function of the hemodynamic environment. For example we have documented that the exposure of canine carotids to ‘hypertensive’ (vs ‘normotensive’) hemodynamics is associated with an increased flux of lipoproteins (LDL) into the intima and luminal media. Alternations in wall compliance, on the other hand, profoundly influence endothelial shape, orientation and cytoskeletal array.     

Peristaltic biofluids flow through vertical porous human vessels using third-grade non-Newtonian fluids model

In this paper, the heat and flow characteristic of third-grade non-Newtonian biofluids flow through a vertical porous human vessel due to peristaltic wall motion are studied. The third-grade model can describe shear thinning (or shear thickening) and normal stress differences, which is acceptable for biofluids modeling. In order to solve the governing equations, the assumption of long-wavelength approximation is utilized. This hypothesis emphasizes that the wavelength of the peristaltic wall motion is large in comparison with the radius of the human vessel, which is widely acceptable in biological investigations. The analytical perturbation method is employed to solve the governing equations. Consequently, analytical expressions for the velocity profile, shear stress, temperature field, and biofluid flow rate are obtained. In addition, the effects of the governing parameters such as the third-grade non-Newtonian parameter, Grashof Number, Eckert number, and porosity, on the results are examined.     

Locations of ectopic beats coincide with spatial gradients of NADH in a regional model of low-flow reperfusion

We studied the origins of ectopic beats during low-flow reperfusion after acute regional ischemia in excised rat hearts. The left anterior descending coronary artery was cannulated. Perfusate was delivered to the cannula using an high-performance liquid chromatography pump. This provided not only precise control of flow rate but also avoided mechanical artifacts associated with vessel occlusion and deocclusion. Optical mapping of epicardial transmembrane potential served to identify activation wavefronts. Imaging of NADH fluorescence was used to quantify local ischemia. Our experiments suggest that low-flow reperfusion of ischemic myocardium leads to a highly heterogeneous ischemic substrate and that the degree of ischemia between adjacent patches of tissue changes in time. In contrast to transient ectopic activity observed during full-flow reperfusion, persistent ectopic arrhythmias were observed during low-flow reperfusion. The origins of ectopic beats were traceable to areas of high spatial gradients of changes in NADH fluorescence caused by low-flow reperfusion.     

Plug Effect of Erythrocytes in Capillary Blood Vessels

As an idealized problem of the motion of blood in small capillary blood vessels, the low Reynolds number flow of plasma (a newtonian fluid) in a circular cylindrical tube involving a series of circular disks is studied. It is assumed in this study that the suspended disks are equally spaced along the axis of the tube, and that their centers remain on the axis of the tube and that their faces are perpendicular to the tube axis. The inertial force of the fluid due to the convective acceleration is neglected on the basis of the smallness of the Reynolds number. The solution of the problem is derived for a quasi-steady flow involving infinitesimally thin disks. The numerical calculation is carried out for a set of different combinations of the interdisk distance and the ratio of the disk radius to the tube radius. The ratio of the velocity of the disk to the average velocity of the fluid is calculated. The different rates of transport of red blood cells and of plasma in capillary blood vessels are discussed. The average pressure gradient along the axis of the tube is computed, and the dependence of the effective viscosity of the blood on the hematocrit and the diameter of the capillary vessel is discussed.     

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.     

The effect of blood viscoelasticity on pulsatile flow in stationary and axially moving tubes

An analytical solution for pulsatile flow of a generalized Maxwell fluid in straight rigid tubes, with and without axial vessel motion, has been used to calculate the effect of blood viscoelasticity on velocity profiles and shear stress in flows representative of those in the large arteries. Measured bulk flow rate Q waveforms were used as starting points in the calculations for the aorta and femoral arteries, from which axial pressure gradient ▿P waves were derived that would reproduce the starting Q waves for viscoelastic flow. The ▿P waves were then used to calculate velocity profiles for both viscoelastic and purely viscous flow. For the coronary artery, published ▿P and axial vessel acceleration waveforms were used in a similar procedure to determine the separate and combined influences of viscoelasticity and vessel motion.Differences in local velocities, comparing viscous flow to viscoelastic flow, were in all cases less than about 2% of the peak local velocity. Differences in peak wall shear stress were less than about 3%.In the coronary artery, wall shear stress differences between viscous and viscoelastic flow were small, regardless of whether axial vessel motion was included. The shape of the wall shear stress waveform and its difference, however, changed dramatically between the stationary and moving vessel cases. The peaks in wall shear stress difference corresponded with large temporal gradients in the combined driving force for the flow.     

Pulsatile blood flow effects on temperature distribution and heat transfer in rigid vessels

The effect of blood velocity pulsations on bioheat transfer is studied. A simple model of a straight rigid blood vessel with unsteady periodic flow is considered. A numerical solution that considers the fully coupled Navier-Stokes and energy equations is used for the simulations. The influence of the pulsation rate on the temperature distribution and energy transport is studied for four typical vessel sizes: aorta, large arteries, terminal arterial branches, and arterioles. The results show that: the pulsating axial velocity produces a pulsating temperature distribution; reversal of flow occurs in the aorta and in large vessels, which produces significant time variation in the temperature profile. Change of the pulsation rate yields a change of the energy transport between the vessel wall and fluid for the large vessels. For the thermally important terminal arteries (0.04-1 mm), velocity pulsations have a small influence on temperature distribution and on the energy transport out of the vessels (8 percent for the Womersley number corresponding to a normal heart rate). Given that there is a small difference between the time-averaged unsteady heat flux due to a pulsating blood velocity and an assumed nonpulsating blood velocity, it is reasonable to assume a nonpulsating blood velocity for the purposes of estimating bioheat transfer.     

Effect of inlet velocity profiles on patient-specific computational fluid dynamics simulations of the carotid bifurcation

Patient-specific computational fluid dynamics (CFD) is a powerful tool for researching the role of blood flow in disease processes. Modern clinical imaging technology such as MRI and CT can provide high resolution information about vessel geometry, but in many situations, patient-specific inlet velocity information is not available. In these situations, a simplified velocity profile must be selected. We studied how idealized inlet velocity profiles (blunt, parabolic, and Womersley flow) affect patient-specific CFD results when compared to simulations employing a "reference standard" of the patient's own measured velocity profile in the carotid bifurcation. To place the magnitude of these effects in context, we also investigated the effect of geometry and the use of subject-specific flow waveform on the CFD results. We quantified these differences by examining the pointwise percent error of the mean wall shear stress (WSS) and the oscillatory shear index (OSI) and by computing the intra-class correlation coefficient (ICC) between axial profiles of the mean WSS and OSI in the internal carotid artery bulb. The parabolic inlet velocity profile produced the most similar mean WSS and OSI to simulations employing the real patient-specific inlet velocity profile. However, anatomic variation in vessel geometry and the use of a nonpatient-specific flow waveform both affected the WSS and OSI results more than did the choice of inlet velocity profile. Although careful selection of boundary conditions is essential for all CFD analysis, accurate patient-specific geometry reconstruction and measurement of vessel flow rate waveform are more important than the choice of velocity profile. A parabolic velocity profile provided results most similar to the patient-specific velocity profile.     

An introduction to biofluid mechanics--basic models and applications

Cardiovascular disease is the primary cause of morbidity and mortality in the western world. Complex hemodynamics play a critical role in the development of atherosclerosis and the processes of aging, as well as many other disease processes. Biofluid mechanics play a major role in the cardiovascular system and it is important to understand the forces and movement of blood cells and whole blood as well as the interaction between blood cells and the vessel wall. Fundamental fluid mechanical, which are important for the understanding of the blood flow in the cardiovascular circulatory system of the human body aspects are presented. Measurement techniques for model studies such as LDA, ultrasound, and MRI studies will be discussed. Viscosity and flow behavior changes specifically the creation of vortices and flow disturbances can be used to show how medication can influence flow behavior. Experiments have shown that hemodynamics may have a strong influence on the creation of aneurysms and varicose veins. Other factors such as vessel wall structure are also important. In preliminary studies, it has been demonstrated that geometry and elasticity of vessel walls help determine flow behavior. High velocity fluctuations indicate flow disturbances that should be avoided. Health care practitioners must understand fluid dynamic factors such as flow rate ratio, pressure and velocity gradients, and flow behavior, velocity distribution, shear stress on the wall and on blood cells. These mechanical factors are largely responsible for the deposit of blood cells and lipids, a leading cause of atherosclerosis. The interaction between blood cells and of the cells with the vessel, leads to the formation of plaques and agglomerations. These deposits are found predominantly at arterial bends and bifurcations where blood flow is disturbed, where a secondary flow is created, and where flow separation regions are found. Experiments on hemodynamic effects in elastic silicon rubber models of the cardiovascular system with flow wire, stents, or patches for vessel surgery will be discussed. These studies can be important in improving diagnostics and therapeutic applications.     

Analytical solutions for the specific activities of a traced substance within the component parts of a simulated kidney-ureter system

Two models for a kidney-ureter system are considered: one model of one vessel in which a traced substance, undergoing exchange between the vessel and an external compartment, is emptying into the ureter; the second model of two approximately parallel, identical vessels in which a traced substance, undergoing exchange between each vessel and an external compartment, is emptying into the ureter. A single impulsive input of label into a vessel is assumed. For mathematical simplicity, the major conditions imposed on each system are: (1) rapid mixing transverse to a vessel axis and no mixing longitudinal to a vessel axis within the plasma; (2) small variation of the specific activity within the plasma in the longitudinal direction to a vessel axis; (3) constant flow rate of urine into the ureter and (4) constant exchange coefficients, tubule flow velocity and traced substance concentrations within individual compartments.     

What measurements are necessary for a comprehensive evaluation of the peripheral arterial circulation?

Several methods are available to detect atherosclerotic lesions with a severe degree of stenosis (>70%), but the diagnosis of atherosclerotic lesions with no stenosis or with a minor degree of stenosis (<20%), is problematic. Hemodynamics associated with stenotic lesions are well described by the relationship of blood pressure and blood flow velocity, both as a function of time and localization (along the length and cross-section of the vessel). The use of this relationship in the clinic is difficult because no precise information is available about the geometry and branching of arteries, blood viscosity, and the velocity distribution over the cross-sectional area of the blood vessel. Besides, the invasiveness of the technique to measure arterial pressure as a function of time and localization does not allow routine application in patients. Because of these limitations, alternative methods have been developed. The degree and extensiveness of atherosclerotic disease can, for instance, be estimated from the changes in maximum blood flow velocity and in velocity profile, i.e., velocity distribution along the cross-section of the vessel. Moreover, the delay between simultaneously recorded arterial blood flow velocity tracings (pulse-wave velocity determination) is used to assess the elastic properties of the vessel. Changes in velocity profile occur at relatively slight degrees of arterial stenosis (around 20%), so that determination of these profiles along diseased arteries may contribute to the early diagnosis of atherosclerotic lesions. In man, transcutaneous information about the maximum and mean blood flow velocities over the cross-sectional area of the artery as an instantaneous function of time as well as the flow pattern can be obtained online with continuous wave Doppler flowmeters, at least when audio spectrum analysis is used as a processing technique. Velocity profiles can be determined with multichannel pulsed Doppler systems if the resolution of the system is adequate and a sufficient number of sample volumes can be obtained, limiting the interpolation between these samples. The on-line recording of velocity profiles can be facilitated by combining the pulsed Doppler device with either a velocity imaging system or a B-mode scan. In systems with a high resolution (sample distance 0.5 mm), one should be able to detect local disturbances in the velocity profile at the site of the lesion (due to local increases in shear stress) and proximal to the lesion (due to reflections), so that lesions with a minor degree of stenosis can be detected. In resistive systems (e.g., internal carotid arteries) in which the relationship between pressure and velocity changes during the cardiac cycle is relatively simple, the elasticity of the arterial wall can be determined by relating the relative diameter changes of the vessel, determined on-line with multichannel pulsed Doppler systems, to the instantaneous velocity pulse. Although the detection of atherosclerotic lesions at an early stage of the disease with sophisticated Doppler devices looks promising, further clinical evaluation is required.