Fluid flow and plaque formation in an aortic bifurcation

Considering steady laminar flow in a two-dimensional symmetric branching channel with local occlusions, a finite element model has been developed to study velocity fields including reverse flow regions, pressure profiles and wall shear stress distributions for different Reynolds numbers, bifurcation angles and lumen reductions. The flow analysis has been extended to include a new submodel for the pseudo-transient formation of plaque at sites and deposition rates defined by the physical characteristics of the flow. Specifically, simulating the onset of atherosclerotic lesions, sinusoidal plaque layers have been placed in areas of critically low wall shear stresses, and simulating the growth of particle depositions, plaque layers have been added in a stepwise fashion in regions of critically high and low shear. Thus two somewhat conflicting hypothetical correlations between critical wall shear stress levels and atheroma have been tested and a solution has been postulated. The validated computer simulation model is a predictive tool for analyzing the effects of local changes in wall curvature due to surgical reconstruction and/or atherosclerotic lesions, and for investigating the design of aortic bifurcations which mitigate plaque formation.     

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.     

Vessel Dimensions in Premature Atheromatous Disease of Aortic Bifurcation

The area ratio of the abdominal aortic bifurcation is defined as the ratio of the sum of the cross-sectional areas of the common iliac arteries to the cross-sectional area of the aorta. Abnormality of the area ratio is associated with initiation of atheroma of the lower aorta. Hypoplasia of the abdominal aorta predisposes to early occlusion as the atheromatous process advances.In five out of six women with thrombosis at the aortic bifurcation the area ratio was low. Hypoplasia of the abdominal aorta was present in two of these patients. All six patients were successfully treated by aortoiliac disobliteration and reconstruction.     

Anti-hypercholesterolemic and anti-atherosclerotic effects of polarized-light therapy in rabbits fed a high-cholesterol diet

The effects of polarized-light therapy (PLT) on high-cholesterol diet (HCD)-induced hypercholesterolemia and atherosclerosis were investigated in comparison with that of lovastatin in rabbits. Hypercholesterolemia was induced by feeding male New Zealand white rabbits with 1% cholesterol in diet for 2 weeks and maintained with 0.5% cholesterol for 6 weeks, followed by normal diet for 2 weeks for recovery. Lovastatin (0.002% in diet) or daily 5-min or 20-min PLT on the outside surface of ears was started 2 weeks after induction of hypercholesterolemia. Hypercholesterolemic rabbits exhibited great increases in serum cholesterol and low-density lipoproteins (LDL) levels, and finally severe atheromatous plaques formation covering 57.5% of the arterial walls. Lovastatin markedly reduced both the cholesterol and LDL, but the reducing effect (47.5%) on atheroma formation was relatively low. By comparison, 5-min PLT preferentially decreased LDL, rather than cholesterol, and thereby potentially reduced the atheroma area to 42.2%. Notably, 20-min PLT was superior to lovastatin in reducing both the cholesterol and LDL levels as well as the atheromatous plaque formation (26.4%). In contrast to the increases in blood alanine transaminase and aspartate transaminase following lovastatin treatment, PLT did not cause hepatotoxicity. In addition, PLT decreased platelets and hematocrit level. The results indicate that PLT attenuates atherosclerosis not only by lowering blood cholesterol and LDL levels, but also by improving blood flow without adverse effects. Therefore, it is suggested that PLT could be a safe alternative therapy for the improvement of hypercholesterolemia and atherosclerosis.     

Measurement of wall motion and wall shear in a compliant arterial cast

Initial measurements of the time-varying wall shear rate at two sites in a compliant cast of a human aortic bifurcation are presented. The shear rates were derived from flow velocities measured by laser Doppler velocimetry (LDV) near the moving walls of the cast. To derive these shear rate values, the distance from the velocimeter sampling volume to the cast wall must be known. The time variation of this distance was obtained from LDV measurements of the velocity of the wall itself.     

The effect of compliance on wall shear in casts of a human aortic bifurcation

