Wavelet compressed, modified Hilbert transform in the space-time discretization of the heat equation

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Harbrecht, Helmut, Schwab, Christoph, Zank, Marco
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866916776307064832
author Harbrecht, Helmut
Schwab, Christoph
Zank, Marco
author_facet Harbrecht, Helmut
Schwab, Christoph
Zank, Marco
contents On a finite time interval $(0,T)$, we consider the multiresolution Galerkin discretization of a modified Hilbert transform $\mathcal H_T$ which arises in the space-time Galerkin discretization of the linear diffusion equation. To this end, we design spline-wavelet systems in $(0,T)$ consisting of piecewise polynomials of degree $\geq 1$ with sufficiently many vanishing moments which constitute Riesz bases in the Sobolev spaces $H^{s}_{0,}(0,T)$ and $H^{s}_{,0}(0,T)$. These bases provide stable multilevel splittings of the temporal discretization spaces into "increment" or "detail" spaces. Furthermore, they allow to optimally compress the nonlocal integrodifferential operators which appear in stable space-time variational formulations of initial-boundary value problems, such as the heat equation and the acoustic wave equation. We then obtain sparse space-time tensor-product spaces via algebraic tensor-products of the temporal multilevel discretizations with standard, hierarchic finite element spaces in the spatial domain (with standard Lagrangian FE bases). Hence, the construction of multiresolutions in the spatial domain is not necessary. An efficient multilevel preconditioner is proposed that solves the linear system of equations resulting from the sparse space-time Galerkin discretization with essentially linear complexity (in work and memory). A substantial reduction in the number of the degrees of freedom and CPU time (compared to time-marching discretizations) is demonstrated in numerical experiments.
format Preprint
id arxiv_https___arxiv_org_abs_2402_10346
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Wavelet compressed, modified Hilbert transform in the space-time discretization of the heat equation
Harbrecht, Helmut
Schwab, Christoph
Zank, Marco
Numerical Analysis
65N30, 65J15
On a finite time interval $(0,T)$, we consider the multiresolution Galerkin discretization of a modified Hilbert transform $\mathcal H_T$ which arises in the space-time Galerkin discretization of the linear diffusion equation. To this end, we design spline-wavelet systems in $(0,T)$ consisting of piecewise polynomials of degree $\geq 1$ with sufficiently many vanishing moments which constitute Riesz bases in the Sobolev spaces $H^{s}_{0,}(0,T)$ and $H^{s}_{,0}(0,T)$. These bases provide stable multilevel splittings of the temporal discretization spaces into "increment" or "detail" spaces. Furthermore, they allow to optimally compress the nonlocal integrodifferential operators which appear in stable space-time variational formulations of initial-boundary value problems, such as the heat equation and the acoustic wave equation. We then obtain sparse space-time tensor-product spaces via algebraic tensor-products of the temporal multilevel discretizations with standard, hierarchic finite element spaces in the spatial domain (with standard Lagrangian FE bases). Hence, the construction of multiresolutions in the spatial domain is not necessary. An efficient multilevel preconditioner is proposed that solves the linear system of equations resulting from the sparse space-time Galerkin discretization with essentially linear complexity (in work and memory). A substantial reduction in the number of the degrees of freedom and CPU time (compared to time-marching discretizations) is demonstrated in numerical experiments.
title Wavelet compressed, modified Hilbert transform in the space-time discretization of the heat equation
topic Numerical Analysis
65N30, 65J15
url https://arxiv.org/abs/2402.10346