Phasic and spatial pressure measurements in a femoral artery branch model for pulsatile flow

Phasic and spatial time-averaged pressure distributions were measured in a 60-deg femoral artery branch model over a large range of branch flow ratios and at physiological Reynolds numbers of about 120 and 700. The results obtained with an in-vivo like flow wave form indicated spatial adverse time average pressure gradients in the branch vicinity which increased in magnitude with branch flow ratio, and the importance of the larger inertial effects at the higher Reynolds numbers. Pressure losses in the branch entrance region were relatively large, and corresponding flow resistances may limit branch flow, particularly at higher Reynolds numbers. The effect of branch flow was to reduce the pressure loss in the main lumen.     

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.     

Fluid dynamic study in a femoral artery branch casting of man with upstream main lumen curvature for steady flow

An in-vitro, steady flow investigation was conducted in a hollow, transparent vascular replica of the profunda femoris branch of man for a range of physiological flow conditions. The replica casting tested was obtained from a human cadaver and indicated some plaque formation along the main lumen and branch. The flow visualization observations and measured pressure distributions indicated the highly three-dimensional flow characteristics with arterial curvature and branching, and the important role of centrifugal effects in fluid transport mechanisms.     

In-vitro pulsatile flow visualization studies in a pulmonary artery model

In-vitro pulsatile flow visualization studies were conducted in an adult-sized pulmonary artery model to observe the effects of valvular pulmonic stenosis on the flow fields of the main, left and right pulmonary arteries. The flow patterns revealed that as the degree of stenosis increased, the jet-type flow created by the valve became narrower, and it impinged on the far (distal) wall of the left pulmonary artery further downstream from the junction of the bifurcation. This in turn led to larger regions of disturbed turbulent flow, as well as helical-type secondary flow motions in the left pulmonary artery, compared to the right pulmonary artery. The flow field in the main pulmonary artery also became more disturbed and turbulent, especially during peak systole and the deceleration phase. The flow visualization observations have been valuable in helping to conduct further quantitative studies such as pressure and velocity field mapping. Such studies are important to understanding the fluid mechanics characteristics of the main pulmonary artery and its two major branches.     

Fluid particle diffusion through high-hematocrit blood flow within a capillary tube

Fluid particle diffusion through blood flow within a capillary tube is an important phenomenon to understand, especially for studies in mass transport in the microcirculation as well as in solving technical issues involved in mixing in biomedical microdevices. In this paper, the spreading of tracer particles through up to 20% hematocrit blood, flowing in a capillary tube, was studied using a confocal micro-PTV system. We tracked hundreds of particles in high-hematocrit blood and measured the radial dispersion coefficient. Results yielded significant enhancement of the particle diffusion, due to a micron-scale flow-field generated by red blood cell motions. By increasing the flow rate, the particle dispersion increased almost linearly under constant hematocrit levels. The particle dispersion also showed near linear dependency on hematocrit up to 20%. A scaling analysis of the results, on the assumption that the tracer trajectories were unbiased random walks, was shown to capture the main features of the results. The dispersion of tracer particles was about 0.7 times that of RBCs. These findings provide good insight into transport phenomena in the microcirculation and in biomedical microdevices.     

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

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

Laminar-to-turbulent transition in pulsatile flow through a stenosis

Laminar-to-turbulent transition in pulsatile flow through a stenosis is studied by means of three-dimensional numerical simulations. The flow transition is associated with the occurrence of a flow instability initiating in the stenosis region. The instability is manifested by a three-dimensional symmetry-breaking and leads to asymmetric separation and intense swirling motion downstream of the stenosis. The above have profound effects on the wall shear stress (WSS). The simulations reveal that the asymmetric separation is extended several radii downstream of the stenosis with substantial WSS fluctuations, in both space and time, occurring in the poststenotic region.     

Numerical simulation of pulsatile flow in a compliant curved tube model of a coronary artery

