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.     

Experimental investigation of steady flow through a model of the human aortic arch

Steady flow through a model of the human aortic arch has been studied with hot-film anemometry. A three sensor hot-film velocity probe was inserted into an acrylic flow chamber fabricated from the in situ casting of a human aorta, and the axial, radial and tangential velocity profiles were determined for steady flows in the region of the aortic arch. These studies demonstrated the presence of a potential core throughout the arch region, with a concomitant boundary layer adjacent to the inner wall of curvature of the arch. Trapped secondary flows in this fluid layer along the inner wall were quantitatively determined. Our steady flow studies in the model human aortic arch suggests that a shear-dependent mass transfer mechanism may play a significant role in the development and propagation of atherosclerotic lesions in this segment of the human cardiovascular system.     

Measurement of wall shear stress distal to a tri-leaflet valve in a rigid model of the aortic arch with branch flows

Flush mounted hot film anemometer probes were used to measure wall shear stress magnitudes on the inside and outside walls of a rigid model of the human aortic arch. The effects of the presence of an Ionescu-Shiley tri-leaflet bioprosthetic heart valve at the entrance of the aortic arch and the side flows through arteries located in the mid-arch region on wall shear stress magnitudes were determined. It was found that the presence of the tri-leaflet valve leads to an elevation of wall shear stress (relative to the same flow without a valve) over the entire aortic arch region by as much as 50 percent. The valve influence extended to about 180 deg from the entrance to the aorta on the inside wall and even further on the outside wall based on extrapolation of available data. Peak wall shear stress magnitudes measured on the outside wall were in the range of 1.5-4.0 N/m2 (15-40 dynes/cm2) over the length of the aortic arch and took on their highest values in the mid-arch region. Inside wall values were of comparable magnitude. It was observed that the presence of the aortic valve and side flow from the top of the aortic arch reduced wall shear stress reversal in the arch region.     

Unsteady and three-dimensional simulation of blood flow in the human aortic arch

A three-dimensional and pulsatile blood flow in a human aortic arch and its three major branches has been studied numerically for a peak Reynolds number of 2500 and a frequency (or Womersley) parameter of 10. The simulation geometry was derived from the three-dimensional reconstruction of a series of two-dimensional slices obtained in vivo using CAT scan imaging on a human aorta. The numerical simulations were obtained using a projection method, and a finite-volume formulation of the Navier-Stokes equations was used on a system of overset grids. Our results demonstrate that the primary flow velocity is skewed towards the inner aortic wall in the ascending aorta, but this skewness shifts to the outer wall in the descending thoracic aorta. Within the arch branches, the flow velocities were skewed to the distal walls with flow reversal along the proximal walls. Extensive secondary flow motion was observed in the aorta, and the structure of these secondary flows was influenced considerably by the presence of the branches. Within the aorta, wall shear stresses were highly dynamic, but were generally high along the outer wall in the vicinity of the branches and low along the inner wall, particularly in the descending thoracic aorta. Within the branches, the shear stresses were considerably higher along the distal walls than along the proximal walls. Wall pressure was low along the inner aortic wall and high around the branches and along the outer wall in the ascending thoracic aorta. Comparison of our numerical results with the localization of early atherosclerotic lesions broadly suggests preferential development of these lesions in regions of extrema (either maxima or minima) in wall shear stress and pressure.     

Experimental study of physiological pulsatile flow past valve prosthesis in a model of human aorta--II. Tilting disc valves and the effect of orientation

In Part II of this two paper sequence, pulsatile flow development past a tilting disc valve in a model human aorta has been studied using quantitative laser Doppler techniques. The valve was mounted in three different orientations with respect to the aortic root in this study. Under pulsatile flow, the region of flow reversal induced near the wall of the minor flow orifice extends to more than one tissue annulus diameter downstream from the valve into the ascending aorta. In a plane perpendicular to the tilt axis, a bi-helical secondary flow is induced distal to the valve. This secondary flow is further compounded by the multiple curvatures in the aorta. Hence the valve orientation affects the velocity profiles as far downstream as the mid-arch region as well as in the brachio-cephalic arterial branch. In the mid-arch region, a flow reversal along the entire cross-section is observed in early diastole for all the three orientations of the disc valve.     

Computational analysis of blood flow in an integrated model of the left ventricle and the aorta

