Analysis of blood flow in the entire coronary arterial tree

A hemodynamic analysis of coronary blood flow must be based on the measured branching pattern and vascular geometry of the coronary vasculature. We recently developed a computer reconstruction of the entire coronary arterial tree of the porcine heart based on previously measured morphometric data. In the present study, we carried out an analysis of blood flow distribution through a network of millions of vessels that includes the entire coronary arterial tree down to the first capillary branch. The pressure and flow are computed throughout the coronary arterial tree based on conservation of mass and momentum and appropriate pressure boundary conditions. We found a power law relationship between the diameter and flow of each vessel branch. The exponent is approximately 2.2, which deviates from Murray's prediction of 3.0. Furthermore, we found the total arterial equivalent resistance to be 0.93, 0.77, and 1.28 mmHg.ml(-1).s(-1).g(-1) for the right coronary artery, left anterior descending coronary artery, and left circumflex artery, respectively. The significance of the present study is that it yields a predictive model that incorporates some of the factors controlling coronary blood flow. The model of normal hearts will serve as a physiological reference state. Pathological states can then be studied in relation to changes in model parameters that alter coronary perfusion.     

Pulsatile blood flow in the entire coronary arterial tree: theory and experiment

The pulsatility of coronary circulation can be accurately simulated on the basis of the measured branching pattern, vascular geometry, and material properties of the coronary vasculature. A Womersley-type mathematical model is developed to analyze pulsatile blood flow in diastole in the absence of vessel tone in the entire coronary arterial tree on the basis of previously measured morphometric data. The model incorporates a constitutive equation of pressure and cross-section area relation based on our previous experimental data. The formulation enables the prediction of the impedance, the pressure distribution, and the pulsatile flow distribution throughout the entire coronary arterial tree. The model is validated by experimental measurements in six diastolic arrested, vasodilated porcine hearts. The agreement between theory and experiment is excellent. Furthermore, the present pulse wave results at low frequency agree very well with previously published steady-state model. Finally, the phase angle of flow is seen to decrease along the trunk of the major coronary artery and primary branches toward the capillary vessels. This study represents the first, most extensive validated analysis of Womersley-type pulse wave transmission in the entire coronary arterial tree down to the first segment of capillaries. The present model will serve to quantitatively test various hypotheses in the coronary circulation under pulsatile flow conditions.     

The role of the Womersley number in pulsatile blood flow a theoretical study of the Casson model.

The purpose of this Note is to clarify the meaning of the Womersley number alpha in pulsatile blood flow in small vessels. In particular. we explain why the use of alpha as aperturbation parameter to obtain approximate solutions of the Casson model (frequently used in the literature) is not appropriate. Using the techniques of dimensional analysis and scaling, we show that alpha is the product of the Reynolds and Strouhal numbers. Since the latter is very small for physiological flows, the result is that alpha < 1 even at relatively high values of the Reynolds number (i.e., for non-negligible inertia) and we validate our perturbation theory results by comparison with a numerical integration of the full model. Although this analysis is based on the Casson model, our method has general validity and may be applied to other models which describe more accurately the rheological properties of blood.     

The inverse Womersley problem for pulsatile flow in straight rigid tubes

In this study a numerical solution for the problem of pulsating flow in rigid tubes is described. The method applies to the case of known flow rate waveform, as opposed to Womersley solution where the pressure gradient was the known quantity. The solution provides the pressure gradient and wall shear stress waveforms as well as the instantaneous velocity profiles. Results show that the method can be used to study the blood flow characteristics in large arteries.     

Validation of a patient-specific one-dimensional model of the systemic arterial tree

The aim of this study is to develop and validate a patient-specific distributed model of the systemic arterial tree. This model is built using geometric and hemodynamic data measured on a specific person and validated with noninvasive measurements of flow and pressure on the same person, providing thus a patient-specific model and validation. The systemic arterial tree geometry was obtained from MR angiographic measurements. A nonlinear viscoelastic constitutive law for the arterial wall is considered. Arterial wall distensibility is based on literature data and adapted to match the wave propagation velocity of the main arteries of the specific subject, which were estimated by pressure waves traveling time. The intimal shear stress is modeled using the Witzig-Womersley theory. Blood pressure is measured using applanation tonometry and flow rate using transcranial ultrasound and phase-contrast-MRI. The model predicts pressure and flow waveforms in good qualitative and quantitative agreement with the in vivo measurements, in terms of wave shape and specific wave features. Comparison with a generic one-dimensional model shows that the patient-specific model better predicts pressure and flow at specific arterial sites. These results obtained let us conclude that a patient-specific one-dimensional model of the arterial tree is able to predict well pressure and flow waveforms in the main systemic circulation, whereas this is not always the case for a generic one-dimensional model.     

