A Reduced Basis Ensemble Kalman Filter for State/parameter Identification in Large-scale Nonlinear Dynamical Systems

The ensemble Kalman filter is nowadays widely employed to solve state and/or parameter identification problems recast in the framework of Bayesian inversion. Unfortunately its cost becomes prohibitive when dealing with systems described by parametrized partial differential equations, because of the cost entailed by each PDE query. This is even worse for nonlinear time-dependent PDEs. In this paper we propose a reduced basis ensemble Kalman filter technique to speed up the numerical solution of Bayesian inverse problems arising from the discretization of nonlinear time dependent PDEs. The reduction stage yields intrinsic approximation errors, whose propagation through the filtering process might affect the accuracy of the identified state/parameters. Since their evaluation is computationally heavy, we equip our reduced basis ensemble Kalman filter with a reduction error model based on ordinary kriging for functional-valued data, to gauge the effect of state reduction on the whole filtering process. The accuracy and efficiency of our method is then verified on two numerical test cases, dealing with the identification of uncertain parameters or fields for a FitzHugh-Nagumo model and a Fisher-Kolmogorov model.

[1]  N. Nguyen,et al.  A general multipurpose interpolation procedure: the magic points , 2008 .

[2]  D. S. McCormick,et al.  Accuracy and stability of filters for dissipative PDEs , 2012, 1203.5845.

[3]  Faming Liang,et al.  Statistical and Computational Inverse Problems , 2006, Technometrics.

[4]  Andrew M. Stuart,et al.  Inverse problems: A Bayesian perspective , 2010, Acta Numerica.

[5]  Danny C. Sorensen,et al.  Nonlinear Model Reduction via Discrete Empirical Interpolation , 2010, SIAM J. Sci. Comput..

[6]  P. Houtekamer,et al.  A Sequential Ensemble Kalman Filter for Atmospheric Data Assimilation , 2001 .

[7]  Oliver G. Ernst,et al.  Bayesian Inverse Problems and Kalman Filters , 2014 .

[8]  Oliver G. Ernst,et al.  Analysis of the Ensemble and Polynomial Chaos Kalman Filters in Bayesian Inverse Problems , 2015, SIAM/ASA J. Uncertain. Quantification.

[9]  T. Başar,et al.  A New Approach to Linear Filtering and Prediction Problems , 2001 .

[10]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[11]  Mario Ohlberger,et al.  Data-driven combined state and parameter reduction for inverse problems , 2014, Advances in Computational Mathematics.

[12]  Dennis McLaughlin,et al.  Efficient Characterization of Uncertain Model Parameters with a Reduced-Order Ensemble Kalman Filter , 2014, SIAM J. Sci. Comput..

[13]  Arnaud Trouvé,et al.  Towards predictive data-driven simulations of wildfire spread - Part II: Ensemble Kalman Filter for the state estimation of a front-tracking simulator of wildfire spread , 2014 .

[14]  R. FitzHugh Impulses and Physiological States in Theoretical Models of Nerve Membrane. , 1961, Biophysical journal.

[15]  Adrian Sandu,et al.  Efficient approximation of Sparse Jacobians for time‐implicit reduced order models , 2014, ArXiv.

[16]  Tiangang Cui,et al.  On dimension reduction in Gaussian filters , 2015, 1508.06452.

[17]  A. Quarteroni,et al.  Reduced Basis Methods for Partial Differential Equations: An Introduction , 2015 .

[18]  Christian Vergara,et al.  Computational Comparison Between Newtonian and Non-Newtonian Blood Rheologies in Stenotic Vessels , 2018 .

[19]  V. M. Tikhomirov,et al.  A Study of the Diffusion Equation with Increase in the Amount of Substance, and its Application to a Biological Problem , 1991 .

[20]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

[21]  J. Kaipio,et al.  Approximation error analysis in nonlinear state estimation with an application to state-space identification , 2007 .

[22]  Pedro Díez,et al.  Real time parameter identification and solution reconstruction from experimental data using the Proper Generalized Decomposition , 2015 .

[23]  Jari P. Kaipio,et al.  Model reduction in state identification problems with an application to determination of thermal parameters , 2009 .

[24]  Karen Willcox,et al.  A Survey of Projection-Based Model Reduction Methods for Parametric Dynamical Systems , 2015, SIAM Rev..

[25]  James D. Murray Mathematical Biology: I. An Introduction , 2007 .

[26]  A. Stuart,et al.  Data Assimilation: A Mathematical Introduction , 2015, 1506.07825.

[27]  Geir Evensen,et al.  Advanced Data Assimilation for Strongly Nonlinear Dynamics , 1997 .

[28]  David Amsallem,et al.  Efficient model reduction of parametrized systems by matrix discrete empirical interpolation , 2015, J. Comput. Phys..

[29]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[30]  A. Stuart,et al.  Ensemble Kalman methods for inverse problems , 2012, 1209.2736.

[31]  Alessandra Menafoglio,et al.  A Universal Kriging predictor for spatially dependent functional data of a Hilbert Space , 2013 .

[32]  J. Mateu,et al.  Ordinary kriging for function-valued spatial data , 2011, Environmental and Ecological Statistics.

[33]  Y. Marzouk,et al.  Large-Scale Inverse Problems and Quantification of Uncertainty , 1994 .

[34]  David Galbally,et al.  Non‐linear model reduction for uncertainty quantification in large‐scale inverse problems , 2009 .

[35]  Andrea Manzoni,et al.  Accurate Solution of Bayesian Inverse Uncertainty Quantification Problems Combining Reduced Basis Methods and Reduction Error Models , 2016, SIAM/ASA J. Uncertain. Quantification.

[36]  Hermann G. Matthies,et al.  A deterministic filter for non-Gaussian Bayesian estimation— Applications to dynamical system estimation with noisy measurements , 2012 .

[37]  A. Stuart,et al.  Well-posedness and accuracy of the ensemble Kalman filter in discrete and continuous time , 2013, 1310.3167.

[38]  B. Haasdonk,et al.  A reduced basis Kalman filter for parametrized partial differential equations , 2016 .

[39]  Tiangang Cui,et al.  Data‐driven model reduction for the Bayesian solution of inverse problems , 2014, 1403.4290.