The reliable and effective assimilation of measurements and numerical simulations in engineering applications involving computational fluid dynamics is an emerging problem as soon as new devices provide more data. In this paper we are mainly driven by hemodynamics applications, a field where the progressive increment of measures and numerical tools makes this problem particularly up-to-date. We adopt a Bayesian approach to the inclusion of noisy data in the incompressible steady Navier-Stokes equations (NSE). The purpose is the quantification of uncertainty affecting velocity and flow related variables of interest, all treated as random variables. The method consists in the solution of an optimization problem where the misfit between data and velocity - in a convenient norm - is minimized under the constraint of the NSE. We derive classical point estimators, namely the maximum a posteriori - MAP - and the maximum likelihood - ML - ones. In addition, we obtain confidence regions for velocity and wall shear stress, a flow related variable of medical relevance. Numerical simulations in 2-dimensional and axisymmetric 3-dimensional domains show the gain yielded by the introduction of a complete statistical knowledge in the assimilation process.
Mots-clés : computational fluid dynamics, optimization, uncertainty quantification, statistical inverse problems, data assimilation, hemodynamics
@article{M2AN_2013__47_4_1037_0, author = {D{\textquoteright}Elia, Marta and Veneziani, Alessandro}, title = {Uncertainty quantification for data assimilation in a steady incompressible {Navier-Stokes} problem}, journal = {ESAIM: Mathematical Modelling and Numerical Analysis }, pages = {1037--1057}, publisher = {EDP-Sciences}, volume = {47}, number = {4}, year = {2013}, doi = {10.1051/m2an/2012056}, mrnumber = {3082288}, zbl = {1271.76062}, language = {en}, url = {http://www.numdam.org/articles/10.1051/m2an/2012056/} }
TY - JOUR AU - D’Elia, Marta AU - Veneziani, Alessandro TI - Uncertainty quantification for data assimilation in a steady incompressible Navier-Stokes problem JO - ESAIM: Mathematical Modelling and Numerical Analysis PY - 2013 SP - 1037 EP - 1057 VL - 47 IS - 4 PB - EDP-Sciences UR - http://www.numdam.org/articles/10.1051/m2an/2012056/ DO - 10.1051/m2an/2012056 LA - en ID - M2AN_2013__47_4_1037_0 ER -
%0 Journal Article %A D’Elia, Marta %A Veneziani, Alessandro %T Uncertainty quantification for data assimilation in a steady incompressible Navier-Stokes problem %J ESAIM: Mathematical Modelling and Numerical Analysis %D 2013 %P 1037-1057 %V 47 %N 4 %I EDP-Sciences %U http://www.numdam.org/articles/10.1051/m2an/2012056/ %R 10.1051/m2an/2012056 %G en %F M2AN_2013__47_4_1037_0
D’Elia, Marta; Veneziani, Alessandro. Uncertainty quantification for data assimilation in a steady incompressible Navier-Stokes problem. ESAIM: Mathematical Modelling and Numerical Analysis , Direct and inverse modeling of the cardiovascular and respiratory systems. Numéro spécial, Volume 47 (2013) no. 4, pp. 1037-1057. doi : 10.1051/m2an/2012056. http://www.numdam.org/articles/10.1051/m2an/2012056/
[1] Sequential parameter estimation for fluid-structure problems. Application to hemodynamics. Inter. J. Numer. Methods Biomed. Eng. 28 (2012) 434-455. RR-7657. | MR
, and ,[2] Data Assimilation for Geophysical Fluids, Handbook of numerical analysis, vol. XIV, chapter 9. Elsevier (2005).
, and ,[3] Introduction to the Theory of Statistics. McGraw-Hill (1974). | Zbl
, and ,[4] Subjective knowledge or objective belief? an oblique look to bayesian methods, in Large-Scale Inverse Problems and Quantification of Uncertainty, edited by G. Biros et al. Wiley Online Library (2011) 33-70. | MR
and ,[5] Ph.D. thesis.
,[6] Applications of Variational Data Assimilation in Computational Hemodynamics, chapter 12. MS & A. Springer (2011) 363-394.
, , , , , and ,[7] A variational Data Assimilation procedure for the incompressible Navier-Stokes equations in hemodynamics. Technical Report TR-2010-19, Department of Mathematics and Computer Science, Emory University, To appear in J. Sci. Comput. Available on www.mathcs.emory.edu (2010). | Zbl
, and ,[8] C. Lamphier, R. Beekman and M. Brummer, Hemodynamic predictors of aortic dilatation in bicuspid aortic valve by velocity-encoded cardiovascular magnetic resonance. J. Cardiovasc. Magn. Reson. 12 (2010) 4.
, , , , , , ,[9] Gmresr: a family of nested gmres methods. Numer. Linear Algebra Appl. 1 (1994) 369-386. | MR | Zbl
and ,[10] Bayesian inference for data assimilation using Least-Squares Finite Element methods, in IOP Conf. Ser. Mat. Sci. Eng., vol. 10. IOP Publishing (2010) 012224.
,[11] SIAM J. Sci. Comput. (2008).
, and .[12] Comput. Methods Appl. Mech. Eng. (2010). | MR
, and .[13] Surrogate and reduced-order modeling: A comparison of approaches for large-scale statistical inverse problems. Large-Scale Inverse Problems and Quantification of Uncertainty (2010) 123-149. | MR
, , and ,[14] Perspectives in flow control and optimization. Society for Industrial Mathematics 5 (2003). | MR | Zbl
,[15] Rank-deficient and discrete ill-posed problems. SIAM Monographs on Mathematical Modeling and Computation. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA (1998). | MR | Zbl
,[16] Weighted least-squares finite elements based on particle imaging velocimetry data. J. Comput. Phys. 229 (2010) 107-118. | MR
, , , , and ,[17] Artificial boundaries and flux pressure conditions for the incompressible navier-stokes equations. Int. J. Numer. Methods Fluids 22 (1996) 325-352. | MR | Zbl
, and ,[18] Applied multivariate statistical analysis. Prentice-Hall, Inc., Upper Saddle River, NJ, USA (1988). | Zbl
and ,[19] Statistical and Computational Inverse Problems. Springer (2005). | MR | Zbl
and ,[20] A new approach to linear filtering and prediction problems. Trans. ASME-J. Basic Eng. 82 (1960) 35-45.
,[21] A preconditioner for the steady-state navier-stokes equations. SIAM J. Sci. Comput. 24 (2002) 237-256. | MR | Zbl
, and ,[22] Sequential identification of boundary support parameters in a fluid-structure vascular model using patient image data. Biomechanics and Modeling in Mechanobiology. Published Online (2012) 1-22.
, , , , , and ,[23] Reduced-order unscented kalman filtering with application to parameter identification in large-dimensional systems. ESAIM: COCV 17 (2011) 380-405. | Numdam | MR | Zbl
and ,[24] Numerical Optimization. Springer (2000). | MR | Zbl
and ,[25] A variational approach for estimating the compliance of the cardiovascular tissue: An inverse fluid-structure interaction problem. SIAM J. Sci. Comput. 33 (2011) 1181-1211. | MR | Zbl
, and ,[26] Certified reduced basis approximation for parametrized partial differential equations and applications. J. Math. Ind. 1 (2011) 3. | MR | Zbl
, and ,[27] Efficient preconditioning of the linearized navier-stokes equations for incompressible flow. J. Comput. Appl. Math. 128 (2001) 261-279. | MR | Zbl
, , and ,[28] Inverse problem theory and methods for model parameter estimation. Society for Industrial Mathematics (2005). | MR | Zbl
,[29] Boundary conditions for blood flow problems, in Proc. of ENUMATH97, edited by R. Rannacher et al., World Sci. Publishing (1998). | MR | Zbl
,[30] Mathematical and Numerical Modeling of Blood flow Problems. Ph.D. thesis, Politecnico di Milano, Italy (1998).
,[31] New insights in gmres-like methods with variable preconditioners. J. Comput. Appl. Math. 61 (1995) 189-204. | MR | Zbl
,Cited by Sources: