Computational hemodynamics in arteries with the one-dimensional augmented fluid-structure interaction system: viscoelastic parameters estimation and comparison with in-vivo data

openBF: an open-source finite volume 1D blood flow solver

Computational simulations are widely adopted in cardiovascular biomechanics because of their capability of producing physiological data otherwise impossible to measure with non-invasive modalities.Objective.This study presents openBF, a computational library for simulating the blood dynamics in the cardiovascular system.Approach.openBF adopts a one-dimensional viscoelastic representation of the arterial system, and is coupled with zero-dimensional windkessel models at the outlets. Equations are solved by means of the finite-volume method and the code is written in Julia. We assess its predictions by performing a multiscale validation study on several domains available from the literature.Main results.At all scales, which range from individual arteries to a population of virtual subjects, openBF's solution show excellent agreement with the solutions from existing software. For reported simulations, openBF requires low computational times.Significance.openBF is easy to install, use, and deploy on multiple platforms and architectures, and gives accurate prediction of blood dynamics in short time-frames. It is actively maintained and available open-source on GitHub, which favours contributions from the biomechanical community.

Rapid estimation of left ventricular contractility with a physics-informed neural network inverse modeling approach

Physics-based computer models based on numerical solutions of the governing equations generally cannot make rapid predictions, which in turn limits their applications in the clinic. To address this issue, we developed a physics-informed neural network (PINN) model that encodes the physics of a closed-loop blood circulation system embedding a left ventricle (LV). The PINN model is trained to satisfy a system of ordinary differential equations (ODEs) associated with a lumped parameter description of the circulatory system. The model predictions have a maximum error of less than 5% when compared to those obtained by solving the ODEs numerically. An inverse modeling approach using the PINN model is also developed to rapidly estimate model parameters (in ∼ 3 min) from single-beat LV pressure and volume waveforms. Using synthetic LV pressure and volume waveforms generated by the PINN model with different model parameter values, we show that the inverse modeling approach can recover the corresponding ground truth values for LV contractility indexed by the end-systolic elastance Ees with a 1% error, which suggests that this parameter is unique. The estimated Ees is about 58% to 284% higher for the data associated with dobutamine compared to those without, which implies that this approach can be used to estimate LV contractility using single-beat measurements. The PINN inverse modeling can potentially be used in the clinic to simultaneously estimate LV contractility and other physiological parameters from single-beat measurements.

Computational framework for the generation of one-dimensional vascular models accounting for uncertainty in networks extracted from medical images

One-dimensional (1D) cardiovascular models offer a non-invasive method to answer medical questions, including predictions of wave-reflection, shear stress, functional flow reserve, vascular resistance and compliance. This model type can predict patient-specific outcomes by solving 1D fluid dynamics equations in geometric networks extracted from medical images. However, the inherent uncertainty in in vivo imaging introduces variability in network size and vessel dimensions, affecting haemodynamic predictions. Understanding the influence of variation in image-derived properties is essential to assess the fidelity of model predictions. Numerous programs exist to render three-dimensional surfaces and construct vessel centrelines. Still, there is no exact way to generate vascular trees from the centrelines while accounting for uncertainty in data. This study introduces an innovative framework employing statistical change point analysis to generate labelled trees that encode vessel dimensions and their associated uncertainty from medical images. To test this framework, we explore the impact of uncertainty in 1D haemodynamic predictions in a systemic and pulmonary arterial network. Simulations explore haemodynamic variations resulting from changes in vessel dimensions and segmentation; the latter is achieved by analysing multiple segmentations of the same images. Results demonstrate the importance of accurately defining vessel radii and lengths when generating high-fidelity patient-specific haemodynamics models. KEY POINTS: This study introduces novel algorithms for generating labelled directed trees from medical images, focusing on accurate junction node placement and radius extraction using change points to provide haemodynamic predictions with uncertainty within expected measurement error. Geometric features, such as vessel dimension (length and radius) and network size, significantly impact pressure and flow predictions in both pulmonary and aortic arterial networks. Standardizing networks to a consistent number of vessels is crucial for meaningful comparisons and decreases haemodynamic uncertainty. Change points are valuable to understanding structural transitions in vascular data, providing an automated and efficient way to detect shifts in vessel characteristics and ensure reliable extraction of representative vessel radii.

Morphology of Abdominal Aortic Aneurysms and Correlation with Biomechanical Tests of Aneurysmal Wall Fragments

BACKGROUND

Evaluate how specific morphologic aspects of abdominal aortic aneurysms (AAAs), including asymmetries, curvatures, tortuosities, and angulations, among others can influence the intrinsic biomechanical properties of the AAA's wall. This study analyzed the correlation of geometric measurements (1-dimensional, 2-dimensional, and 3-dimensional) of preoperative tomographic images of AAA with uniaxial biomechanical tests of the arterial wall fragments of these AAA obtained in open surgical repair of aneurysms.

METHODS

It was a multicenter, experimental, and observational study, and initially 54 individuals were selected who underwent open surgical of AAA, with valid biomechanical tests of the anterior wall of the AAA. Seven individuals were excluded because they had poor preoperative quality computed tomography scans and/or artifacts that impeded image segmentation and extraction of AAA geometric indices. The aortic fragments were subjected to uniaxial biomechanical destructive tests to obtain the following data: maximum load, failure stress, failure tension, failure strain energy, strain, and fragment thickness. In the same patients, preoperative computed tomography scans were performed with the extraction of 26 geometric indices, subdivided into 9 1-dimensional indices, 6 2-dimensional indices, and 11 3-dimensional indices. Data were subjected to statistical analysis using SPSS version 28.

RESULTS

