Numerical simulation of local blood flow in the carotid and cerebral arteries under altered gravity

A computational fluid dynamics (CFD) approach was presented to model the blood flows in the carotid bifurcation and the brain arteries under altered gravity. Physical models required for CFD simulation were introduced including a model for arterial wall motion due to fluid-wall interactions, a shear thinning fluid model of blood, a vascular bed model for outflow boundary conditions, and a model for autoregulation mechanism. The three-dimensional unsteady incompressible Navier-Stokes equations coupled with these models were solved iteratively using the pseudocompressibility method and dual time stepping. Gravity source terms were added to the Navier-Stokes equations to take the effect of gravity into account. For the treatment of complex geometry, a chimera overset grid technique was adopted to obtain connectivity between arterial branches. For code validation, computed results were compared with experimental data for both steady-state and time-dependent flows. This computational approach was then applied to blood flows through a realistic carotid bifurcation and two Circle of Willis models, one using an idealized geometry and the other using an anatomical data set. A three-dimensional Circle of Willis configuration was reconstructed from subject-specific magnetic resonance images using an image segmentation method. Through the numerical simulation of blood flow in two model problems, namely, the carotid bifurcation and the brain arteries, it was observed that the altered gravity has considerable effects on arterial contraction/dilatation and consequent changes in flow conditions.     

Validation of a fluid–structure interaction numerical model for predicting flow transients in arteries

Fluid–structure interaction (FSI) numerical models are now widely used in predicting blood flow transients. This is because of the importance of the interaction between the flowing blood and the deforming arterial wall to blood flow behaviour. Unfortunately, most of these FSI models lack rigorous validation and, thus, cannot guarantee the accuracy of their predictions. This paper presents the comprehensive validation of a two-way coupled FSI numerical model, developed to predict flow transients in compliant conduits such as arteries. The model is validated using analytical solutions and experiments conducted on polyurethane mock artery. Flow parameters such as pressure and axial stress (and precursor) wave speeds, wall deformations and oscillating frequency, fluid velocity and Poisson coupling effects, were used as the basis of this validation. Results show very good comparison between numerical predictions, analytical solutions and experimental data. The agreement between the three approaches is generally over 95%. The model also shows accurate prediction of Poisson coupling effects in unsteady flows through flexible pipes, which up to this stage have only being predicted analytically. Therefore, this numerical model can accurately predict flow transients in compliant vessels such as arteries.     

Abnormal arterial flows by a distributed model of the fetal circulation

Modeling the propagation of blood pressure and flow along the fetoplacental arterial tree may improve interpretation of abnormal flow velocity waveforms in fetuses. The current models, however, either do not include a wide range of gestational ages or do not account for variation in anatomical, vascular, or rheological parameters. We developed a mathematical model of the pulsating fetoumbilical arterial circulation using Womersley's oscillatory flow theory and viscoelastic arterial wall properties. Arterial flow waves are calculated at different arterial locations from which the pulsatility index (PI) can be determined. We varied blood viscosity, placental and brain resistances, placental compliance, heart rate, stiffness of the arterial wall, and length of the umbilical arteries. The PI increases in the umbilical artery and decreases in the cerebral arteries, as a result of increasing placental resistance or decreasing brain resistance. Both changes in resistance decrease the flow through the placenta. An increased arterial stiffness increases the PIs in the entire fetoplacental circulation. Blood viscosity and peripheral bed compliance have limited influence on the flow profiles. Bradycardia and tachycardia increase and decrease the PI in all arteries, respectively. Umbilical arterial length has limited influence on the PI but affects the mean arterial pressure at the placental cord insertion. The model may improve the interpretation of arterial flow pulsations and thus may advance both the understanding of pathophysiological processes and clinical management.     

Waveform dependence of pulsatile flow in a stenosed channel

