Effects of arterial wall models and measurement uncertainties on cardiovascular model predictions

We developed a methodology to assess and compare the prediction quality of cardiovascular models for patient-specific simulations calibrated with uncertainty-hampered measurements. The methodology was applied in a one-dimensional blood flow model to estimate the impact of measurement uncertainty in wall model parameters on the predictions of pressure and flow in an arterial network. We assessed the prediction quality of three wall models that have been widely used in one-dimensional blood flow simulations. A 37-artery network, previously used in one experimental and several simulation studies, was adapted to patient-specific conditions with a set of three clinically measurable inputs: carotid–femoral wave speed, mean arterial pressure and area in the brachial artery. We quantified the uncertainty of the predicted pressure and flow waves in eight locations in the network and assessed the sensitivity of the model prediction with respect to the measurements of wave speed, pressure and cross-sectional area. Furthermore, we developed novel time-averaged sensitivity indices to assess the contribution of model parameters to the uncertainty of time-varying quantities (e.g., pressure and flow). The results from our patient-specific network model demonstrated that our novel indices allowed for a more accurate sensitivity analysis of time-varying quantities compared to conventional Sobol sensitivity indices.     

Pulse wave propagation in a model human arterial network: assessment of 1-D numerical simulations against in vitro measurements

A numerical model based on the nonlinear, one-dimensional (1-D) equations of pressure and flow wave propagation in conduit arteries is tested against a well-defined experimental 1:1 replica of the human arterial tree. The tree consists of 37 silicone branches representing the largest central systemic arteries in the human, including the aorta, carotid arteries and arteries that perfuse the upper and lower limbs and the main abdominal organs. The set-up is mounted horizontally and connected to a pulsatile pump delivering a periodic output similar to the aortic flow. Terminal branches end in simple resistance models, consisting of stiff capillary tubes leading to an overflow reservoir that reflects a constant venous pressure. The parameters required by the numerical algorithm are directly measured in the in vitro set-up and no data fitting is involved. Comparison of experimental and numerical pressure and flow waveforms shows the ability of the 1-D time-domain formulation to capture the main features of pulse wave propagation measured throughout the system test. As a consequence of the simple resistive boundary conditions used to reduce the uncertainty of the parameters involved in the simulation, the experimental set-up generates waveforms at terminal branches with additional non-physiological oscillations. The frequencies of these oscillations are well captured by the 1-D model, even though amplitudes are overestimated. Adding energy losses in bifurcations and including fluid inertia and compliance to the purely resistive terminal models does not reduce the underdamped effect, suggesting that wall visco-elasticity might play an important role in the experimental results. Nevertheless, average relative root-mean-square errors between simulations and experimental waveforms are smaller than 4% for pressure and 19% for the flow at all 70 locations studied.     

Toward a noninvasive subject-specific estimation of abdominal aortic pressure

A method for estimation of central arterial pressure based on linear one-dimensional wave propagation theory is presented in this paper. The equations are applied to a distributed model of the arterial tree, truncated by three-element windkessels. To reflect individual differences in the properties of the arterial trees, we pose a minimization problem from which individual parameters are identified. The idea is to take a measured waveform in a peripheral artery and use it as input to the model. The model subsequently predicts the corresponding waveform in another peripheral artery in which a measurement has also been made, and the arterial tree model is then calibrated in such a way that the computed waveform matches its measured counterpart. For the purpose of validation, invasively recorded abdominal aortic, brachial, and femoral pressures in nine healthy subjects are used. The results show that the proposed method estimates the abdominal aortic pressure wave with good accuracy. The root mean square error (RMSE) of the estimated waveforms was 1.61 +/- 0.73 mmHg, whereas the errors in systolic and pulse pressure were 2.32 +/- 1.74 and 3.73 +/- 2.04 mmHg, respectively. These results are compared with another recently proposed method based on a signal processing technique, and it is shown that our method yields a significantly (P < 0.01) lower RMSE. With more extensive validation, the method may eventually be used in clinical practice to provide detailed, almost individual, specific information as a valuable basis for decision making.     

Pulse wave propagation in a model human arterial network: Assessment of 1-D visco-elastic simulations against in vitro measurements

