Adaptive Galerkin approximation algorithms for Kolmogorov equations in infinite dimensions

Space-time variational formulations and adaptive Wiener–Hermite polynomial chaos Galerkin discretizations of Kolmogorov equations in infinite dimensions, such as Fokker–Planck and Ornstein–Uhlenbeck equations for functions defined on an infinite-dimensional separable Hilbert space H, are developed....

Full description

Bibliographic Details
Main Authors: Schwab, C, Süli, E
Format: Journal article
Language:English
Published: Springer-Verlag 2013
_version_ 1797073293164937216
author Schwab, C
Süli, E
author_facet Schwab, C
Süli, E
author_sort Schwab, C
collection OXFORD
description Space-time variational formulations and adaptive Wiener–Hermite polynomial chaos Galerkin discretizations of Kolmogorov equations in infinite dimensions, such as Fokker–Planck and Ornstein–Uhlenbeck equations for functions defined on an infinite-dimensional separable Hilbert space H, are developed. The wellposedness of these equations in the Hilbert space L2(H, μ) of functions on the infinite-dimensional domain H, which are square-integrable with respect to a Gaussian measure μ with trace class covariance operator Q on H, is proved. Specifically, for the infinite-dimensional Fokker–Planck equation, adaptive space-time Galerkin discretizations, based on a wavelet polynomial chaos Riesz basis obtained by tensorization of biorthogonal piecewise polynomial wavelet bases in time with a spatial Wiener–Hermite polynomial chaos arising from the Wiener–Itô decomposition of L2(H, μ), are introduced. The resulting space-time adaptive Wiener–Hermite polynomial Galerkin discretization algorithms of the infinite-dimensional PDE are proved to converge quasioptimally in the sense that they produce sequences of finite-dimensional approximations that attain the best possible convergence rates afforded by best N-term approximations of the solution from tensor-products of multiresolution (wavelet) time-discretizations and theWiener–Hermite polynomial chaos in L2(H, μ). As a consequence, the proposed adaptive Galerkin solution algorithms exhibit dimension-independent performance, which is optimal with respect to the algebraic best N-term rate afforded by the solution and the polynomial degree and regularity of the multiresolution (wavelet) time-discretizations in the finite-dimensional case, in particular. All constants in our error and complexity bounds are shown to be independent of the number of “active” coordinates identified by the proposed adaptive Galerkin approximation algorithms. The computational work and memory required by the proposed algorithms scale linearly with the support size of the coefficient vectors that arise in the approximations, with dimension-independent constants.
first_indexed 2024-03-06T23:19:59Z
format Journal article
id oxford-uuid:6869e853-cef3-44bf-8cad-4db3a29b0fe0
institution University of Oxford
language English
last_indexed 2024-03-06T23:19:59Z
publishDate 2013
publisher Springer-Verlag
record_format dspace
spelling oxford-uuid:6869e853-cef3-44bf-8cad-4db3a29b0fe02022-03-26T18:44:37ZAdaptive Galerkin approximation algorithms for Kolmogorov equations in infinite dimensionsJournal articlehttp://purl.org/coar/resource_type/c_dcae04bcuuid:6869e853-cef3-44bf-8cad-4db3a29b0fe0EnglishSymplectic Elements at OxfordSpringer-Verlag2013Schwab, CSüli, ESpace-time variational formulations and adaptive Wiener–Hermite polynomial chaos Galerkin discretizations of Kolmogorov equations in infinite dimensions, such as Fokker–Planck and Ornstein–Uhlenbeck equations for functions defined on an infinite-dimensional separable Hilbert space H, are developed. The wellposedness of these equations in the Hilbert space L2(H, μ) of functions on the infinite-dimensional domain H, which are square-integrable with respect to a Gaussian measure μ with trace class covariance operator Q on H, is proved. Specifically, for the infinite-dimensional Fokker–Planck equation, adaptive space-time Galerkin discretizations, based on a wavelet polynomial chaos Riesz basis obtained by tensorization of biorthogonal piecewise polynomial wavelet bases in time with a spatial Wiener–Hermite polynomial chaos arising from the Wiener–Itô decomposition of L2(H, μ), are introduced. The resulting space-time adaptive Wiener–Hermite polynomial Galerkin discretization algorithms of the infinite-dimensional PDE are proved to converge quasioptimally in the sense that they produce sequences of finite-dimensional approximations that attain the best possible convergence rates afforded by best N-term approximations of the solution from tensor-products of multiresolution (wavelet) time-discretizations and theWiener–Hermite polynomial chaos in L2(H, μ). As a consequence, the proposed adaptive Galerkin solution algorithms exhibit dimension-independent performance, which is optimal with respect to the algebraic best N-term rate afforded by the solution and the polynomial degree and regularity of the multiresolution (wavelet) time-discretizations in the finite-dimensional case, in particular. All constants in our error and complexity bounds are shown to be independent of the number of “active” coordinates identified by the proposed adaptive Galerkin approximation algorithms. The computational work and memory required by the proposed algorithms scale linearly with the support size of the coefficient vectors that arise in the approximations, with dimension-independent constants.
spellingShingle Schwab, C
Süli, E
Adaptive Galerkin approximation algorithms for Kolmogorov equations in infinite dimensions
title Adaptive Galerkin approximation algorithms for Kolmogorov equations in infinite dimensions
title_full Adaptive Galerkin approximation algorithms for Kolmogorov equations in infinite dimensions
title_fullStr Adaptive Galerkin approximation algorithms for Kolmogorov equations in infinite dimensions
title_full_unstemmed Adaptive Galerkin approximation algorithms for Kolmogorov equations in infinite dimensions
title_short Adaptive Galerkin approximation algorithms for Kolmogorov equations in infinite dimensions
title_sort adaptive galerkin approximation algorithms for kolmogorov equations in infinite dimensions
work_keys_str_mv AT schwabc adaptivegalerkinapproximationalgorithmsforkolmogorovequationsininfinitedimensions
AT sulie adaptivegalerkinapproximationalgorithmsforkolmogorovequationsininfinitedimensions