Identification of artery wall stiffness: in vitro validation and in vivo results of a data assimilation procedure applied to a 3D fluid-structure interaction model

Guiding patient-specific cardiac simulations through data-assimilation of soft tissue kinematics from dynamic CT scan

Fluid-structure interaction (FSI) can be key in the generation of accurate digital replica of cardiovascular systems. To personalize these models, however, several patient-specific parameters need to be measured, which can be challenging to accomplish in a non-invasive manner. Alternatively, the cardiac kinematics of the patient can be extracted from imaging data and then directly imposed as a dynamic boundary condition in the computational model, also incorporating temporal and spatial measurement errors. A more advanced method combines FSI with kinematic driven simulations using data-assimilation. Despite its potential, the application of this technique to complex multi-physics cardiovascular simulations remains limited. In this study, we develop an FSI model of a patient's left ventricle (LV) and aorta, personalized with dynamic imaging data using a Nudging algorithm-a data assimilation technique-which is tailored to each cardiac chamber. In particular, for the LV, which embeds small-scale and irregular endocardial structures (higher measurement errors), the active contraction of the patient is replicated primarily using integral measurements (ventricular volume and surface area). On the other hand, the passive motion of the aorta is guided in the simulation relying directly on the local tissue positions from CT scan. The algorithm's simplicity and zero additional computational cost make it particularly suitable for multi-physics problems. Our results show that the assimilation procedure must be tuned to guide the system toward the measurements within the uncertainty range of the in-vivo data.

Fluid-Structure Interaction Within Models of Patient-Specific Arteries: Computational Simulations and Experimental Validations

Cardiovascular disease (CVD) is the leading cause of mortality worldwide and its incidence is rising due to an aging population. The development and progression of CVD is directly linked to adverse vascular hemodynamics and biomechanics, whose in-vivo measurement remains challenging but can be simulated numerically and experimentally. The ability to evaluate these parameters in patient-specific CVD cases is crucial to better predict future disease progression, risk of adverse events, and treatment efficacy. While significant progress has been made toward patient-specific hemodynamic simulations, blood vessels are often assumed to be rigid, which does not consider the compliant mechanical properties of vessels whose malfunction is implicated in disease. In an effort to simulate the biomechanics of flexible vessels, fluid-structure interaction (FSI) simulations have emerged as promising tools for the characterization of hemodynamics within patient-specific cardiovascular anatomies. Since FSI simulations combine the blood's fluid domain with the arterial structural domain, they pose novel challenges for their experimental validation. This paper reviews the scientific work related to FSI simulations for patient-specific arterial geometries and the current standard of FSI model validation including the use of compliant arterial phantoms, which offer novel potential for the experimental validation of FSI results.

Uncertainty Quantification in the In Vivo Image-Based Estimation of Local Elastic Properties of Vascular Walls

Introduction: Patient-specific computational models are a powerful tool for planning cardiovascular interventions. However, the in vivo patient-specific mechanical properties of vessels represent a major source of uncertainty. In this study, we investigated the effect of uncertainty in the elastic module (E) on a Fluid–Structure Interaction (FSI) model of a patient-specific aorta. Methods: The image-based χ-method was used to compute the initial E value of the vascular wall. The uncertainty quantification was carried out using the generalized Polynomial Chaos (gPC) expansion technique. The stochastic analysis was based on four deterministic simulations considering four quadrature points. A deviation of about ±20% on the estimation of the E value was assumed. Results: The influence of the uncertain E parameter was evaluated along the cardiac cycle on area and flow variations extracted from five cross-sections of the aortic FSI model. Results of stochastic analysis showed the impact of E in the ascending aorta while an insignificant effect was observed in the descending tract. Conclusions: This study demonstrated the importance of the image-based methodology for inferring E, highlighting the feasibility of retrieving useful additional data and enhancing the reliability of in silico models in clinical practice.

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.

Introduction of a Novel Image-Based and Non-Invasive Method for the Estimation of Local Elastic Properties of Great Vessels