The accuracy of the nonlinear one-dimensional (1-D) equations of pressure and flow wave propagation in Voigt-type visco-elastic arteries was tested against measurements in a well-defined experimental 1:1 replica of the 37 largest conduit arteries in the human systemic circulation. The parameters required by the numerical algorithm were directly measured in the in vitro setup and no data fitting was involved. The inclusion of wall visco-elasticity in the numerical model reduced the underdamped high-frequency oscillations obtained using a purely elastic tube law, especially in peripheral vessels, which was previously reported in this paper [Matthys et al., 2007. Pulse wave propagation in a model human arterial network: Assessment of 1-D numerical simulations against in vitro measurements. J. Biomech. 40, 3476-3486]. In comparison to the purely elastic model, visco-elasticity significantly reduced the average relative root-mean-square errors between numerical and experimental waveforms over the 70 locations measured in the in vitro model: from 3.0% to 2.5% (p<0.012) for pressure and from 15.7% to 10.8% (p<0.002) for the flow rate. In the frequency domain, average relative errors between numerical and experimental amplitudes from the 5th to the 20th harmonic decreased from 0.7% to 0.5% (p<0.107) for pressure and from 7.0% to 3.3% (p<10(-6)) for the flow rate. These results provide additional support for the use of 1-D reduced modelling to accurately simulate clinically relevant problems at a reasonable computational cost.     

Modelling pulse wave propagation in the rabbit systemic circulation to assess the effects of altered nitric oxide synthesis

Pulse wave propagation in the mature rabbit systemic circulation was simulated using the one-dimensional equations of blood flow in compliant vessels. A corrosion cast of the rabbit circulation was manufactured to obtain arterial lengths and diameters. Pulse wave speeds and inflow and outflow boundary conditions were derived from in vivo data. Numerical results captured the main features of in vivo pressure and velocity pulse waveforms in the aorta, brachiocephalic artery and central ear artery. This model was used to elucidate haemodynamic mechanisms underlying changes in peripheral pulse waveforms observed in vivo after administering drugs that alter nitric oxide synthesis in the endothelial cells lining blood vessels. According to our model, these changes can be explained by single or combined alterations of blood viscosity, peripheral resistance and compliance, and the elasticity of conduit arteries.     

One-dimensional model for propagation of a pressure wave in a model of the human arterial network: comparison of theoretical and experimental results

Pulse wave evaluation is an effective method for arteriosclerosis screening. In a previous study, we verified that pulse waveforms change markedly due to arterial stiffness. However, a pulse wave consists of two components, the incident wave and multireflected waves. Clarification of the complicated propagation of these waves is necessary to gain an understanding of the nature of pulse waves in vivo. In this study, we built a one-dimensional theoretical model of a pressure wave propagating in a flexible tube. To evaluate the applicability of the model, we compared theoretical estimations with measured data obtained from basic tube models and a simple arterial model. We constructed different viscoelastic tube set-ups: two straight tubes; one tube connected to two tubes of different elasticity; a single bifurcation tube; and a simple arterial network with four bifurcations. Soft polyurethane tubes were used and the configuration was based on a realistic human arterial network. The tensile modulus of the material was similar to the elasticity of arteries. A pulsatile flow with ejection time 0.3 s was applied using a controlled pump. Inner pressure waves and flow velocity were then measured using a pressure sensor and an ultrasonic diagnostic system. We formulated a 1D model derived from the Navier-Stokes equations and a continuity equation to characterize pressure propagation in flexible tubes. The theoretical model includes nonlinearity and attenuation terms due to the tube wall, and flow viscosity derived from a steady Hagen-Poiseuille profile. Under the same configuration as for experiments, the governing equations were computed using the MacCormack scheme. The theoretical pressure waves for each case showed a good fit to the experimental waves. The square sum of residuals (difference between theoretical and experimental wave-forms) for each case was <10.0%. A possible explanation for the increase in the square sum of residuals is the approximation error for flow viscosity. However, the comparatively small values prove the validity of the approach and indicate the usefulness of the model for understanding pressure propagation in the human arterial network.     

On the choice of outlet boundary conditions for patient-specific analysis of aortic flow using computational fluid dynamics

