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.     

Comparative study of flow in right-sided and left-sided aortas: numerical simulations in patient-based models

A right-sided aorta is a rare malformation which may be associated with other various types of congenital heart disease. We utilised haemodynamic, echocardiographic measurements, computerised tomography and image reconstruction software packages that were integrated in a computational fluid dynamics model to determine blood flow patterns in patient-based aortas. In the left-sided aorta, a systolic clockwise rotational component was present, while helical flow was depicted in the aortic arch that was converted in the descending aorta as counter-rotating vortices with accompanying retrograde flow. The right-sided configuration has not altered the orientation of the three-dimensional vortex, but intensification of polymorphic flow patterns, alterations in wall shear stress distribution and development of a lateral pressure gradient at the area of an aneurysmal anomaly was observed. Moreover, increments of Reynolds, Womersley and Dean numbers were evident. These phenomena along with the formation of the aneurysm might influence cardiovascular risk in patients with right-sided aortas.     

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.     

Pulsatile flow visualization in the abdominal aorta under differing physiologic conditions: implications for increased susceptibility to atherosclerosis.

The infrarenal abdominal aorta is a common site for clinically significant atherosclerosis. As has been shown in other susceptible locations, vessel geometry, flow division rates, and pulsatility may result in hemodynamic conditions which influence the preferential localization of disease in the abdominal aorta segment. Pulsatile flow visualization was performed in a glass model of the aorta constructed from measurements of angiograms and cadaver aortas. Flow rates and pulsatile waveforms were varied to reflect typical physiological conditions. Under normal resting conditions, the flow patterns in the infrarenal aorta were more complex than those in the suprarenal location. Time varying vortex patterns appeared at the level of the renal arteries and propagated through the infrarenal aorta into the common iliac arteries. A region of oscillating velocity direction extended from the renal arteries to the aortic bifurcation along the posterior wall. Dye became trapped along the posterior wall, requiring several cardiac cycles for clearance. In contrast, there was rapid clearance of the dye in the anterior aorta. Under postprandial conditions, the flow patterns in the aorta were basically unchanged. Simulated exercise conditions created laminar hemodynamic features very different from the resting conditions, including a decrease in dye residence time. This study reveals significant time-dependent variations in the hemodynamics of the abdominal aorta under differing physiologic conditions. Hemodynamic factors such as low wall shear stress, oscillating shear direction, and high particle residence time may be related to the clinically seen preferential plaque localization in the infrarenal aorta.     

Velocity distribution along an elastic model of human arterial tree

An experimental investigation of an elastic model of the human arterial tree, has been performed for physiological type flow by pulsed Doppler ultrasonic velocimetry. The arterial tree model, fabricated in clear polyurethane, includes the aortic arch, with a Starr-Edwards ball valve mounted in the root of the aorta, the descending aorta and the iliac bifurcation. Our study showed that the velocity profile, a few centimeters beyond the valve, is skewed, with higher velocities towards the top and the inner wall (anatomically the posterior and left lateral wall). An inward shift of the maximum velocity and reverse flow are denoted along the inner wall of the aortic arch. The velocity profiles in the descending aorta are blunted. Downstream from the vertex of the iliac bifurcation, there is vorticity creation, but the branching effect is quickly damped by the pulsatility of the flow and the elasticity of the wall.     

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.     

Two-component laser velocimeter measurements downstream of heart valve prostheses in pulsatile flow

Elevated turbulent shear stresses resulting from disturbed blood flow through prosthetic heart valves can cause damage to red blood cells and platelets. The purpose of this study was to measure the turbulent shear stresses occurring downstream of aortic prosthetic valves during in-vitro pulsatile flow. By matching the indices of refraction of the blood analog fluid and model aorta, correlated, simultaneous two-component laser velocimeter measurements of the axial and radial velocity components were made immediately downstream of two aortic prosthetic valves. Velocity data were ensemble averaged over 200 or more cycles for a 15-ms window opened at peak systolic flow. The systolic duration for cardiac flows of 8.4 L/min was 200 ms. Ensemble-averaged total shear stress levels of 2820 dynes/cm2 and 2070 dynes/cm2 were found downstream of a trileaflet valve and a tilting disk valve, respectively. These shear stress levels decreased with axial distance downstream much faster for the tilting disk valve than for the trileaflet valve.     

