Steady Flow of Natural Gas in Pipeline Networks via Solution of a Nonlinear Differential-Algebraic System of Equations

TL;DR

This paper develops a method to analyze steady gas flow in pipeline networks considering gravity and inertia, using differential-algebraic equations and Newton-Raphson, revealing gravity's importance but inertial effects' insignificance.

Contribution

It introduces a novel approach to solve coupled nonlinear differential-algebraic equations for pipeline flow, incorporating gravity and inertial effects, and demonstrates the method's effectiveness on real data.

Findings

Gravity significantly affects flow dynamics.

Inertial effects are negligible in studied cases.

Method applicable to various pipeline networks.

Abstract

In the consideration of steady-state flow of gas in pipeline networks, the exclusion of gravity and nonlinear inertial effects (convective acceleration) leads to a fortuitous simplification in the governing equations to yield a system of nonlinear algebraic equations. Consequently, there are no studies that quantify the effect of gravity and inertial effects on the flow of gas in pipeline networks or delineate regimes of flow conditions wherein the effects are significant or negligible. In addressing this need, we consider the steady-state flow equations in pipeline networks without neglecting the gravitational and inertial terms and in place of a system of algebraic equations (one for each pipe), this approach results in a nonlinear system of first-order ordinary differential equations (ODEs) which are coupled through algebraic equations that appear in the form of boundary conditions on the pressure and balance of mass flows at either end. One of our main contributions in this article is to demonstrate how the Newton-Raphson algorithm can still be used to solve the coupled nonlinear differential-algebraic system by utilizing the appropriate forward sensitivity ODEs to evaluate the Jacobian terms arising in the iterative scheme. We also propose a variable transformation that alleviates the poor scaling of the ODE, and we introduce a two-point collocation scheme as a coarse approximation of the system from which to find initial guesses for the Newton iterations. Simulation studies were conducted for a single pipe as well as a large-scale pipeline network with real data. From these studies, we concluded that while the effect of gravity is important, the inertial effect was negligible in all cases. The proposed methodology is applicable to a wide class of pipeline and thermal-fluid networks beyond natural gas, including liquid pipelines and hydrogen transport.