Rigid and compliant casts of a human aortic bifurcation were subjected to physiologically realistic pulsatile fluid flows. At a number of sites near the wall in the approximate median plane of the bifurcation of these models, fluid velocity was measured with a laser Doppler velocimeter, and wall motion (in the case of the compliant cast) was determined with a Reticon linescan camera. The velocity and wall motion data were combined to estimate the instantaneous shear rates at the cast wall. Analysis showed that at the outer walls the cast compliance reduced shear rates, while at the walls of the flow divider the shear rate was increased.     

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.     

Some features of oscillatory flow in a model bifurcation

Oscillatory flow in the lung is studied using an order-of-magnitude analysis and flow visualization experiments in a single bifurcation with lung-like geometry. The results are used to obtain a classification scheme that identifies three major flow regimes, distinguished on the basis of whether the flow is dominated by unsteadiness, viscous effects, or the effects of convective acceleration. The unsteady regime is found to exist for values of a dimensionless stroke length (L/a, i.e., stroke volume/local cross-sectional area) less than or equal to 3 and for values of a dimensionless frequency (alpha 2 = alpha 2 omega/nu, where alpha is airway radius, omega the oscillatory frequency, and nu the kinematic viscosity) less than or equal to 10 in basic agreement with previous studies. The viscous regime is found when alpha 2(L/a)(a/R)1/2 less than 10 and alpha 2 less than 10 where R is the local radius of curvature in the bifurcation; the convective regime is found when alpha 2(L/a)(a/R)1/2 greater than 10 and L/a greater than 3. This same approach yields scaling laws for the magnitude of secondary flow velocities and shows that the ratio of secondary-to-axial velocity is small everywhere outside of the convective regime where it scales with (a/R)1/2. Comparison of these results to related simple flows shows that many of the features observed can be attributed to the effects of curvature, suggesting that the influence of the flow divider and of area change may be of lesser importance than previously thought.     

Development of a general method for designing microvascular networks using distribution of wall shear stress

In the present study, theoretical formulations for calculation of optimal bifurcation angle and relationship between the diameters of mother and daughter vessels using the power law model for non-Newtonian fluids are developed. The method is based on the distribution of wall shear stress in the mother and daughter vessels. Also, the effect of distribution of wall shear stress on the minimization of energy loss and flow resistance is considered. It is shown that constant wall shear stress in the mother and daughter vessels provides the minimum flow resistance and energy loss of biological flows. Moreover, the effects of different wall shear stresses in the mother and daughter branches, different lengths of daughter branches in the asymmetric bifurcations and non-Newtonian effect of biological fluid flows on the bifurcation angle and the relationship between the diameters of mother and daughter branches are considered. Using numerical simulations for non-Newtonian models such as power law and Carreau models, the effects of optimal bifurcation angle on the pressure drop and flow resistance of blood flow in the symmetric bifurcation are investigated. Numerical simulations show that optimal bifurcation angle decreases the pressure drop and flow resistance especially for bifurcations at large Reynolds number.     

The effects of non-Newtonian viscoelasticity and wall elasticity on flow at a 90 degrees bifurcation

To study the fundamentals of hemodynamics in arteries, the flow parameters: pulsatility, elasticity and non-Newtonian viscoelasticity were considered in detail in a 90 degrees-T-bifurcation of a rigid and elastic model. The velocity distribution 2.5 mm behind the bifurcation in the straight tube was measured with a laser-Doppler-anemometer. The fluid used was an aqueous glycerine solution and a viscoelastic Separan mixture. Flow visualization studies were done with a sheet of laser light in the plane of the bifurcation. The velocity distribution was measured for both steady and pulsatile flows with a laser-Doppler-anemometer in a backward scattered way. From the velocity measurements the shear gradients were calculated. Substantial differences were found in the flow behavior of Newtonian and non-Newtonian fluids, especially behind the bifurcation in the main tube, where secondary flows and flow separation started. Also, differences due to the elastic and rigid wall could be seen. Very high shear gradients were found in the flow between main flow and the separation zone which can lead to a damage of the blood cells.     

Turbulence modeling in three-dimensional stenosed arterial bifurcations