Bloodflow in arteries often shows a rich variety of vortical flows, which are dominated by the complex geometry of blood vessels, the dynamic pulsation of blood flow, and the complicated boundary conditions. With a two-dimensional model of unsteady flow in a stenosed channel, the pulsatile influence on such vortical fluid dynamics has been numerically studied in terms of waveform dependence on physiological pulsation. Results are presented for unsteady flows downstream of the stenosed portion with variation in the wavefiorms of systole and diastole. Overall, a train of propagating vortex waves is observed for all the cases, but it shows great sensitivity to the waveforms. The generation and development of the vortex waves may be linked to the presence of an adverse pressure gradient within a specific interval between two points of inflection of the systolic waveform. The adverse pressure gradient consists of a global pressure gradient that is found to be closely related to the dynamnics of' the pulsation, and a local pressure gradient, which is obsented to be dominated by the nonlinear vortex dynamics.     

Three-dimensional hemodynamics analysis of the circle of Willis in the patient-specific nonintegral arterial structures

The hemodynamic alteration in the cerebral circulation caused by the geometric variations in the cerebral circulation arterial network of the circle of Wills (CoW) can lead to fatal ischemic attacks in the brain. The geometric variations due to impairment in the arterial network result in incomplete cerebral arterial structure of CoW and inadequate blood supply to the brain. Therefore, it is of great importance to understand the hemodynamics of the CoW, for efficiently and precisely evaluating the status of blood supply to the brain. In this paper, three-dimensional computational fluid dynamics of the main CoW vasculature coupled with zero-dimensional lumped parameter model boundary condition for the CoW outflow boundaries is developed for analysis of the blood flow distribution in the incomplete CoW cerebral arterial structures. The geometric models in our study cover the arterial segments from the aorta to the cerebral arteries, which can allow us to take into account the innate patient-specific resistance of the arterial trees. Numerical simulations of the governing fluid mechanics are performed to determine the CoW arterial structural hemodynamics, for illustrating the redistribution of the blood flow in CoW due to the structural variations. We have evaluated our coupling methodology in five patient-specific cases that were diagnosed with the absence of efferent vessels or impairment in the connective arteries in their CoWs. The velocity profiles calculated by our approach in the segments of the patient-specific arterial structures are found to be very close to the Doppler ultrasound measurements. The accuracy and consistency of our hemodynamic results have been improved (to \(16.1 \pm 18.5\) %) compared to that of the pure-resistance boundary conditions (of 43.5 \(\pm \) 28 %). Based on our grouping of the five cases according to the occurrence of unilateral occlusion in vertebral arteries, the inter-comparison has shown that (i) the flow reduction in posterior cerebral arteries is the consequence of the unilateral vertebral arterial occlusion, and (ii) the flow rate in the anterior cerebral arteries is correlated with the posterior structural variations. This study shows that our coupling approach is capable of providing comprehensive information of the hemodynamic alterations in the pathological CoW arterial structures. The information generated by our methodology can enable evaluation of both the functional and structural status of the clinically significant symptoms, for assisting the treatment decision-making.     

A hybrid one-dimensional/Womersley model of pulsatile blood flow in the entire coronary arterial tree

Using a frequency-domain Womersley-type model, we previously simulated pulsatile blood flow throughout the coronary arterial tree. Although this model represents a good approximation for the smaller vessels, it does not take into account the nonlinear convective energy losses in larger vessels. Here, using Womersley's theory, we present a hybrid model that considers the nonlinear effects for the larger epicardial arteries while simulating the distal vessels (down to the 1st capillary segments) with the use of Womersley's Theory. The main trunk and primary branches were discretized and modeled with one-dimensional Navier-Stokes equations, while the smaller-diameter vessels were treated as Womersley-type vessels. Energy losses associated with vessel bifurcations were incorporated in the present analysis. The formulation enables prediction of impedance and pressure and pulsatile flow distribution throughout the entire coronary arterial tree down to the first capillary segments in the arrested, vasodilated state. We found that the nonlinear convective term is negligible and the loss of energy at a bifurcation is small in the larger epicardial vessels of an arrested heart. Furthermore, we found that the flow waves along the trunk or at the primary branches tend to scale (normalized with respect to their mean values) to a single curve, except for a small phase angle difference. Finally, the model predictions for the inlet pressure and flow waves are in excellent agreement with previously published experimental results. This hybrid one-dimensional/Womersley model is an efficient approach that captures the essence of the hemodynamics of a complex large-scale vascular network. The present model has numerous applications to understanding the dynamics of coronary circulation.     

