Data Assimilation in Reduced Modeling
Abstract
We consider the problem of optimal recovery of an element of a Hilbert space from measurements obtained through known linear functionals on . Problems of this type are well studied \cite{MRW} under an assumption that belongs to a prescribed model class, e.g. a known compact subset of . Motivated by reduced modeling for parametric partial differential equations, this paper considers another setting where the additional information about is in the form of how well can be approximated by a certain known subspace of of dimension , or more generally, how well can be approximated by each -dimensional subspace of a sequence of nested subspaces . A recovery algorithm for the one-space formulation, proposed in \cite{MPPY}, is proven here to be optimal and to have a simple formulation, if certain favorable bases are chosen to represent and the measurements. The major contribution of the present paper is to analyze the multi-space case for which it is shown that the set of all satisfying the given information can be described as the intersection of a family of known ellipsoids in . It follows that a near optimal recovery algorithm in the multi-space problem is to identify any point in this intersection which can provide a much better accuracy than in the one-space problem. Two iterative algorithms based on alternating projections are proposed for recovery in the multi-space problem. A detailed analysis of one of them provides a posteriori performance estimates for the iterates, stopping criteria, and convergence rates. Since the limit of the algorithm is a point in the intersection of the aforementioned ellipsoids, it provides a near optimal recovery for .
Cite
@article{arxiv.1506.04770,
title = {Data Assimilation in Reduced Modeling},
author = {Peter Binev and Albert Cohen and Wolfgang Dahmen and Ronald DeVore and Guergana Petrova and Przemyslaw Wojtaszczyk},
journal= {arXiv preprint arXiv:1506.04770},
year = {2015}
}
Comments
27 pages