Comparing ruptured and unruptured AAA, no statistical difference was observed between the biomechanical and geometric parameters. The fragment thickness of the ruptured AAA was lower than that of the unruptured AAA (P < 0.05). By comparing tomographic geometric indices and biomechanical parameters of the aortic fragments using Pearson's coefficient, positive and linear correlations (P < 0.05) were observed between the geometric variable maximum diameter (Dmax) of the AAA with maximum load (r = 0.408), failure tension (r = 0.372), and failure stress (r = 0.360). Positive and linear correlations were also observed between the variable diameter/height ratio (DHr) and the maximum load (r = 0.360), failure tension (r = 0.354), and failure stress (r = 0.289). The geometric variable DHr was dependent and correlated with Dmax. Simple regression analysis showed that R2 varied between 8.3% and 16.7%, and all models were significant (P < 0.05).

CONCLUSIONS

Dmax and DHr were linearly and positively correlated with the resistance parameters (maximum load, failure tension, and failure stress) of the AAA fragments. The DHr variable is dependent and correlated with Dmax. There was no correlation between the other geometric indices and the biomechanical parameters of the AAA wall. The asymmetries did not globally influence the biomechanics of AAA wall; however, they may influence regionally. Larger AAAs were stronger than smaller ones. Therefore, such findings may point toward Dmax is still the main geometric parameter, which influences the anterior wall, and possibly globally in the AAA.

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.

The effect of cardiac properties on arterial pulse waves: An in-silico study

This study investigated the effects of cardiac properties variability on arterial pulse wave morphology using blood flow modelling and pulse wave analysis. A lumped-parameter model of the left part of the heart was coupled to a one-dimensional model of the arterial network and validated using reference pulse waveforms in turn verified by comparison with in vivo measurements. A sensitivity analysis was performed to assess the effects of variations in cardiac parameters on central and peripheral pulse waveforms. Results showed that left ventricle contractility, stroke volume, cardiac cycle duration, and heart valves impairment are determinants of central waveforms morphology, pulse pressure and its amplification. Contractility of the left atrium has negligible effects on arterial pulse waves. Results also suggested that it might be possible to infer left ventricular dysfunction by analysing the timing of the dicrotic notch and cardiac function by analysing PPG signals. This study has identified cardiac properties that may be extracted from in vivo central and peripheral pulse waves to assess cardiac function.

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.

Determining Clinically-Viable Biomarkers for Ischaemic Stroke Through a Mechanistic and Machine Learning Approach

Assessment of distal cerebral perfusion after ischaemic stroke is currently only possible through expensive and time-consuming imaging procedures which require the injection of a contrast medium. Alternative approaches that could indicate earlier the impact of blood flow occlusion on distal cerebral perfusion are currently lacking. The aim of this study was to identify novel biomarkers suitable for clinical implementation using less invasive diagnostic techniques such as Transcranial Doppler (TCD). We used 1D modelling to simulate pre- and post-stroke velocity and flow wave propagation in a typical arterial network, and Sobol’s sensitivity analysis, supported by the use of Gaussian process emulators, to identify biomarkers linked to cerebral perfusion. We showed that values of pulsatility index of the right anterior cerebral artery > 1.6 are associated with poor perfusion and may require immediate intervention. Three additional biomarkers with similar behaviour, all related to pulsatility indices, were identified. These results suggest that flow pulsatility measured at specific locations could be used to effectively estimate distal cerebral perfusion rates, and ultimately improve clinical diagnosis and management of ischaemic stroke.

A method of parameter estimation for cardiovascular hemodynamics based on deep learning and its application to personalize a reduced-order model

Precise model personalization is a key step towards the application of cardiovascular physical models. In this manuscript, we propose to use deep learning (DL) to solve the parameter estimation problem in cardiovascular hemodynamics. Based on the convolutional neural network (CNN) and fully connected neural network (FCNN), a multi-input deep neural network (DNN) model is developed to map the nonlinear relationship between measurements and the parameters to be estimated. In this model, two separate network structures are designed to extract the features of two types of measurement data, including pressure waveforms and a vector composed of heart rate (HR) and pulse transit time (PTT), and a shared structure is used to extract their combined dependencies on the parameters. Besides, we try to use the transfer learning (TL) technology to further strengthen the personalized characteristics of a trained-well network. For assessing the proposed method, we conducted the parameter estimation using synthetic data and in vitro data respectively, and in the test with synthetic data, we evaluated the performance of the TL algorithm through two individuals with different characteristics. A series of estimation results show that the estimated parameters are in good agreement with the true values. Furthermore, it is also found that the estimation accuracy can be significantly improved by a multicycle combination strategy. Therefore, we think that the proposed method has the potential to be used for parameter estimation in cardiovascular hemodynamics, which can provide an immediate, accurate, and sustainable personalization process, and deserves more attention in the future.

On the Periodicity of Cardiovascular Fluid Dynamics Simulations

Three-dimensional cardiovascular fluid dynamics simulations typically require computation of several cardiac cycles before they reach a periodic solution, rendering them computationally expensive. Furthermore, there is currently no standardized method to determine whether a simulation has yet reached that periodic state. In this work, we propose the use of an asymptotic error measurement to quantify the difference between simulation results and their ideal periodic state using open-loop lumped-parameter modeling. We further show that initial conditions are crucial in reducing computational time and develop an automated framework to generate appropriate initial conditions from a one-dimensional model of blood flow. We demonstrate the performance of our initialization method using six patient-specific models from the Vascular Model Repository. In our examples, our initialization protocol achieves periodic convergence within one or two cardiac cycles, leading to a significant reduction in computational cost compared to standard methods. All computational tools used in this work are implemented in the open-source software platform SimVascular. Automatically generated initial conditions have the potential to significantly reduce computation time in cardiovascular fluid dynamics simulations.