Data Assimilation in Reduced Modeling

Abstract

We consider the problem of optimal recovery of an element u of a Hilbert space H from m measurements obtained through known linear functionals on H. Problems of this type are well studied MRW under an assumption that u belongs to a prescribed model class, e.g. a known compact subset of H. Motivated by reduced modeling for parametric partial differential equations, this paper considers another setting where the additional information about u is in the form of how well u can be approximated by a certain known subspace Vn of H of dimension n, or more generally, how well u can be approximated by each k-dimensional subspace Vk of a sequence of nested subspaces V0⊂ V1·s⊂ Vn. A recovery algorithm for the one-space formulation, proposed in MPPY, is proven here to be optimal and to have a simple formulation, if certain favorable bases are chosen to represent Vn 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 u satisfying the given information can be described as the intersection of a family of known ellipsoids in H. 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 u.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…