Given the increasing importance of multi-energy carrier system modeling, this paper focuses on modeling the dynamics that occur in natural gas transmission lines. In high-pressure pipelines, these dynamics are governed by a set of hyperbolic partial differential equations. Dissipative finite volume discretization schemes are proposed that honor the discrete maximum principle. The convergence of the proposed discretization schemes is investigated and compared with existing methods suggested in the literature for natural gas transmission such as the implicit cell-centered method. High order time discretization methods are also tested and their suitability with respect to the maximum principle is discussed.