The influence of boundary conditions on wall shear stress distribution in patients specific coronary trees

Patient specific geometrical data on human coronary arteries can be reliably obtained multislice computer tomography (MSCT) imaging. MSCT cannot provide hemodynamic variables, and the outflow through the side branches must be estimated. The impact of two different models to determine flow through the side branches on the wall shear stress (WSS) distribution in patient specific geometries is evaluated. Murray's law predicts that the flow ratio through the side branches scales with the ratio of the diameter of the side branches to the third power. The empirical model is based on flow measurements performed by Doriot et al. (2000) in angiographically normal coronary arteries. The fit based on these measurements showed that the flow ratio through the side branches can best be described with a power of 2.27. The experimental data imply that Murray's law underestimates the flow through the side branches. We applied the two models to study the WSS distribution in 6 coronary artery trees. Under steady flow conditions, the average WSS between the side branches differed significantly for the two models: the average WSS was 8% higher for Murray's law and the relative difference ranged from -5% to +27%. These differences scale with the difference in flow rate. Near the bifurcations, the differences in WSS were more pronounced: the size of the low WSS regions was significantly larger when applying the empirical model (13%), ranging from -12% to +68%. Predicting outflow based on Murray's law underestimates the flow through the side branches. Especially near side branches, the regions where atherosclerotic plaques preferentially develop, the differences are significant and application of Murray's law underestimates the size of the low WSS region.     

PIV-measured versus CFD-predicted flow dynamics in anatomically realistic cerebral aneurysm models

Computational fluid dynamics (CFD) modeling of nominally patient-specific cerebral aneurysms is increasingly being used as a research tool to further understand the development, prognosis, and treatment of brain aneurysms. We have previously developed virtual angiography to indirectly validate CFD-predicted gross flow dynamics against the routinely acquired digital subtraction angiograms. Toward a more direct validation, here we compare detailed, CFD-predicted velocity fields against those measured using particle imaging velocimetry (PIV). Two anatomically realistic flow-through phantoms, one a giant internal carotid artery (ICA) aneurysm and the other a basilar artery (BA) tip aneurysm, were constructed of a clear silicone elastomer. The phantoms were placed within a computer-controlled flow loop, programed with representative flow rate waveforms. PIV images were collected on several anterior-posterior (AP) and lateral (LAT) planes. CFD simulations were then carried out using a well-validated, in-house solver, based on micro-CT reconstructions of the geometries of the flow-through phantoms and inlet/outlet boundary conditions derived from flow rates measured during the PIV experiments. PIV and CFD results from the central AP plane of the ICA aneurysm showed a large stable vortex throughout the cardiac cycle. Complex vortex dynamics, captured by PIV and CFD, persisted throughout the cardiac cycle on the central LAT plane. Velocity vector fields showed good overall agreement. For the BA, aneurysm agreement was more compelling, with both PIV and CFD similarly resolving the dynamics of counter-rotating vortices on both AP and LAT planes. Despite the imposition of periodic flow boundary conditions for the CFD simulations, cycle-to-cycle fluctuations were evident in the BA aneurysm simulations, which agreed well, in terms of both amplitudes and spatial distributions, with cycle-to-cycle fluctuations measured by PIV in the same geometry. The overall good agreement between PIV and CFD suggests that CFD can reliably predict the details of the intra-aneurysmal flow dynamics observed in anatomically realistic in vitro models. Nevertheless, given the various modeling assumptions, this does not prove that they are mimicking the actual in vivo hemodynamics, and so validations against in vivo data are encouraged whenever possible.     

The evaluation of Murray's law in Psilotum nudum (Psilotaceae), an analogue of ancestral vascular plants

Previous work has shown that the xylem of seed plants follows Murray's law when conduits do not provide structural support to the plant. Here, compliance with Murray's law was tested in the stem photosynthesizer Psilotum nudum, a seedless vascular plant. Psilotum nudum was chosen because the central stele does not provide structural support, which means that Murray's law is applicable, and because its simple shoot structure resembles the earliest vascular plants. Murray's law predicts that the sum of the conduit radii cubed (Σr(3)) should decrease in direct proportion with the volume flow rate (Q) to maximize the hydraulic conductance per unit vascular investment. Agreement with Murray's law was assessed by estimating the transpiration rate distal to a cross-section, which should determine Q under steady state conditions, and comparing that with the Σr(3) of that cross-section. As predicted, regressions between the Σr(3) of the cross-section and Q resulted in a linear relationship with a y-intercept that was not different from zero. Two more rigorous statistical tests were also unable to reject Murray's law. Psilotum nudum plants also increased their conductance per investment by having more conduits distally than proximally, which is more efficient hydraulically than equal or declining conduit numbers distally.     