The arterial circulation of the heart

Blood flows into the aorta and its branches during left ventricular systole. Most of the arterial walls in the body stretch during systole in accordance with their elastic properties (Roston, 1962a, b). During diastole, the rebound of the distended walls supplies an additional propulsive force pushing the blood forward. Since the metabolic exchange between most of the tissues in the body and their blood vessels is ordinarily the same throughout the cardiac cycle, it makes little difference whether or not the blood flow occurs during systole or diastole. The circulation in the coronary arteries behaves in a quite different way. Because the muscle fibers of the heart contract during systole and relax during diastole, different conditions for blood flow and metabolic exchange exist during the phases of the cardiac cycle. As a result, specification of whether blood flows in the coronary arteries during systole or diastole may be important. Such specification complicates the study of the coronary artery circulation. For example, because of the arterial elasticity, some of the blood which enters the coronary arteries during diastole comes in contact with the muscle fibers during systole. The present work contains a theoretical study of the coronary artery circulation.     

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.     

Flow in a mechanical bileaflet heart valve at laminar and near-peak systole flow rates: CFD simulations and experiments

Time-accurate, fully 3D numerical simulations and particle image velocity laboratory experiments are carried out for flow through a fully open bileaflet mechanical heart valve under steady (nonpulsatile) inflow conditions. Flows at two different Reynolds numbers, one in the laminar regime and the other turbulent (near-peak systole flow rate), are investigated. A direct numerical simulation is carried out for the laminar flow case while the turbulent flow is investigated with two different unsteady statistical turbulence modeling approaches, unsteady Reynolds-averaged Navier-Stokes (URANS) and detached-eddy simulation (DES) approach. For both the laminar and turbulent cases the computed mean velocity profiles are in good overall agreement with the measurements. For the turbulent simulations, however, the comparisons with the measurements demonstrate clearly the superiority of the DES approach and underscore its potential as a powerful modeling tool of cardiovascular flows at physiological conditions. The study reveals numerous previously unknown features of the flow.     

Dynamic viscous flow in distensible vessels of skeletal muscle microcirculation: application to pressure and flow transients

Blood flow in the microcirculation of the rat skeletal muscle during transient changes of arterial pressure is analyzed theoretically. Although flow in such small vessels is quasi-steady and has a very low Reynolds number, time-dependent nonuniform flows along the length of the blood vessels can be observed due to vessel distensibility. The governing equations for a single microvessel are derived using previously measured microvessel elasticity, and several solutions to different inflow and outflow pressures and flow conditions are investigated. The results indicate that when such distensible microvessels are subjected to a step increase of arterial pressure, the arterial flow shows a rapid overshoot followed by a progressive decay to steady-state. An arterial step flow induces a different response which takes the form of a monotonically increasing pressure. Pressure and flows are nonuniform along the vessel length during such transients. In-vitro whole organ pressure-flow data are presented in the dilated rat gracilis muscle which qualitatively agree with the theoretical predictions.     

Non-Newtonian Behavior of Blood in Oscillatory Flow

Sinusoidal oscillatory flow of blood and of aqueous glycerol solutions was produced in rigid cylindrical tubes. For aqueous glycerol, the amplitude of the measured pressure gradient wave form conformed closely to that predicted by Womersley's theory of oscillatory flow, up to Reynolds numbers approaching 2000. Blood differed significantly from aqueous glycerol solutions of comparable viscosity, especially at low frequencies and high hematocrits. As frequency increased, the hydraulic impedance of blood decreased to a minimum at a frequency of about 1-2 CPS, increasing monotonically at higher frequencies. The dynamic apparent viscosity of blood, calculated from Womersley's theory, decreased with increasing flow amplitude. The reactive component of the hydraulic impedance increased with frequency as predicted by theory; the resistive component decreased with increasing frequency, differing from the resistance of a Newtonian fluid which increased with frequency.     

