In this work, various high-accuracy numerical schemes for transport problems in fractured media are further developed and compared. Specifically, to capture sharp gradients and abrupt changes in time, schemes with low order of accuracy are not always sufficient. To this end, discontinuous Galerkin up to order two, Streamline Upwind Petrov-Galerkin, and finite differences, are formulated. The resulting schemes are solved with sparse direct numerical solvers. Moreover, time discontinuous Galerkin methods of order one and two are solved monolithically and in a decoupled fashion, respectively, employing finite elements in space on locally refined meshes. Our algorithmic developments are substantiated with one regular fracture network and several further configurations in fractured media with large parameter contrasts on small length scales. Therein, the evaluation of the numerical schemes and implementations focuses on three key aspects, namely accuracy, monotonicity, and computational costs.
