An interactive simulation tool for patient-specific clinical decision support in single-ventricle physiology

Patient-specific closed-loop model of the fontan circulation: Calibration and validation

The Fontan circulation, designed for managing patients with a single functional ventricle, presents challenges in long-term outcomes. Computational methods offer potential solutions, yet their application in cardiology practice remains largely unexplored. Our aim was to assess the ability of a patient-specific, closed-loop, reduced-order blood flow model to simulate pulsatile blood flow in the Fontan circulation. Using one-dimensional models, we simulated the aorta, superior and inferior venae cavae, and right and left pulmonary arteries, while lumping heart chambers and remaining vessels into zero-dimensional models. The model was calibrated with patient-specific haemodynamic data from combined cardiac catheterisation and magnetic resonance exams, using a novel physics-based stepwise methodology involving simpler open-loop models. Testing on a 10-year-old, anesthetised patient, demonstrated the model's capability to replicate pulsatile pressure and flow in the larger vessels and ventricular pressure. Average relative errors in mean pressure and flow were 2.9 % and 3.6 %, with average relative point-to-point errors (RPPE) in pressure and flow at 5.2 % and 16.0 %. Comparing simulation results to measurements, mean aortic pressure and flow values were 50.7 vs. 50.4 mmHg and 41.6 vs. 41.9 ml/s, respectively, while ventricular pressure values were 28.7 vs. 27.4 mmHg. The model accurately described time-varying ventricular volume with a RPPE of 2.9 %, with mean, minimum, and maximum ventricular volume values for simulation results vs. measurements at 59.2 vs. 58.2 ml, 38.0 vs. 37.6 ml, and 76.0 vs. 74.4 ml, respectively. It provided physiologically realistic predictions of haemodynamic changes from pulmonary vasodilation and atrial fenestration opening. The new model and calibration methodology are freely available, offering a platform to virtually investigate the Fontan circulation's response to clinical interventions and explore potential mechanisms of Fontan failure. Future efforts will concentrate on broadening the model's applicability to a wider range of patient populations and clinical scenarios, as well as testing its effectiveness.

Computational approaches for mechanobiology in cardiovascular development and diseases

The cardiovascular development in vertebrates evolves in response to genetic and mechanical cues. The dynamic interplay among mechanics, cell biology, and anatomy continually shapes the hydraulic networks, characterized by complex, non-linear changes in anatomical structure and blood flow dynamics. To better understand this interplay, a diverse set of molecular and computational tools has been used to comprehensively study cardiovascular mechanobiology. With the continual advancement of computational capacity and numerical techniques, cardiovascular simulation is increasingly vital in both basic science research for understanding developmental mechanisms and disease etiologies, as well as in clinical studies aimed at enhancing treatment outcomes. This review provides an overview of computational cardiovascular modeling. Beginning with the fundamental concepts of computational cardiovascular modeling, it navigates through the applications of computational modeling in investigating mechanobiology during cardiac development. Second, the article illustrates the utility of computational hemodynamic modeling in the context of treatment planning for congenital heart diseases. It then delves into the predictive potential of computational models for elucidating tissue growth and remodeling processes. In closing, we outline prevailing challenges and future prospects, underscoring the transformative impact of computational cardiovascular modeling in reshaping cardiovascular science and clinical practice.

Passive performance evaluation and validation of a viscous impeller pump for subpulmonary fontan circulatory support

Patients with single ventricle defects undergoing the Fontan procedure eventually face Fontan failure. Long-term cavopulmonary assist devices using rotary pump technologies are currently being developed as a subpulmonary power source to prevent and treat Fontan failure. Low hydraulic resistance is a critical safety requirement in the event of pump failure (0 RPM) as a modest 2 mmHg cavopulmonary pressure drop can compromise patient hemodynamics. The goal of this study is therefore to assess the passive performance of a viscous impeller pump (VIP) we are developing for Fontan patients, and validate flow simulations against in-vitro data. Two different blade heights (1.09 mm vs 1.62 mm) and a blank housing model were tested using a mock circulatory loop (MCL) with cardiac output ranging from 3 to 11 L/min. Three-dimensional flow simulations were performed and compared against MCL data. In-silico and MCL results demonstrated a pressure drop of < 2 mmHg at a cardiac output of 7 L/min for both blade heights. There was good agreement between simulation and MCL results for pressure loss (mean difference - 0.23 mmHg 95% CI [0.24-0.71]). Compared to the blank housing model, low wall shear stress area and oscillatory shear index on the pump surface were low, and mean washout times were within 2 s. This study demonstrated the low resistance characteristic of current VIP designs in the failed condition that results in clinically acceptable minimal pressure loss without increased washout time as compared to a blank housing model under normal cardiac output in Fontan patients.

