Generalized finite difference schemes with higher order Whitney forms

• Published in: ESAIM. Mathematical modelling and numerical analysis 2021-07, Vol.55 (4), p.1439-1460
• Main Authors: Kettunen, Lauri, Lohi, Jonni, Räbinä, Jukka, Mönkölä, Sanna, Rossi, Tuomo
• Format: Article
• Language: English
• Subjects: Constitutive relationships; Finite difference method; Flow control; Interpolation; Mathematical analysis; Partial differential equations; Wave propagation
• ISSN: 0764-583X; 1290-3841
• DOI: 10.1051/m2an/2021026

Finite difference kind of schemes are popular in approximating wave propagation problems in finite dimensional spaces. While Yee’s original paper on the finite difference method is already from the sixties, mathematically there still remains questions which are not yet satisfactorily covered. In this paper, we address two issues of this kind. Firstly, in the literature Yee’s scheme is constructed separately for each particular type of wave problem. Here, we explicitly generalize the Yee scheme to a class of wave problems that covers at large physics field theories. For this we introduce Yee’s scheme for all problems of a class characterised on a Minkowski manifold by (i) a pair of first order partial differential equations and by (ii) a constitutive relation that couple the differential equations with a Hodge relation. In addition, we introduce a strategy to systematically exploit higher order Whitney elements in Yee-like approaches. This makes higher order interpolation possible both in time and space. For this, we show that Yee-like schemes preserve the local character of the Hodge relation, which is to say, the constitutive laws become imposed on a finite set of points instead of on all ordinary points of space. As a result, the usage of higher order Whitney forms does not compel to change the actual solution process at all. This is demonstrated with a simple example.