A 4D flow MRI evaluation of the impact of shear-dependent fluid viscosity on in vitro Fontan circulation flow

In vitro hemodynamic performance of a blood pump for self-powered venous assist in univentricular hearts

Univentricular heart anomalies represent a group of severe congenital heart defects necessitating early surgical intervention in infancy. The Fontan procedure, the final stage of single-ventricle palliation, establishes a serial connection between systemic and pulmonary circulation by channeling venous return to the lungs. The absence of the subpulmonary ventricle in this peculiar circulation progressively eventuates in failure, primarily due to chronic elevation in inferior vena cava (IVC) pressure. This study experimentally validates the effectiveness of an intracorporeally-powered venous ejector pump (VEP) in reducing IVC pressure in Fontan patients. The VEP exploits a fraction of aortic flow to create a jet-venturi effect for the IVC, negating the external power requirement and driveline infections. An invitro Fontan mock-up circulation loop is developed and the impact of VEP design parameters and physiological conditions is assessed using both idealized and patient-specific total cavopulmonary connection (TCPC) phantoms. The VEP performance in reducing IVC pressure exhibited an inverse relationship with the cardiac output and extra-cardiac conduit (ECC) size and a proportional relationship with the transpulmonary pressure gradient (TPG) and mean arterial pressure (MAP). The ideal VEP with fail-safe features provided an IVC pressure drop of 1.82 ± 0.49, 2.45 ± 0.54, and 3.12 ± 0.43 mm Hg for TPG values of 6, 8, and 10 mm Hg, respectively, averaged over all ECC sizes and cardiac outputs. Furthermore, the arterial oxygen saturation was consistently maintained above 85% during full-assist mode. These results emphasize the potential utility of the VEP to mitigate elevated venous pressure in Fontan patients.

Hemi-Fontan and bidirectional Glenn operations result in flow-mediated viscous energy loss at the time of stage II palliation

Background