The endothelial cells (ECs) lining a blood vessel wall are exposed to both the wall shear stress (WSS) of blood flow and the circumferential strain (CS) of pulsing artery wall motion. These two forces and their interaction are believed to play a role in determining remodeling of the vessel wall and development of arterial disease (atherosclerosis). This study focused on the WSS and CS dynamic behavior in a compliant model of a coronary artery taking into account the curvature of the bending artery and physiological radial wall motion. A three-dimensional finite element model with transient flow and moving boundaries was set up to simulate pulsatile flow with physiological pressure and flow wave forms characteristic of the coronary arteries. The characteristic coronary artery curvature and flow conditions applied to the simulation were: aspect ratio (lambda) = 10, diameter variation (DV) = 6 percent, mean Reynolds number (Re) = 150, and unsteadiness parameter (alpha) = 3. The results show that mean WSS is about 50 percent lower on the inside wall than the outside wall while WSS oscillation is stronger on the inside wall. The stress phase angle (SPA) between CS and WSS, which characterizes the dynamics of the mechanical force pattern applied to the endothelial cell layer, shows that CS and WSS are more out of phase in the coronaries than in any other region of the circulation (-220 deg on the outside wall, -250 deg on the inside wall). This suggests that in addition to WSS, SPA may play a role in localization of coronary atherosclerosis.     

Simple pulsatile flow in an artery with a constriction

A smooth isolated, axisymmetric occlusion in a straight vascular tube is a tractable problem for pulsatile flow calculations via finite-difference approximations to the Navier-Stokes equation. Steady flow depends on the Reynolds number and two geometric parameters which describe the stenosis. The mere addition of a simple harmonic to the mean flow adds two more parameters. One is the reduced frequency, or Strokes number, and the other epsilon, the ratio of unsteady to steady flux. After describing steady stenosis flow examples, the dynamic patterns of pulsatile flow are illustrated indicating the inadequacy of basing hypotheses of atherosclerosis on mean (steady) flow.     

A fluid--structure interaction finite element analysis of pulsatile blood flow through a compliant stenotic artery.

A new model is used to analyze the fully coupled problem of pulsatile blood flow through a compliant, axisymmetric stenotic artery using the finite element method. The model uses large displacement and large strain theory for the solid, and the full Navier-Stokes equations for the fluid. The effect of increasing area reduction on fluid dynamic and structural stresses is presented. Results show that pressure drop, peak wall shear stress, and maximum principal stress in the lesion all increase dramatically as the area reduction in the stenosis is increased from 51 to 89 percent. Further reductions in stenosis cross-sectional area, however, produce relatively little additional change in these parameters due to a concomitant reduction in flow rate caused by the losses in the constriction. Inner wall hoop stretch amplitude just distal to the stenosis also increases with increasing stenosis severity, as downstream pressures are reduced to a physiological minimum. The contraction of the artery distal to the stenosis generates a significant compressive stress on the downstream shoulder of the lesion. Dynamic narrowing of the stenosis is also seen, further augmenting area constriction at times of peak flow. Pressure drop results are found to compare well to an experimentally based theoretical curve, despite the assumption of laminar flow.     

Multiphase hemodynamic simulation of pulsatile flow in a coronary artery

A multiphase transient non-Newtonian three-dimensional (3-D) computational fluid dynamics (CFD) simulation has been performed for pulsatile hemodynamics in an idealized curved section of a human coronary artery. We present the first prediction, to the authors' knowledge, of particulate buildup on the inside curvature using the multiphase theory of dense suspension hemodynamics. In this study, the particulates are red blood cells (RBCs). The location of RBC buildup on the inside curvature correlates with lower wall shear stress (WSS) relative to the outside curvature. These predictions provide insight into how blood-borne particulates interact with artery walls and hence, have relevance for understanding atherogenesis since clinical observations show that atherosclerotic plaques generally form on the inside curvatures of arteries. The buildup of RBCs on the inside curvature is driven by the secondary flow and higher residence times. The higher viscosity in the central portion of the curved vessel tends to block their flow, causing them to migrate preferentially through the boundary layer. The reason for this is the nearly neutrally buoyant nature of the dense two-phase hemodynamic flow. The two-phase non-Newtonian viscosity model predicts greater shear thinning than the single-phase non-Newtonian model. Consequently, the secondary flow induced in the curvature is weaker. The waveforms for computed hemodynamic parameters, such as hematocrit, WSS, and viscosity, follow the prescribed inlet velocity waveforms. The lower oscillatory WSS produced on the inside curvature has implications for understanding thickening of the intimal layer.     

Vortex generation in pulsatile flow through arterial bifurcation models including the human carotid artery

