Automated generation of 0D and 1D reduced-order models of patient-specific blood flow.

Three-dimensional (3D) cardiovascular fluid dynamics simulations typically require hours to days of computing time on a high-performance computing cluster. One-dimensional (1D) and lumped-parameter zero-dimensional (0D) models show great promise for accurately predicting blood bulk flow and pressure waveforms with only a fraction of the cost. They can also accelerate uncertainty quantification, optimization, and design parameterization studies. Despite several prior studies generating 1D and 0D models and comparing them to 3D solutions, these were typically limited to either 1D or 0D and a singular category of vascular anatomies. This work proposes a fully automated and openly available framework to generate and simulate 1D and 0D models from 3D patient-specific geometries, automatically detecting vessel junctions and stenosis segments. Our only input is the 3D geometry; we do not use any prior knowledge from 3D simulations. All computational tools presented in this work are implemented in the open-source software platform SimVascular. We demonstrate the reduced-order approximation quality against rigid-wall 3D solutions in a comprehensive comparison with N = 72 publicly available models from various anatomies, vessel types, and disease conditions. Relative average approximation errors of flows and pressures typically ranged from 1% to 10% for both 1D and 0D models, measured at the outlets of terminal vessel branches. In general, 0D model errors were only slightly higher than 1D model errors despite requiring only a third of the 1D runtime. Automatically generated ROMs can significantly speed up model development and shift the computational load from high-performance machines to personal computers.

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.

Hemodynamic performance of tissue-engineered vascular grafts in Fontan patients.

In the field of congenital heart surgery, tissue-engineered vascular grafts (TEVGs) are a promising alternative to traditionally used synthetic grafts. Our group has pioneered the use of TEVGs as a conduit between the inferior vena cava and the pulmonary arteries in the Fontan operation. The natural history of graft remodeling and its effect on hemodynamic performance has not been well characterized. In this study, we provide a detailed analysis of the first U.S. clinical trial evaluating TEVGs in the treatment of congenital heart disease. We show two distinct phases of graft remodeling: an early phase distinguished by rapid changes in graft geometry and a second phase of sustained growth and decreased graft stiffness. Using clinically informed and patient-specific computational fluid dynamics (CFD) simulations, we demonstrate how changes to TEVG geometry, thickness, and stiffness affect patient hemodynamics. We show that metrics of patient hemodynamics remain within normal ranges despite clinically observed levels of graft narrowing. These insights strengthen the continued clinical evaluation of this technology while supporting recent indications that reversible graft narrowing can be well tolerated, thus suggesting caution before intervening clinically.

Patient-Specific Multiscale Modeling of the Assisted Bidirectional Glenn.

BACKGROUND: First-stage palliation of neonates with single-ventricle physiology is associated with poor outcomes and challenging clinical management. Prior computational modeling and in vitro experiments introduced the assisted bidirectional Glenn (ABG), which increased pulmonary flow and oxygenation over the bidirectional Glenn (BDG) and the systemic-to-pulmonary shunt in idealized models. In this study, we demonstrate that the ABG achieves similar performance in patient-specific models and assess the influence of varying shunt geometry. METHODS: In a small cohort of single-ventricle prestage 2 patients, we constructed three-dimensional in silico models and tuned lumped parameter networks to match clinical measurements. Each model was modified to produce virtual BDG and ABG surgeries. We simulated the hemodynamics of the stage 1 procedure, BDG, and ABG by using multiscale computational modeling, coupling a finite-element flow solver to the lumped parameter network. Two levels of pulmonary vascular resistances (PVRs) were investigated: baseline (low) PVR of the patients and doubled (high) PVR. The shunt nozzle diameter, anastomosis location, and shape were also manipulated. RESULTS: The ABG increased the pulmonary flow rate and pressure by 15% to 20%, which was accompanied by a rise in superior vena caval pressure (2 to 3 mm Hg) at both PVR values. Pulmonary flow rate and superior vena caval pressures were most sensitive to the shunt nozzle diameter. CONCLUSIONS: Patient-specific ABG performance was similar to prior idealized simulations and experiments, with good performance at lower PVR values in the range of measured clinical data. Larger shunt outlet diameters and lower PVR led to improved ABG performance.

Optimization of the Assisted Bidirectional Glenn Procedure for First Stage Single Ventricle Repair.

BACKGROUND: First-stage single-ventricle palliation is challenging to manage, and significant interstage morbidity and mortality remain. Prior computational and in vitro studies of the assisted bidirectional Glenn (ABG), a novel first-stage procedure that has shown potential for early conversion to a more stable augmented Glenn physiology, demonstrated increased pulmonary flow and oxygen delivery while decreasing cardiac work, as compared to conventional stage-1 alternatives. This study aims to identify optimal shunt designs for the ABG to improve pulmonary flow while maintaining or decreasing superior vena caval (SVC) pressure. METHODS: A representative three-dimensional model of a neonatal bidirectional Glenn (BDG) was created, with a shunt connecting the innominate artery to the SVC. The shunt design was studied as a six-parameter constrained shape optimization problem. We simulated hemodynamics for each candidate designs using a multiscale finite element flow solver and compared performance against designs with taper-less shunts, the standalone BDG, and a simplified control volume model. Three values of pulmonary vascular resistance (PVR) of 2.3, 4.3, and 7.1 WUm2 were studied. RESULTS: Increases in pulmonary flow were generally accompanied by increases in SVC pressure, except at low PVR (2.3 WUm2), where the optimal shunt geometry achieved a 13% increase in pulmonary flow without incurring any increase in SVC pressure. Shunt outlet area was the most influential design parameter, while others had minimal effect. CONCLUSION: Assisted bidirectional Glenn performance is sensitive to PVR and shunt outlet diameter. An increase in pulmonary flow without a corresponding increase in SVC pressure is possible only when PVR is low.