Discrete Vector Calculus and Helmholtz Hodge Decomposition for Classical Finite Difference Summation by Parts Operators

Abstract

In this article, discrete variants of several results from vector calculus are studied for classical finite difference summation by parts operators in two and three space dimensions. It is shown that existence theorems for scalar/vector potentials of irrotational/solenoidal vector fields cannot hold discretely because of grid oscillations, which are characterised explicitly. This results in a non-vanishing remainder associated to grid oscillations in the discrete Helmholtz Hodge decomposition. Nevertheless, iterative numerical methods based on an interpretation of the Helmholtz Hodge decomposition via orthogonal projections are proposed and applied successfully. In numerical experiments, the discrete remainder vanishes and the potentials converge with the same order of accuracy as usual in other first order partial differential equations. Motivated by the successful application of the Helmholtz Hodge decomposition in theoretical plasma physics, applications to the discrete analysis of magnetohydrodynamic (MHD) wave modes are presented and discussed.

References in corpus (10)

• Review of Summation-by-parts schemes for initial-boundary-value problems
• Velocity Field of Compressible MHD Turbulence: Wavelet Decomposition and Mode Scalings
• On discretely entropy conservative and entropy stable discontinuous Galerkin methods
• Summation-By-Parts Operators and High-Order Quadrature
• Summation-by-parts operators for correction procedure via reconstruction
• Shallow water equations: Split-form, entropy stable, well-balanced, and positivity preserving numerical methods
• Extended Skew-Symmetric Form for Summation-by-Parts Operators and Varying Jacobians
• Mimetic Properties of Difference Operators: Product and Chain Rules as for Functions of Bounded Variation and Entropy Stability of Second Derivatives
• Some Notes on Summation by Parts Time Integration Methods
• Numerical generation of vector potentials from specified magnetic fields

Cited by in corpus (3)

• Summation-by-parts operators for general function spaces: The second derivative
• A new variational discretization technique for initial value problems bypassing governing equations
• A New Class of $A$ Stable Summation by Parts Time Integration Schemes with Strong Initial Conditions