Under normal healthy conditions, blood flow in the carotid artery bifurcation is laminar. However, in the presence of a stenosis, the flow can become turbulent at the higher Reynolds numbers during systole. There is growing consensus that the transitional k-omega model is the best suited Reynolds averaged turbulence model for such flows. Further confirmation of this opinion is presented here by a comparison with the RNG k-epsilon model for the flow through a straight, nonbifurcating tube. Unlike similar validation studies elsewhere, no assumptions are made about the inlet profile since the full length of the experimental tube is simulated. Additionally, variations in the inflow turbulence quantities are shown to have no noticeable affect on downstream turbulence intensity, turbulent viscosity, or velocity in the k-epsilon model, whereas the velocity profiles in the transitional k-omega model show some differences due to large variations in the downstream turbulence quantities. Following this validation study, the transitional k-omega model is applied in a three-dimensional parametrically defined computer model of the carotid artery bifurcation in which the sinus bulb is manipulated to produce mild, moderate, and severe stenosis. The parametric geometry definition facilitates a powerful means for investigating the effect of local shape variation while keeping the global shape fixed. While turbulence levels are generally low in all cases considered, the mild stenosis model produces higher levels of turbulent viscosity and this is linked to relatively high values of turbulent kinetic energy and low values of the specific dissipation rate. The severe stenosis model displays stronger recirculation in the flow field with higher values of vorticity, helicity, and negative wall shear stress. The mild and moderate stenosis configurations produce similar lower levels of vorticity and helicity.     

A numerical study on hemodynamics in the left coronary bifurcation with normal and hypertension conditions

In this study, a three-dimensional analysis of the non-Newtonian blood flow was carried out in the left coronary bifurcation. The Casson model and hyperelastic and rigid models were used as the constitutive equation for blood flow and vessel wall model, respectively. Physiological conditions were considered first normal and then compliant with hypertension disease with the aim of evaluating hemodynamic parameters and a better understanding of the onset and progression of atherosclerosis plaques in the coronary artery bifurcation. Two-way fluid–structure interaction method applying a fully implicit second-order backward Euler differencing scheme has been used which is performed in the commercial code ANSYS and ANSYS CFX (version 15.0). When artery deformations and blood pressure are associated, arbitrary Lagrangian–Eulerian formulation is employed to calculate the artery domain response using the temporal blood response. As a result of bifurcation, noticeable velocity reduction and backflow formation decrease shear stress and made it oscillatory at the starting point of the LCx branch which caused the shear stress to be less than 1 and 2 Pa in the LCx and the LAD branches, respectively. Oscillatory shear index (OSI) as a hemodynamic parameter represents the increase in residence time and oscillatory wall shear stress. Because of using the ideal 3D geometry and realistic physiological conditions, the values obtained for shear stress are more accurate than the previous studies. Comparing the results of this study with previous clinical investigations shows that the regions with low wall shear stress less than 1.20 Pa and with high OSI value more than 0.3 are in more potential risk to the atherosclerosis plaque development, especially in the posterior after the bifurcation.     

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

The non-Newtonian fluid flow in a bifurcation model with a non-planar daughter branch is investigated by using finite element method to solve the three-dimensional Navier–Stokes equations coupled with a non-Newtonian constitutive model, in which the shear thinning behavior of the blood fluid is incorporated by the Carreau–Yasuda model. The objective of this study is to investigate the influence of the non-Newtonian property of fluid as well as of curvature and out-of-plane geometry in the non-planar daughter vessel on wall shear stress (WSS) and flow phenomena. In the non-planar daughter vessel, the flows are typified by the skewing of the velocity profile towards the outer wall, creating a relatively low WSS at the inner wall. In the downstream of the bifurcation, the velocity profiles are shifted towards the flow divider. The low WSS is found at the inner walls of the curvature and the lateral walls of the bifurcation. Secondary flow patterns that swirl fluid from the inner wall of curvature to the outer wall in the middle of the vessel are also well documented for the curved and bifurcating vessels. The numerical results for the non-Newtonian fluid and the Newtonian fluid with original Reynolds number and the corresponding rescaled Reynolds number are presented. Significant difference between the non-Newtonian flow and the Newtonian flow is revealed; however, reasonable agreement between the non-Newtonian flow and the rescaled Newtonian flow is found. Results of this study 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.     

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.     