Visualization experiments were performed to elucidate the complicated flow pattern in pulsatile flow through arterial bifurcations. Human common carotid arteries, which were made transparent, and glass-models simulating Y- and T-shaped bifurcations were used. Pulsatile flow with wave forms similar to those of arterial flow was generated with a piston pump, elastic tube, airchamber, and valves controlling the outflow resistance. Helically recirculating flow with a pattern similar to that of the horseshoe vortex produced around wall-based protuberances in circular tubes was observed in pulsatile flow through all the bifurcations used in the present study. This flow type, which we shall refer to as the horseshoe vortex, has also been demonstrated to occur at the human common carotid bifurcation in steady flow with Reynolds numbers above 100. Time-varying flows also produced the horseshoe vortex mostly during the decelerating phase. Fluid particles of dye solution approaching the bifurcation apex diverged, divided into two directions perpendicularly, and then showed helical motion representing the horseshoe vortex formation. While this helical flow was produced, the stagnation points appeared on the wall upstream of the apex. Their position was dependent upon the flow distribution ratio between the branches in the individual arteries. The region affected by the horseshoe vortex was smaller during pulsatile flow than during steady flow. Lowering the Reynolds number together with the Womersley number weakened the intensity of helical flow. A separation bubble, resulting from the divergence or wall roughness, was observed at the outer or inner wall of the branch vessels and made the flow more complicated.     

Experimental study of pulsatile and steady flow through a smooth tube and an atherosclerotic coronary artery casting of man

In vitro investigation of pulsatile and steady flows through a smooth, straight circular tube and a diseased human coronary artery cast was conducted with sugar-water solutions simulating the viscosity of blood. Time averaged pressure drops for pulsatile flows measured in the circular tube over a Reynolds number ranging from 50 to 1,000 were found to be identical to those for steady flows in the same tube, both of which were in excellent agreement with the Poiseuille flow prediction. For the polyurethane case (# 124) made from a human main coronary with significant but 'non obstructive' diffuse atherosclerotic disease, pressure drops for steady flows were found to be greater than Poiseuille flow predictions by a factor of 3-8 in the physiological Reynolds number range from about 100 to 400. Pulsatile flows in the same artery cast resulted in additional 30% increases in time averaged pressure drops, and thus flow resistance, compared to the steady flow data. Steady and pulsatile flow data measured in a straight, axisymmetric model of cast # 124 showed considerably smaller increases in flow resistance than those observed in # 124 casting.     

Combined effects of pulsatile flow and dynamic curvature on wall shear stress in a coronary artery bifurcation model

A three-dimensional model with simplified geometry for the branched coronary artery is presented. The bifurcation is defined by an analytical intersection of two cylindrical tubes lying on a sphere that represents an idealized heart surface. The model takes into account the repetitive variation of curvature and motion to which the vessel is subject during each cardiac cycle, and also includes the phase difference between arterial motion and blood flowrate, which may be nonzero for patients with pathologies such as aortic regurgitation. An arbitrary Lagrangian Eulerian (ALE) formulation of the unsteady, incompressible, three-dimensional Navier-Stokes equations is employed to solve for the flow field, and numerical simulations are performed using the spectral/hp element method. The results indicate that the combined effect of pulsatile inflow and dynamic geometry depends strongly on the aforementioned phase difference. Specifically, the main findings of this work show that the time-variation of flowrate ratio between the two branches is minimal (less than 5%) for the simulation with phase difference angle equal to 90 degrees, and maximal (51%) for 270 degrees. In two flow pulsatile simulation cases for fixed geometry and dynamic geometry with phase angle 270 degrees, there is a local minimum of the normalized wall shear rate amplitude in the vicinity of the bifurcation, while in other simulations a local maximum is observed.     

Characteristics of secondary flow in steady and pulsatile flows through a symmetrical bifurcation

Steady and pulsatile flow in a glass model simulating an arterial bifurcation was investigated by flow visualization techniques. Secondary flow generated at the bifurcation has a similar pattern to a vortex, called the horseshoe vortex, produced around a wall-based protuberance in a circular tube. The same flow disturbance was clearly observed during the decelerating phase of pulsatile flow. The vortex produces a stagnation point on the top and bottom wall just upstream from the bifurcation apex. When aluminium dust was suspended in the test fluid perfusing the blood vessel model, particles deposited over an area spreading from the stagnation point to the lateral corners of the bifurcation. Comparison between the present results and topographical patterns of atherosclerosis reported in the literature suggests that it is in such low shear regions that lipid deposition tends to occur most.     

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.