Topographical mapping of sites of enhanced HRP permeability in the normal rabbit aorta.

The spatial distribution of sites of enhanced permeability to the macromolecule horseradish peroxidase (HRP) in the normal rabbit aorta after one min circulation was studied using image analysis. These sites, referred to as "HRP spots," exhibit a nonuniform distribution that is qualitatively similar in all rabbits studied. The density of HRP spots is highest in the aortic arch, decreases distally, reaches a minimum in the lower descending thoracic aorta, and then increases again in the abdominal aorta. The region of highest spot density follows a clockwise helical pattern in the aortic arch and outside the arch occurs in streaks largely oriented in the bulk flow direction. The streaks in the abdominal aorta localize along the anatomical right lateral wall and occasionally along the left lateral wall proximal to the celiac artery and along the ventral wall between the celiac and superior mesenteric arteries. The density of spots is high in the immediate vicinity of aortic ostia with the most elevated density being distal to ostia in most cases. At a short distance from the ostium edge of the celiac and superior mesenteric branches the proximal density is comparably high, and no preferred spot orientation is observed around the brachiocephalic vessel. These results are consistent with an influence of localizing factors such as detailed hemodynamic phenomena and/or arterial wall structural and/or functional variations.     

Some flow visualization and laser-Doppler-velocity measurements in a true-to-scale elastic model of a human aortic arch--a new model technique.

Flow studies were done in an elastic true-to-scale silicone rubber model of an aortic arch to study further hemodynamic influences on atherosclerosis. The model was prepared from a cast of a young woman. A revised model technique was used. The model had a compliance similar to that of the human aortic arch. Velocity measurements were done in the model with a two component laser-Doppler-anemometer in steady and pulsatile flow using a calcium chloride solution with a viscosity of eta = 3.18 mPas and density of rho = 1.28 kg/m3 at 20 degrees C. The time average Reynolds numbers over a whole cycle in the ascending aorta was Re = 1350. The Womersley parameter for pulsatile flow was a = 20. The pulse wave velocity in the ascending aorta was about c = 5.4 m/sec. The secondary flow behavior was discussed for steady and pulsatile flow. Reverse flows were found, especially along the inner radius of the aortic arch in the descending aorta in steady and pulsatile flow and also in small areas of the ascending aorta and at the branches of the aortic arch. The formation of atherosclerotic plaques at preferred local flow regions is discussed.     

Comparison of physiological and post-endovascular aneurysm repair infrarenal blood flow

Endovascular aneurysm repair (EVAR) of abdominal aortic aneurysms results in redirection of blood through the deployed endograft (EG). Even though EVAR is clinically effective, the absolute flow restoration is not warranted. Our purpose was to compare the physiological with the post-EVAR infrarenal flow conditions. We developed patient-specific models based on computed tomography data of five healthy volunteers and ten patients treated with the Endurant^® stent-graft system. Wall shear stress (WSS), helicity, pressure and velocity fields were calculated using computational fluid dynamics. The results showed a decrease of peak WSS on the part of the EG that resides in the iliac arteries, compared to the physiological value (p = 0.01). At the abdominal part, the average helicity seems to increase after EVAR, while at the iliac arteries part, the intensity of helical flow seems physiological. Pressure drop and peak velocity in the iliac arteries part are lower than the physiological values (p = 0.04). The comparison revealed that most hemodynamic properties converge to normal levels at the abdominal part whereas statistically significant variations were observed in the iliac arteries part. The delineation of the differences between physiological and postoperative flow data could pave the way for the improvement of EG designs.     

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

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

The effect of blood viscoelasticity on pulsatile flow in stationary and axially moving tubes

