Non-invasive estimation of the parameters of a three-element windkessel model of aortic arch arteries in patients undergoing thoracic endovascular aortic repair

Background: Image-based computational hemodynamic modeling and simulations are important for personalized diagnosis and treatment of cardiovascular diseases. However, the required patient-specific boundary conditions are often not available and need to be estimated. Methods: We propose a pipeline for estimating the parameters of the popular three-element Windkessel (WK3) models (a proximal resistor in series with a parallel combination of a distal resistor and a capacitor) of the aortic arch arteries in patients receiving thoracic endovascular aortic repair of aneurysms. Pre-operative and post-operative 1-week duplex ultrasound scans were performed to obtain blood flow rates, and intra-operative pressure measurements were also performed invasively using a pressure transducer pre- and post-stent graft deployment in arch arteries. The patient-specific WK3 model parameters were derived from the flow rate and pressure waveforms using an optimization algorithm reducing the error between simulated and measured pressure data. The resistors were normalized by total resistance, and the capacitor was normalized by total resistance and heart rate. The normalized WK3 parameters can be combined with readily available vessel diameter, brachial blood pressure, and heart rate data to estimate WK3 parameters of other patients non-invasively. Results: Ten patients were studied. The medians (interquartile range) of the normalized proximal resistor, distal resistor, and capacitor parameters are 0.10 (0.07-0.15), 0.90 (0.84-0.93), and 0.46 (0.33-0.58), respectively, for common carotid artery; 0.03 (0.02-0.04), 0.97 (0.96-0.98), and 1.91 (1.63-2.26) for subclavian artery; 0.18 (0.08-0.41), 0.82 (0.59-0.92), and 0.47 (0.32-0.85) for vertebral artery. The estimated pressure showed fairly high tolerance to patient-specific inlet flow rate waveforms using the WK3 parameters estimated from the medians of the normalized parameters. Conclusion: When patient-specific outflow boundary conditions are not available, our proposed pipeline can be used to estimate the WK3 parameters of arch arteries.

Image-Based Computational Hemodynamics Analysis of Systolic Obstruction in Hypertrophic Cardiomyopathy

Hypertrophic Cardiomyopathy (HCM) is a pathological condition characterized by an abnormal thickening of the myocardium. When affecting the medio-basal portion of the septum, it is named Hypertrophic Obstructive Cardiomyopathy (HOCM) because it induces a flow obstruction in the left ventricular outflow tract. In any type of HCM, the myocardial function can become compromised, possibly resulting in cardiac death. In this study, we investigated with computational analysis the hemodynamics of patients with different types of HCM. The aim was quantifying the effects of this pathology on the intraventricular blood flow and pressure gradients, and providing information potentially useful to guide the indication and the modality of the surgical treatment (septal myectomy). We employed an image-based computational approach, integrating fluid dynamics simulations with geometric and functional data, reconstructed from standard cardiac cine-MRI acquisitions. We showed that with our approach we can better understand the patho-physiological behavior of intraventricular blood flow dynamics due to the abnormal morphological and functional aspect of the left ventricle. The main results of our investigation are: (a) a detailed patient-specific analysis of the blood velocity, pressure and stress distribution associated to HCM; (b) a computation-based classification of patients affected by HCM that can complement the current clinical guidelines for the diagnosis and treatment of HOCM.

Identification of vortex structures in a cohort of 204 intracranial aneurysms

An intracranial aneurysm (IA) is a cerebrovascular pathology that can lead to death or disability if ruptured. Abnormal wall shear stress (WSS) has been associated with IA growth and rupture, but little is known about the underlying flow physics related to rupture-prone IAs. Previous studies, based on analysis of a few aneurysms or partial views of three-dimensional vortex structures, suggest that rupture is associated with complex vortical flow inside IAs. To further elucidate the relevance of vortical flow in aneurysm pathophysiology, we studied 204 patient IAs (56 ruptured and 148 unruptured). Using objective quantities to identify three-dimensional vortex structures, we investigated the characteristics associated with aneurysm rupture and if these features correlate with previously proposed WSS and morphological characteristics indicative of IA rupture. Based on the Q-criterion definition of a vortex, we quantified the degree of the aneurysmal region occupied by vortex structures using the volume vortex fraction (vVF) and the surface vortex fraction (sVF). Computational fluid dynamics simulations showed that the sVF, but not the vVF, discriminated ruptured from unruptured aneurysms. Furthermore, we found that the near-wall vortex structures co-localized with regions of inflow jet breakdown, and significantly correlated to previously proposed haemodynamic and morphologic characteristics of ruptured IAs.