Multilevel ensemble Kalman filtering for spatio-temporal processes

Hakon Andreas Hoel, Fabio Nobile
2021
Journal paper

We design and analyse the performance of a multilevel ensemble Kalman filter method (MLEnKF) for filtering settings where the underlying state-space model is an infinite-dimensional spatio-temporal process. We consider underlying models that needs to be simulated by numerical methods, with discretization in both space and time. The multilevel Monte Carlo sampling strategy, achieving variance reduction through pairwise coupling of ensemble particles on neighboring resolutions, is used in the sample-moment step of MLEnKF to produce an efficent hierarchical filtering method for spatio-temporal models. Under sufficent regularity, MLEnKF is proven to be more efficient for weak approximations than EnKF, asymptotically in the large-ensemble and fine-numerical-resolution limit. Numerical examples support our theoretical findings.

About this result

This page is automatically generated and may contain information that is not correct, complete, up-to-date, or relevant to your search query. The same applies to every other page on this website. Please make sure to verify the information with EPFL's official sources.

Related concepts

Related publications

Hakon Andreas Hoel, Fabio Nobile
2021
Journal paper

Abstract

Official source

https://infoscience.epfl.ch/record/283159?ln=en

About this result

Related concepts (7)

Related concepts

Related publications (3)

Related publications

Finite element heterogeneous multiscale method for the wave equation

Assyr Abdulle

A finite element heterogeneous multiscale method is proposed for the wave equation with highly oscillatory coefficients. It is based on a finite element discretization of an effective wave equation at the macro scale, whose a priori unknown effective coefficients are computed on sampling domains at the micro scale within each macro finite element. Hence the computational work involved is independent of the highly heterogeneous nature of the medium at the smallest scale. Optimal error estimates in the energy norm and the L-2 norm and convergence to the homogenized solution are proved, when both the macro and the micro scales are refined simultaneously. Numerical experiments corroborate the theoretical convergence rates and illustrate the behavior of the numerical method for periodic and heterogeneous media.

2011

Fully discrete analysis of the heterogeneous multiscale method for elliptic problems with multiple scales

Assyr Abdulle, Yun Bai

A fully discrete analysis of the finite element heterogeneous multiscale method (FE-HMM) for elliptic problems with N+1 well-separated scales is discussed. The FE-HMM is a numerical homogenization method that relies on a macroscopic scheme (macro FEM) for the approximation of the effective solution corresponding to the multiscale problem. The effective data are recovered from micro scale computation (micro FEM) on sampling domains located at appropriate quadrature points of the macroscopic mesh. At the macroscopic level, the numerical method can be seen as an FEM with numerical quadrature for a modified effective problem, hence variational crimes are made when designing this method. Up to now, the method has been analysed for two scales and the micro FEM was assumed to be conforming. For more than two scales, variational crimes are committed also at the intermediate (meso) scales and the effective data of the macroscopic scheme are obtained from a cascade of FEMs with numerical integration, which require a careful analysis. Numerical experiments for three-scale problems illustrate the theoretical convergence rates.

Oxford Univ Press2015

Reduced basis finite element heterogeneous multiscale method for high-order discretizations of elliptic homogenization problems

Assyr Abdulle, Yun Bai

A new finite element method for the efficient discretization of elliptic homogenization problems is proposed. These problems, characterized by data varying over a wide range of scales cannot be easily solved by classical numerical methods that need mesh resolution down to the finest scales and multiscale methods capable of capturing the large scale components of the solution on macroscopic meshes are needed. Recently, the finite element heterogeneous multiscale method (FE-HMM) has been proposed for such problems, based on a macroscopic solver with effective data recovered from the solution of micro problems on sampling domains at quadrature points of a macroscopic mesh. Departing from the approach used in the FE-HMM, we show that interpolation techniques based on the reduced basis methodology (an offline-online strategy) allow one to design an efficient numerical method relying only on a small number of accurately computed micro solutions. This new method, called the reduced basis finite element heterogeneous multiscale method (RB-FE-HMM) is significantly more efficient than the FE-HMM for high order macroscopic discretizations and for three-dimensional problems, when the repeated computation of micro problems over the whole computational domain is expensive. A priori error estimates of the RB-FE-HMM are derived. Numerical computations for two and three dimensional problems illustrate the applicability and efficiency of the numerical method. (C) 2012 Elsevier Inc. All rights reserved.

Elsevier2012

Numerical analysis

Numerical analysis is the study of algorithms that use numerical approximation (as opposed to symbolic manipulations) for the problems of mathematical analysis (as distinguished from discrete mathema

Space

Space is a three-dimensional continuum containing positions and directions. In classical physics, physical space is often conceived in three linear dimensions. Modern physicists usually consider it,

Monte Carlo method

Monte Carlo methods, or Monte Carlo experiments, are a broad class of computational algorithms that rely on repeated random sampling to obtain numerical results. The underlying concept is to use rando