Integral Representations and Quadrature Schemes for the Modified Hilbert Transformation

TL;DR

This paper develops quadrature schemes for accurately computing matrices involving the modified Hilbert transformation, crucial in finite element discretizations of parabolic and hyperbolic problems, achieving exponential convergence.

Contribution

It introduces weakly singular integral representations and quadrature schemes that handle singularities effectively, enabling machine-precision calculations for arbitrary polynomial degrees and non-uniform meshes.

Findings

Quadrature schemes achieve exponential convergence.

Numerical results confirm high accuracy and efficiency.

Applicable to various polynomial degrees and mesh types.

Abstract

We present quadrature schemes to calculate matrices, where the so-called modified Hilbert transformation is involved. These matrices occur as temporal parts of Galerkin finite element discretizations of parabolic or hyperbolic problems when the modified Hilbert transformation is used for the variational setting. This work provides the calculation of these matrices to machine precision for arbitrary polynomial degrees and non-uniform meshes. The proposed quadrature schemes are based on weakly singular integral representations of the modified Hilbert transformation. First, these weakly singular integral representations of the modified Hilbert transformation are proven. Second, using these integral representations, we derive quadrature schemes, which treat the occurring singularities appropriately. Thus, exponential convergence with respect to the number of quadrature nodes for the proposed quadrature schemes is achieved. Numerical results, where this exponential convergence is observed, conclude this work.