An edge-based scheme on polyhedral meshes for vector advection-reaction equations

Abstract

We devise and analyze an edge-based scheme on polyhedral meshes to approximate a vector advection-reaction problem. The well-posedness of the discrete problem is analyzed first under the classical positivity hypothesis of Friedrichs’ systems that requires a lower bound on the lowest eigenvalue of some tensor depending on the model parameters. We also prove stability when the lowest eigenvalue is null or even slightly negative if the mesh size is small enough.A priorierror estimates are established for solutions inW1,q(Ω) withq ∈ ((3/2),2]. Numerical results are presented on three-dimensional polyhedral meshes.