An analytical solution for pulsatile flow of a generalized Maxwell fluid in straight rigid tubes, with and without axial vessel motion, has been used to calculate the effect of blood viscoelasticity on velocity profiles and shear stress in flows representative of those in the large arteries. Measured bulk flow rate Q waveforms were used as starting points in the calculations for the aorta and femoral arteries, from which axial pressure gradient ▿P waves were derived that would reproduce the starting Q waves for viscoelastic flow. The ▿P waves were then used to calculate velocity profiles for both viscoelastic and purely viscous flow. For the coronary artery, published ▿P and axial vessel acceleration waveforms were used in a similar procedure to determine the separate and combined influences of viscoelasticity and vessel motion.Differences in local velocities, comparing viscous flow to viscoelastic flow, were in all cases less than about 2% of the peak local velocity. Differences in peak wall shear stress were less than about 3%.In the coronary artery, wall shear stress differences between viscous and viscoelastic flow were small, regardless of whether axial vessel motion was included. The shape of the wall shear stress waveform and its difference, however, changed dramatically between the stationary and moving vessel cases. The peaks in wall shear stress difference corresponded with large temporal gradients in the combined driving force for the flow.     

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.     

Wave intensity in the ascending aorta: effects of arterial occlusion

We examine the effects of arterial occlusion on the pressure, velocity and the reflected waves in the ascending aorta using wave intensity analysis. In 11 anaesthetised, open-chested dogs, snares were used to produce total arterial occlusion at 4 sites: the upper descending aorta at the level of the aortic valve (thoracic); the lower thoracic aorta at the level of the diaphragm (diaphragm); the abdominal aorta between the renal arteries (abdominal) and the left iliac artery, 2 cm downstream from the aorta iliac bifurcation (iliac). Pressure and flow in the ascending aorta were measured, and data were collected before and during the occlusion. During thoracic and diaphragm occlusions a significant increase in mean aortic pressure (46% and 23%) and in wave speed (25% and 10%) was observed, while mean flow rate decreased significantly (23% and 17%). Also, the reflected compression wave arrived significantly earlier (45% and 15%) and its peak intensity was significantly greater (257% and 125%), all compared with control. Aortic occlusion distal to the renal arteries, however, caused an indiscernible change in the pressure and velocity waveforms, and in the intensities and timing of the waves in the forward and backward directions. The measured pressure and velocity waveforms are the result of the interaction between the heart and the arterial system. The separated pressure, velocity and wave intensity are required to provide information about arterial hemodynamic such as the timing and magnitude of the forward and backward waves. The net wave intensity is simpler to calculate but provides information only about the predominant direction of the waves and can be misleading when forward and backward waves of comparable magnitudes are present simultaneously.     

Three patterns of distribution characterize the organization of endothelial microfilaments at aortic flow dividers

Since actin microfilaments are essential in the maintenance of endothelial integrity and in the repair of injured endothelium, we have carried out a detailed study of the distribution of microfilaments in the immediate vicinity of aortic branches. Branches are of major interest because there is a predilection for atherosclerotic lesions near branch ostia. We made an extensive, systematic examination of branches of the aorta and iliac arteries using in situ staining of perfusion-fixed arteries. Microfilaments were localized using rhodamine phalloidin. Three patterns of staining were observed. Some endothelial cells showed prominent central stress fibers. Others had few central stress fibers but prominent peripheral fibers. Still others showed an intermediate pattern with some central and some peripheral fibers present. At small branch sites, the lip of the divider was more blunt, and there were more cells with peripheral actin. At large branches, cells with peripheral actin were confined mainly to the lip, while there were many more cells with prominent central fibers. We also found that major differences can occur over very small distances, so adjacent cells may have strikingly different patterns of microfilament distribution. These patterns appear to reflect the geometry of the flow divider and local variations in hemodynamic shear stress. The differences in microfilament distribution may reflect differences in endothelial functions which are essential in maintaining endothelial integrity.     

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.     

Pulsatile non-Newtonian haemodynamics in a 3D bifurcating abdominal aortic aneurysm model

Numerical prediction of non-Newtonian blood flow in a 3D abdominal aortic aneurysm bifurcating model is carried out. The non-Newtonian Carreau model is used to characterise the shear thinning behaviour of the human blood. A physical inlet velocity waveform incorporating a radial velocity distribution reasonably representative of a practical case configuration is employed. Case studies subject to both equal and unequal outlet pressures at iliac bifurcations are presented to display convincingly the downstream pressure influences on the flow behaviour within the aneurysm. Simulations indicate that the non-Newtonian aspects of the blood cannot at all be neglected or given a cursory treatment. The wall shear stress (WSS) is found to change significantly at both the proximal and distal ends of the aneurysm. At the peak systole, the WSS is peak around the bifurcation point, whereas the WSS becomes zero in the bifurcation point. Differential downstream pressure fields display significant effects regarding the flow evolution in the iliac arteries, whereas little or no effects are observed directly on the flow details in the aneurysm.     