Monitoring variability in parameter estimates for lumped parameter models of the systemic circulation using longitudinal hemodynamic measurements

BACKGROUND

Physics-based cardiovascular models are only recently being considered for disease diagnosis or prognosis in clinical settings. These models depend on parameters representing the physical and physiological properties of the modeled system. Personalizing these parameters may give insight into the specific state of the individual and etiology of disease. We applied a relatively fast model optimization scheme based on common local optimization methods to two model formulations of the left ventricle and systemic circulation. One closed-loop model and one open-loop model were applied. Intermittently collected hemodynamic data from an exercise motivation study were used to personalize these models for data from 25 participants. The hemodynamic data were collected for each participant at the start, middle and end of the trial. We constructed two data sets for the participants, both consisting of systolic and diastolic brachial pressure, stroke volume, and left-ventricular outflow tract velocity traces paired with either the finger arterial pressure waveform or the carotid pressure waveform.

RESULTS

We examined the feasibility of separating parameter estimates for the individual from population estimates by assessing the variability of estimates using the interquartile range. We found that the estimated parameter values were similar for the two model formulations, but that the systemic arterial compliance was significantly different ([Formula: see text]) depending on choice of pressure waveform. The estimates of systemic arterial compliance were on average higher when using the finger artery pressure waveform as compared to the carotid waveform.

CONCLUSIONS

We found that for the majority of participants, the variability of parameter estimates for a given participant on any measurement day was lower than the variability both across all measurement days combined for one participant, and for the population. This indicates that it is possible to identify individuals from the population, and that we can distinguish different measurement days for the individual participant by parameter values using the presented optimization method.

svMorph: Interactive Geometry-Editing Tools for Virtual Patient-Specific Vascular Anatomies

We propose svMorph, a framework for interactive virtual sculpting of patient-specific vascular anatomic models. Our framework includes three tools for the creation of tortuosity, aneurysms, and stenoses in tubular vascular geometries. These shape edits are performed via geometric operations on the surface mesh and vessel centerline curves of the input model. The tortuosity tool also uses the physics-based Oriented Particles method, coupled with linear blend skinning, to achieve smooth, elastic-like deformations. Our tools can be applied separately or in combination to produce simulation-suitable morphed models. They are also compatible with popular vascular modeling software, such as simvascular. To illustrate our tools, we morph several image-based, patient-specific models to create a range of shape changes and simulate the resulting hemodynamics via three-dimensional, computational fluid dynamics. We also demonstrate the ability to quickly estimate the hemodynamic effects of the shape changes via the automated generation of associated zero-dimensional lumped-parameter models.

Optimal Fenestration of the Fontan Circulation

In this paper, we develop a pulsatile compartmental model of the Fontan circulation and use it to explore the effects of a fenestration added to this physiology. A fenestration is a shunt between the systemic and pulmonary veins that is added either at the time of Fontan conversion or at a later time for the treatment of complications. This shunt increases cardiac output and decreases systemic venous pressure. However, these hemodynamic benefits are achieved at the expense of a decrease in the arterial oxygen saturation. The model developed in this paper incorporates fenestration size as a parameter and describes both blood flow and oxygen transport. It is calibrated to clinical data from Fontan patients, and we use it to study the impact of a fenestration on several hemodynamic variables, including systemic oxygen availability, effective oxygen availability, and systemic venous pressure. In certain scenarios corresponding to high-risk Fontan physiology, we demonstrate the existence of a range of fenestration sizes in which the systemic oxygen availability remains relatively constant while the systemic venous pressure decreases.

Interactive, Visual Simulation of a Spatio-Temporal Model of Gas Exchange in the Human Alveolus

In interdisciplinary fields such as systems biology, good communication between experimentalists and theorists is crucial for the success of a project. Theoretical modeling in physiology usually describes complex systems with many interdependencies. On one hand, these models have to be grounded on experimental data. On the other hand, experimenters must be able to understand the interdependent complexities of the theoretical model in order to interpret the model’s results in the physiological context. We promote interactive, visual simulations as an engaging way to present theoretical models in physiology and to make complex processes tangible. Based on a requirements analysis, we developed a new model for gas exchange in the human alveolus in combination with an interactive simulation software named Alvin. Alvin exceeds the current standard with its spatio-temporal resolution and a combination of visual and quantitative feedback. In Alvin, the course of the simulation can be traced in a three-dimensional rendering of an alveolus and dynamic plots. The user can interact by configuring essential model parameters. Alvin allows to run and compare multiple simulation instances simultaneously. We exemplified the use of Alvin for research by identifying unknown dependencies in published experimental data. Employing a detailed questionnaire, we showed the benefits of Alvin for education. We postulate that interactive, visual simulation of theoretical models, as we have implemented with Alvin on respiratory processes in the alveolus, can be of great help for communication between specialists and thereby advancing research.