On connecting large vessels to small. The meaning of Murray's law

A large part of the branching vasculature of the mammalian circulatory and respiratory systems obeys Murray's law, which states that the cube of the radius of a parent vessel equals the sum of the cubes of the radii of the daughters. Where this law is obeyed, a functional relationship exists between vessel radius and volumetric flow, average linear velocity of flow, velocity profile, vessel-wall shear stress, Reynolds number, and pressure gradient in individual vessels. In homogeneous, full-flow sets of vessels, a relation is also established between vessel radius and the conductance, resistance, and cross- sectional area of a full-flow set.     

On the choice of outlet boundary conditions for patient-specific analysis of aortic flow using computational fluid dynamics

Boundary conditions (BCs) are an essential part in computational fluid dynamics (CFD) simulations of blood flow in large arteries. Although several studies have investigated the influence of BCs on predicted flow patterns and hemodynamic wall parameters in various arterial models, there is a lack of comprehensive assessment of outlet BCs for patient-specific analysis of aortic flow. In this study, five different sets of outlet BCs were tested and compared using a subject-specific model of a normal aorta. Phase-contrast magnetic resonance imaging (PC-MRI) was performed on the same subject and velocity profiles extracted from the in vivo measurements were used as the inlet boundary condition. Computational results obtained with different outlet BCs were assessed in terms of their agreement with the PC-MRI velocity data and key hemodynamic parameters, such as pressure and flow waveforms and wall shear stress related indices. Our results showed that the best overall performance was achieved by using a well-tuned three-element Windkessel model at all model outlets, which not only gave a good agreement with in vivo flow data, but also produced physiological pressure waveforms and values. On the other hand, opening outlet BCs with zero pressure at multiple outlets failed to reproduce any physiologically relevant flow and pressure features.     

Investigation of hemodynamics in the development of dissecting aneurysm within patient-specific dissecting aneurismal aortas using computational fluid dynamics (CFD) simulations

Aortic dissecting aneurysm is one of the most catastrophic cardiovascular emergencies that carries high mortality. It was pointed out from clinical observations that the aneurysm development is likely to be related to the hemodynamics condition of the dissected aorta. In order to gain more insight on the formation and progression of dissecting aneurysm, hemodynamic parameters including flow pattern, velocity distribution, aortic wall pressure and shear stress, which are difficult to measure in vivo, are evaluated using numerical simulations. Pulsatile blood flow in patient-specific dissecting aneurismal aortas before and after the formation of lumenal aneurysm (pre-aneurysm and post-aneurysm) is investigated by computational fluid dynamics (CFD) simulations. Realistic time-dependent boundary conditions are prescribed at various arteries of the complete aorta models. This study suggests the helical development of false lumen around true lumen may be related to the helical nature of hemodynamic flow in aorta. Narrowing of the aorta is responsible for the massive recirculation in the poststenosis region in the lumenal aneurysm development. High pressure difference of 0.21 kPa between true and false lumens in the pre-aneurismal aorta infers the possible lumenal aneurysm site in the descending aorta. It is also found that relatively high time-averaged wall shear stress (in the range of 4-8 kPa) may be associated with tear initiation and propagation. CFD modeling assists in medical planning by providing blood flow patterns, wall pressure and wall shear stress. This helps to understand various phenomena in the development of dissecting aneurysm.     

On the relative importance of rheology for image-based CFD models of the carotid bifurcation