Background: In the context of a growing demand for the use of in silico models to meet clinical requests, image-based methods play a crucial role. In this study, we present a parametric equation able to estimate the elasticity of vessel walls, non-invasively and indirectly, from information uniquely retrievable from imaging. Methods: A custom equation was iteratively refined and tuned from the simulations of a wide range of different vessel models, leading to the definition of an indirect method able to estimate the elastic modulus E of a vessel wall. To test the effectiveness of the predictive capability to infer the E value, two models with increasing complexity were used: a U-shaped vessel and a patient-specific aorta. Results: The original formulation was demonstrated to deviate from the ground truth, with a difference of 89.6%. However, the adoption of our proposed equation was found to significantly increase the reliability of the estimated E value for a vessel wall, with a mean percentage error of 9.3% with respect to the reference values. Conclusion: This study provides a strong basis for the definition of a method able to estimate local mechanical information of vessels from data easily retrievable from imaging, thus potentially increasing the reliability of in silico cardiovascular models.

Patient-Specific Inverse Modeling of In Vivo Cardiovascular Mechanics with Medical Image-Derived Kinematics as Input Data: Concepts, Methods, and Applications

Inverse modeling approaches in cardiovascular medicine are a collection of methodologies that can provide non-invasive patient-specific estimations of tissue properties, mechanical loads, and other mechanics-based risk factors using medical imaging as inputs. Its incorporation into clinical practice has the potential to improve diagnosis and treatment planning with low associated risks and costs. These methods have become available for medical applications mainly due to the continuing development of image-based kinematic techniques, the maturity of the associated theories describing cardiovascular function, and recent progress in computer science, modeling, and simulation engineering. Inverse method applications are multidisciplinary, requiring tailored solutions to the available clinical data, pathology of interest, and available computational resources. Herein, we review biomechanical modeling and simulation principles, methods of solving inverse problems, and techniques for image-based kinematic analysis. In the final section, the major advances in inverse modeling of human cardiovascular mechanics since its early development in the early 2000s are reviewed with emphasis on method-specific descriptions, results, and conclusions. We draw selected studies on healthy and diseased hearts, aortas, and pulmonary arteries achieved through the incorporation of tissue mechanics, hemodynamics, and fluid-structure interaction methods paired with patient-specific data acquired with medical imaging in inverse modeling approaches.

Uncertainty quantification implementations in human hemodynamic flows

BACKGROUND AND OBJECTIVE

Human hemodynamic modeling is usually influenced by uncertainties occurring from a considerable unavailability of information linked to the boundary conditions and the physical properties used in the numerical models. Calculating the effect of these uncertainties on the numerical findings along the cardiovascular system is a demanding process due to the complexity of the morphology of the body and the area dynamics. To cope with all these difficulties, Uncertainty Quantification (UQ) methods seem to be an ideal tool.

RESULTS

This study focuses on analyzing and summarizing some of the recent research efforts and directions of implementing UQ in human hemodynamic flows by analyzing 139 research papers. Initially, the suitability of applying this approach is analyzed and demonstrated. Then, an overview of the most significant research work in various fields of biomedical hemodynamic engineering is presented. Finally, it is attempted to identify any possible forthcoming directions for research and methodological progress of UQ in biomedical sciences.

CONCLUSION

This review concludes that by finding the best statistical methods and parameters to represent the propagated uncertainties, while achieving a good interpretation of the interaction between input-output, is crucial for implementing UQ in biomedical sciences.

A surrogate model for plaque modeling in carotids based on Robin conditions calibrated by cine MRI data

We propose a surrogate model for the fluid-structure interaction (FSI) problem for the study of blood dynamics in carotid arteries in presence of plaque. This is based on the integration of a numerical model with subject-specific data and clinical imaging. We propose to model the plaque as part of the tissues surrounding the vessel wall through the application of an elastic support boundary condition. In order to characterize the plaque and other surrounding tissues, such as the close-by jugular vein, the elastic parameters of the boundary condition were spatially differentiated and their values were estimated by minimizing the discrepancies between computed vessel displacements and reference values obtained from CINE Magnetic Resonance Imaging data. We applied the model to three subjects with a degree of stenosis greater than 70%. We found that accounting for both plaque and jugular vein in the estimation of the elastic parameters increases the accuracy. In particular, in all patients, mismatches between computed and in vivo measured wall displacements were one to two orders of magnitude lower than the spatial resolution of the original MRI data. These results confirmed the validity of the proposed surrogate plaque model. We also compared fluid-dynamics results with those obtained in a fixed wall setting and in a full FSI model, used as gold standard, highlighting the better accordance of our results in comparison to the rigid ones.