Engineering Perspective on Cardiovascular Simulations of Fontan Hemodynamics: Where Do We Stand with a Look Towards Clinical Application

BACKGROUND

Cardiovascular simulations for patients with single ventricles undergoing the Fontan procedure can assess patient-specific hemodynamics, explore surgical advances, and develop personalized strategies for surgery and patient care. These simulations have not yet been broadly accepted as a routine clinical tool owing to a number of limitations. Numerous approaches have been explored to seek innovative solutions for improving methodologies and eliminating these limitations.

PURPOSE

This article first reviews the current state of cardiovascular simulations of Fontan hemodynamics. Then, it will discuss the technical progress of Fontan simulations with the emphasis of its clinical impact, noting that substantial improvements have been made in the considerations of patient-specific anatomy, flow, and blood rheology. The article concludes with insights into potential future directions involving clinical validation, uncertainty quantification, and computational efficiency. The advancements in these aspects could promote the clinical usage of Fontan simulations, facilitating its integration into routine clinical practice.

A Tribute to Ajit Yoganathan's Cardiovascular Fluid Mechanics Lab: A Survey of Its Contributions to Our Understanding of the Physiology and Management of Single-Ventricle Patients

INTRODUCTION

Among patients with congenital heart disease, those born with only a single working ventricle represent a particularly complex sub-population, typically requiring multiple surgeries and suffering from high levels of mortality and morbidity. Their cardiac care is complex and has evolved considerably since surgical palliation was first introduced more than 50 years ago. Improvements in treatment have been driven both by growing clinical experience and by knowledge gained through experimental and computational studies of blood flow in these patients. The Cardiovascular Fluid Mechanics Lab at the Georgia Institute of Technology, founded 30 years ago by Dr. Ajit Yoganathan, has pioneered work in this field.

METHODS

In this review, key contributions of Dr. Yoganathan's Cardiovascular Fluid Dynamics Lab are surveyed, including experimental flow loop studies as well as computational fluid dynamics analyses that address many of the critical challenges that cardiologists and surgeons face in treating these patients, including how to reconstruct cardiovascular anatomy to minimize power loss, balance blood flow distribution at key bifurcation points, and avoid other unfavorable hemodynamic conditions.

CONCLUSIONS

Among many contributions in this field, work from the Cardiovascular Fluid Mechanics Lab has led to novel medical devices and patient-specific computational modeling workflows and software tools. These key contributions from this group have enhanced our understanding of the physiology and management of single-ventricle patients.

Computational Modeling to Support Surgical Decision Making in Single Ventricle Physiology

Many of the advances in congenital heart surgery were built upon lessons and insights gained from model simulations. While animal and mock-circuit models have historically been the main arena to test new operative techniques and concepts, the recognition that complex cardiovascular anatomy and circulation can be modeled mathematically ushered a new era of collaboration between surgeons and engineers. In 1996, the computational age in congenital heart surgery began when investigators in London and Milan tapped the power of the computer to simulate the Fontan procedure and introduced operative improvements. Since then, computational modeling has led to numerous contributions in congenial heart surgery as continuing sophistication and advances in numerical and imaging methods furthered the ability to refine anatomic and physiologic details. Idealized generic models have given way to precise patient-specific simulations of the 3-dimensional anatomy, reconstructed circulation, affected hemodynamics, and altered physiology. Tools to perform virtual surgery, and predict flow dynamic and circulatory results, have been developed for some of the most complex defects, such as those requiring single ventricle palliation. In today's quest for personalized medicine and precision care, computational modeling's role to assist surgical planning in complex congenital heart surgery will continue to grow and evolve. With ever closer collaboration between surgeons and engineers, and clear understanding of modeling limitations, computational simulations can be a valuable adjunct to support preoperative surgical decision making.

A personalized computational model of edema formation in myocarditis based on long-axis biventricular MRI images

Background

