- Poster presentation
- Open Access
Quantification of wall shear stress using finite-element interpolations in multidimensional phase contrast MR data of the thoracic aorta
© Sotelo et al.; licensee BioMed Central Ltd. 2014
- Published: 16 January 2014
- Wall Shear Stress
- Wall Shear Stress Distribution
- Average Wall Shear Stress
- Shear Stress Tensor
- Local Wall Shear Stress
Different methods have been proposed to estimate Wall Shear Stress (WSS). Morgan and co-workers (J Thorac Cardiovasc Surg and Ann Biomed Eng 1998) used a finite difference scheme to quantify the WSS tensor. However, it is well known that finite-difference methods cannot effectively handle complex geometries, as those found in the cardiovascular system. To account for arbitrary cross-section shapes, Stalder et al., MRM 2008, used B-spline (BS) interpolations to smoothly describe the lumen contours. In this work, we propose and validate a new method for calculating the WSS distribution based on Finite-Element (FE) interpolations.
The velocity field obtained at discrete locations from 2D and 3D CINE PC-MRI was interpolated using linear triangular FE. The shear stress tensor over the entire section cut was obtained from a global least-squares stress-projection method, from which the axial WSS vector was obtained. The proposed approach was benchmarked against a modified Poiseuille flow profile, and the robustness of the method was assessed by changing the level of resolution and noise. Additionally, we computed the WSS distribution in different aortic sections from a pulsatile aortic phantom, and from 5 healthy volunteers. In the aortic phantom and volunteers, 3D CINE PC-MRI flow data was acquired in 2D cutting planes in 5 with a spatial resolution 0.8 mm2 and temporal resolution 37 ms. We have also compared our framework with a BS based method previously reported in the literature.
In conclusion, we have developed a novel methodology to calculate WSS based on FE interpolations, which provides an excellent approximation of local WSS values, stability when subjected to noise and remarkable convergence properties as the pixel size is decreased.
VRI # 44/2011 (Pontificia Universidad Católica de Chile), Anillo ACT 079 and FONDECYT #11100427 and #11121224. JS thanks CONICYT for scholarship for doctoral studies.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.