BACKGROUND: Patient-specific computational fluid dynamics (CFD) models derived from medical images often require simplifying assumptions to render the simulations conceptually or computationally tractable. In this study, we investigated the sensitivity of image-based CFD models of the carotid bifurcation to assumptions regarding the blood rheology. METHOD OF APPROACH: CFD simulations of three different patient-specific models were carried out assuming: a reference high-shear Newtonian viscosity, two different non-Newtonian (shear-thinning) rheology models, and Newtonian viscosities based on characteristic shear rates or, equivalently, assumed hematocrits. Sensitivity of wall shear stress (WSS) and oscillatory shear index (OSI) were contextualized with respect to the reproducibility of the reconstructed geometry, and to assumptions regarding the inlet boundary conditions. RESULTS: Sensitivity of WSS to the various rheological assumptions was roughly 1.0 dyn/cm(2) or 8%, nearly seven times less than that due to geometric uncertainty (6.7 dyn/cm(2) or 47%), and on the order of that due to inlet boundary condition assumptions. Similar trends were observed regarding OSI sensitivity. Rescaling the Newtonian viscosity based on time-averaged inlet shear rate served to approximate reasonably, if overestimate slightly, non-Newtonian behavior. CONCLUSIONS: For image-based CFD simulations of the normal carotid bifurcation, the assumption of constant viscosity at a nominal hematocrit is reasonable in light of currently available levels of geometric precision, thus serving to obviate the need to acquire patient-specific rheological data.     

Approximating hemodynamics of cerebral aneurysms with steady flow simulations

Computational fluid dynamics (CFD) simulations can be employed to gain a better understanding of hemodynamics in cerebral aneurysms and improve diagnosis and treatment. However, introduction of CFD techniques into clinical practice would require faster simulation times. The aim of this study was to evaluate the use of computationally inexpensive steady flow simulations to approximate the aneurysm's wall shear stress (WSS) field. Two experiments were conducted. Experiment 1 compared for two cases the time-averaged (TA), peak systole (PS) and end diastole (ED) WSS field between steady and pulsatile flow simulations. The flow rate waveform imposed at the inlet was varied to account for variations in heart rate, pulsatility index, and TA flow rate. Consistently across all flow rate waveforms, steady flow simulations accurately approximated the TA, but not the PS and ED, WSS field. Following up on experiment 1, experiment 2 tested the result for the TA WSS field in a larger population of 20 cases covering a wide range of aneurysm volumes and shapes. Steady flow simulations approximated the space-averaged WSS with a mean error of 4.3%. WSS fields were locally compared by calculating the absolute error per node of the surface mesh. The coefficient of variation of the root-mean-square error over these nodes was on average 7.1%. In conclusion, steady flow simulations can accurately approximate the TA WSS field of an aneurysm. The fast computation time of 6 min per simulation (on 64 processors) could help facilitate the introduction of CFD into clinical practice.     

Investigating Murray's law in the chick embryo.

According to the optimization principle known as Murray's law, the blood vessel geometry at a bifurcation satisfies the relation alpha = (D3(1) + D3(2))/D3(0) = 1, where D0, D1, and D2 are the diameters of the parent and two daughter vessels, respectively. Previous investigations have shown that mature blood vessels adhere to this law fairly closely. The purpose of this study was to test Murray's law in the developing extraembryonic blood vessels of 2-4 day-old chick embryos. Vessel diameters were measured manually using image analysis software. The measurements for the group of all vessels at all studied stages (n = 449) gave alpha = 1.01+/-0.34 (mean +/- SD), and the value of alpha is similar at all stages. These results indicate that Murray's law holds in the chick embryo, even before medial smooth muscle becomes functional, suggesting that blood vessels follow the same basic morphogenetic rules throughout life.     

The effect of inlet waveforms on computational hemodynamics of patient-specific intracranial aneurysms

Due to the lack of patient-specific inlet flow waveform measurements, most computational fluid dynamics (CFD) simulations of intracranial aneurysms usually employ waveforms that are not patient-specific as inlet boundary conditions for the computational model. The current study examined how this assumption affects the predicted hemodynamics in patient-specific aneurysm geometries. We examined wall shear stress (WSS) and oscillatory shear index (OSI), the two most widely studied hemodynamic quantities that have been shown to predict aneurysm rupture, as well as maximal WSS (MWSS), energy loss (EL) and pressure loss coefficient (PLc). Sixteen pulsatile CFD simulations were carried out on four typical saccular aneurysms using 4 different waveforms and an identical inflow rate as inlet boundary conditions. Our results demonstrated that under the same mean inflow rate, different waveforms produced almost identical WSS distributions and WSS magnitudes, similar OSI distributions but drastically different OSI magnitudes. The OSI magnitude is correlated with the pulsatility index of the waveform. Furthermore, there is a linear relationship between aneurysm-averaged OSI values calculated from one waveform and those calculated from another waveform. In addition, different waveforms produced similar MWSS, EL and PLc in each aneurysm. In conclusion, inlet waveform has minimal effects on WSS, OSI distribution, MWSS, EL and PLc and a strong effect on OSI magnitude, but aneurysm-averaged OSI from different waveforms has a strong linear correlation with each other across different aneurysms, indicating that for the same aneurysm cohort, different waveforms can consistently stratify (rank) OSI of aneurysms.     

