Purpose: Pediatric and adult patients with sickle cell anemia (SCA) are at increased risk of stroke and cerebrovascular accident. In the general adult population, there is a relationship between arterial hemodynamics and pathology; however, this relationship in SCA patients remains to be elucidated. The aim of this work was to characterize circle of Willis hemodynamics in patients with SCA and quantify the impact of viscosity choice on pathophysiologically-relevant hemodynamics measures. Methods: Based on measured vascular geometries, time-varying flow rates, and blood parameters, detailed patient-specific simulations of the circle of Willis were conducted for SCA patients (n = 6). Simulations quantified the impact of patient-specific and standard blood viscosities on wall shear stress (WSS). Results: These results demonstrated that use of a standard blood viscosity introduces large errors into the estimation of pathophysiologically-relevant hemodynamic parameters. Standard viscosity models overpredicted peak WSS by 55% and 49% for steady and pulsatile flow, respectively. Moreover, these results demonstrated non-uniform, spatial patterns of positive and negative WSS errors related to viscosity, and standard viscosity simulations overpredicted the time-averaged WSS by 32% (standard deviation = 7.1%). Finally, differences in shear rate demonstrated that the viscosity choice alters the simulated near-wall flow field, impacting hemodynamics measures. Conclusions: This work presents simulations of circle of Willis arterial flow in SCA patients and demonstrates the importance and feasibility of using a patient-specific viscosity in these simulations. Accurately characterizing cerebrovascular hemodynamics in SCA populations has potential for elucidating the pathophysiology of large-vessel occlusion, aneurysms, and tissue damage in these patients.
This study aims to investigate the capability of smoothed particle hydrodynamics (SPH), a fully Lagrangian mesh-free method, to simulate the bulk blood flow dynamics in two realistic left ventricular (LV) models. Three dimensional geometries and motion of the LV, proximal left atrium and aortic root are extracted from cardiac magnetic resonance imaging and multi-slice computed tomography imaging data. SPH simulation results are analyzed and compared with those obtained using a traditional finite volume-based numerical method, and to in vivo phase contrast magnetic resonance imaging and echocardiography data, in terms of the large-scale blood flow phenomena usually clinically measured. A quantitative comparison of the velocity fields and global flow parameters between the in silico models and the in vivo data shows a reasonable agreement, given the inherent uncertainties and limitations in the modeling and imaging techniques. The results indicate the capability of SPH as a promising tool for predicting clinically relevant large-scale LV flow information.
Total cavopulmonary connection is the result of a series of palliative surgical repairs performed on patients with single ventricle heart defects. The resulting anatomy has complex and unsteady hemodynamics characterized by flow mixing and flow separation. Although varying degrees of flow pulsatility have been observed in vivo, non-pulsatile (time-averaged) boundary conditions have traditionally been assumed in hemodynamic modeling, and only recently have pulsatile conditions been incorporated without completely characterizing their effect or importance. In this study, 3D numerical simulations with both pulsatile and non-pulsatile boundary conditions were performed for 24 patients with different anatomies and flow boundary conditions from Georgia Tech database. Flow structures, energy dissipation rates and pressure drops were compared under rest and simulated exercise conditions. It was found that flow pulsatility is the primary factor in determining the appropriate choice of boundary conditions, whereas the anatomic configuration and cardiac output had secondary effects. Results show that the hemodynamics can be strongly influenced by the presence of pulsatile flow. However, there was a minimum pulsatility threshold, identified by defining a weighted pulsatility index (wPI), above which the influence was significant. It was shown that when wPI < 30%, the relative error in hemodynamic predictions using time-averaged boundary conditions was less than 10% compared to pulsatile simulations. In addition, when wPI < 50, the relative error was less than 20%. A correlation was introduced to relate wPI to the relative error in predicting the flow metrics with non-pulsatile flow conditions.
Cardiovascular simulations have great potential as a clinical tool for planning and evaluating patient-specific treatment strategies for those suffering from congenital heart diseases, specifically Fontan patients. However, several bottlenecks have delayed wider deployment of the simulations for clinical use; the main obstacle is simulation cost. Currently, time-averaged clinical flow measurements are utilized as numerical boundary conditions (BCs) in order to reduce the computational power and time needed to offer surgical planning within a clinical time frame. Nevertheless, pulsatile blood flow is observed in vivo, and its significant impact on numerical simulations has been demonstrated. Therefore, it is imperative to carry out a comprehensive study analyzing the sensitivity of using time-averaged BCs. In this study, sensitivity is evaluated based on the discrepancies between hemodynamic metrics calculated using time-averaged and pulsatile BCs; smaller discrepancies indicate less sensitivity. The current study incorporates a comparison between 3D patient-specific CFD simulations using both the time-averaged and pulsatile BCs for 101 Fontan patients. The sensitivity analysis involves two clinically important hemodynamic metrics: hepatic flow distribution (HFD) and indexed power loss (iPL). Paired demographic group comparisons revealed that HFD sensitivity is significantly different between single and bilateral superior vena cava cohorts but no other demographic discrepancies were observed for HFD or iPL. Multivariate regression analyses show that the best predictors for sensitivity involve flow pulsatilities, time-averaged flow rates, and geometric characteristics of the Fontan connection. These predictors provide patient-specific guidelines to determine the effectiveness of analyzing patient-specific surgical options with time-averaged BCs within a clinical time frame.
The ultimate goal of Fontan surgical planning is to provide additional insights into the clinical decision-making process. In its current state, surgical planning offers an accurate hemodynamic assessment of the pre-operative condition, provides anatomical constraints for potential surgical options, and produces decent post-operative predictions if boundary conditions are similar enough between the pre-operative and post-operative states. Moving forward, validation with post-operative data is a necessary step in order to assess the accuracy of surgical planning and determine which methodological improvements are needed. Future efforts to automate the surgical planning process will reduce the individual expertise needed and encourage use in the clinic by clinicians. As post-operative physiologic predictions improve, Fontan surgical planning will become an more effective tool to accurately model patient-specific hemodynamics.
Edge-to-edge repair for mitral valve regurgitation is being increasingly performed in high-surgical risk patients using minimally invasive mitral clipping devices. Known procedural complications include chordal rupture and mitral leaflet perforation. Hence, it is important to quantitatively evaluate the effect of edge-to-edge repair on chordal integrity. in this study, we employ a computational mitral valve model to simulate functional mitral regurgitation (FMR) by creating papillary muscle displacement. Edge-to-edge repair is then modeled by simulated coaptation of the mid portion of the mitral leaflets. in the setting of simulated FMR, edge-to-edge repair was shown to sustain low regurgitant orifice area, until a two fold increase in the inter-papillary muscle distance as compared to the normal mitral valve. Strain in the chordae was evaluated near the papillary muscles and the leaflets. Following edge-to-edge repair, strain near the papillary muscles did not significantly change relative to the unrepaired valve, while strain near the leaflets increased significantly relative to the unrepaired valve. These data demonstrate the potential for computational simulations to aid in the pre-procedural evaluation of possible complications such as chordal rupture and leaflet perforation following percutaneous edge-to-edge repair.
Background
Astronauts undergoing long-duration spaceflight are exposed to numerous health risks, including Spaceflight-Associated Neuro-Ocular Syndrome (SANS), a spectrum of ophthalmic changes that can result in permanent loss of visual acuity. The etiology of SANS is not well understood but is thought to involve changes in cerebrovascular flow dynamics in response to microgravity. There is a paucity of knowledge in this area; in particular, cerebrospinal fluid (CSF) flow dynamics have not been well characterized under microgravity conditions. Our study was designed to determine the effect of simulated microgravity (head-down tilt [HDT]) on cerebrovascular flow dynamics. We hypothesized that microgravity conditions simulated by acute HDT would result in increases in CSF pulsatile flow.
Methods
In a prospective cohort study, we measured flow in major cerebral arteries, veins, and CSF spaces in fifteen healthy volunteers using phase contrast magnetic resonance (PCMR) before and during 15° HDT.
Results
We found a decrease in all CSF flow variables [systolic peak flow (p = 0.009), and peak-to-peak pulse amplitude (p = 0.001)]. Cerebral arterial average flow (p = 0.04), systolic peak flow (p = 0.04), and peak-to-peak pulse amplitude (p = 0.02) all also significantly decreased. We additionally found a decrease in average cerebral arterial flow (p = 0.040). Finally, a significant increase in cerebral venous cross-sectional area under HDT (p = 0.005) was also observed.
Conclusions
These results collectively demonstrate that acute application of −15° HDT caused a reduction in CSF flow variables (systolic peak flow and peak-to-peak pulse amplitude) which, when coupled with a decrease in average cerebral arterial flow, systolic peak flow, and peak-to-peak pulse amplitude, is consistent with a decrease in cardiac-related pulsatile CSF flow. These results suggest that decreases in cerebral arterial inflow were the principal drivers of decreases in CSF pulsatile flow.
Computational fluid dynamics (CFD) tools have been extensively applied to study the hemodynamics in the total cavopulmonary connection (TCPC) in patients with only a single functioning ventricle. Without the contraction of a sub-pulmonary ventricle, pulsatility of flow through this connection is low and variable across patients, which is usually neglected in most numerical modeling studies. Recent studies suggest that such pulsatility can be non-negligible and can be important in hemodynamic predictions. The goal of this work is to compare the results of an in-house numerical methodology for simulating pulsatile TCPC flow with experimental results. Digital particle image velocimetry (DPIV) was acquired on TCPC in vitro models to evaluate the capability of the CFD tool in predicting pulsatile TCPC flow fields. In vitro hemodynamic measurements were used to compare the numerical prediction of power loss across the connection. The results demonstrated the complexity of the pulsatile TCPC flow fields and the validity of the numerical approach in simulating pulsatile TCPC flow dynamics in both idealized and complex patient specific models.