Pressure and Flow Wave Propagation in Patient-Specific Models of the Arterial Tree

Blood flow in the arterial circulation induces hemodynamic forces that play an important role in various forms of vascular diseases. Temporal variation of the wall shear stress seems to play a significant role in atherogenesis and plaque stability. Flow induced wall shear stress has been linked to growth and possibly rupture of the aneurysm wall. Hemodynamic forces are patient-specific and difficult to assess in the clinic. At present, there is no in vivo measurement technique that enables measurement of hemodynamic forces to the degree of precision needed. However, when imaging modalities used frequently in clinical routine re-create high-definition, patient geometric quantification of the blood vessel, they can be employed as a base for creating predictive hemodynamic models. Which in the case of understanding healthy vs. pathologic blood flow within the cerebral or systemic circulation, renders this an interesting approach. First, we developed a "generic 1D" distributed model of the human arterial tree including the primary systemic arteries and coupled this to a heart model. The fluid mechanics equations were solved numerically to obtain pressure and flow throughout the arterial tree. A nonlinear viscoelastic constitutive law for the arterial wall was considered while distal vessels were terminated with a three-element Windkessel model. The coronary arteries were modeled assuming a systolic flow impediment proportional to ventricular varying elastance. The model predictions were validated with noninvasive measurements of pressure and flow performed in young volunteers. Flow in the large arteries was visualized with magnetic resonance imaging, cerebral flow detected with ultrasound Doppler and blood pressure measured with applanation tonometry. Model predictions at different arterial locations compared well to measured flow and pressure waves at the same anatomical points. Thus, the generic 1D model reflected the flow and pressure measurements of the "average subject" of our volunteer population. Following the same approach as the generic 1D model, we built and validated a patient-specific model. In this case, geometric data, flow and pressure measurements were obtained for one person. The model predicted pressure and flow waveforms in good agreement with the in-vivo measurements with regards to wave shape and features. Comparison with a generic 1-D model has shown that the patient-specific model better predicted pressure and flow at specific arterial sites. Overall, the patient-specific 1-D model was able to predict pressure and flow waveforms in the main systemic circulation, whereas this was not always the case for a generic 1-D model. The inherent underestimation of energy losses of the 1-D wave propagation model, due to bifurcations, non-planarity and complex geometry, were examined. The 1-D model was compared to a rigid wall 3-D computational fluid dynamic model. Newtonian and non-Newtonian blood properties were studied and the longitudinal pressure distribution along the arteries was compared with the 1-D patient-specific model mean pressure prediction. The results indicated that pressure drop is significant only in small diameter vessels such as the precerebral and cerebral arteries. In these vessels the 1-D model in comparison to 3-D models consistently underestimated pressure drop. The complex flow patterns resulting from asymmetry and bifurcation yield shear stresses in the 3-D model that were greater than the 1-D model. A 3-D unsteady fluid structure interaction simulation in a patient-specific model was performed to simultaneously capture the flow details, given by the 3-D model, and wave propagation phenomena, provided by the 1-D model. The 3-D unsteady fluid structure interaction approach is the most computationally intense and cumbersome, but it allows physiological simulations with a high level of detail and accuracy. For instance, this approach could be relevant to obtain blood flow details in regions that are prone to atherosclerotic plaques or development of aneurysms. The 3-D fluid structure interaction simulation was performed for a patient-specific aorta. Important clinical parameters such as wall shear stress were quantified and significant differences were found in comparison to the rigid wall 3-D simulation indicating the relevance of a fluid structure interaction approach. A comparison of the fluid structure interaction to an equivalent 1-D model resulted in good reproduction of the pressure and flow waveforms. The effect of a decreased compliance of the arterial tree on hemodynamical parameters has been assessed with the use of a 1-D model. Local, proximal aorta and global stiffening of the arterial tree were modeled and led to two different mechanisms that contribute to the increase in central pulse pressure. They probably both contribute to systolic hypertension and their relative contribution depends on the topology of arterial stiffening and geometrical alterations taking place in aging or in disease. All these patient-specific models are about to being in use in a clinical environment and will be useful for providing better diagnostics and treatment planning in a near future.