To study the effects of intraventricular flow dynamics on the aortic flow, we created an integrated model of the left ventricle and aorta and conducted a computer simulation of diastolic and systolic blood flow within this model. The results demonstrated that the velocity profile at the aortic annulus changed dynamically, and was influenced by the intraventricular flow dynamics. The profile was almost flat in early systole but became nonuniform as systole progressed, and was skewed toward the posterior side in midsystole and toward the anterior side in later systole. At a distance from the aortic annulus, a different velocity profile was induced by the twisting and torsion of the aorta. In the ascending aorta, the fastest flow was initially located in the posteromedial sector, and it moved to the posterior section along the circumference as systole progressed. The nonuniformity of the aortic inflow gave rise to a complex wall shear stress (WSS) distribution in the aorta. A comparison of the WSS distribution obtained in this integrated analysis with that obtained in flow calculations using an isolated aorta model with Poiseuille and flat inlet conditions showed that intraventricular flow affected the WSS distribution in the ascending aorta. These results address the importance of an integrated analysis of flow in the left ventricle and aorta.     

Systolic fluid–structure interaction model of the congenitally bicuspid aortic valve: assessment of modelling requirements

A transient fluid–structure interaction (FSI) model of a congenitally bicuspid aortic valve has been developed which allows simultaneous calculation of fluid flow and structural deformation. The valve is modelled during the systolic phase (the stage when blood pressure is elevated within the heart to pump blood to the body). The geometry was simplified to represent the bicuspid aortic valve in two dimensions. A congenital bicuspid valve is compared within the aortic root only and within the aortic arch. Symmetric and asymmetric cusps were simulated, along with differences in mechanical properties. A moving arbitrary Lagrange–Euler mesh was used to allow FSI. The FSI model requires blood flow to induce valve opening and induced strains in the region of 10%. It was determined that bicuspid aortic valve simulations required the inclusion of the ascending aorta and aortic arch. The flow patterns developed were sensitive to cusp asymmetry and differences in mechanical properties. Stiffening of the valve amplified peak velocities, and recirculation which developed in the ascending aorta. Model predictions demonstrate the need to take into account the category, including any existing cusp asymmetry, of a congenital bicuspid aortic valve when simulating its fluid flow and mechanics.     

In vitro hemodynamic investigation of the embryonic aortic arch at late gestation

This study focuses on the dynamic flow through the fetal aortic arch driven by the concurrent action of right and left ventricles. We created a parametric pulsatile computational fluid dynamics (CFD) model of the fetal aortic junction with physiologic vessel geometries. To gain a better biophysical understanding, an in vitro experimental fetal flow loop for flow visualization was constructed for identical CFD conditions. CFD and in vitro experimental results were comparable. Swirling flow during the acceleration phase of the cardiac cycle and unidirectional flow following mid-deceleration phase were observed in pulmonary arteries (PA), head-neck vessels, and descending aorta. Right-to-left (oxygenated) blood flowed through the ductus arteriosus (DA) posterior relative to the antegrade left ventricular outflow tract (LVOT) stream and resembled jet flow. LVOT and right ventricular outflow tract flow mixing had not completed until approximately 3.5 descending aorta diameters downstream of the DA insertion into the aortic arch. Normal arch model flow patterns were then compared to flow patterns of four common congenital heart malformations that include aortic arch anomalies. Weak oscillatory reversing flow through the DA junction was observed only for the Tetralogy of Fallot configuration. PA and hypoplastic left heart syndrome configurations demonstrated complex, abnormal flow patterns in the PAs and head-neck vessels. Aortic coarctation resulted in large-scale recirculating flow in the aortic arch proximal to the DA. Intravascular flow patterns spatially correlated with abnormal vascular structures consistent with the paradigm that abnormal intravascular flow patterns associated with congenital heart disease influence vascular growth and function.     

Wall shear stress distribution in a model human aortic arch: assessment by an electrochemical technique

Wall shear stress (WSS) distribution in a human aortic arch model is studied using 130 cathode electrodes flush-mounted on the model walls. Flow visualizations are made in a transparent geometry model to identify the regions of fluid mechanical interests, e.g. regions of flow separation, eddy formation and flow stagnancy. The 130 electrodes are strategically positioned in the arch based on information obtained from the flow visualizations. The measured data indicate that the aortic arch may be categorized into eight regions: three along the inner wall of the arch (A,B,C); and five near the outer wall (D,E,F,G,H). (1) The regions of low WSS are distributed along the inner wall of the ascending aorta A; the inner wall of the descending aorta C; and the upstream inner wall of the innominate and the common carotid branchings F. (2) The high WSS regions are distributed along the outer wall of the arch E; and the inner wall in the arch opposite to the left subclavian branching B. (3) In certain regions, high and low WSS may be found next to each other (e.g. G and H) without a definable boundary in between; and (4) as the Reynolds number increases, the areas of low WSS decrease, while the high WSS areas increase with no obvious change in magnitude of the stress along the inner wall of the arch. At the branchings, the WSS distribution is not affected by the Reynolds number within the range of observations. The measured WSS distribution is compared with Rodkiewicz's map of early atherosclerotic lesions in the aortic arch of cholesterol fed rabbits.     

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

