“…This method is a hybrid concept between finite elements and finite volume discretizations that features some desirable properties including the ability to deal with unstructured meshes on arbitrarily shaped domains, the conservativity of inter-element fluxes, and the feasibility of error estimates in L 2 and H 1 norms. FVE methods have historically been applied for flow equations [11,12,34,46] and recently also for several applicative time-dependent convection-diffusion problems [8,9,17,18,35,45]. Other numerical methods for similar problems include finite differences [20,53,59], finite volume methods [1,5,51], and finite elements [4,21,60].…”