Natalia L. Rubio

Luca Pegolotti, Martin R. Pfaller, Natalia L. Rubio et al.

Reduced-order models based on physics are a popular choice in cardiovascular modeling due to their efficiency, but they may experience reduced 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 2% and 3% for pressure and flow rate, respectively, 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.

Natalia L. Rubio, Eric F. Darve, Alison L. Marsden · stanford

Physics-based 0D reduced-order models provide computationally lightweight predictions of cardiovascular flows, resolving bulk hemodynamics in fractions of a second that would take days to solve using traditional 3D finite-element techniques. However, the accuracy of 0D models is limited as a result of the dramatic simplifications made in their derivations. In this work, we use 0D parameters learned from high-fidelity 3D data to improve 0D model accuracy without sacrificing its low computational cost or interpretability. We use the resistor-quadratic resistor-inductor (RRI) model to predict pressure drops over 0D vessels and bifurcations, where the resistances and inductance (0D parameters) are predicted from the bifurcation or vessel geometry using neural networks trained on high-fidelity 3D simulations. We validate the hybrid physics-based data-driven framework in three types of patient-specific vasculature - aortic, aortofemoral, and pulmonary anatomies. Use of learned 0D parameters reduces error by at least 50% compared to baseline 0D parameters across all anatomical cohorts. The improvements are especially marked for the more complex pulmonary anatomies, where 0D models with learned parameters reduced error from 30% to 7%. Exclusion of the quadratic resistor in the RRI model improved convergence compared to using the full RRI model. The resulting hybrid model presents a means of real-time (personal laptop runtime of <2 seconds for the most complex pulmonary anatomies), interpretable, and accurate cardiovascular flow modeling, enabling digital twins that support clinical decision-making as well as cardiovascular science and engineering research.

Natalia L. Rubio, Eric F. Darve, Alison L. Marsden

Three-dimensional (3D) finite-element simulations of cardiovascular flows provide high-fidelity predictions to support cardiovascular medicine, but their high computational cost limits clinical practicality. Reduced-order models (ROMs) offer computationally efficient alternatives but suffer reduced accuracy, particularly at vessel bifurcations where complex flow physics are inadequately captured by standard Poiseuille flow assumptions. We present an enhanced numerical framework that integrates machine learning-predicted bifurcation coefficients into zero-dimensional (0D) hemodynamic ROMs to improve accuracy while maintaining computational efficiency. We develop a resistor-resistor-inductor (RRI) model that uses neural networks to predict pressure-flow relationships from bifurcation geometry, incorporating linear and quadratic resistances along with inductive effects. The method employs non-dimensionalization to reduce training data requirements and apriori flow split prediction for improved bifurcation characterization. We incorporate the RRI model into a 0D model using an optimization-based solution strategy. We validate the approach in isolated bifurcations and vascular trees, across Reynolds numbers from 0 to 5,500, defining ROM accuracy by comparison to 3D finite element simulation. Results demonstrate substantial accuracy improvements: averaged across all trees and Reynolds numbers, the RRI method reduces inlet pressure errors from 54 mmHg (45%) for standard 0D models to 25 mmHg (17%), while a simplified resistor-inductor (RI) variant achieves 31 mmHg (26%) error. The enhanced 0D models show particular effectiveness at high Reynolds numbers and in extensive vascular networks. This hybrid numerical approach enables accurate, real-time hemodynamic modeling for clinical decision support, uncertainty quantification, and digital twins in cardiovascular biomedical engineering.