Low Reynolds number turbulence modeling of blood flow in arterial stenoses.

Moderate and severe arterial stenoses can produce highly disturbed flow regions with transitional and or turbulent flow characteristics. Neither laminar flow modeling nor standard two-equation models such as the kappa-epsilon turbulence ones are suitable for this kind of blood flow. In order to analyze the transitional or turbulent flow distal to an arterial stenosis, authors of this study have used the Wilcox low-Re turbulence model. Flow simulations were carried out on stenoses with 50, 75 and 86% reductions in cross-sectional area over a range of physiologically relevant Reynolds numbers. The results obtained with this low-Re turbulence model were compared with experimental measurements and with the results obtained by the standard kappa-epsilon model in terms of velocity profile, vortex length, wall shear stress, wall static pressure, and turbulence intensity. The comparisons show that results predicted by the low-Re model are in good agreement with the experimental measurements. This model accurately predicts the critical Reynolds number at which blood flow becomes transitional or turbulent distal an arterial stenosis. Most interestingly, over the Re range of laminar flow, the vortex length calculated with the low-Re model also closely matches the vortex length predicted by laminar flow modeling. In conclusion, the study strongly suggests that the proposed model is suitable for blood flow studies in certain areas of the arterial tree where both laminar and transitional/turbulent flows coexist.     

Effect of alterations in femoral artery flow on abdominal vessel hemodynamics in swine

In support of an in vivo investigation in swine of the influence of changes in fluid dynamic wall shear on arterial macromolecular permeability, a procedure has been developed to alter the flows in the porcine posterior arterial vasculature by opening and closing a reversible arteriovenous shunt placed on one of the femoral arteries. Laparoscopic techniques were used to place appropriately modified Transonic Systems ultrasonic flow probes on both external and circumflex iliac arteries, and on the terminal aorta. Flow measurements were made prior to shunt placement, and with the shunt open and closed, to measure the influence of altered external iliac artery flow on the distribution to the infrarenal abdominal vessels. Similar experiments were carried out to relate the flow rates in the external iliac arteries to those in the femoral arteries, which are more accessible. Based on the relationships among the measured flow rates, rules have been developed to estimate the major infrarenal flows in the pig, at baseline and with the shunt opened and closed, from only the flow rates measured at the two femoral arteries.     

Fluid-solid interaction in arteries incorporating the autoregulation concept in boundary conditions

In pre-surgery decisions in hospital emergency cases, fast and reliable results of the solid and fluid mechanics problems are of great interest to clinicians. In the current investigation, an iterative process based on a pressure-type boundary condition is proposed in order to reduce the computational costs of blood flow simulations in arteries, without losing control of the important clinical parameters. The incorporation of cardiovascular autoregulation, together with the well-known impedance boundary condition, forms the basis of the proposed methodology. With autoregulation, the instabilities associated with conventional pressure-type or impedance boundary conditions are avoided without an excessive increase in computational costs. The general behaviour of pulsatile blood flow in arteries, which is important from the clinical point of view, is well reproduced through this new methodology. In addition, the interaction between the blood and the arterial walls occurs via a modified weak coupling, which makes the simulation more stable and computationally efficient. Based on in vitro experiments, the hyperelastic behaviour of the wall is characterised and modelled. The applications and benefits of the proposed pressure-type boundary condition are shown in a model of an idealised aortic arch with and without an ascending aorta dissection, which is a common cardiovascular disorder.     

Development of wall surface tangent DPIV measurement techniques for arterial branch models.

