  Natural gas ranks second in consumption among primary energy sources in the
United States. The majority of production sites are in remote locations, hence
natural gas needs to be transported through a pipeline network equipped with a
variety of physical components such as compressors, valves, etc. Thus, from the
point of view of both economics and reliability, it is desirable to achieve
optimal transportation of natural gas using these pipeline networks. The
physics that governs the flow of natural gas through various components in a
pipeline network is governed by nonlinear and non-convex equality and
inequality constraints and the most general steady-flow operations problem
takes the form of a Mixed Integer Nonlinear Program. In this paper, we consider
one example of steady-flow operations -- the Optimal Gas Flow (OGF) problem for
a natural gas pipeline network that minimizes the production cost subject to
the physics of steady-flow of natural gas. The ability to quickly determine
global optimal solution and a lower bound to the objective value of the OGF for
different demand profiles plays a key role in efficient day-to-day operations.
One strategy to accomplish this relies on tight relaxations to the nonlinear
constraints of the OGF. Currently, many nonlinear constraints that arise due to
modeling the non-ideal equation of state either do not have relaxations or have
relaxations that scale poorly for realistic network sizes. In this work, we
combine recent advancements in the development of polyhedral relaxations for
univariate functions to obtain tight relaxations that can be solved within a
few seconds on a standard laptop. We demonstrate the quality of these
relaxations through extensive numerical experiments on very large scale test
networks available in the literature and find that the proposed relaxation is
able to prove optimality in 92% of the instances.
