Uncertainty Quantification in a Patient-Specific One-Dimensional Arterial Network Model: EnKF-Based Inflow Estimator

Reducing Geometric Uncertainty in Computational Hemodynamics by Deep Learning-Assisted Parallel-Chain MCMC

Computational hemodynamic modeling has been widely used in cardiovascular research and healthcare. However, the reliability of model predictions is largely dependent on the uncertainties of modeling parameters and boundary conditions, which should be carefully quantified and further reduced with available measurements. In this work, we focus on propagating and reducing the uncertainty of vascular geometries within a Bayesian framework. A novel deep learning (DL)-assisted parallel Markov chain Monte Carlo (MCMC) method is presented to enable efficient Bayesian posterior sampling and geometric uncertainty reduction. A DL model is built to approximate the geometry-to-hemodynamic map, which is trained actively using online data collected from parallel MCMC chains and utilized for early rejection of unlikely proposals to facilitate convergence with less expensive full-order model evaluations. Numerical studies on two-dimensional aortic flows are conducted to demonstrate the effectiveness and merit of the proposed method.

Inverse problems in blood flow modeling: A review

Mathematical and computational modeling of the cardiovascular system is increasingly providing non-invasive alternatives to traditional invasive clinical procedures. Moreover, it has the potential for generating additional diagnostic markers. In blood flow computations, the personalization of spatially distributed (i.e., 3D) models is a key step which relies on the formulation and numerical solution of inverse problems using clinical data, typically medical images for measuring both anatomy and function of the vasculature. In the last years, the development and application of inverse methods has rapidly expanded most likely due to the increased availability of data in clinical centers and the growing interest of modelers and clinicians in collaborating. Therefore, this work aims to provide a wide and comparative overview of literature within the last decade. We review the current state of the art of inverse problems in blood flows, focusing on studies considering fully dimensional fluid and fluid-solid models. The relevant physical models and hemodynamic measurement techniques are introduced, followed by a survey of mathematical data assimilation approaches used to solve different kinds of inverse problems, namely state and parameter estimation. An exhaustive discussion of the literature of the last decade is presented, structured by types of problems, models and available data.

Data Assimilation by Stochastic Ensemble Kalman Filtering to Enhance Turbulent Cardiovascular Flow Data From Under-Resolved Observations

We propose a data assimilation methodology that can be used to enhance the spatial and temporal resolution of voxel-based data as it may be obtained from biomedical imaging modalities. It can be used to improve the assessment of turbulent blood flow in large vessels by combining observed data with a computational fluid dynamics solver. The methodology is based on a Stochastic Ensemble Kalman Filter (SEnKF) approach and geared toward pulsatile and turbulent flow configurations. We describe the observed flow fields by a mean value and its covariance. These flow fields are combined with forecasts obtained from a direct numerical simulation of the flow field. The method is validated against canonical pulsatile and turbulent flows. Finally, it is applied to a clinically relevant configuration, namely the flow downstream of a bioprosthetic valve in an aorta phantom. It is demonstrated how the 4D flow field obtained from experimental observations can be enhanced by the data assimilation algorithm. Results show that the presented method is promising for future use with in vivo data from 4D Flow Magnetic Resonance Imaging (4D Flow MRI). 4D Flow MRI returns spatially and temporally averaged flow fields that are limited by the spatial and the temporal resolution of the tool. These averaged flow fields and the associated uncertainty might be used as observation data in the context of the proposed methodology.

Data-driven cardiovascular flow modelling: examples and opportunities

High-fidelity blood flow modelling is crucial for enhancing our understanding of cardiovascular disease. Despite significant advances in computational and experimental characterization of blood flow, the knowledge that we can acquire from such investigations remains limited by the presence of uncertainty in parameters, low resolution, and measurement noise. Additionally, extracting useful information from these datasets is challenging. Data-driven modelling techniques have the potential to overcome these challenges and transform cardiovascular flow modelling. Here, we review several data-driven modelling techniques, highlight the common ideas and principles that emerge across numerous such techniques, and provide illustrative examples of how they could be used in the context of cardiovascular fluid mechanics. In particular, we discuss principal component analysis (PCA), robust PCA, compressed sensing, the Kalman filter for data assimilation, low-rank data recovery, and several additional methods for reduced-order modelling of cardiovascular flows, including the dynamic mode decomposition and the sparse identification of nonlinear dynamics. All techniques are presented in the context of cardiovascular flows with simple examples. These data-driven modelling techniques have the potential to transform computational and experimental cardiovascular research, and we discuss challenges and opportunities in applying these techniques in the field, looking ultimately towards data-driven patient-specific blood flow modelling.

Data-Augmented Modeling of Intracranial Pressure

Precise management of patients with cerebral diseases often requires intracranial pressure (ICP) monitoring, which is highly invasive and requires a specialized ICU setting. The ability to noninvasively estimate ICP is highly compelling as an alternative to, or screening for, invasive ICP measurement. Most existing approaches for noninvasive ICP estimation aim to build a regression function that maps noninvasive measurements to an ICP estimate using statistical learning techniques. These data-based approaches have met limited success, likely because the amount of training data needed is onerous for this complex applications. In this work, we discuss an alternative strategy that aims to better utilize noninvasive measurement data by leveraging mechanistic understanding of physiology. Specifically, we developed a Bayesian framework that combines a multiscale model of intracranial physiology with noninvasive measurements of cerebral blood flow using transcranial Doppler. Virtual experiments with synthetic data are conducted to verify and analyze the proposed framework. A preliminary clinical application study on two patients is also performed in which we demonstrate the ability of this method to improve ICP prediction.

Practical identifiability and uncertainty quantification of a pulsatile cardiovascular model

Mathematical models are essential tools to study how the cardiovascular system maintains homeostasis. The utility of such models is limited by the accuracy of their predictions, which can be determined by uncertainty quantification (UQ). A challenge associated with the use of UQ is that many published methods assume that the underlying model is identifiable (e.g. that a one-to-one mapping exists from the parameter space to the model output). In this study we present a novel workflow to calibrate a lumped-parameter model to left ventricular pressure and volume time series data. Key steps include using (1) literature and available data to determine nominal parameter values; (2) sensitivity analysis and subset selection to determine a set of identifiable parameters; (3) optimization to find a point estimate for identifiable parameters; and (4) frequentist and Bayesian UQ calculations to assess the predictive capability of the model. Our results show that it is possible to determine 5 identifiable model parameters that can be estimated to our experimental data from three rats, and that computed UQ intervals capture the measurement and model error.