Boundary conditions (BCs) are an essential part in computational fluid dynamics (CFD) simulations of blood flow in large arteries. Although several studies have investigated the influence of BCs on predicted flow patterns and hemodynamic wall parameters in various arterial models, there is a lack of comprehensive assessment of outlet BCs for patient-specific analysis of aortic flow. In this study, five different sets of outlet BCs were tested and compared using a subject-specific model of a normal aorta. Phase-contrast magnetic resonance imaging (PC-MRI) was performed on the same subject and velocity profiles extracted from the in vivo measurements were used as the inlet boundary condition. Computational results obtained with different outlet BCs were assessed in terms of their agreement with the PC-MRI velocity data and key hemodynamic parameters, such as pressure and flow waveforms and wall shear stress related indices. Our results showed that the best overall performance was achieved by using a well-tuned three-element Windkessel model at all model outlets, which not only gave a good agreement with in vivo flow data, but also produced physiological pressure waveforms and values. On the other hand, opening outlet BCs with zero pressure at multiple outlets failed to reproduce any physiologically relevant flow and pressure features.     

Experimental validation of a time-domain-based wave propagation model of blood flow in viscoelastic vessels

Time-domain-based one-dimensional wave propagation models of the arterial system are preferable over one-dimensional wave propagation models in the frequency domain since the latter neglect the non-linear convection forces present in the physiological situation, especially when the vessel is tapered. Moreover, one-dimensional wave propagation models of the arterial system can be used to provide boundary conditions for fully three-dimensional fluid-structure interaction computations that are usually defined in the time domain. In this study, a time-domain-based one-dimensional wave propagation model in a cross-sectional area, flow and pressure (A,q,p)-formulation is developed. Using this formulation, a constitutive law that includes viscoelasticity based on the mechanical behaviour of a Kelvin body, is introduced. The resulting pressure and flow waves travelling through a straight and tapered vessel are compared to experimental data obtained from measurements in an in vitro setup. The model presented shows to be well suited to predict wave propagation through these straight and tapered vessels with viscoelastic wall properties and hereto can serve as a time-domain-based method to model wave propagation in the human arterial system.     

Wave propagation in a model of the arterial circulation

The propagation of the arterial pulse wave in the large systemic arteries has been calculated using a linearised method of characteristics analysis to follow the waves generated by the heart. The model includes anatomical and physiological data for the 55 largest arteries adjusted so that the bifurcating tree of arteries is well matched for forward travelling waves. The peripheral arteries in the model are terminated by resistance elements which are adjusted to produce a physiologically reasonable distribution of mean blood flow. In the model, the pressure and velocity wave generated by the contraction of the left ventricle propagates to the periphery where it is reflected. These reflected waves are re-reflected by each of the bifurcations that they encounter and a very complex pattern of waves is generated. The results of the calculations exhibit many of the features of the systemic arteries, including the increase of the pulse pressure with distance away from the heart as well as the initial decrease and then the large increase in the magnitude of back flow during late systole going from the ascending aorta to the abdominal aorta to the arteries of the leg. The model is then used to study the effects of the reflection or absorption of waves by the heart and the mechanisms leading to the incisura are investigated. Calculations are carried out with the total occlusion of different arterial segments in order to model experiments in which the effects of the occlusion of different arteries on pressure and flow in the ascending aorta were measured. Finally, the effects of changes in peripheral resistance on pressure and velocity waveforms are also studied. We conclude from these calculations that the complex pattern of wave propagation in the large arteries may be the most important determinant of arterial haemodynamics.     

Inlet boundary conditions for blood flow simulations in truncated arterial networks