Background  

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

Steady flow development past valve prostheses in a model human aorta—II. Tilting disc valves

The flow development in the model human aorta with uniform entry as well as with centrally occuluding valves mounted at the root of the aorta was described in Part I of this two-paper sequence. Part II deals with the flow development in the model aorta with tilting disc valves mounted at the root of the aorta. Bjork-Shiley and Hall-Kaster tilting disc valves were mounted in three different orientations with respect to the root of the aorta. The velocity profiles and turbulent stresses were measured with laser-Doppler anemometry. Our results under steady flow conditions in the model human aorta show quantitatively that the flow development in the ascending aorta as well as in the brachio-cephalic artery are strongly dependent on the orientation of the tilting disc valves. With the valves tilting towards the outer wall of curvature, our results suggest a tendency for flow separation at the flow divider region of the brachio-cephalic artery.     

Pulsatile flow during cardiopulmonary bypass preserves postoperative microcirculatory perfusion irrespective of systemic hemodynamics

The onset of nonpulsatile cardiopulmonary bypass is known to deteriorate microcirculatory perfusion, but it has never been investigated whether this may be prevented by restoration of pulsatility during extracorporeal circulation. We therefore investigated the distinct effects of nonpulsatile and pulsatile flow on microcirculatory perfusion during on-pump cardiac surgery. Patients undergoing coronary artery bypass graft surgery were randomized into a nonpulsatile (n = 17) or pulsatile (n = 16) cardiopulmonary bypass group. Sublingual mucosal microvascular perfusion was measured at distinct perioperative time intervals using sidestream dark field imaging, and quantified as the level of perfused small vessel density and microvascular flow index (vessel diameter < 20 μm). Microcirculation measurements were paralleled by hemodynamic and free hemoglobin analyses. The pulse wave during pulsatile bypass estimated 58 ± 17% of the baseline blood pressure waveform. The observed reduction in perfused vessel density during aorta cross-clamping was only restored in the pulsatile flow group and increased from 15.5 ± 2.4 to 20.3 ± 3.7 mm/mm(2) upon intensive care admission (P < 0.01). The median postoperative microvascular flow index was higher in the pulsatile group [2.6 (2.5-2.9)] than in the nonpulsatile group [2.1 (1.7-2.5); P = 0.001]. Pulsatile flow was not associated with augmentation of free hemoglobin production and was paralleled by improved oxygen consumption from 70 ± 14 to 82 ± 16 ml·min(-1)·m(-2) (P = 0.01) at the end of aortic cross-clamping. In conclusion, pulsatile cardiopulmonary bypass preserves microcirculatory perfusion throughout the early postoperative period, irrespective of systemic hemodynamics. This observation is paralleled by an increase in oxygen consumption during pulsatile flow, which may hint toward decreased microcirculatory heterogeneity during extracorporeal circulation and preservation of microcirculatory perfusion throughout the perioperative period.     

Comparison of hinge microflow fields of bileaflet mechanical heart valves implanted in different sinus shape and downstream geometry

The characterization of the bileaflet mechanical heart valves (BMHVs) hinge microflow fields is a crucial step in heart valve engineering. Earlier in vitro studies of BMHV hinge flow at the aorta position in idealized straight pipes have shown that the aortic sinus shapes and sizes may have a direct impact on hinge microflow fields. In this paper, we used a numerical study to look at how different aortic sinus shapes, the downstream aortic arch geometry, and the location of the hinge recess can influence the flow fields in the hinge regions. Two geometric models for sinus were investigated: a simplified axisymmetric sinus and an idealized three-sinus aortic root model, with two different downstream geometries: a straight pipe and a simplified curved aortic arch. The flow fields of a 29-mm St Jude Medical BMHV with its four hinges were investigated. The simulations were performed throughout the entire cardiac cycle. At peak systole, recirculating flows were observed in curved downsteam aortic arch unlike in straight downstream pipe. Highly complex three-dimensional leakage flow through the hinge gap was observed in the simulation results during early diastole with the highest velocity at 4.7 m/s, whose intensity decreased toward late diastole. Also, elevated wall shear stresses were observed in the ventricular regions of the hinge recess with the highest recorded at 1.65 kPa. Different flow patterns were observed between the hinge regions in straight pipe and curved aortic arch models. We compared the four hinge regions at peak systole in an aortic arch downstream model and found that each individual hinge did not vary much in terms of the leakage flow rate through the valves.     