Myocarditis is defined as the inflammation of the myocardium, i.e. the cardiac muscle. Among the reasons that lead to this disease, we may include infections caused by a virus, bacteria, protozoa, fungus, and others. One of the signs of the inflammation is the formation of edema, which may be a consequence of the interaction between interstitial fluid dynamics and immune response. This complex physiological process was mathematically modeled using a nonlinear system of partial differential equations (PDE) based on porous media approach. By combing a model based on Biot’s poroelasticity theory with a model for the immune response we developed a new hydro-mechanical model for inflammatory edema. To verify this new computational model, T2 parametric mapping obtained by Magnetic Resonance (MR) imaging was used to identify the region of edema in a patient diagnosed with unspecific myocarditis.

Results

A patient-specific geometrical model was created using MRI images from the patient with myocarditis. With this model, edema formation was simulated using the proposed hydro-mechanical mathematical model in a two-dimensional domain. The computer simulations allowed us to correlate spatiotemporal dynamics of representative cells of the immune systems, such as leucocytes and the pathogen, with fluid accumulation and cardiac tissue deformation.

Conclusions

This study demonstrates that the proposed mathematical model is a very promising tool to better understand edema formation in myocarditis. Simulations obtained from a patient-specific model reproduced important aspects related to the formation of cardiac edema, its area, position, and shape, and how these features are related to immune response.

Preload Sensitivity with TORVAD Counterpulse Support Prevents Suction and Overpumping

PURPOSE

This study compares preload sensitivity of continuous flow (CF) VAD support to counterpulsation using the Windmill toroidal VAD (TORVAD). The TORVAD is a two-piston rotary pump that ejects 30 mL in early diastole, which increases cardiac output while preserving aortic valve flow.

METHODS

Preload sensitivity was compared for CF vs. TORVAD counterpulse support using two lumped parameter models of the cardiovascular system: (1) an open-loop model of the systemic circulation was used to obtain ventricular function curves by isolating the systemic circulation and prescribing preload and afterload boundary conditions, and (2) a closed-loop model was used to test the physiological response to changes in pulmonary vascular resistance, systemic vascular resistance, heart rate, inotropic state, and blood volume. In the open-loop model, ventricular function curves (cardiac output vs left ventricular preload) are used to assess preload sensitivity. In the closed-loop model, left ventricular end systolic volume is used to assess the risk of left ventricular suction.

RESULTS

At low preloads of 5 mmHg, CF support overpumps the circulation compared to TORVAD counterpulse support (cardiac output of 3.3 L/min for the healthy heart, 4.7 with CF support, and 3.5 with TORVAD counterpulse support) and has much less sensitivity than counterpulse support (0.342 L/min/mmHg for the healthy heart, 0.092 with CF support, and 0.306 with TORVAD counterpulse support). In the closed-loop model, when PVR is increased beyond 0.035 mmHg s/mL, CF support overpumps the circulation and causes ventricular suction events, but TORVAD counterpulse support maintains sufficient ventricular volume and does not cause suction.

CONCLUSIONS

Counterpulse support with the TORVAD preserves aortic valve flow and provides physiological sensitivity across all preload conditions. This should prevent overpumping and minimize the risk of suction.

3D bioprinting of vascular conduits for pediatric congenital heart repairs

In children with congenital heart defects, surgical correction often involves the use of valves, patches or vascular conduits to establish anatomic continuity. Due to the differences between the pediatric and adult populations, tissue reconstruction in pediatric patients requires a substantially different approach from those in adults. Cardiovascular anatomy of children with congenital heart defect vary, which requires tailored surgical operations for each patient. Since grafts used in these palliative surgeries are sensitive to the local hemodynamic environments, their geometries need to be precisely designed to ensure long-term performance. Tissue engineered vascular grafts (TEVGs) have made tremendous progress over the past decade, but it remains difficult to fabricate patient- and operation-specific vascular grafts. This review summarizes historical milestones of TEVG development for repairing pediatric congenital defects and current clinical outcomes. We also highlight ongoing works on 3D bioprinting of TEVGs with complex geometries and address the current limitations of each technique. Although 3D bioprinted vascular grafts with appropriate functions are yet to be developed, some of the current researches are promising to create better patient specific tissue engineered vascular grafts in the future.

Detailed structural assessment of healthy interventricular septum in the presence of remodeling infarct in the free wall - A finite element model

Purpose

Computational modelling may improve the fundamental understanding of various mechanisms of diseases more particularly related to clinical challenges. In this study the effect of remodeling infarct presence in the left ventricle on the interventricular septal wall is studied using the finite element methods.

Methods