Reducing the impact of geometric errors in flow computations using velocity measurements

Numerical blood flow simulations are typically set up from anatomical medical images and calibrated using velocity measurements. However, the accuracy of the computational geometry itself is limited by the resolution of the anatomical image. We first show that applying standard no-slip boundary conditions on inaccurately extracted boundaries can cause large errors in the results, in particular the pressure gradient. In this work, we therefore propose to augment the flow model calibration by slip/transpiration boundary conditions, whose parameters are then estimated using velocity measurements. Numerical experiments show that this methodology can considerably improve the accuracy of the estimated pressure gradients and 3D velocity fields when the vessel geometry is uncertain.

Mechanical Characterization of the Vessel Wall by Data Assimilation of Intravascular Ultrasound Studies

Atherosclerotic plaque rupture and erosion are the most important mechanisms underlying the sudden plaque growth, responsible for acute coronary syndromes and even fatal cardiac events. Advances in the understanding of the culprit plaque structure and composition are already reported in the literature, however, there is still much work to be done toward in-vivo plaque visualization and mechanical characterization to assess plaque stability, patient risk, diagnosis and treatment prognosis. In this work, a methodology for the mechanical characterization of the vessel wall plaque and tissues is proposed based on the combination of intravascular ultrasound (IVUS) imaging processing, data assimilation and continuum mechanics models within a high performance computing (HPC) environment. Initially, the IVUS study is gated to obtain volumes of image sequences corresponding to the vessel of interest at different cardiac phases. These sequences are registered against the sequence of the end-diastolic phase to remove transversal and longitudinal rigid motions prescribed by the moving environment due to the heartbeat. Then, optical flow between the image sequences is computed to obtain the displacement fields of the vessel (each associated to a certain pressure level). The obtained displacement fields are regarded as observations within a data assimilation paradigm, which aims to estimate the material parameters of the tissues within the vessel wall. Specifically, a reduced order unscented Kalman filter is employed, endowed with a forward operator which amounts to address the solution of a hyperelastic solid mechanics model in the finite strain regime taking into account the axially stretched state of the vessel, as well as the effect of internal and external forces acting on the arterial wall. Due to the computational burden, a HPC approach is mandatory. Hence, the data assimilation and computational solid mechanics computations are parallelized at three levels: (i) a Kalman filter level; (ii) a cardiac phase level; and (iii) a mesh partitioning level. To illustrate the capabilities of this novel methodology toward the in-vivo analysis of patient-specific vessel constituents, mechanical material parameters are estimated using in-silico and in-vivo data retrieved from IVUS studies. Limitations and potentials of this approach are exposed and discussed.

Relative pressure estimation from velocity measurements in blood flows: State-of-the-art and new approaches

The relative pressure difference across stenotic blood vessels serves as an important clinical index for the diagnosis of many cardiovascular diseases. While the clinical gold standard for relative pressure difference measurements is invasive catheterization, Phase-Contrast Magnetic Resonance Imaging has emerged as a promising tool for enabling a noninvasive quantification, by linking highly spatially resolved velocity measurements with relative pressures via the incompressible Navier-Stokes equations. In this work, we provide a review and analysis of current methods for relative pressure estimation and propose 3 additional techniques. Methods are compared using synthetic data from numerical examples, and sensitivity to subsampling and noise was explored. Through our analysis, we verify that the newly proposed approaches are more robust with respect to spatial subsampling and less sensitive to noise and therefore provide improved means for estimating relative pressure differences noninvasively.

Coupling schemes for the FSI forward prediction challenge: Comparative study and validation

This paper presents a numerical study in which several partitioned solution procedures for incompressible fluid-structure interaction are compared and validated against the results of an experimental fluid-structure interaction benchmark. The numerical methods discussed cover the three main families of coupling schemes: strongly coupled, semi-implicit, and loosely coupled. Very good agreement is observed between the numerical and experimental results. The comparisons confirm that strong coupling can be efficiently avoided, via semi-implicit and loosely coupled schemes, without compromising stability and accuracy. Copyright © 2016 John Wiley & Sons, Ltd.

The simulation of magnetic resonance elastography through atherosclerosis