In the context of patient-specific cardiovascular applications, hemodynamics models (going from 3D to 0D) are often limited to a part of the arterial tree. This restriction implies the set up of artificial interfaces with the remaining parts of the cardiovascular system. In particular, the inlet boundary condition is crucial: it supplies the impulsion to the system and receives the reflected backward waves created by the distal network. Some aspects of this boundary condition need to be properly defined such as the treatment of backward waves (reflected or absorbed) and the value of the imposed hemodynamic wave (total or forward component). Most authors prescribe as inlet boundary condition (BC) the total measured variable (pressure, velocity or flow rate) in a reflective way. We show that with this type of inlet boundary condition, the model does not produce physiological waveforms. We suggest instead to prescribe only the forward component of the prescribed variable in an absorbing way. In this way, the computed reflected waves superpose with the prescribed forward waves to produce the total wave at the inlet. In this work, different inlet boundary conditions are implemented and compared for a 1D blood flow model. We test our boundary conditions on a truncated arterial model presented in the literature as well as on a patient-specific lower-limb model of a femoral bypass. We show that with this new boundary condition, a much better fitting is observed on the shape and intensity of the simulated pressure and velocity waves.     

Pressure and flow in the systemic arterial system

The dynamic characteristics of the proximal arterial system are studied by solving the nonlinear momentum and mass conservation equations for pressure and flow. The equations are solved for a model systemic arterial system that includes the aorta, common iliacs, and the internal and external iliac arteries. The model includes geometric and elastic taper of the aorta, nonlinearly elastic arteries, side flows, and a complex distal impedance. The model pressure wave shape, inlet and outlet impedance, wave travel, and apparent wave velocity compare favorably with the values measured on humans. Calculations indicate that: (i) reflections are the major factor determining the shape and distal amplification of the pressure wave in the arterial tree; (ii) although important in attenuating the proximal transmission of reflecting waves, geometric taper is not the major cause of the distal pressure wave amplification; (iii) the dicrotic wave is a result of peripheral reflection and is not due to the sudden change in flow at the end of systole; (iv) the elastic taper and nonlinearity of the wall elasticity are of minor significance in determining the flow and pressure profiles; and (v) in spite of numerous nonlinearities, the system behaves in a somewhat linear fashion for the lower frequency components.     

Computer modeling of the human systemic arterial tree

A model of the human systemic arterial tree has been devised, based on a lumped-parameter-circuit approximate form. This model has been set up and studied on an analog computer. A feature of this simulation is the division of the arterial system into sections whose lengths are inversely proportional (approximately) to their cross-sectional area-or what is termed ‘equal-volume’ modeling.

Great care was exercised in the determination of the model parameters, using expressions for these parameters from a recent paper by Rideout and Dick on fluid flow in distensible tubes, with numerical values based on measurements reported in the medical literature.

The simulated pressure and flow waveforms obtained with the model compare favorably with data recorded from the normal adult human, and exhibit such well-known features as distal delay and peaking of pressure pulses. The aortic input impedance vs. frequency curve checks well against measurements on the human. The model also provides a simple means for determination of cardiac output, cardiac work and cardiac power under various assumed conditions such as variation of heart rate.     

Heterogeneous mechanics of the mouse pulmonary arterial network

Individualized modeling and simulation of blood flow mechanics find applications in both animal research and patient care. Individual animal or patient models for blood vessel mechanics are based on combining measured vascular geometry with a fluid structure model coupling formulations describing dynamics of the fluid and mechanics of the wall. For example, one-dimensional fluid flow modeling requires a constitutive law relating vessel cross-sectional deformation to pressure in the lumen. To investigate means of identifying appropriate constitutive relationships, an automated segmentation algorithm was applied to micro-computerized tomography images from a mouse lung obtained at four different static pressures to identify the static pressure–radius relationship for four generations of vessels in the pulmonary arterial network. A shape-fitting function was parameterized for each vessel in the network to characterize the nonlinear and heterogeneous nature of vessel distensibility in the pulmonary arteries. These data on morphometric and mechanical properties were used to simulate pressure and flow velocity propagation in the network using one-dimensional representations of fluid and vessel wall mechanics. Moreover, wave intensity analysis was used to study effects of wall mechanics on generation and propagation of pressure wave reflections. Simulations were conducted to investigate the role of linear versus nonlinear formulations of wall elasticity and homogeneous versus heterogeneous treatments of vessel wall properties. Accounting for heterogeneity, by parameterizing the pressure/distention equation of state individually for each vessel segment, was found to have little effect on the predicted pressure profiles and wave propagation compared to a homogeneous parameterization based on average behavior. However, substantially different results were obtained using a linear elastic thin-shell model than were obtained using a nonlinear model that has a more physiologically realistic pressure versus radius relationship.     

