TY - JOUR
T1 - Convergence analysis of pressure reconstruction methods from discrete velocities
AU - Araya, Rodolfo
AU - Bertoglio, Cristobal
AU - Carcamo, Cristian
AU - Nolte, David
AU - Uribe, Sergio
N1 - Publisher Copyright:
© The authors. Published by EDP Sciences, SMAI 2023.
PY - 2023/5/1
Y1 - 2023/5/1
N2 - Magnetic resonance imaging allows the measurement of the three-dimensional velocity field in blood flows. Therefore, several methods have been proposed to reconstruct the pressure field from such measurements using the incompressible Navier- Stokes equations, thereby avoiding the use of invasive technologies. However, those measurements are obtained at limited spatial resolution given by the voxel sizes in the image. In this paper, we propose a strategy for the convergence analysis of state-of-the-art pressure reconstruction methods. The methods analyzed are the so-called Pressure Poisson Estimator (PPE) and Stokes Estimator (STE). In both methods, the right-hand side corresponds to the terms that involving the field velocity in the Navier- Stokes equations, with a piecewise linear interpolation of the exact velocity. In the theoretical error analysis, we show that many terms of different order of convergence appear. These are certainly dominated by the lowest-order term, which in most cases stems from the interpolation of the velocity field. However, the numerical results in academic examples indicate that only the PPE may profit of increasing the polynomial order, and that the STE presents a higher accuracy than the PPE, but the interpolation order of the velocity field always prevails. Furthermore, we compare the pressure estimation methods on real MRI data, assessing the impact of different spatial resolutions and polynomial degrees on each method. Here, the results are consistent with the academic test cases in terms of sensitivity to polynomial order as well as the STE showing to be potentially more accurate when compared to reference pressure measurements.
AB - Magnetic resonance imaging allows the measurement of the three-dimensional velocity field in blood flows. Therefore, several methods have been proposed to reconstruct the pressure field from such measurements using the incompressible Navier- Stokes equations, thereby avoiding the use of invasive technologies. However, those measurements are obtained at limited spatial resolution given by the voxel sizes in the image. In this paper, we propose a strategy for the convergence analysis of state-of-the-art pressure reconstruction methods. The methods analyzed are the so-called Pressure Poisson Estimator (PPE) and Stokes Estimator (STE). In both methods, the right-hand side corresponds to the terms that involving the field velocity in the Navier- Stokes equations, with a piecewise linear interpolation of the exact velocity. In the theoretical error analysis, we show that many terms of different order of convergence appear. These are certainly dominated by the lowest-order term, which in most cases stems from the interpolation of the velocity field. However, the numerical results in academic examples indicate that only the PPE may profit of increasing the polynomial order, and that the STE presents a higher accuracy than the PPE, but the interpolation order of the velocity field always prevails. Furthermore, we compare the pressure estimation methods on real MRI data, assessing the impact of different spatial resolutions and polynomial degrees on each method. Here, the results are consistent with the academic test cases in terms of sensitivity to polynomial order as well as the STE showing to be potentially more accurate when compared to reference pressure measurements.
KW - 4D flow MRI
KW - Error analysis
KW - Navier- Stokes
KW - Poisson
KW - Stokes
UR - http://www.scopus.com/inward/record.url?scp=85163189140&partnerID=8YFLogxK
U2 - 10.1051/m2an/2023021
DO - 10.1051/m2an/2023021
M3 - Article
AN - SCOPUS:85163189140
SN - 2822-7840
VL - 57
SP - 1839
EP - 1861
JO - ESAIM: Mathematical Modelling and Numerical Analysis
JF - ESAIM: Mathematical Modelling and Numerical Analysis
IS - 3
ER -