In this study, two rat heart (one model with healthy myocardium and one model with remodeling free wall and healthy septal wall) with magnetic resonance imaging data was gathered to reconstruct three-dimensional (3D) rat heart models. 3D data points from Segment® were imported into SolidEdge® for creation of 3D rat heart models. Abaqus® was used for finite element modeling.

Results

The strain in the healthy interventricular septum of the infarcted left ventricle wall increased when compared to the healthy interventricular septum in the healthy left ventricle. Similarly, the average stress in the healthy left ventricle was observed to have increased on the healthy the interventricular septum where the free wall is subjected to remodeling infarct. When comparing the infarcted models to the healthy model, it was found that the average strain had greatly increased by up to 50.0 %.

Conclusions

The remodeling infarct in the left ventricle has an impact on the healthy interventricular septal wall. Even though the interventricular septal wall was modelled as healthy, it was observed that it has undergone considerable changes in stresses and strains in circumferential and longitudinal direction. The observed changes in myocardial stresses and strains may result in poor global functioning of the heart.

A Hybrid Experimental-Computational Modeling Framework for Cardiovascular Device Testing

Significant advances in biomedical science often leverage powerful computational and experimental modeling platforms. We present a framework named physiology simulation coupled experiment ("PSCOPE") that can capitalize on the strengths of both types of platforms in a single hybrid model. PSCOPE uses an iterative method to couple an in vitro mock circuit to a lumped-parameter numerical simulation of physiology, obtaining closed-loop feedback between the two. We first compared the results of Fontan graft obstruction scenarios modeled using both PSCOPE and an established multiscale computational fluid dynamics method; the normalized root-mean-square error values of important physiologic parameters were between 0.1% and 2.1%, confirming the fidelity of the PSCOPE framework. Next, we demonstrate an example application of PSCOPE to model a scenario beyond the current capabilities of multiscale computational methods-the implantation of a Jarvik 2000 blood pump for cavopulmonary support in the single-ventricle circulation; we found that the commercial Jarvik 2000 controller can be modified to produce a suitable rotor speed for augmenting cardiac output by approximately 20% while maintaining blood pressures within safe ranges. The unified modeling framework enables a testing environment which simultaneously operates a medical device and performs computational simulations of the resulting physiology, providing a tool for physically testing medical devices with simulated physiologic feedback.

An overview of mechanical circulatory support in single-ventricle patients

The population of people with a single-ventricle is continually increasing due to improvements across the spectrum of medical care. Unfortunately, a proportion of these patients will develop heart failure. Often, for these patients, mechanical circulatory support (MCS) represents the only available treatment option. While single-ventricle patients currently represent a small proportion of the total number of patients who receive MCS, as the single-ventricle patient population increases, this number will increase as well. Outcomes for these complex single-ventricle patients who require MCS has begun to be evaluated. When considering the entire population, survival to hospital discharge is 30-50%, though this must be considered with the significant heterogeneity of the single-ventricle patient population. Patients with a single-ventricle have unique anatomy, mechanisms of failure, indications for MCS and the type of support utilized. This has made the interpretation and the generalizability of the limited available data difficult. It is likely that some subsets will have a significantly worse prognosis and others a better one. Unfortunately, with these limited data, indications of a favorable or poor outcome have not yet been elucidated. Though currently, a database has been constructed to address this issue. While the outcomes for these complex patients is unclear, at least in some situations, they are poor. However, significant advances may provide improvements going forward, including new devices, computer simulations and 3D printed models. The most important factor, however, will be the increased experience gained by the heart failure team to improve patient selection, timing, device and configuration selection and operative approach.

Toward a patient-specific tissue engineered vascular graft

Integrating three-dimensional printing with the creation of tissue-engineered vascular grafts could provide a readily available, patient-specific, autologous tissue source that could significantly improve outcomes in newborns with congenital heart disease. Here, we present the recent case of a candidate for our tissue-engineered vascular graft clinical trial deemed ineligible due to complex anatomical requirements and consider the application of three-dimensional printing technologies for a patient-specific graft. We 3D-printed a closed-disposable seeding device and validated that it performed equivalently to the traditional open seeding technique using ovine bone marrow–derived mononuclear cells. Next, our candidate’s preoperative imaging was reviewed to propose a patient-specific graft. A seeding apparatus was then designed to accommodate the custom graft and 3D-printed on a commodity fused deposition modeler. This exploratory feasibility study represents an important proof of concept advancing progress toward a rationally designed patient-specific tissue-engineered vascular graft for clinical application.