Analysis of pulsatile blood flow: A carotid siphon model

Numerical results for axial and secondary flow velocity and pressure in a three-dimensional model of the human carotid siphon have been calculated; the investigations were carried out under physiologically relevant pulsatile flow conditoins. Time-dependent, three-dimensional Navier-Stokes equations were solved numerically by using a special finite element method. The results of the computer simulation presented here concentrate on the secondary motion effect during the pulsatile flow cycle in multiple three-dimensional curvatures.     

A model of blood flow in the mesenteric arterial system

Background  

There are some early clinical indicators of cardiac ischemia, most notably a change in a person's electrocardiogram. Less well understood, but potentially just as dangerous, is ischemia that develops in the gastrointestinal system. Such ischemia is difficult to diagnose without angiography (an invasive and time-consuming procedure) mainly due to the highly unspecific nature of the disease.     

A novel,FFT-based one-dimensional blood flow solution method for arterial network

In the present work, we propose an FFT-based method for solving blood flow equations in an arterial network with variable properties and geometrical changes. An essential advantage of this approach is in correctly accounting for the vessel skin friction through the use of Womersley solution. To incorporate nonlinear effects, a novel approximation method is proposed to enable calculation of nonlinear corrections. Unlike similar methods available in the literature, the set of algebraic equations required for every harmonic is constructed automatically. The result is a generalized, robust and fast method to accurately capture the increasing pulse wave velocity downstream as well as steepening of the pulse front. The proposed method is shown to be appropriate for incorporating correct convection and diffusion coefficients. We show that the proposed method is fast and accurate and it can be an effective tool for 1D modelling of blood flow in human arterial networks.

     

A turbulence model for pulsatile arterial flows

Difficulties in predicting the behavior of some high Reynolds number flows in the circulatory system stem in part from the severe requirements placed on the turbulence model chosen to close the time-averaged equations of fluid motion. In particular, the successful turbulence model is required to (a) correctly capture the "nonequilibrium" effects wrought by the interactions of the organized mean-flow unsteadiness with the random turbulence, (b) correctly reproduce the effects of the laminar-turbulent transitional behavior that occurs at various phases of the cardiac cycle, and (c) yield good predictions of the near-wall flow behavior in conditions where the universal logarithmic law of the wall is known to be not valid. These requirements are not immediately met by standard models of turbulence that have been developed largely with reference to data from steady, fully turbulent flows in approximate local equilibrium. The purpose of this paper is to report on the development of a turbulence model suited for use in arterial flows. The model is of the two-equation eddy-viscosity variety with dependent variables that are zero-valued at a solid wall and vary linearly with distance from it. The effects of transition are introduced by coupling this model to the local value of the intermittency and obtaining the latter from the solution of a modeled transport equation. Comparisons with measurements obtained in oscillatory transitional flows in circular tubes show that the model produces substantial improvements over existing closures. Further pulsatile-flow predictions, driven by a mean-flow wave form obtained in a diseased human carotid artery, indicate that the intermittency-modified model yields much reduced levels of wall shear stress compared to the original, unmodified model. This result, which is attributed to the rapid growth in the thickness of the viscous sublayer arising from the severe acceleration of systole, argues in favor of the use of the model for the prediction of arterial flows.     

Comparison of steady and pulsatile flow in a double branching arterial model