The clinical diagnosis of atherosclerosis via the measurement of stenosis size is widely acknowledged as an imperfect criterion. The vulnerability of an atherosclerotic plaque to rupture is associated with its mechanical properties. The potential to image these mechanical properties using magnetic resonance elastography (MRE) was investigated through synthetic datasets. An image of the steady state wave propagation, equivalent to the first harmonic, can be extracted directly from finite element analysis. Inversion of this displacement data yields a map of the shear modulus, known as an elastogram. The variation of plaque composition, stenosis size, Gaussian noise, filter thresholds and excitation frequency were explored. A decreasing mean shear modulus with an increasing lipid composition was identified through all stenosis sizes. However the inversion algorithm showed sensitivity to parameter variation leading to artefacts which disrupted both the elastograms and quantitative trends. As noise was increased up to a realistic level, the contrast was maintained between the fully fibrous and lipid plaques but lost between the interim compositions. Although incorporating a Butterworth filter improved the performance of the algorithm, restrictive filter thresholds resulted in a reduction of the sensitivity of the algorithm to composition and noise variation. Increasing the excitation frequency improved the techniques ability to image the magnitude of the shear modulus and identify a contrast between compositions. In conclusion, whilst the technique has the potential to image the shear modulus of atherosclerotic plaques, future research will require the integration of a heterogeneous inversion algorithm.

Computational fluid dynamics modelling in cardiovascular medicine

This paper reviews the methods, benefits and challenges associated with the adoption and translation of computational fluid dynamics (CFD) modelling within cardiovascular medicine. CFD, a specialist area of mathematics and a branch of fluid mechanics, is used routinely in a diverse range of safety-critical engineering systems, which increasingly is being applied to the cardiovascular system. By facilitating rapid, economical, low-risk prototyping, CFD modelling has already revolutionised research and development of devices such as stents, valve prostheses, and ventricular assist devices. Combined with cardiovascular imaging, CFD simulation enables detailed characterisation of complex physiological pressure and flow fields and the computation of metrics which cannot be directly measured, for example, wall shear stress. CFD models are now being translated into clinical tools for physicians to use across the spectrum of coronary, valvular, congenital, myocardial and peripheral vascular diseases. CFD modelling is apposite for minimally-invasive patient assessment. Patient-specific (incorporating data unique to the individual) and multi-scale (combining models of different length- and time-scales) modelling enables individualised risk prediction and virtual treatment planning. This represents a significant departure from traditional dependence upon registry-based, population-averaged data. Model integration is progressively moving towards ‘digital patient’ or ‘virtual physiological human’ representations. When combined with population-scale numerical models, these models have the potential to reduce the cost, time and risk associated with clinical trials. The adoption of CFD modelling signals a new era in cardiovascular medicine. While potentially highly beneficial, a number of academic and commercial groups are addressing the associated methodological, regulatory, education- and service-related challenges.

A methodological paradigm for patient-specific multi-scale CFD simulations: from clinical measurements to parameter estimates for individual analysis

A new framework for estimation of lumped (for instance, Windkessel) model parameters from uncertain clinical measurements is presented. The ultimate aim is to perform patient-specific haemodynamic analysis. This framework is based on sensitivity analysis tools and the sequential estimation approach of the unscented Kalman filter. Sensitivity analysis and parameter estimation are performed in lumped parameter models, which act as reduced order surrogates of the 3D domain for haemodynamic analysis. While the goal of sensitivity analysis is to assess potential identifiability problems, the unscented Kalman filter estimation leads to parameter estimates based on clinical measurements and modelling assumptions. An application of such analysis and parameter estimation methodology is demonstrated for synthetic and real data. Equality constraints on various physiological parameters are enforced. Since the accuracy of the Windkessel parameter estimates depends on the lumped parameter representativeness, the latter is iteratively improved by running few 3D simulations while simultaneously improving the former. Such a method is applied on a patient-specific aortic coarctation case. Less than 3% and 9% errors between the clinically measured quantities and 3D simulation results for rest and stress are obtained, respectively. Knowledge on how these Windkessel parameters change from rest to stress can thus be learned by such an approach. Lastly, it is demonstrated that the proposed approach is capable of dealing with a wide variety of measurements and cases where the pressure and flow clinical measurements are not taken simultaneously.