Estimation of wall shear stress using 4D flow cardiovascular MRI and computational fluid dynamics

Electronic version of an article published as Journal of mechanics in medicine and biology, 0, 1750046 (2016), 16 pages. DOI:10.1142/S0219519417500464 © World Scientific Publishing CompanyIn the last few years, wall shear stress (WSS) has arisen as a new diagnostic indicator in patients with arterial disease. There is a substantial evidence that the WSS plays a significant role, together with hemodynamic indicators, in initiation and progression of the vascular diseases. Estimation of WSS values, therefore, may be of clinical significance and the methods employed for its measurement are crucial for clinical community. Recently, four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) has been widely used in a number of applications for visualization and quantification of blood flow, and although the sensitivity to blood flow measurement has increased, it is not yet able to provide an accurate three-dimensional (3D) WSS distribution. The aim of this work is to evaluate the aortic blood flow features and the associated WSS by the combination of 4D flow cardiovascular magnetic resonance (4D CMR) and computational fluid dynamics technique. In particular, in this work, we used the 4D CMR to obtain the spatial domain and the boundary conditions needed to estimate the WSS within the entire thoracic aorta using computational fluid dynamics. Similar WSS distributions were found for cases simulated. A sensitivity analysis was done to check the accuracy of the method. 4D CMR begins to be a reliable tool to estimate the WSS within the entire thoracic aorta using computational fluid dynamics. The combination of both techniques may provide the ideal tool to help tackle these and other problems related to wall shear estimation.Peer ReviewedPostprint (author's final draft

On the numerical treatment of viscous and convective effects in relative pressure reconstruction methods

The mechanism of many cardiovascular diseases can be understood by studying the pressure distribution in blood vessels. Direct pressure measurements, however, require invasive probing and provide only single‐point data. Alternatively, relative pressure fields can be reconstructed from imaging‐based velocity measurements by considering viscous and inertial forces. Both contributions can be potential troublemakers in pressure reconstruction: the former due to its higher‐order derivatives, and the latter because of the quadratic nonlinearity in the convective acceleration. Viscous and convective terms can be treated in various forms, which, although equivalent for ideal measurements, can perform differently in practice. In fact, multiple versions are often used in literature, with no apparent consensus on the more suitable variants. In this context, the present work investigates the performance of different versions of relative pressure estimators. For viscous effects, in particular, two new modified estimators are presented to circumvent second‐order differentiation without requiring surface integrals. In‐silico and in‐vitro data in the typical regime of cerebrovascular flows are considered, allowing a systematic noise sensitivity study. Convective terms are shown to be the main source of error, even for flows with pronounced viscous component. Moreover, the conservation (often integrated) form of convection exhibits higher noise sensitivity than the standard convective description, in all three families of estimators considered here. For the classical pressure Poisson estimator, the present modified version of the viscous term tends to yield better accuracy than the (recently introduced) integrated form, although this effect is in most cases negligible when compared to convection‐related errors

Validation of 4D Flow based relative pressure maps in aortic flows

While the clinical gold standard for pressure difference measurements is invasive catheterization, 4D Flow MRI is a promising tool for enabling a non-invasive quantification, by linking highly spatially resolved velocity measurements with pressure differences via the incompressible Navier–Stokes equations. In this work we provide a validation and comparison with phantom and clinical patient data of pressure difference maps estimators. We compare the classical Pressure Poisson Estimator (PPE) and the new Stokes Estimator (STE) against catheter pressure measurements under a variety of stenosis severities and flow intensities. Specifically, we use several 4D Flow data sets of realistic aortic phantoms with different anatomic and hemodynamic severities and two patients with aortic coarctation. The phantom data sets are enriched by subsampling to lower resolutions, modification of the segmentation and addition of synthetic noise, in order to study the sensitivity of the pressure difference estimators to these factors. Overall, the STE method yields more accurate results than the PPE method compared to catheterization data. The superiority of the STE becomes more evident at increasing Reynolds numbers with a better capacity of capturing pressure gradients in strongly convective flow regimes. The results indicate an improved robustness of the STE method with respect to variation in lumen segmentation. However, with heuristic removal of the wall-voxels, the PPE can reach a comparable accuracy for lower Reynolds’ numbers

Computational fluid dynamics simulations of blood flow regularized by 3D phase contrast MRI

Abstract Background Phase contrast magnetic resonance imaging (PC-MRI) is used clinically for quantitative assessment of cardiovascular flow and function, as it is capable of providing directly-measured 3D velocity maps. Alternatively, vascular flow can be estimated from model-based computation fluid dynamics (CFD) calculations. CFD provides arbitrarily high resolution, but its accuracy hinges on model assumptions, while velocity fields measured with PC-MRI generally do not satisfy the equations of fluid dynamics, provide limited resolution, and suffer from partial volume effects. The purpose of this study is to develop a proof-of-concept numerical procedure for constructing a simulated flow field that is influenced by both direct PC-MRI measurements and a fluid physics model, thereby taking advantage of both the accuracy of PC-MRI and the high spatial resolution of CFD. The use of the proposed approach in regularizing 3D flow fields is evaluated. Methods The proposed algorithm incorporates both a Newtonian fluid physics model and a linear PC-MRI signal model. The model equations are solved numerically using a modified CFD algorithm. The numerical solution corresponds to the optimal solution of a generalized Tikhonov regularization, which provides a flow field that satisfies the flow physics equations, while being close enough to the measured PC-MRI velocity profile. The feasibility of the proposed approach is demonstrated on data from the carotid bifurcation of one healthy volunteer, and also from a pulsatile carotid flow phantom. Results The proposed solver produces flow fields that are in better agreement with direct PC-MRI measurements than CFD alone, and converges faster, while closely satisfying the fluid dynamics equations. For the implementation that provided the best results, the signal-to-error ratio (with respect to the PC-MRI measurements) in the phantom experiment was 6.56 dB higher than that of conventional CFD; in the in vivo experiment, it was 2.15 dB higher. Conclusions The proposed approach allows partial or complete measurements to be incorporated into a modified CFD solver, for improving the accuracy of the resulting flow fields estimates. This can be used for reducing scan time, increasing the spatial resolution, and/or denoising the PC-MRI measurements.http://deepblue.lib.umich.edu/bitstream/2027.42/116061/1/12938_2015_Article_104.pd

Hemodynamic Measurement Using Four-Dimensional Phase-Contrast MRI: Quantification of Hemodynamic Parameters and Clinical Applications

Recent improvements have been made to the use of time-resolved, three-dimensional phase-contrast (PC) magnetic resonance imaging (MRI), which is also named four-dimensional (4D) PC-MRI or 4D flow MRI, in the investigation of spatial and temporal variations in hemodynamic features in cardiovascular blood flow. The present article reviews the principle and analytical procedures of 4D PC-MRI. Various fluid dynamic biomarkers for possible clinical usage are also described, including wall shear stress, turbulent kinetic energy, and relative pressure. Lastly, this article provides an overview of the clinical applications of 4D PC-MRI in various cardiovascular regions.113Ysciescopuskc

Vessel segmentation for time resolved phase contrast MRI

Quantication of cardiovascular flow and blood volumes are useful tools in diagnosing cardiovascular disease such as congenital heart defects and different kinds of valve leakage. Medical imaging techniques enables non-invasive analysis of anatomy and physiology. In order to perform flow quantification from medical image sequences, the boundaries of the vessels of interest are usually delineated manually by medical professionals. This is a time consuming process and the result depends, to a high degree, on user experience. This thesis presents an automated vessel segmentation method for the main vessels around the heart from velocity encoded Magnetic Resonance Imaging sequences. The proposed method only require a manual delineation in one image. The algorithm is based on an active contour, using the Euler-Lagrange equation together with internal and external forces designed from a set of fundamental assumptions regarding vessel shape and behaviour. More specically, constraints were applied to the geometrical shape and elasticity of a vessel. Validation of the method was performed by comparing the detected stroke volume with manual delineation, and also by measuring the segmentation overlapping of the two methods. On a test set of 20 patients, 19 resulted in excellent segmentation agreement with manual delineations, with a mean Dice coefficient over 0.8. However, performance instability was observed when changing the values of two algorithm parameters, and one of the patients in the set constantly resulted in segmentation failure for all tested parameter combinations. The relative variability in stroke volume between the proposed algorithm and manual delineations was, at best, 6 +- 3.6%. This is comparable to the interobserver variability from a previous physiological study 1 of 3 +- 4%, which indicates the potential of the suggested method if improvements in robustness and stability is implemented

Modelling the Human Cardiac Fluid Mechanics. 4th ed

With the Karlsruhe Heart Model (KaHMo) we aim to share our vision of integrated computational simulation across multiple disciplines of cardiovascular research, and emphasis yet again the importance of Modelling the Human Cardiac Fluid Mechanics within the framework of the international STICH study. The focus of this work is on integrated cardiovascular fluid mechanics, and the potential benefits to future cardiovascular research and the wider bio-medical community

Estimation of wall shear stress using 4D flow cardiovascular MRI and computational fluid dynamics

Company In the last few years, wall shear stress (WSS) has arisen as a new diagnostic indicator in patients with arterial disease. There is a substantial evidence that the WSS plays a significant role, together with hemodynamic indicators, in initiation and progression of the vascular diseases. Estimation of WSS values, therefore, may be of clinical significance and the methods employed for its measurement are crucial for clinical community. Recently, four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) has been widely used in a number of applications for visualization and quantification of blood flow, and although the sensitivity to blood flow measurement has increased, it is not yet able to provide an accurate three-dimensional (3D) WSS distribution. The aim of this work is to evaluate the aortic blood flow features and the associated WSS by the combination of 4D flow cardiovascular magnetic resonance (4D CMR) and computational fluid dynamics technique. In particular, in this work, we used the 4D CMR to obtain the spatial domain and the boundary conditions needed to estimate the WSS within the entire thoracic aorta using computational fluid dynamics. Similar WSS distributions were found for cases simulated. A sensitivity analysis was done to check the accuracy of the method. 4D CMR begins to be a reliable tool to estimate the WSS within the entire thoracic aorta using computational fluid dynamics. The combination of both techniques may provide the ideal tool to help tackle these and other problems related to wall shear estimatio