Comparison of Immersed Boundary Simulations of Heart Valve Hemodynamics Against In Vitro 4D Flow MRI Data.

The immersed boundary (IB) method is a mathematical framework for fluid-structure interaction problems (FSI) that was originally developed to simulate flows around heart valves. Direct comparison of FSI simulations around heart valves against experimental data is challenging, however, due to the difficulty of performing robust and effective simulations, the complications of modeling a specific physical experiment, and the need to acquire experimental data that is directly comparable to simulation data. Such comparators are a necessary precursor for further formal validation studies of FSI simulations involving heart valves. In this work, we performed physical experiments of flow through a pulmonary valve in an in vitro pulse duplicator, and measured the corresponding velocity field using 4D flow MRI (4-dimensional flow magnetic resonance imaging). We constructed a computer model of this pulmonary artery setup, including modeling valve geometry and material properties via a technique called design-based elasticity, and simulated flow through it with the IB method. The simulated flow fields showed excellent qualitative agreement with experiments, excellent agreement on integral metrics, and reasonable relative error in the entire flow domain and on slices of interest. These results illustrate how to construct a computational model of a physical experiment for use as a comparator.

A Mechanistic Lumped Parameter Model of the Berlin Heart EXCOR to Analyze Device Performance and Physiologic Interactions.

PURPOSE: The Berlin Heart EXCOR (BH) is the only FDA-approved, extracorporeal pulsatile ventricular assist device (VAD) for infants and children with heart failure. Clinicians control four settings on the device-systolic and diastolic drive pressures, device pump rate, and systolic time as a percentage of the pump cycle. However, interactions between BH pneumatics and the native circulation remain poorly understood. Thus, establishing appropriate device size and settings can be challenging on a patient-to-patient basis. METHODS: In this study we develop a novel lumped parameter network based on simplified device mechanics. We perform parametric studies to characterize device behavior, study interactions between the left ventricle (LV) and BH across different device settings, and develop patient-specific simulations. We then simulate the impact of changing device parameters for each of three patients. RESULTS: Increasing systolic pressure and systolic time increased device output. We identified previously unobserved cycle-to-cycle variations in LV-BH interactions that may impact patient health. Patient-specific simulations demonstrated the model's ability to replicate BH performance, captured trends in LV behavior after device implantation, and emphasized the importance of device rate and volume in optimizing BH efficiency. CONCLUSION: We present a novel, mechanistic model that can be readily adjusted to study a wide range of device settings and clinical scenarios. Physiologic interactions between the BH and the native LV produced large variability in cardiac loading. Our findings showed that operating the BH at a device rate greater than the patient's native heart decreases variability in physiological interactions between the BH and LV, increasing cardiac offloading while maintaining cardiac output. Device rates that are close to the resting heart rate may result in unfavorable cardiac loading conditions. Our work demonstrates the utility of our model to investigate BH performance for patient-specific physiologies.

In Vitro Assessment of Right Ventricular Outflow Tract Anatomy and Valve Orientation Effects on Bioprosthetic Pulmonary Valve Hemodynamics.

PURPOSE: The congenital heart defect Tetralogy of Fallot (ToF) affects 1 in 2500 newborns annually in the US and typically requires surgical repair of the right ventricular outflow tract (RVOT) early in life, with variations in surgical technique leading to large disparities in RVOT anatomy among patients. Subsequently, often in adolescence or early adulthood, patients usually require surgical placement of a xenograft or allograft pulmonary valve prosthesis. Valve longevity is highly variable for reasons that remain poorly understood. METHODS: This work aims to assess the performance of bioprosthetic pulmonary valves in vitro using two 3D printed geometries: an idealized case based on healthy subjects aged 11 to 13 years and a diseased case with a 150% dilation in vessel diameter downstream of the valve. Each geometry was studied with two valve orientations: one with a valve leaflet opening posterior, which is the native pulmonary valve position, and one with a valve leaflet opening anterior. RESULTS: Full three-dimensional, three-component, phase-averaged velocity fields were obtained in the physiological models using 4D flow MRI. Flow features, particularly vortex formation and reversed flow regions, differed significantly between the RVOT geometries and valve orientations. Pronounced asymmetry in streamwise velocity was present in all cases, while the diseased geometry produced additional asymmetry in radial flows. Quantitative integral metrics demonstrated increased secondary flow strength and recirculation in the rotated orientation for the diseased geometry. CONCLUSIONS: The compound effects of geometry and orientation on bioprosthetic valve hemodynamics illustrated in this study could have a crucial impact on long-term valve performance.