Data are presented to compare fluid flow parameters for steady flow with those for time-varying flow in a simplified two branch model which simulates the region of the abdominal aorta near the celiac and superior mesenteric branches of the dog. Measurements in the model included laser doppler anemometry velocity profiles during steady flow, sinusoidal flow with a superimposed mean flow (referred to as simple oscillatory flow) and arterial pulsatile flow. Shear rate measurements were made by an electrochemical technique during steady flow. Flow visualization studies were done during steady and pulsatile flow. Fluid flow effects in the simplified model during steady flow showed many similarities to the results from previous steady flow studies in a canine aortic cast. Shear rates in the region of the proximal (first, or celiac) branch were independent of flow rates in the distal (second, or mesenteric) branch, but the shear pattern within the proximal branch changed significantly as flow in the proximal branch increased. Shear rates on the proximal flow divider (leading edge into the distal branch) depended primarily on the flow rate to the proximal branch, but not on flow to the distal branch. At certain daughter branch flow ratios (approximately 2:1, proximal to distal), flow separation was promoted at the outer wall of the second branch, but flow separation did not occur in the first branch. In contrast to the canine aortic case results, flow separation was never detected on the distal (mesenteric) flow divider of the simplified model. This observation reflects the subtle effects of geometry on flow since the mesenteric flow divider in the canine cast protrudes into the main flow whereas the distal flow divider in the simplified model does not. There were distinct differences in the flow phenomena between steady, simple oscillatory and arterial pulsatile flow. Peak shear rates during pulsatile flow were as much as 10--100 times greater than steady flow shear rates at comparable mean flow rates. Particularly noteworthy for the pulsatile flow with a Womersley parameter of sixteen were very blunt velocity profiles throughout systole, and the absence of flow separation or reversal in those regions of the model that exhibited flow separation during steady flow. The shape of the waveform influences the nature of the flow during time-varying flows. Future studies of fluid dynamics in model systems must consider the pulsatile nature of the flow if a true interpretation of arterial flow phenomena is to be made.     

Analysis of pulmonary arterial pressure profile after occlusion of pulsatile blood flow

In isolated canine lung lobes perfused with a pulsatile pump, arterial occlusions were performed and the postocclusion arterial pressure profiles were analyzed to estimate the pulmonary capillary pressure. A solenoid valve interposed between the pump and the lobar artery was used to perform arterial occlusions at several instants equally distributed within a pressure cycle. Double occlusions were also accomplished by simultaneously activating the solenoid valve and clamping the venous outflow of the lung lobe. To analyze an arterial occlusion pressure profile, we computed the best monoexponential fit of the pressure decay over a short period of time after the occlusion maneuvers. Two estimates of the capillary pressure were derived from this analysis: 1) the extrapolation of the exponential fit to the instant of occlusion, and 2) the point at which the recorded pressure decay curve merges with the exponential fit. The pressures thus determined were compared with the double occlusion pressure that provided an independent estimate of the pulmonary capillary pressure. Our results show that, under a wide range of conditions, the estimates of the capillary pressure obtained from the arterial occlusion data are nearly equal to the double occlusion pressures. Additionally, we estimated the capillary pressure variations within a pressure cycle by examining the occlusion pressures sampled at different instants of the cycle. The pulsatility of the pulmonary microvascular pressure varied with the pump frequency as well as the state of arterial and venous vasoaction. These variations are consistent with the representation of the lung vasculature as a low-pass filter.     

A simple coronary blood flow model to study the collateral flow index

Biomechanics and Modeling in Mechanobiology - In this work, we present a novel modeling framework to investigate the effects of collateral circulation into the coronary blood flow physiology. A...     

Computational study of pulsatile blood flow in prototype vessel geometries of coronary segments

The spatial and temporal distributions of wall shear stress (WSS) in prototype vessel geometries of coronary segments are investigated via numerical simulation, and the potential association with vascular disease and specifically atherosclerosis and plaque rupture is discussed. In particular, simulation results of WSS spatio-temporal distributions are presented for pulsatile, non-Newtonian blood flow conditions for: (a) curved pipes with different curvatures, and (b) bifurcating pipes with different branching angles and flow division. The effects of non-Newtonian flow on WSS (compared to Newtonian flow) are found to be small at Reynolds numbers representative of blood flow in coronary arteries. Specific preferential sites of average low WSS (and likely atherogenesis) were found at the outer regions of the bifurcating branches just after the bifurcation, and at the outer-entry and inner-exit flow regions of the curved vessel segment. The drop in WSS was more dramatic at the bifurcating vessel sites (less than 5% of the pre-bifurcation value). These sites were also near rapid gradients of WSS changes in space and time – a fact that increases the risk of rupture of plaque likely to develop at these sites. The time variation of the WSS spatial distributions was very rapid around the start and end of the systolic phase of the cardiac cycle, when strong fluctuations of intravascular pressure were also observed. These rapid and strong changes of WSS and pressure coincide temporally with the greatest flexion and mechanical stresses induced in the vessel wall by myocardial motion (ventricular contraction). The combination of these factors may increase the risk of plaque rupture and thrombus formation at these sites.     

Microcontinuum model for pulsatile blood flow through a stenosed tube