Descending aorta subject-specific one-dimensional model validated against in vivo data

The aorta plays a major role in the cardiovascular system and its function and structure are primarily affected by aging, eating habits, life style and other cardiovascular risk factors, inducing increased stiffness which is associated with cardiovascular and cerebral morbi-mortality. Our objective was to develop and validate a robust subject-specific one-dimensional wave propagation numerical model of the descending aorta. This model with a cross-sectional area, velocity and pressure formulation is built using geometric and hemodynamic data measured on a specific person and is validated against in vivo data acquired on the same subject at three distinct anatomical locations along the thoracic aorta. We studied seven healthy volunteers, who underwent carotid applanation tonometry and aortic cardiovascular magnetic resonance (CMR). Responses of our model in terms of changes in central pressure waveform with arterial alterations were consistent with previously described physiological knowledge. Quantitative validation averaged over the three descending aortic locations and the seven subjects provided low rms errors (given in percentage of the maximal clinical value) between simulated and CMR data, i.e. area: 10±6%, velocity: 11±3%, flow rate: 9±3%. Finally, we also found low rms (5±2%) when comparing simulated pressure in the proximal aortic location against tonometric carotid pressure curves. In conclusion, this simple model performs similar to more complex models of the entire systemic arterial tree at a fraction of the cost, and could be of major usefulness in the non-invasive and local estimation of proximal biomechanical and hemodynamic indices.     

Pulsatile blood flow in the entire coronary arterial tree: theory and experiment

The pulsatility of coronary circulation can be accurately simulated on the basis of the measured branching pattern, vascular geometry, and material properties of the coronary vasculature. A Womersley-type mathematical model is developed to analyze pulsatile blood flow in diastole in the absence of vessel tone in the entire coronary arterial tree on the basis of previously measured morphometric data. The model incorporates a constitutive equation of pressure and cross-section area relation based on our previous experimental data. The formulation enables the prediction of the impedance, the pressure distribution, and the pulsatile flow distribution throughout the entire coronary arterial tree. The model is validated by experimental measurements in six diastolic arrested, vasodilated porcine hearts. The agreement between theory and experiment is excellent. Furthermore, the present pulse wave results at low frequency agree very well with previously published steady-state model. Finally, the phase angle of flow is seen to decrease along the trunk of the major coronary artery and primary branches toward the capillary vessels. This study represents the first, most extensive validated analysis of Womersley-type pulse wave transmission in the entire coronary arterial tree down to the first segment of capillaries. The present model will serve to quantitatively test various hypotheses in the coronary circulation under pulsatile flow conditions.     

A novel porous media-based approach to outflow boundary resistances of 1D arterial blood flow models

In this paper we introduce a novel method for prescribing terminal boundary conditions in one-dimensional arterial flow networks. This is carried out by coupling the terminal arterial vessel with a poro-elastic tube, representing the flow resistance offered by microcirculation. The performance of the proposed porous media-based model has been investigated through several different numerical examples. First, we investigate model parameters that have a profound influence on the flow and pressure distributions of the system. The simulation results have been compared against the waveforms generated by three elements (RCR) Windkessel model. The proposed model is also integrated into a realistic arterial tree, and the results obtained have been compared against experimental data at different locations of the network. The accuracy and simplicity of the proposed model demonstrates that it can be an excellent alternative for the existing models.

     

The effect of inlet waveforms on computational hemodynamics of patient-specific intracranial aneurysms

