6533b7d6fe1ef96bd126700c
RESEARCH PRODUCT
Generalized finite difference schemes with higher order Whitney forms
Tuomo RossiJukka RäbinäSanna MönköläLauri KettunenJonni Lohisubject
Differential equationDifferential formsähkömagnetismiFirst-order partial differential equationdifferential formselectromagnetism010103 numerical & computational mathematics01 natural sciencesdifferentiaaligeometriaMinkowski spaceApplied mathematicsdifferential geometry0101 mathematicsFinite setfinite difference methodMathematicsNumerical AnalysisSpacetimeApplied MathematicsFinite difference methodFinite differencevector-valued formswhitney forms010101 applied mathematicsComputational MathematicsModeling and Simulationelasticityco-vector valued formsAnalysisdescription
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.
year | journal | country | edition | language |
---|---|---|---|---|
2021-07-01 | ESAIM: Mathematical Modelling and Numerical Analysis |