The effects of polar nature of blood and pulsatility on flow through a stenosed tube have been analysed by assuming blood as a micropolar fluid. Linearized solutions of basic equations are obtained through consecutive applications of finite Hankel and Laplace transforms. The analytical expressions for axial and particle angular velocities, wall shear stress, resistance to flow and apparent viscosity have been obtained. The axial velocity profiles for Newtonian and micropolar fluids have been compared. The interesting observation of this analysis is velocity, in certain parts of cycle, for micropolar fluid is higher than Newtonain fluid. Variation of apparent viscosity eta a with tube radius shows both inverse Fahraeus-Lindqvist and Fahraeus-Lindqvist effects. Finally, the resistance to flow and wall shear stress for normal and diseased blood have been computed and compared.     

Numerical simulations of pulsatile blood flow using a new constitutive model

In the present paper we use a new constitutive equation for whole human blood [R.G. Owens, A new microstructure-based constitutive model for human blood, J. Non-Newtonian Fluid Mech. (2006), to appear] to investigate the steady, oscillatory and pulsatile flow of blood in a straight, rigid walled tube at modest Womersley numbers. Comparisons are made with the experimental results of Thurston [Elastic effects in pulsatile blood flow, Microvasc. Res. 9 (1975), 145-157] for the pressure drop per unit length against volume flow rate and oscillatory flow rate amplitude. Agreement in all cases is very good. In the presentation of the numerical and experimental results we discuss the microstructural changes in the blood that account for its rheological behaviour in this simple class of flows. In this context, the concept of an apparent complex viscosity proves to be useful.     

Balance between myogenic,flow-dependent,and metabolic flow control in coronary arterial tree: a model study

Myogenic response, flow-dependent dilation, and direct metabolic control are important mechanisms controlling coronary flow. A model was developed to study how these control mechanisms interact at different locations in the arteriolar tree and to evaluate their contribution to autoregulatory and metabolic flow control. The model consists of 10 resistance compartments in series, each representing parallel vessel units, with their diameters determined by tone depending on either flow and pressure [flow-dependent tone reduction factor (TRF(flow)) x Tone(myo)] or directly on metabolic factors (Tone(meta)). The pressure-Tone(myo) and flow-TRF(flow) relations depend on the vessel size obtained from interpolation of data on isolated vessels. Flow-dependent dilation diminishes autoregulatory properties compared with pressure-flow lines obtained from vessels solely influenced by Tone(myo). By applying Tone(meta) to the four distal compartments, the autoregulatory properties are restored and tone is equally distributed over the compartments. Also, metabolic control and blockage of nitric oxide are simulated. We conclude that a balance is required between the flow-dependent properties upstream and the constrictive metabolic properties downstream. Myogenic response contributes significantly to flow regulation.     

Stress analysis in a layered aortic arch model under pulsatile blood flow

Background  

Many cardiovascular diseases, such as aortic dissection, frequently occur on the aortic arch and fluid-structure interactions play an important role in the cardiovascular system. Mechanical stress is crucial in the functioning of the cardiovascular system; therefore, stress analysis is a useful tool for understanding vascular pathophysiology. The present study is concerned with the stress distribution in a layered aortic arch model with interaction between pulsatile flow and the wall of the blood vessel.     

Developed pulsatile flow in a deployed coronary stent

The patho-physiologic process of restenosis and tissue growth may not be completely eliminated and is the primary concern of clinicians performing angioplasty and stent implantation procedures. Recent evidence suggests that the restenosis process is influenced by several factors: (1) geometry and size of vessel; (2) stent design; and (3) it's location that alter hemodynamic parameters, including local wall shear stress (WSS) distributions. The present three-dimensional (3D) analysis of pulsatile flow in a deployed coronary stent: (1) shows complex 3D variation of hemodynamic parameters; and (2) quantifies the changes in local WSS distributions for developed flow and compares with recently published WSS data for developing flow. Higher order of magnitude of WSS of 290 dyn/cm(2) is observed on the surface of cross-link intersections at the entrance of the stent for developed flow, which is about half of that for developing flow. Low WSS of 0.8 dyn/cm(2) and negative WSS of -8 dyn/cm(2) are seen at the immediate upstream and downstream regions of strut intersections. Persistent recirculation is observed at the downstream region of each strut cross-link and the regions of low and negative WSS may lead to patho-physiologic conditions near the stented region. The key finding of this study is that the location of stent in the coronary artery determines the developing or developed nature of the flow, which in turn, results in varied level of WSS.