Due to the lack of patient-specific inlet flow waveform measurements, most computational fluid dynamics (CFD) simulations of intracranial aneurysms usually employ waveforms that are not patient-specific as inlet boundary conditions for the computational model. The current study examined how this assumption affects the predicted hemodynamics in patient-specific aneurysm geometries. We examined wall shear stress (WSS) and oscillatory shear index (OSI), the two most widely studied hemodynamic quantities that have been shown to predict aneurysm rupture, as well as maximal WSS (MWSS), energy loss (EL) and pressure loss coefficient (PLc). Sixteen pulsatile CFD simulations were carried out on four typical saccular aneurysms using 4 different waveforms and an identical inflow rate as inlet boundary conditions. Our results demonstrated that under the same mean inflow rate, different waveforms produced almost identical WSS distributions and WSS magnitudes, similar OSI distributions but drastically different OSI magnitudes. The OSI magnitude is correlated with the pulsatility index of the waveform. Furthermore, there is a linear relationship between aneurysm-averaged OSI values calculated from one waveform and those calculated from another waveform. In addition, different waveforms produced similar MWSS, EL and PLc in each aneurysm. In conclusion, inlet waveform has minimal effects on WSS, OSI distribution, MWSS, EL and PLc and a strong effect on OSI magnitude, but aneurysm-averaged OSI from different waveforms has a strong linear correlation with each other across different aneurysms, indicating that for the same aneurysm cohort, different waveforms can consistently stratify (rank) OSI of aneurysms.     

Abnormal arterial flows by a distributed model of the fetal circulation

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

Beyond the aorta: partial transmission of reflected waves from aortic coarctation into supra-aortic branches modulates cerebral hemodynamics and left ventricular load

Wave reflection from the site of aortic coarctation produces a reflected backward compression wave (BCW) that raises left ventricular (LV) afterload. However, not all reflected wave power will propagate back to the LV. This study investigated the hypothesis that the BCW is partially transmitted into supra-aortic vessels as a forward wave and explored the consequences of this phenomenon for cerebral and LV haemodynamic load. In eight sheep, high fidelity pressure and flow were measured in the aortic trunk (AoT) and brachiocephalic trunk (BCT, the single supra-aortic vessel present in sheep) at baseline and during two levels of proximal descending aortic constriction. Wave power analysis showed that aortic constriction produced not only a BCW in the AoT, but also a second forward compression wave (\(\mathrm{FCW}_{2})\) in the BCT that augmented pressure and flow after the initial forward compression wave (\(\mathrm{FCW}_{1})\). Mathematical analysis and a one-dimensional model of the human systemic arteries and aortic coarctation suggested that the relative transmission of waves into supra-aortic vessels versus the aorta was determined by the relative admittances of these vessels. Reducing supra-aortic admittance (1) increased pressure and flow pulsatility in cerebral arteries, (2) produced carotid and middle cerebral arterial flow waveforms with an older adult phenotype, (3) promoted transmission of reflected wave power towards the LV and (4) substantially increased mid- to late-systolic myocardial stress, which may promote LV hypertrophy. These findings suggest that wave transmission into supra-aortic branches has an important impact on both cerebral hemodynamics and LV load in aortic coarctation.     

Exercise reduces arterial pressure augmentation through vasodilation of muscular arteries in humans

Exercise markedly influences pulse wave morphology, but the mechanism is unknown. We investigated whether effects of exercise on the arterial pulse result from alterations in stroke volume or pulse wave velocity (PWV)/large artery stiffness or reduction of pressure wave reflection. Healthy subjects (n = 25) performed bicycle ergometry. with workload increasing from 25 to 150 W for 12 min. Digital arterial pressure waveforms were recorded using a servo-controlled finger cuff. Radial arterial pressure waveforms and carotid-femoral PWV were determined by applanation tonometry. Stroke volume was measured by echocardiography, and brachial and femoral artery blood flows and diameters were measured by ultrasound. Digital waveforms were recorded continuously. Other measurements were made before and after exercise. Exercise markedly reduced late systolic and diastolic augmentation of the peripheral pressure pulse. At 15 min into recovery, stroke volume and PWV were similar to baseline values, but changes in pulse wave morphology persisted. Late systolic augmentation index (radial pulse) was reduced from 54 +/- 3.9% at baseline to 42 +/- 3.7% (P < 0.01), and diastolic augmentation index (radial pulse) was reduced from 37 +/- 1.8% to 25 +/- 2.9% (P < 0.001). These changes were accompanied by an increase in femoral blood flow (from 409 +/- 44 to 773 +/- 48 ml/min, P < 0.05) and an increase in femoral artery diameter (from 8.2 +/- 0.4 to 8.6 +/- 0.4 mm, P < 0.05). In conclusion, exercise dilates muscular arteries and reduces arterial pressure augmentation, an effect that will enhance ventricular-vascular coupling and reduce load on the left ventricle.