Experimental techniques for measuring unsteady flow in a glass arterial bifurcation model have been developed to aid in quantifying three-dimensional wall shear fluctuations associated with arterial disease. The unique feature of the current technique is the use of a "curved" laser sheet, which was everywhere tangent to the inner wall of a daughter tube in an arterial bifurcation model. Surface tangent velocity vector field measurements were made to demonstrate the potential of this technique. Ensemble-averaged data showing weak secondary flows as well as statistical distributions of flow angles are presented. Measurements of this type may be used to estimate mean and instantaneous wall shear magnitude and direction, data that are necessary for understanding the importance of circumferential motions on arterial disease.     

Comparison of hemodynamic endpoints between normal subject and tetralogy patient using Womersley velocity profile and MR based flow measurements

Right ventricular (RV) enlargement and pulmonary valve insufficiency (PI) are well-known, unavoidable long term sequelae encountered by patients who undergo tetralogy of Fallot (TOF) surgery. Despite their lifelong need for cardiac surveillance and occasional re-intervention, there is a paucity of numerical data characterizing blood flows in their pulmonary arteries (PA). Specifically, although PA regurgitation is well-known to be ubiquitously present in adult repaired TOF (rTOF) patients yet, there have been only limited numerical studies to fully characterize this process. The few studies available have utilized idealized, simplistic geometric models or overly simplistic boundary conditions that fail to account for flow reversals near the arterial walls as observed in in-vitro and MRI based in-vivo studies. The objective of this study was to establish and validate a numerical methodology of PA blood flow using actual patient specific geometry and flow measurements obtained using phase-contrast MRI, employing Womersley type velocity profiles that model flow reversals near walls. The results from computation were validated with measurements. For the normal subject, the time averaged right PA pressure from computation (13.8 mmHg) and experiment (14.6 mmHg) differed by 6%. The time-averaged main PA pressure from computation (16.5 mmHg) and experiment (16.3 mmHg) differed by 1%. The numerically computed left PA regurgitant fraction was 89% compared to measured 77.5%, while the same for the rTOF was 43% (computation), compared to 39.6% (measured). We conclude that the use of numerical computations using the Womersley boundary condition allows reliable modeling of the pathophysiology of PA flow in rTOF.     

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 role of the surrounding tissue in the propagation of waves through the arterial system.

A theoretical analysis of the flow in arteries is presented, taking into consideration the role played by the surrounding tissues in determining the speed of propagatoion and the damping of the blood pressure pulse. This study was undertaken (a) to exhibit a method of computing the flow properties with a more nearly accurate model, (b) to see if the displacement on the skin would be related to the arterial wall displacement, and hence to pressure, velocity and flow rate of blood in the artery, and if it is likely to be measurable. It was found that the pressure of the 'viscous' part in the surrounding tissue increases the pulse velocity and the damping of the wave over the values found by other models which considered only thick-walled elastic tubes with no surrounding tissue. This study also shows that measurements on the skin can provide information about changes in arterial circulation due to diseases such as: edema, arteriosclerosis and others where the Young's modulus for either the arterial wall or the surrounding tissues is altered.     

An immersed boundary lattice Boltzmann approach to simulate deformable liquid capsules and its application to microscopic blood flows

In this paper, we develop an immersed boundary lattice Boltzmann approach to simulate deformable capsules in flows. The lattice Boltzmann method is utilized to solve the incompressible flow field over a regular Eulerian grid, while the immersed boundary method is employed to incorporate the fluid-membrane interaction with a Lagrangian representation of the capsule membrane. This algorithm was validated for the Laplace relationship, the dispersion relationship for interfacial waves and the drag coefficient for cylinders; excellent agreement with theoretical results was observed. Furthermore, simulations of single and multiple red blood cells in shear and channel flows were performed. Several characteristic hemodynamic and hemorheological features were successfully reproduced, including the tank-treading motions, cell migration from the vessel wall, slipper-shaped cell deformation, cell-free layers, blunt velocity profiles and the Fahraeus effect. These simulations therefore demonstrate the potential usefulness of this computational model for microscopic biofluidic systems. However, extension of this algorithm to three-dimensional situations is necessary for more realistic simulations.