Pulsatile non-Newtonian haemodynamics in a 3D bifurcating abdominal aortic aneurysm model

Numerical prediction of non-Newtonian blood flow in a 3D abdominal aortic aneurysm bifurcating model is carried out. The non-Newtonian Carreau model is used to characterise the shear thinning behaviour of the human blood. A physical inlet velocity waveform incorporating a radial velocity distribution reasonably representative of a practical case configuration is employed. Case studies subject to both equal and unequal outlet pressures at iliac bifurcations are presented to display convincingly the downstream pressure influences on the flow behaviour within the aneurysm. Simulations indicate that the non-Newtonian aspects of the blood cannot at all be neglected or given a cursory treatment. The wall shear stress (WSS) is found to change significantly at both the proximal and distal ends of the aneurysm. At the peak systole, the WSS is peak around the bifurcation point, whereas the WSS becomes zero in the bifurcation point. Differential downstream pressure fields display significant effects regarding the flow evolution in the iliac arteries, whereas little or no effects are observed directly on the flow details in the aneurysm.     

Numerical simulation of steady flow fields in a model of abdominal aorta with its peripheral branches

In the present study, a numerical calculation procedure based on a finite volume method was developed to simulate steady flow fields in a model of abdominal aorta with its peripheral branches. The study focused on the steady baseline flow fields and the wall shear stress (WSS) distribution as well as the localization of the reversed flow regions and results were compared to those obtained by other investigators. In the case of resting conditions, the existence of a region of reversed flow of about one to two diameters in size and next to the renal arteries and along the posterior wall as observed by other researchers was confirmed. However, under the exercise conditions this region could be wiped out. The flow reversal along the lateral walls proximal to the bifurcation persisted in both rest and exercise conditions. The WSS distribution and the wall shear stress gradient distribution were obtained. The lowest WSS occurred near the ostia of the renal arteries and the lateral walls of the iliac arteries. And the highest is always at the turn to the branch. The results were generally consistent with those obtained experimentally and numerically by other investigators. It was also shown that the steady flow might be used to depict the averaged behavior of pulsatile flow. The present computer code provides a platform for the future more realistic simulations.     

Hemodynamic analysis of a novel stent graft design with slit perforations in thoracic aortic aneurysm

Thoracic endovascular aortic repair (TEVAR) has been introduced as a less invasive approach to the treatment of thoracic aortic aneurysm (TAA). However, the effectiveness of TEVAR in the treatment of TAA is often limited due to the complex anatomy of aortic arch. Flow preservation at the three supra-aortic branches further increases the overall technical difficulty. This study proposes a novel stent graft design with slit perforations that can positively alter the hemodynamics at the aortic arch while maintaining blood flow to supra-aortic branches. We carried out a computational fluid dynamic (CFD) analysis to evaluate flow characteristics near stented aortic arch in simplified TAA models, followed by in-vitro experiments using particle image velocimetry (PIV) in a mock circulatory loop. The hemodynamics result was studied in terms of time-averaged wall shear stress (TAWSS), oscillating shear index (OSI), and endothelial cell action potential (ECAP). The results showed that the stent graft with slit perforations can reduce the disturbed flow region considerably. Furthermore, the effect of the slits on flow preservation to the supra-aortic branches was simulated and compared with experimental results. The effectiveness of the stent graft with slit perforations in preserving flow to the branches was demonstrated by both simulated and experimental results. Low TAWSS and elevated ECAP were observed in the aortic arch aneurysm after the placement of the stent graft with slits, implying the potential of thrombus formation in the aneurysm. On the other hand, the effects of the stent grafts with full-slit design and half-slit design on the shear stress did not differ significantly. The present analysis indicated that not only could the stent graft with slit perforations shield the aneurysm from rupture, but also it resulted in a favorable environment for thrombus that can contribute to the shrinkage of the aneurysm.