Space-time virtual elements for the heat equation

Sergio Gómez,Lorenzo Mascotto,Andrea Moiola,Ilaria Perugia

DOI: https://doi.org/10.1137/22M154140X

2023-09-04

Abstract:We propose and analyze a space-time virtual element method for the discretization of the heat equation in a space-time cylinder, based on a standard Petrov-Galerkin formulation. Local discrete functions are solutions to a heat equation problem with polynomial data. Global virtual element spaces are nonconforming in space, so that the analysis and the design of the method are independent of the spatial dimension. The information between time slabs is transmitted by means of upwind terms involving polynomial projections of the discrete functions. We prove well posedness and optimal error estimates for the scheme, and validate them with several numerical tests.

Numerical Analysis

What problem does this paper attempt to address?

This paper attempts to solve the problem of numerical discretization of the heat conduction equation in the space - time domain. Specifically, the author proposes and analyzes a space - time virtual element method (VEM) based on the standard Petrov - Galerkin variational formula for solving the heat conduction equation.

### Main Objectives of the Paper

1. **Propose a New Space - Time Virtual Element Method**: The author designs a new space - time VEM, which is applicable to spaces of any dimension and can handle complex polygonal meshes.

2. **Prove the Well - Posedness and Error Estimation of the Method**: The author proves that the proposed discrete scheme is well - posed and gives the optimal a priori error estimate.

3. **Numerical Verification**: The effectiveness of the theoretical results is verified through multiple numerical experiments.

### Characteristics of the Method

- **Local Virtual Element Space**: The local VE space on each space - time element consists of functions that satisfy the heat conduction equation with polynomial data, which makes the method suitable for further extension, for example, to its Trefftz variant.

- **Non - Tensor - Product - Type Mesh**: Although a time - tensor - product (prism) mesh is used, the VE space is not of the tensor - product type. Even for prism elements with a simple base, the proposed VE space is not the same as the standard tensor - product finite - element space.

- **Continuity Constraints across Mesh Faces**: The global VE space involves approximate continuity constraints across mesh faces. Specifically, non - conforming conditions similar to those in the Poisson problem are imposed on time - like faces and discontinuous functions in time are allowed; while information is transmitted through upwind terms on space - like faces. These upwind terms are defined by polynomial projection.

- **Simplified Time - Step Processing**: Since the discrete space is discontinuous in time, the global (expensive) problem can be decomposed into a series of local (cheap) problems for solution.

### Main Contributions

- A novel space - time VEM for the heat conduction equation in spaces of arbitrary dimension is proposed.
- The well - posedness and optimal a priori error estimate of this method are proved.
- The theoretical results are verified by numerical experiments.

The advantage of this method lies in its flexibility and efficiency, especially in dealing with complex geometries and high - order accuracy. In addition, the non - conforming setting enables this method to be designed independently of the spatial dimension, thus simplifying the analysis and implementation process.