Short abstract A new scheme using scalar-valued edge dofs is presented and analyzed to approximate the vector advection-reaction problem. The discrete analysis is done under general assumptions on the parameters, including e.g. divergence-free fields. Numerical results are presented on 3D polyhedral meshes.