Steady flow development past valve prostheses in a model human aorta—I. Centrally occluding valves

In this paper, laser-Doppler anemometry measurement of steady flow development in a model human aorta has been reported. Studies were made with uniform entry flow at the root of the aorta and our measurements showed the establishment of a pair of Dean vortices in the mid-arch region. Subsequently, the nature of flow development past centrally occluding caged ball valves in the model aorta was investigated. Our studies showed that in the ascending aorta, an asymmetric velocity profile is obtained with larger velocity gradients towards the inner wall of tertiary curvature (anatomically the left lateral wall) with centrally occluding valves. The peripheral flow past these valves prevented the development of Dean vortices in the mid-arch region. The caged ball valves at the root of the aorta had no discernible effect on the velocity profiles in the brachio-cephalic artery.     

Pulsatile flow past aortic valve bioprostheses in a model human aorta

Pulsatile flow development past tissue valve prostheses in a model human aorta has been studied using qualitative flow visualization and quantitative laser-Doppler techniques. Experiments were conducted both in steady and physiological pulsatile flow situations and the measurements included the pressure drop across the valve, the instantaneous flow rate as well as the velocity profiles and turbulent stresses downstream to the valves. Our study shows that the velocity profiles with pericardial valves are closer to those measured past natural aortic valves. The porcine valves with a smaller valve opening area produce a narrower and stronger jet downstream from the valve with relatively larger turbulent axial stresses in the boundary of the jet. Our study suggests that the pericardial valves with turbulent stresses comparable to those of caged ball and tilting disc valves are preferable from a hemodynamic point of view.     

Steady flow visualization in a rigid canine aortic cast

Steady flow studies were conducted in a transparent canine aortic cast. The cast segment stretched from the aortic valve to beyond the renal arteries and included all major branches. Flow was visualized by analysis of dye streaklines. Flow rates for basal and exercising cardiovascular states were simulated. The Reynolds numbers in the ascending aorta for basal and exercising conditions were 900 and 1587 respectively. Aortic core flow was laminar in basal simulations. Disturbed flow commenced in the upper descending aorta with exercising flow rates. Separation zones existed along the inner curvature of the aortic arch and the proximal walls of the brachiocephalic, left subclavian, and coeliac arteries. Such zones may exist over a portion of the cardiac cycle. If either renal artery was occluded, then a vortex formed. This vortex is associated with high shear regions which correlate well with sites where sudanophilic lesions have been reported in cholesterol-fed nephrectomized rabbits.     

Local hemodynamics affect monocytic cell adhesion to a three-dimensional flow model coated with E-selectin.

Monocyte adhesion to the endothelium depends on concentrations of receptors/ligands, local concentrations of chemoattractants, monocyte transport to the endothelial surface and hemodynamic forces. Monocyte adhesion to the inert surface of a three-dimensional perfusion model was shown to correlate inversely with wall shear stress, but was also affected by flow patterns which influenced the near-wall cell availability. We hypothesized that (a) under the same flow conditions, insolubilized E-selectin on the model's surface may mediate adhesive interactions at higher wall shear stresses, compared to an uncoated model, and (b) pulsatile flow may modify the adhesion profile obtained under steady flow. An axisymmetric flow model with a stenosis and a sudden expansion produced a range of wall shear stresses and a separated flow with recirculation and reattachment. Pre-activated U937 cells were perfused through the model under either steady (Re = 100, 140) or pulsatile (Remean = 107) flow. The velocity field was characterized through computational fluid dynamics and validated by inert particle tracking. Surface E-selectin greatly increased cell adhesion in all regions at Re = 100 and 140, compared to an uncoated model under the same flow conditions. In regions where the cells near the wall were abundant (taper and stenosis), adhesion to E-selectin correlated with the reciprocal of local wall shear stress when flow was steady. Pulsatile flow distributed the adherent cells more evenly throughout the coated model. Hence, characterizing both the local hemodynamics and the biological activity on the vessel wall is important in leukocyte adhesion.     

Estimating absolute aortic pressure using MRI and a one-dimensional model

Aortic blood pressure is a strong indicator to cardiovascular diseases and morbidity. Clinically, pressure measurements are done by inserting a catheter in the aorta. However, imaging techniques have been used to avoid the invasive procedure of catheterization. In this paper, we combined MRI measurements to a one-dimensional model in order to simulate blood flow in an aortic segment. Absolute pressure was estimated in the aorta by using MRI measured flow as boundary conditions and MRI measured compliance as a pressure law for solving the model. Model computed pressure was compared to catheter measured pressure in an aortic phantom. Furthermore, aortic pressure was estimated in vivo in three healthy volunteers.