Learning reduced-order models for cardiovascular simulations with graph neural networks.

Reduced-order models based on physics are a popular choice in cardiovascular modeling due to their efficiency, but they may experience loss in accuracy when working with anatomies that contain numerous junctions or pathological conditions. We develop one-dimensional reduced-order models that simulate blood flow dynamics using a graph neural network trained on three-dimensional hemodynamic simulation data. Given the initial condition of the system, the network iteratively predicts the pressure and flow rate at the vessel centerline nodes. Our numerical results demonstrate the accuracy and generalizability of our method in physiological geometries comprising a variety of anatomies and boundary conditions. Our findings demonstrate that our approach can achieve errors below 3% for pressure and flow rate, provided there is adequate training data. As a result, our method exhibits superior performance compared to physics-based one-dimensional models while maintaining high efficiency at inference time.

Beyond CFD: Emerging methodologies for predictive simulation in cardiovascular health and disease.

Physics-based computational models of the cardiovascular system are increasingly used to simulate hemodynamics, tissue mechanics, and physiology in evolving healthy and diseased states. While predictive models using computational fluid dynamics (CFD) originated primarily for use in surgical planning, their application now extends well beyond this purpose. In this review, we describe an increasingly wide range of modeling applications aimed at uncovering fundamental mechanisms of disease progression and development, performing model-guided design, and generating testable hypotheses to drive targeted experiments. Increasingly, models are incorporating multiple physical processes spanning a wide range of time and length scales in the heart and vasculature. With these expanded capabilities, clinical adoption of patient-specific modeling in congenital and acquired cardiovascular disease is also increasing, impacting clinical care and treatment decisions in complex congenital heart disease, coronary artery disease, vascular surgery, pulmonary artery disease, and medical device design. In support of these efforts, we discuss recent advances in modeling methodology, which are most impactful when driven by clinical needs. We describe pivotal recent developments in image processing, fluid-structure interaction, modeling under uncertainty, and reduced order modeling to enable simulations in clinically relevant timeframes. In all these areas, we argue that traditional CFD alone is insufficient to tackle increasingly complex clinical and biological problems across scales and systems. Rather, CFD should be coupled with appropriate multiscale biological, physical, and physiological models needed to produce comprehensive, impactful models of mechanobiological systems and complex clinical scenarios. With this perspective, we finally outline open problems and future challenges in the field.

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.

Deep Neural Network to Accurately Predict Left Ventricular Systolic Function Under Mechanical Assistance.

Characterizing left ventricle (LV) systolic function in the presence of an LV assist device (LVAD) is extremely challenging. We developed a framework comprising a deep neural network (DNN) and a 0D model of the cardiovascular system to predict parameters of LV systolic function. DNN input data were systemic and pulmonary arterial pressure signals, and rotation speeds of the device. Output data were parameters of LV systolic function, including end-systolic maximal elastance (E max,lv ), a variable essential for adequate hemodynamic assessment of the LV. A 0D model of the cardiovascular system, including a wide range of LVAD settings and incorporating the whole spectrum of heart failure, was used to generate data for the training procedure of the DNN. The DNN predicted E max,lv with a mean relative error of 10.1%, and all other parameters of LV function with a mean relative error of <13%. The framework was then able to retrieve a number of LV physiological variables (i.e., pressures, volumes, and ejection fraction) with a mean relative error of <5%. Our method provides an innovative tool to assess LV hemodynamics under device assistance, which could be helpful for a better understanding of LV-LVAD interactions, and for therapeutic optimization.

Model order reduction of flow based on a modular geometrical approximation of blood vessels.

We are interested in a reduced order method for the efficient simulation of blood flow in arteries. The blood dynamics is modeled by means of the incompressible Navier-Stokes equations. Our algorithm is based on an approximated domain-decomposition of the target geometry into a number of subdomains obtained from the parametrized deformation of geometrical building blocks (e.g., straight tubes and model bifurcations). On each of these building blocks, we build a set of spectral functions by Proper Orthogonal Decomposition of a large number of snapshots of finite element solutions (offline phase). The global solution of the Navier-Stokes equations on a target geometry is then found by coupling linear combinations of these local basis functions by means of spectral Lagrange multipliers (online phase). Being that the number of reduced degrees of freedom is considerably smaller than their finite element counterpart, this approach allows us to significantly decrease the size of the linear system to be solved in each iteration of the Newton-Raphson algorithm. We achieve large speedups with respect to the full order simulation (in our numerical experiments, the gain is at least of one order of magnitude and grows inversely with respect to the reduced basis size), whilst still retaining satisfactory accuracy for most cardiovascular simulations.

Implementation and Calibration of a Deep Neural Network to Predict Parameters of Left Ventricular Systolic Function Based on Pulmonary and Systemic Arterial Pressure Signals.

The evaluation of cardiac contractility by the assessment of the ventricular systolic elastance function is clinically challenging and cannot be easily obtained at the bedside. In this work, we present a framework characterizing left ventricular systolic function from clinically readily available data, including systemic and pulmonary arterial pressure signals. We implemented and calibrated a deep neural network (DNN) consisting of a multi-layer perceptron with 4 fully connected hidden layers and with 16 neurons per layer, which was trained with data obtained from a lumped model of the cardiovascular system modeling different levels of cardiac function. The lumped model included a function of circulatory autoregulation from carotid baroreceptors in pulsatile conditions. Inputs for the DNN were systemic and pulmonary arterial pressure curves. Outputs from the DNN were parameters of the lumped model characterizing left ventricular systolic function, especially end-systolic elastance. The DNN adequately performed and accurately recovered the relevant hemodynamic parameters with a mean relative error of less than 2%. Therefore, our framework can easily provide complex physiological parameters of cardiac contractility, which could lead to the development of invaluable tools for the clinical evaluation of patients with severe cardiac dysfunction.