Algebraic flux correction schemes preserving the eigenvalue range of symmetric tensor fields

Abstract

This work extends the algebraic flux correction (AFC) paradigm to finite element discretizations of conservation laws for symmetric tensor fields. The proposed algorithms are designed to enforce discrete maximum principles and preserve the eigenvalue range of evolving tensors. To that end, a continuous Galerkin approximation is modified by adding a linear artificial diffusion operator and a nonlinear antidiffusive correction. The latter is decomposed into edge-based fluxes and constrained to prevent violations of local bounds for the minimal and maximal eigenvalues. In contrast to the flux-corrected transport (FCT) algorithm developed previously by the author and existing slope limiting techniques for stress tensors, the admissible eigenvalue range is defined implicitly and the limited antidiffusive terms are incorporated into the residual of the nonlinear system. In addition to scalar limiters that use a common correction factor for all components of a tensor-valued antidiffusive flux, tensor limiters are designed using spectral decompositions. The new limiter functions are analyzed using tensorial extensions of the existing AFC theory for scalar convection-diffusion equations. The proposed methodology is backed by rigorous proofs of eigenvalue range preservation and Lipschitz continuity. Convergence of pseudo time-stepping methods to stationary solutions is demonstrated in numerical studies.