Superior cavopulmonary connection (SCPC) for stage II palliation of hypoplastic left heart syndrome (HLHS) is achieved most frequently by either a bidirectional Glenn (BDG) or hemi-Fontan (HF) operation. The comparison of flow hemodynamic efficiency at the region of surgical reconstruction and in proximal pulmonary arteries has been evaluated primarily using computational modeling techniques with conflicting reports. The purpose of this descriptive study was to compare flow hemodynamics following stage II (BDG vs HF) using 4-dimensional flow magnetic resonance imaging (4D-Flow MRI) with particular focus on flow-mediated viscous energy loss (EL') under matched hemodynamic conditions.

Methods

Patients with hypoplastic left heart syndrome (HLHS) who underwent either HF or BDG as part of stage II palliation underwent pre-Fontan 4D-Flow MRI. Patients were matched by the pulmonary vascular resistance index, net superior vena cava (SVC) flow, right pulmonary artery (RPA) and left pulmonary artery (LPA) size, and age. Maximum EL' throughout the cardiac cycle was calculated along the SVC-RPA and SVC-LPA tracts.

Results

Eight patients who underwent HF as part of their stage II single ventricle palliation were matched with 8 patients who underwent BDG. There were no differences between the 2 groups in median volumetric indices, including end-diastolic volume (P = .278) and end-systolic volume (P = .213). Moreover, no differences were observed in ejection fraction (P = .091) and cardiac index (P = .324). There also were no differences in peak EL' measured along the SVC-RPA tract (median, 0.05 mW for HF vs 0.04 mW for BDG; P = .365) or along the SVC-LPA tract (median, 0.05 mW vs 0.04 mW; P = .741).

Conclusions

The second stage of surgical palliation of HLHS using either HF or BDG results in similar flow-mediated viscous energy loss throughout the SCPC junction. 4D-Flow MRI and computational methods should be applied together to investigate flow hemodynamic patterns throughout the Fontan palliation and overall efficiency of the Fontan circuit.

Central Venous Waveform Patterns in the Fontan Circulation Independently Contribute to the Prediction of Composite Survival

It is well appreciated that the Fontan circulation perturbs central venous hemodynamics, with elevated pressure being the clearest change associated with Fontan comorbidities, such as Fontan-associated liver disease (FALD) and protein-losing enteropathy (PLE). Our group has better quantity of these venous perturbations through single- and multi-location analyses of flow waveforms obtained from magnetic resonance imaging of Fontan patients. Here, we determine if such analyses, which yield principal components (PC) that describe flow features, are associated with Fontan survival. Patients with a Fontan circulation (N = 140) that underwent free-breathing and mechanically ventilated cardiac MRI were included in this study. Standard volumetric and functional hemodynamics, as well as flow analysis principal components, were subjected to univariate and bivariate Cox regression analyses to determine composite clinical outcome, including plastic bronchitis, PLE, and referral and receipt of transplant. Unsurprisingly, ventricular function measures of ejection fraction (EF; HR = 0.88, p < 0.0001), indexed end-systolic volume (ESVi; HR 1.02, p < 0.0001), and indexed end-diastolic volume (EDVi; HR = 1.02, p = 0.0007) were found as specific predictors of clinical events, with specificities uniformly > 0.75. Additionally a feature of IVC flow (PC2) indicating increased flow in systole was found as a highly sensitive predictor (HR = 0.851, p = 0.027, sensitivity 0.93). In bivariate prediction, combinations of ventricular function (EF, ESVi, EDVi) with this IVC flow feature yielded best overall prediction of composite outcome. This suggests that central venous waveform analysis relays additional information about Fontan patient survival and that coupling sensitive and specific measures in bivariate analysis is a useful approach for obtaining superior prediction of survival.

Numerical simulation of the cavopulmonary connection flow with conduit stenoses of varying configurations

The circulation in the total cavopulmonary connection (TCPC) is a low-energy system which operation and efficiency are subjected to multiple factors. Some retrospective studies report that the abnormal narrowing of vessels in the system, i.e. stenosis, is one of the most dangerous geometric factors which can result in heart failure. In the present study, the effect of varying extracardiac conduit (ECC) stenosis on the hemodynamics in a surrogate TCPC model is investigated using high-fidelity numerical simulations. The efficiency of the surrogate TCPC model was quantified according to the power loss, relative perfusion in lungs and the percentage of conduit surface area with abnormally low and high wall shear stress for venous flow. Additionally, the impact of respiration and asymmetry in the stenosis geometry to the system was examined. The results show that the flow in the TCPC model exhibits pronounced unsteadiness even under the steady initial boundary conditions, while the uneven pulmonary flow distribution and the presence of the ECC stenosis amplify the chaotic nature of the flow. Energy efficiency of the system is shown to strongly correlate with amount of vortical structures in the model and their range of scales. Finally, the study demonstrates that the presence of respiration in the model adds to perturbations in the flow which causes increase in the power loss. Results obtained in the study provide valuable insights on how the ECC stenosis effect the flow in the surrogate TCPC model under different flow conditions.

In Silico Evaluation of a Self-powered Venous Ejector Pump for Fontan Patients

PURPOSE

The Fontan circulation carries a dismal prognosis in the long term due to its peculiar physiology and lack of a subpulmonic ventricle. Although it is multifactorial, elevated IVC pressure is accepted to be the primary cause of Fontan's high mortality and morbidity. This study presents a self-powered venous ejector pump (VEP) that can be used to lower the high IVC venous pressure in single-ventricle patients.

METHODS

A self-powered venous assist device that exploits the high-energy aortic flow to lower IVC pressure is designed. The proposed design is clinically feasible, simple in structure, and is powered intracorporeally. The device's performance in reducing IVC pressure is assessed by conducting comprehensive computational fluid dynamics simulations in idealized total cavopulmonary connections with different offsets. The device was finally applied to complex 3D reconstructed patient-specific TCPC models to validate its performance.

RESULTS

The assist device provided a significant IVC pressure drop of more than 3.2 mm Hg in both idealized and patient-specific geometries, while maintaining a high systemic oxygen saturation of more than 90%. The simulations revealed no significant caval pressure rise (< 0.1 mm Hg) and sufficient systemic oxygen saturation (> 84%) in the event of device failure, demonstrating its fail-safe feature.

CONCLUSIONS

A self-powered venous assist with promising in silico performance in improving Fontan hemodynamics is proposed. Due to its passive nature, the device has the potential to provide palliation for the growing population of patients with failing Fontan.

Hemodynamic characteristics in a cerebral aneurysm model using non-Newtonian blood analogues

This study aims to develop an experimentally validated computational fluid dynamics (CFD) model to estimate hemodynamic characteristics in cerebral aneurysms (CAs) using non-Newtonian blood analogues. Blood viscosities varying with shear rates were measured under four temperatures first, which serves as the reference for the generation of blood analogues. Using the blood analogue, particle image velocimetry (PIV) measurements were conducted to quantify flow characteristics in a CA model. Then, using the identical blood properties in the experiment, CFD simulations were executed to quantify the flow patterns, which were used to compare with the PIV counterpart. Additionally, hemodynamic characteristics in the simplified Newtonian and non-Newtonian models were quantified and compared using the experimentally validated CFD model. Results showed the proposed non-Newtonian viscosity model can predict blood shear-thinning properties accurately under varying temperatures and shear rates. Another developed viscosity model based on the blood analogue can well represent blood rheological properties. The comparisons in flow characteristics show good agreements between PIV and CFD, demonstrating the developed CFD model is qualified to investigate hemodynamic factors within CAs. Furthermore, results show the differences of absolute values were insignificant between Newtonian and non-Newtonian fluids in the distributions of wall shear stress (WSS) and oscillatory shear index (OSI) on arterial walls. However, not only does the simplified Newtonian model underestimate WSS and OSI in most regions of the aneurysmal sac, but it also makes mistakes in identifying the high OSI regions on the sac surface, which may mislead the hemodynamic assessment on the pathophysiology of CAs.

Patient-Specific Inverse Modeling of In Vivo Cardiovascular Mechanics with Medical Image-Derived Kinematics as Input Data: Concepts, Methods, and Applications

Inverse modeling approaches in cardiovascular medicine are a collection of methodologies that can provide non-invasive patient-specific estimations of tissue properties, mechanical loads, and other mechanics-based risk factors using medical imaging as inputs. Its incorporation into clinical practice has the potential to improve diagnosis and treatment planning with low associated risks and costs. These methods have become available for medical applications mainly due to the continuing development of image-based kinematic techniques, the maturity of the associated theories describing cardiovascular function, and recent progress in computer science, modeling, and simulation engineering. Inverse method applications are multidisciplinary, requiring tailored solutions to the available clinical data, pathology of interest, and available computational resources. Herein, we review biomechanical modeling and simulation principles, methods of solving inverse problems, and techniques for image-based kinematic analysis. In the final section, the major advances in inverse modeling of human cardiovascular mechanics since its early development in the early 2000s are reviewed with emphasis on method-specific descriptions, results, and conclusions. We draw selected studies on healthy and diseased hearts, aortas, and pulmonary arteries achieved through the incorporation of tissue mechanics, hemodynamics, and fluid-structure interaction methods paired with patient-specific data acquired with medical imaging in inverse modeling approaches.

Assessment of turbulent blood flow and wall shear stress in aortic coarctation using image-based simulations

In this study, we analyzed turbulent flows through a phantom (a 180[Formula: see text] bend with narrowing) at peak systole and a patient-specific coarctation of the aorta (CoA), with a pulsating flow, using magnetic resonance imaging (MRI) and computational fluid dynamics (CFD). For MRI, a 4D-flow MRI is performed using a 3T scanner. For CFD, the standard [Formula: see text], shear stress transport [Formula: see text], and Reynolds stress (RSM) models are applied. A good agreement between measured and simulated velocity is obtained for the phantom, especially for CFD with RSM. The wall shear stress (WSS) shows significant differences between CFD and MRI in absolute values, due to the limited near-wall resolution of MRI. However, normalized WSS shows qualitatively very similar distributions of the local values between MRI and CFD. Finally, a direct comparison between in vivo 4D-flow MRI and CFD with the RSM turbulence model is performed in the CoA. MRI can properly identify regions with locally elevated or suppressed WSS. If the exact values of the WSS are necessary, CFD is the preferred method. For future applications, we recommend the use of the combined MRI/CFD method for analysis and evaluation of the local flow patterns and WSS in the aorta.

The Rheology of the Carotid Sinus: A Path Toward Bioinspired Intervention

The association between blood viscosity and pathological conditions involving a number of organ systems is well known. However, how the body measures and maintains appropriate blood viscosity is not well-described. The literature endorsing the function of the carotid sinus as a site of baroreception can be traced back to some of the earliest descriptions of digital pressure on the neck producing a drop in blood delivery to the brain. For the last 30 years, improved computational fluid dynamic (CFD) simulations of blood flow within the carotid sinus have demonstrated a more nuanced understanding of the changes in the region as it relates to changes in conventional metrics of cardiovascular function, including blood pressure. We suggest that the unique flow patterns within the carotid sinus may make it an ideal site to transduce flow data that can, in turn, enable real-time measurement of blood viscosity. The recent characterization of the PIEZO receptor family in the sinus vessel wall may provide a biological basis for this characterization. When coupled with other biomarkers of cardiovascular performance and descriptions of the blood rheology unique to the sinus region, this represents a novel venue for bioinspired design that may enable end-users to manipulate and optimize blood flow.

Diagnosis and Treatment of Right Heart Failure in Pulmonary Vascular Diseases: A National Heart, Lung, and Blood Institute Workshop

Right ventricular dysfunction is a hallmark of advanced pulmonary vascular, lung parenchymal, and left heart disease, yet the underlying mechanisms that govern (mal)adaptation remain incompletely characterized. Owing to the knowledge gaps in our understanding of the right ventricle (RV) in health and disease, the National Heart, Lung, and Blood Institute (NHLBI) commissioned a working group to identify current challenges in the field. These included a need to define and standardize normal RV structure and function in populations; access to RV tissue for research purposes and the development of complex experimental platforms that recapitulate the in vivo environment; and the advancement of imaging and invasive methodologies to study the RV within basic, translational, and clinical research programs. Specific recommendations were provided, including a call to incorporate precision medicine and innovations in prognosis, diagnosis, and novel RV therapeutics for patients with pulmonary vascular disease.

Is Doppler Echocardiography Adequate for Surgical Planning of Single Ventricle Patients?

PURPOSE

Surgical planning has shown great potential for optimizing outcomes for patients affected by single ventricle (SV) malformations. Phase-contrast magnetic resonance imaging (PC-MRI) is the routine technique used for flow acquisition in the surgical planning paradigm. However, PC-MRI may suffer from possible artifacts in certain cases; furthermore, this technology may not be readily available for patients in low and lower-middle-income countries. Therefore, this study aims to investigate the effectiveness of using Doppler echocardiography (echo-Doppler) for flow acquisitions of SV surgical planning.

METHODS

This study included eight patients whose blood flow data was acquired by both PC-MRI and echo-Doppler. A virtual surgery platform was used to generate two surgical options for each patient: (1) a traditional Fontan conduit and (2) a Y-graft. Computational fluid dynamics (CFD) simulations were conducted using the two flow acquisitions to assess clinically relevant hemodynamic metrics: indexed power loss (iPL) and hepatic flow distribution (HFD).

RESULTS

Differences exist in flow data acquired by PC-MRI and echo-Doppler, but no statistical significance was obtained. Flow fields, therefore, exhibit discrepancies between simulations using flow acquisitions by PC-MRI and echo-Doppler. In virtual surgery, the two surgical options were ranked based on these metrics. No difference was observed in the ranking of surgical options between using different flow acquisitions.

CONCLUSION

Doppler echocardiography is an adequate alternative approach to acquire flow data for SV surgical planning. This finding encourages broader usage of SV surgical planning with echo-Doppler when MRI may present artifacts or is not available, especially in low and lower-middle-income countries.

Hemodynamic interplay of vorticity, viscous energy loss, and kinetic energy from 4D Flow MRI and link to cardiac function in healthy subjects and Fontan patients

The purpose of this study was to directly assess (patho)physiology of intraventricular hemodynamic interplay between four-dimensional flow cardiovascular magnetic resonance imaging (4D Flow MRI)-derived vorticity with kinetic energy (KE) and viscous energy loss (EL) over the cardiac cycle and their association to ejection fraction (EF) and stroke volume (SV). Fifteen healthy subjects and thirty Fontan patients underwent whole heart 4D Flow MRI. Ventricular vorticity, KE, and EL were computed over systole (vorticity_vol_{avg systole}, KE_{avg systole}, and EL_{avg systole}) and diastole (vorticity_vol_{avg diastole}, KE_{avg diastole}, and EL_{avg diastole}). The association between vorticity_vol and KE and EL was tested by Spearman correlation. Fontan patients were grouped to normal and impaired EF groups. A significant correlation was found between SV and vorticity in healthy subjects (systolic: ρ = 0.84, P < 0.001; diastolic: ρ = 0.81, P < 0.001) and in Fontan patients (systolic: ρ = 0.61, P < 0.001; diastolic: ρ = 0.54, P = 0.002). Healthy subjects showed positive correlation between vorticity_vol versus KE (systole: ρ = 0.96, P < 0.001; diastole: ρ = 0.90, P < 0.001) and EL (systole: ρ = 0.85, P < 0.001; diastole: ρ = 0.84, P < 0.001). Fontan patients showed significantly elevated vorticity_vol compared with healthy subjects (vorticity_vol_{avg systole}: 3.1 [2.3-3.9] vs. 1.7 [1.3-2.4] L/s, P < 0.001; vorticity_vol_{avg diastole}: 3.1 [2.0-3.7] vs. 2.1 [1.6-2.8] L/s, P = 0.002). This elevated vorticity in Fontan patients showed strong association with KE (systole: ρ = 0.91, P < 0.001; diastole: ρ = 0.85, P < 0.001) and EL (systole: ρ = 0.82, P < 0.001; diastole: ρ = 0.89, P < 0.001). Fontan patients with normal EF showed significantly higher vorticity_vol_{avg systole} and EL_{avg systole}, but significantly decreased KE _{avg diastole}, in the presence of normal SV, compared with healthy subjects. Healthy subjects show strong physiological hemodynamic interplay between vorticity with KE and EL. Fontan patients demonstrate a pathophysiological hemodynamic interplay characterized by correlation of elevated vorticity with KE and EL in the presence of maintained normal stroke volume. Altered vorticity and energetic hemodynamics are found in the presence of normal EF in Fontan patients.NEW & NOTEWORTHY Physiologic intraventricular hemodynamic interplay/coupling is present in the healthy left ventricle between vorticity versus viscous energy loss and kinetic energy from four-dimensional flow cardiovascular magnetic resonance imaging (4D Flow MRI). Conversely, Fontan patients present compensatory pathophysiologic hemodynamic coupling by an increase in intraventricular vorticity that positively correlates to viscous energy loss and kinetic energy levels in the presence of maintained normal stroke volume. Altered vorticity and energetics are found in the presence of normal ejection fraction in Fontan patients.