Toward translating near-infrared spectroscopy oxygen saturation data for the non-invasive prediction of spatial and temporal hemodynamics during exercise

Image-based computational fluid dynamics (CFD) studies conducted at rest have shown that atherosclerotic plaque in the thoracic aorta (TA) correlates with adverse wall shear stress (WSS), but there is a paucity of such data under elevated flow conditions. We developed a pedaling exercise protocol to obtain phase contrast magnetic resonance imaging (PC-MRI) blood flow measurements in the TA and brachiocephalic arteries during three-tiered supine pedaling at 130, 150, and 170 % of resting heart rate (HR), and relate these measurements to non-invasive tissue oxygen saturation \((\hbox {StO}_{2})\) acquired by near-infrared spectroscopy (NIRS) while conducting the same protocol. Local quantification of WSS indices by CFD revealed low time-averaged WSS on the outer curvature of the ascending aorta and the inner curvature of the descending aorta (dAo) that progressively increased with exercise, but that remained low on the anterior surface of brachiocephalic arteries. High oscillatory WSS observed on the inner curvature of the aorta persisted during exercise as well. Results suggest locally continuous exposure to potentially deleterious indices of WSS despite benefits of exercise. Linear relationships between flow distributions and tissue oxygen extraction calculated from \(\hbox {StO}_{2}\) were found between the left common carotid versus cerebral tissue \((r^{2}=0.96)\) and the dAo versus leg tissue \((r^{2}=0.87)\). A resulting six-step procedure is presented to use NIRS data as a surrogate for exercise PC-MRI when setting boundary conditions for future CFD studies of the TA under simulated exercise conditions. Relationships and ensemble-averaged PC-MRI inflow waveforms are provided in an online repository for this purpose.     

Volume flow estimation in the precapillary mesenteric microvasculature in vivo and the principle of constant pressure gradient

Volume flow was estimated from axial erythrocyte velocity measurements in 30 mesenteric microvessels of 6 rabbits and was compared to Murray's law predictions. The diameters of capillaries and precapillary arterioles ranged between 5.6 and 12 microm. The significant pulsating flow component existing in these microvessels was taken into account by measuring instantaneous axial blood velocity throughout the course of a cardiac cycle and then averaging over the period. In addition, the effect of the velocity profile variation with diameter was taken into account, for the first time, by using a profile factor function. According to Murray's law, the relation between blood volume flow and diameter is governed by a 'cube' law. Curve fitting to volume flow and diameter data pairs for rabbits, showed a dependence of volume flow on diameter raised to the 4th power (with a correlation coefficient equal to 0.97). The above result supports the hypothesis that, in the precapillary part of microvasculature, the principle of constant longitudinal pressure gradient rather than the principle of minimum work may be valid.     

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.     

Computing the ankle-brachial index with parallel computational fluid dynamics

The ankle-brachial index (ABI), a ratio of arterial blood pressure in the ankles and upper arms, is used to diagnose and monitor circulatory conditions such as coarctation of the aorta and peripheral artery disease. Computational simulations of the ABI can potentially determine the parameters that produce an ABI indicative of ischemia or other abnormalities in blood flow. However, 0- and 1-D computational methods are limited in describing a 3-D patient-derived geometry. Thus, we present a massively parallel framework for computational fluid dynamics (CFD) simulations in the full arterial system. Using the lattice Boltzmann method to solve the Navier–Stokes equations, we employ highly parallelized and scalable methods to generate the simulation domain and efficiently distribute the computational load among processors. For the first time, we compute an ABI with 3-D CFD. In this proof-of-concept study, we investigate the dependence of ABI on the presence of stenoses, or narrowed regions of the arteries, by directly modifying the arterial geometry. As a result, our framework enables the computation a hemodynamic factor characterizing flow at the scale of the full arterial system, in a manner that is extensible to patient-specific imaging data and holds potential for treatment planning.     

On the principles of the vascular network branching

