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

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 priori error estimates are established for solutions in W^1,q(Ω) with q ∈ (3/2),2). Numerical results are presented on three-dimensional polyhedral meshes.