Numerical and experimental models of post-operative realistic flows in stenosed coronary bypasses.

By means of both experimental and finite element methods, we simulated three-dimensional unsteady flows through coronary bypass anastomosis. The host artery includes a stenosis shape located at two different distances of grafting. The inflow rates are issued from in vivo measurements in patients who had undergone coronary bypass surgery a few days before. We provide a comparison between experimental and numerical velocity profiles coupled with the numerical analysis of spatial and temporal wall shear stress evolution. The interaction between the graft and coronary flows has been demonstrated. The phase inflow difference can partly be responsible for specific flow phenomena: jet deflection towards a preferential wall or feedback phenomenon that causes the flapping of the post-stenotic jet during the cardiac cycle. In conclusion, we showed the sensitivity of these typical flows to distance of grafting, inflows waveforms but also to their phase difference.     

An experimental analysis of the flow field in a three-dimensional model of the human carotid artery bifurcation

Steady flow measurements were carried out in a rigid three-dimensional model of the human carotid artery bifurcation at a Reynolds number of 640 and a flow division ratio of 50/50. Both axial and secondary velocities were measured with a laser-Doppler anemometer. In the bulb opposite to the flow divider a zone with negative axial velocities was found with a maximal diameter of about 60% of the local diameter of the branch and a cross-sectional extent of about 25% of the local cross-sectional area. In the bulb the maximum axial velocity shifted towards the divider wall and at the end of the bulb an axial velocity plateau arose near the non-divider wall. Halfway through the bulb, secondary flow showed a vortex through which fluid flowed towards the divider wall near the bifurcation plane and back towards the non-divider wall near the upper walls.     

Influence of pulsatility on the development of intracardiac jets: an in vitro laser Doppler study.

So far, it has been hypothesized that numerical data obtained in steady flow conditions apply to pulsatile flows. In order to study the modifications of the velocity fields due to pulsatility, jets were produced by 8 orifices (with a diameter "D" of 4.4 to 11.3 mm) included in a chamber of 50 mm. The velocity was measured using laser Doppler anemometry with a pulsatile flow ("pf") and compared to the values obtained in steady ("sf"): at maximum velocity, the longitudinal velocity profile is qualitatively similar to this observed in steady flow: it is made of a plateau followed by an hyperbolic velocity decay in the turbulent area. The length of the core ("Lpf") is strongly related to "D" (Lpf = 3.72 D + 5.49, r = .99) and the velocity decay depends on the ratio between the distance "x" from the orifice and "D" (V/Vo = 2.83D/x + 3.46, r = .85, where V is the velocity at "x" and Vo the initial velocity). During the acceleration and the deceleration, the laminar core is disturbed by turbulences. The comparison of "pf" data with "sf" data demonstrated similar diameters at the origin of the jets (Dpf = 0.96 Dsf + .12, r = .99), but significant (p less than .0001) differences both for "L" and "V/Vo": Lpf = .91Lsf + 6.58, r = .97, V/Vopf = .63 V/Vosf + .34, r = .76. Thus, pulsatility modifies velocity fields and the results obtained in steady flow conditions do not apply to pulsatile jets.     

Blood flow downstream of a two-dimensional bifurcation

Many cardiovascular lesions such as aneurysms, intimal cushions, and atherosclerotic plaques tend to occur near bifurcations. This suggests that hemodynamic factors may be involved. Since measuring devices (such as anemometers) are still too large to allow local measurements of flow disturbances, we have attempted to predict the nature of these factors mathematically. Biological variables include pulsatile flow of a nonNewtonian fluid in distensible branching vessels with different angles and flow rates. Our initial analysis considers the flow in a two-dimensional bifurcation with a symmetrical flow divider perfused with steady flow at variable Reynolds numbers. At all flows, high shear forces develop on either side of the flow divider (i.e. at the apex of the bifurcation). With high flows, regions of sluggish or reverse flow develop near the outer walls of the bifurcation. The analysis confirms that the flow at the apex is quite different from that at the outer angles and that the latter varies more with flow rate than the former.