We propose an explanation of Murray's law without applying the minimality principles. The model deals with a "delivering" artery system of an organ that is characterized, first, by the space-filling embedding into the organ tissue and, second, by the uniform distribution of the blood pressure drop over it. The latter assumption is justified using the available physiological data and the idea about conditions needed for perfect self-regulation. Based on the two statements we get Murray's law, and so, demonstrate that it can be also regarded as a direct consequence of the organism's capacity for controlling finely the blood flow redistribution over peripheral vascular networks.     

What is needed to make low-density lipoprotein transport in human aorta computational models suitable to explore links to atherosclerosis? Impact of initial and inflow boundary conditions

Personalized computational hemodynamics (CH) is a promising tool to clarify/predict the link between low density lipoproteins (LDL) transport in aorta, disturbed shear and atherogenesis. However, CH uses simplifying assumptions that represent sources of uncertainty. In particular, modelling blood-side to wall LDL transfer is challenged by the cumbersomeness of protocols needed to obtain reliable LDL concentration profile estimations. This paucity of data is limiting the establishment of rigorous CH protocols able to balance the trade-offs among the variety of in vivo data to be acquired, and the accuracy required by biological/clinical applications.In this study, we analyze the impact of LDL concentration initialization (initial conditions, ICs) and inflow boundary conditions (BCs) on CH models of LDL blood-to-wall transfer in aorta. Technically, in an image-based model of human aorta, two different inflow BCs are generated imposing subject-specific inflow 3D PC-MRI measured or idealized (flat) velocity profiles. For each simulated BC, four different ICs for LDL concentration are applied, imposing as IC the LDL distribution resulting from steady-state simulations with average conditions, or constant LDL concentration values.Based on CH results, we conclude that: (1) the imposition of realistic 3D velocity profiles as inflow BC reduces the uncertainty affecting the representation of LDL transfer; (2) different LDL concentration ICs lead to markedly different patterns of LDL transfer.Given that it is not possible to verify in vivo the proper LDL concentration initialization to be applied, we suggest to carefully set and unambiguously declare the imposed BCs and LDL concentration IC when modelling LDL transfer in aorta, in order to obtain reproducible and ultimately comparable results among different laboratories.     

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.     

Flow mixing during peripheral veno-arterial extra corporeal membrane oxygenation – A simulation study

Peripheral veno-arterial extra-corporeal membrane oxygenation (ECMO) is an artificial circulation that supports patients with severe cardiac and respiratory failure. Differential hypoxia during ECMO support has been reported, and it has been suggested that it is due to the mixing of well-perfused retrograde ECMO flow and poorly-perfused antegrade left ventricle (LV) flow in the aorta. This study aims to quantify the relationship between ECMO support level and location of the mixing zone (MZ) of the ECMO and LV flows. Steady-state and transient computational fluid dynamics (CFD) simulations were performed using a patient-specific geometrical model of the aorta. A range of ECMO support levels (from 5% to 95% of total cardiac output) were evaluated. For ECMO support levels above 70%, the MZ was located in the aortic arch, resulting in perfusion of the arch branches with poorly perfused LV flow. The MZ location was stable over the cardiac cycle for high ECMO flows (>70%), but moved 5 cm between systole and diastole for ECMO support level of 60%. This CFD approach has potential to improve individual patient care and ECMO design.     

Which diameter and angle rule provides optimal flow patterns in a coronary bifurcation?

The branching angle and diameter ratio in epicardial coronary artery bifurcations are two important determinants of atherogenesis. Murray's cubed diameter law and bifurcation angle have been assumed to yield optimal flows through a bifurcation. In contrast, we have recently shown a 7/3 diameter law (HK diameter model), based on minimum energy hypothesis in an entire tree structure. Here, we derive a bifurcation angle rule corresponding to the HK diameter model and critically evaluate the streamline flow through HK and Murray-type bifurcations. The bifurcations from coronary casts were found to obey the HK diameter model and angle rule much more than Murray's model. A finite element model was used to investigate flow patterns for coronary artery bifurcations of various types. The inlet velocity and pressure boundary conditions were measured by ComboWire. Y-bifurcation of Murray type decreased wall shear stress-WSS (10%-40%) and created an increased oscillatory shear index-OSI in atherosclerosis-prone regions as compared with HK-type bifurcations. The HK-type bifurcations were found to have more optimal flow patterns (i.e., higher WSS and lower OSI) than Murray-type bifurcations which have been traditionally believed to be optimized. This study has implications for changes in bifurcation angles and diameters in percutaneous coronary intervention.