Paper detail

On numerical methods for hyperbolic PDE with curl involutions

In this paper we present three different numerical approaches to account for curl-type involution constraints in hyperbolic partial differential equations for continuum physics. All approaches have a direct analogy to existing and well-known divergence-preserving schemes for the Maxwell and MHD equations. The first method consists in a generalization of the Godunov-Powell terms, which means adding suitable multiples of the involution constraints to the PDE system in order to achieve the symmetric Godunov form. The second method is an extension of the generalized Lagrangian multiplier (GLM) approach of Munz et al., where the numerical errors in the involution constraint are propagated away via an augmented PDE system. The last method is an exactly involution preserving discretization, similar to the exactly divergence-free schemes for the Maxwell and MHD equations, making use of appropriately staggered meshes. We present some numerical results that allow to compare all three approaches with each other.