Steady State Gas Flow in Pipeline Networks: Existence and Uniqueness of Solution
Abstract: In this paper we discuss the uniqueness and existence of solution to a real gas flow network by employing graph theory. A directed graph is an efficient way to represent a gas network. We consider steady state real gas flow network that includes pipelines, compressors, and the connectors. The pipelines and compressors are represented as edges of the graph and the interconnecting points are represented as nodes of the graph representing the network. We show that a unique solution of such a system exists. We use monotonicity property of a mapping to proof uniqueness, and the contraction mapping theorem is used to prove existence.

1. Introduction

In this paper we investigate the existence and uniqueness of solution to a real gas flow pipeline network by representing the network by graph. A directed graph is an efficient way to represent a gas network. The pipelines and compressors are represented as edges of the graph and the interconnecting points are represented as nodes of the graph representing the network. In a directed graph representation of a gas network, each edge has been assigned a flow direction. Existence and Uniqueness of Solutions to the Generalized Riemann Problem for Isentropic Flow is discussed in . In that paper a gas flow in networks of pipelines is considered. Their models are based on the generalized Riemann problem formulation, where the flow in each connected pipe section is described by the hyperbolic conservation law supplemented by initial conditions within each section of the flow network. In our case, we consider a steady state and isothermal flow in the pipelines, which enables us to simplify the three conservation laws (mass, momentum, and energy) into a single equation for each pipeline, see (3.1). The flow equation of the gas in a compressor is given in (3.2). As we are employing graph theory to study the gas flow network, reviews of some preliminary concepts from graph theory are presented in the first section of this paper. In the second section we present the mathematical model for steady state real gas flow in a pipeline network. Section three contains the discussion about existence and uniqueness of solution for a real gas flow network.

2. Preliminary

Graph Theory

A graph $G=\left(V,E\right)$ is a structure consisting of a finite set V of elements called vertices or nodes and a set E whose elements are called edges . A directed graph or digraph is a graph in which each edge has a direction from one node to another. A walk of a graph G is defined as a finite alternating sequence of vertices and edges, beginning and ending with vertices, such that each edge is incident with the two vertices immediately preceding and following it. A walk in which no vertex appears more than once is called a path. A path beginning and ending with the same vertex is called a cycle.

A graph G is said to connected if there is at least one path between every pair of vertices in G. A tree is a connected graph with no cycles. A spanning tree T of G is a tree consisting of all vertices in G. For a given spanning tree T of G, any edge in G which is not in the tree T is called a chord. A basic result from graph theory states that adding a chord to a spanning tree T will create exactly one cycle. Such a cycle formed by adding a chord to a spanning tree is called a fundamental cycle.

Theorem 1: A tree with n vertices has $n-1$ edges.

Proof: The theorem will be proved by induction. on the number of vertices. It is easy to see that the theorem is true for $n=1,2,3$. Assume that the theorem holds true for all trees with fewer than n vertices. Let us now consider a tree T with n vertices. In T let ${e}_{k}$ be an edge with end vertices ${v}_{i}$ and ${v}_{j}$. Since there is one and only one path between every pair of vertices in a tree, there is no other path between ${v}_{i}$ and ${v}_{j}$ other than ${e}_{k}$. Therefore, deletion of ${e}_{k}$ from T will disconnect the graph. Furthermore, $T-{e}_{k}$ consists of exactly two components, and there are no cycles in T, each of these components is a tree. Both of these trees ahve fewer than n vertices each, and therefore, by the induction assumption, each contains one less edge than the vertices in it. Thus $T-{e}_{k}$ has $n-2$ edges. Hence T has exactly $n-1$ vertices.

Theorem 2: Let n and e be the number of vertices and edges, respectively, in a connected graph G. Let T be a spanning tree of G. Then

1) The number of edges in T is $n-1$ and the number of chords corresponding to the spanning tree T is $e-n+1$ chords.

2) The number of the fundamental cycles corresponding to the spanning tree T is exactly the number of chords, $e-n+1$.

Proof:

1) Since a spanning tree is a tree, theorem 1 can be used to show the first part of 1). Since a chord is an edge in G which is not in the given spanning tree T, the number of chords corresponding to T is the total number of edges e in the graph G minus those $n-1$ edges which belong to the spanning tree T, i.e. $e-n+1$.

2) From part 1), we know that we have only $e-n+1$ chords. Each chord creates a fundamental cycle, and hence the number of the fundamental cycles with respect to the given spanning tree is $e-n+1$.

Incidence matrix of a digraph

Let G be a directed graph with n vertices and e edges, with no circuits (cycles). The node-edge incident matrix A of G is an n × e matrix, defined as

${a}_{ij}=\left(\begin{array}{l}1,\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\text{edge}\text{\hspace{0.17em}}j\text{\hspace{0.17em}}\text{is}\text{\hspace{0.17em}}\text{coming}\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}\text{to}\text{\hspace{0.17em}}\text{node}\\ -1,\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\text{edge}\text{\hspace{0.17em}}j\text{\hspace{0.17em}}\text{is}\text{\hspace{0.17em}}\text{going}\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}\text{to}\text{\hspace{0.17em}}\text{node}\\ 0,\text{\hspace{0.17em}}\text{otherwise}.\end{array}$

Theorem 3: If A is the incident matrix of a connected graph with n nodes, the rank of A is $n-1.$

Proof: Since there are exactly one −1 and one 1 in every column of the incidence matrix A, the sum of all these row vectors is 0. Thus the n vectors are not linearly independent. Therefore, the rank of A is less than n. Since the graph G is connected, it cannot be partitioned like

$A=\left(\begin{array}{cc}{A}_{1}& 0\\ 0& {A}_{2}\end{array}\right)$

such that ${A}_{1}$ is with m rows and ${A}_{2}$ with $n-m$ rows. In other words, no m by m sub matrix of A can be found, for $m\le n-1$, such that the sum of those m rows is equal to zero. Since there are only three constants −1, 0, and 1 in this field, the additions of all vectors taken m at a time for $m=1,2,\cdots ,n-1$ exhausts all possible linear combinations of $n-1$ row vectors. Thus we have just shown that no linear combination of m row vectors of A, for $m\le n-1$ can be equal to zero. Therefore, the rank of A must be at least $n-1$. Since the rank of A is no more than $n-1$ and no less than $n-1$, it must be exactly equal to $n-1$.

A matrix ${A}_{f}$, which is $n-1$ by e obtained from A by deleting one of its rows is called reduced incident matrix.

Theorem 4: The reduced incidence matrix of a tree is invertible.

Proof: A is the incidence matrix of a tree means, by theorem 2, it has n rows and $n-1$ columns. Hence, ${A}_{f}$ is an $n-1$ by $n-1$ matrix. By theorem 3, A has rank $n-1$ implies that ${A}_{f}$ is of rank $n-1$. Since ${A}_{f}$ is a square matrix of order $n-1$ and has rank $n-1$, ${A}_{f}$ is invertible.

Cycle (circuit) matrix of a digraph

Let ${G}^{\prime }$ be the un-directed version of a digraph G, i.e., ${G}^{\prime }$ is the graph G with out considering the directions. Each cycle in ${G}^{\prime }$, after being assigned an arbitrary orientation, can be represented by a vector whose components are are −1, 1, or 0 according to whether and how the edge is included in the cycle. A cycle matrix B is a matrix where each row corresponds to a cycle vector, and is defined by

${b}_{ij}=\left(\begin{array}{l}1,\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\text{cycle}\text{\hspace{0.17em}}i\text{\hspace{0.17em}}\text{contains}\text{\hspace{0.17em}}\text{edge}\text{\hspace{0.17em}}j\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{their}\text{\hspace{0.17em}}\text{orientation}\text{\hspace{0.17em}}\text{coincides}\\ -1,\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\text{cycle}\text{\hspace{0.17em}}i\text{\hspace{0.17em}}\text{contains}\text{\hspace{0.17em}}\text{edge}\text{\hspace{0.17em}}j\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{their}\text{\hspace{0.17em}}\text{orientation}\text{\hspace{0.17em}}\text{are}\text{\hspace{0.17em}}\text{opposite}\\ 0,\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\text{the}\text{\hspace{0.17em}}i\text{th}\text{\hspace{0.17em}}\text{cycle}\text{\hspace{0.17em}}\text{does}\text{\hspace{0.17em}}\text{not}\text{\hspace{0.17em}}\text{include}\text{\hspace{0.17em}}\text{edge}\text{\hspace{0.17em}}j\end{array}$

A cycle is also called circuit.

A set of fundamental cycles with respect to any spanning tree in a connected digraph are the only independent cycles in the digraph, since each contains an edge not in any of the others. The rest of cycles can obtained as linear combinations of these fundamental cycles. A sub matrix ${B}_{f}$ of the matrix B in which all rows correspond to the fundamental cycles is called Fundamental Cycle matrix. If n is the number of vertices and e the number of edges in a connected digraph, then ${B}_{f}$ is an $e-n+1$ by e matrix, because the number of fundamental cycles is $e-n+1$, each fundamental cycle being produced by one chord. Since the fundamental cycles are linearly independent the rank of ${B}_{f}$ is $e-n+1$.

Theorem 5: Let B and A be, respectively, the circuit matrix and incidence matrix of a digraph such that the columns in B and A are arranged using the same order of edges. Then

$A{B}^{\text{T}}={B}^{\text{T}}A=0$

Proof: Consider the mth row in B and the kth row in A. If the cycle m does not include any edge incident on vertex K, the product of the two rows is clearly zero. If, on the other hand, vertex k in cycle m, there are exactly two edges (say x and y) incident on K that are also in cycle m. This situation can occur in only four different ways. The possible entries in row k of A and row m of B in column positions x and y are tabulated in Table 1 for each of these four cases. In each case, the dot product is zero.

3. Mathematical Model of a Steady State Gas Flow in Pipeline Networks

A typical gas network consists of one or more gas sources, one or more gas deliveries, pipelines, compressors, and other devices, such as valves and regulators  . The compressors are installed in the network to increase the gas pressure so that the gas can flow through the pipeline to the locations where it is consumed. Valves and regulators provide control of the gas flow rate, prevent excessive growth of pressure in the network.

Table 1. Possible entries in row k of A and row m of Bin column positions x and y.

We consider the three basic components of a gas network to model the gas flow in pipeline networks:

1) pipelines,

2) compressors, and

3) nodes (interconnection points).

A directed graph is an efficient way to represent a gas network. The pipelines and compressors are represented as edges of the graph and the interconnecting points are represented as nodes of the graph representing the network. In a directed graph representation of a gas network, each edge has been assigned a flow direction. If the actual direction of the gas flow coincides with the direction of the edge, the flow rate is positive, otherwise the flow rate is negative.

3.1. Gas Flow in Pipes

Gas flow in a pipe is described by the three conservation laws (mass, momentum, and Energy) closed by an equation of state. If we consider the flow to be steady state and isothermal, the gas flow can be described by the following equations

$\left(\begin{array}{l}\frac{\partial \left(\rho {u}^{2}+p\right)}{\partial x}=-\frac{fu|u|}{2D}\\ P=P\left(\rho ,T\right)\end{array}$ (1)

where $\rho$ is the gas density, u is the average flow velocity, P is the pressure, T is the temperature, f is friction factor, x is position, and D is the diameter of the pipe. If we use the volumetric flow rate Q and pressure P as flow variables, and

use the substitution $c=\frac{P}{\rho }$ the above equations can be written as:

$\frac{\partial \left(\frac{{Q}^{2}c}{P}+p\right)}{\partial x}=-\frac{fQ|Q|}{2D}$ (2)

If we consider that c is the average of its values at the two end points of the pipe, we can integrate the above equation to get

$\frac{{P}_{u}^{2}}{2}-c{Q}^{2}log{P}_{u}-\frac{{P}_{d}^{2}}{2}+c{Q}^{2}log{P}_{d}+\frac{fQ|Q|L}{2D}=0$ (3)

The logarithmic terms are because of the inertia term included in eqn which is neglected by many authors.

3.2. Compressor Equations

Compressors are installed in the gas network to transport gas and compensate for the loss of energy due to frictional resistance which results in a loss of pressure at the downstream of the pipe. The flow rate through the compressor and pressures at the upstream and downstream of the compressor are related by the following equation.

$\eta \ast W={z}_{cu}RT{Q}_{c}\frac{\gamma }{\gamma -1}\left[{\left(\frac{{P}_{cd}}{{P}_{cu}}\right)}^{\frac{\gamma -1}{\gamma }}-1\right]\text{ }\text{for}\text{\hspace{0.17em}}\text{ }c=1:t$ (4)

where:

${P}_{cu}$ is inlet pressure.

${P}_{cd}$ is outlet pressure.

${Q}_{c}$ is flow rate through the compressor.

W is power of the compressor.

$\eta$ is compressor efficiency.

$\gamma$ is ratio of specific heats.

${z}_{cu}$ is compressibility factor at entry.

3.3. Node Equations

We can categorize the nodes in the network as source nodes (where gas is supplied to the system), junction nodes (where two or more edges are connected), and sink nodes (where gas is delivered for consumption). At a source nodes either the pressure or the flow rate is specified. While at junction or a sink node a flow rate balance is made.

3.3.1. Pressure Node

If the pressure value ${S}_{p}$ is specified at a node of the network, then we get the equation

$P={S}_{p}$ (5)

at this node.

3.3.2. Flow Node

The node at which we balance the flow rates coming into and going out of the is called flow node. Suppose a flow node connects N edges (pipes or compressors) and has a nodal flow rate ${S}_{q}$. ${S}_{q}$ is positive if gas is added in to the network, negative if gas is tapped out of the network, and zero if there is no gas entering or leaving the network at this node.

Then balancing the flow rate at this node gives us

$\underset{j\in I{n}_{i}}{\sum }{Q}_{j}-\underset{j\in Ou{t}_{i}}{\sum }{Q}_{j}={S}_{q}\text{ }\text{for}\text{ }\text{\hspace{0.17em}}i=1:{N}_{q}$ (6)

$I{n}_{i}$ is the set of incoming flows in to the ith Flow Node.

$Ou{t}_{i}$ is the set of outgoing flows from the ith Flow Node.

Therefore, the steady state, isothermal gas flow in pipeline networks is described by a non-linear simultaneous equation containing (3), (4), (5), and (6).

This system of equations can be written as

$F\left(x\right)=y,\text{ }x=\left(P,{Q}_{p},{Q}_{c}\right)$ (7)

4. Existence and Uniqueness of the Solution to a Network Which Can Be Represented by a Tree

In the present section, proofs for existence and uniqueness of solutions to a gas flow network will be given. First, we present a proof of uniqueness of solutions which is based on a monotonicity property of a mapping. Existence is then proved with the help of the contraction mapping theorem (4.2). Assume the gas flow network of n nodes can be represented by a tree. Suppose the pressure value ${S}_{p}$, of one node, say node 1 is given. By theorem (1) this tree has $n-1$ edges. Let ${A}_{l}$ and A be, respectively, the pipe-node incidence matrix and the edge-node incidence matrix of the tree. Let ${A}_{f}$ be the reduced incidence matrix. Now the network equations can be written as:

$\left(\begin{array}{l}{A}_{l}^{\text{T}}\left(\frac{{P}^{2}}{2}-c{Q}_{p}^{2}\otimes \mathrm{log}P\right)=\varphi \left({Q}_{p}\right),\text{\hspace{0.17em}}\text{ }m\text{\hspace{0.17em}}\text{equations}\\ \eta \ast W=b{Q}_{c}\frac{\gamma }{\gamma -1}\left[{\left(\frac{{P}_{cd}}{{P}_{cu}}\right)}^{\frac{\gamma -1}{\gamma }}-1\right],\text{\hspace{0.17em}}\text{\hspace{0.17em}}t=n-m-1\text{\hspace{0.17em}}\text{equations}\\ {P}_{1}={S}_{p}\\ {A}_{f}Q={S}_{q}\end{array}$ (8)

where ${Q}_{p}$ and ${Q}_{c}$ are, respectively, pipe flows and compressor flows. Since ${A}_{f}$ is obtained from A by deleting the row corresponding to the reference node, it is $n-1$ by $n-1$ matrix and is invertible. Hence, there is a unique solution for the last Equation in (8). If we have t compressors involves 2 pressures, there fore the t pressures can be determined if the values of the other t pressures are known. Assuming the t pressures from the compressor equations and one pressure from the reference node are known, we have only $n-1-t=m$ unknown pressures. By substituting the known flow rates and the known pressures in to the pipe equations we get $n-1-t$ equations for $n-1-c$ unknown pressures. Since the rows of ${A}_{l}^{\text{T}}$ are rows of the reduced incidence matrix ${A}_{f}$, they are linearly independent. This implies ${A}_{l}^{\text{T}}$ is invertible and we get the unique solution $\frac{{P}^{2}}{2}-c{Q}_{p}^{2}\otimes \mathrm{log}P$ for the first Equation of (8). This gives us unique P as $\frac{{P}^{2}}{2}-c{Q}_{p}^{2}\otimes \mathrm{log}P$ is one-to-one.

4.1. Uniqueness of the Solution to a Gas Flow Network

In this section, let us assume we have a network of pipes, with no compressor, which can be represented by directed graph. Suppose the gas flow network has n nodes. Assume the pressure ${S}_{p}$ is given at a reference node, say at node 1. Then our network system of Equations (7) takes the form:

$\left(\begin{array}{l}{A}^{\text{T}}\left(\frac{{P}^{2}}{2}-c{Q}^{2}\otimes \mathrm{log}P\right)=\varphi \left(Q\right)\\ {P}_{1}={S}_{p}\\ {A}_{f}Q={S}_{q}\end{array}$ (9)

Let ${B}_{f}$ be the reduced cycle matrix with respect to some spanning tree. The size of ${B}_{f}$ is $e-n+1$ by e, and its rank is $e-n+1$. The above system is equivalent to

$\left(\begin{array}{l}{A}^{\text{T}}\left(\frac{{P}^{2}}{2}-c{Q}^{2}\otimes \mathrm{log}P\right)=\varphi \left(Q\right)\\ {P}_{1}={S}_{p}\\ {B}_{f}\varphi \left(Q\right)=0\\ {A}_{f}Q={S}_{q}\end{array}$ (10)

Now let us consider the system

$\left(\begin{array}{l}{B}_{f}\varphi \left(Q\right)=0\\ {A}_{f}Q={S}_{q}\end{array}$ (11)

The system (11) contains only the flow variable Q. The first equation contains $e-n+1$ equations and the second contains $n-1$ equations. Hence, system (11) consists of e equations for e unknown flwo rate variables. First, we show that the solution of the system (11) is unique, and then it can be shown that the first two equations of system (10)

$\left(\begin{array}{l}{A}^{\text{T}}\left(\frac{{P}^{2}}{2}-c{Q}^{2}\otimes \mathrm{log}P\right)=\varphi \left(Q\right)\\ {P}_{1}={S}_{p}\end{array}$ (12)

give unique P. System (12) consists of $e+1$ equations for n unknowns. But, we can eliminate $e-n+1$ equations as they are linar combinations of the others, since $e-n+1$ ${\varphi }_{i}\left(Q\right)$ are expressed in terms of the others in the equation ${B}_{f}\varphi \left(Q\right)=0$.

Now let us show that the solution of system (11) is unique.

Definition: A mapping $\Phi :H\to H$ is said to be strictly monotonic if for every $x,y\in H$ we have

$\left(\Phi \left(x\right)-\Phi \left(y\right),x-y\right)\ge 0,$

and equality holds if and only if $x=y$.

Theorem 6: Let $\Phi :{ℝ}^{d}\to {ℝ}^{d}$, for some positive integer d, is given by

$\Phi \left(x\right)={\left({\Phi }_{1}\left({x}_{1}\right),{\Phi }_{2}\left({x}_{2}\right),{\Phi }_{3}\left({x}_{3}\right),\cdots ,{\Phi }_{d}\left({x}_{d}\right)\right)}^{\text{T}}$

where

${\Phi }_{i}\left({x}_{j}\right)={c}_{j}{x}_{j}|{x}_{j}|,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}1\le j\le d$

with ${c}_{j}>0$. Then $\Phi$ is strictly monotonic.

Proof: For every $x={\left[{x}_{1},{x}_{2},{x}_{3},\cdots ,{x}_{d}\right)}^{\text{T}},y={\left[{y}_{1},{y}_{2},{y}_{3},\cdots ,{y}_{d}\right)}^{\text{T}}\in {ℝ}^{d}$,

$\left(\Phi \left(x\right)-\Phi \left(y\right),x-y\right)=\underset{j=1}{\overset{d}{\sum }}\text{ }\text{ }{c}_{j}\left({x}_{j}|{x}_{j}|-{y}_{j}|{y}_{j}|\right)\left({x}_{j}-{y}_{j}\right)$

The function $h\left(s\right)=s|s|$ is a strictly increasing function for all s. Hence, each term on the right-hand side is non-negative. Thus,

$\left(\Phi \left(x\right)-\Phi \left(y\right),x-y\right)\ge 0.$

Equality holds if and only if every term at the right hand side is zero, so ${x}_{j}={y}_{j}$, for every j. There fore, $\Phi$ is strictly monotonic.

Definition: Let $r>0,t>0$, be two integers, and $d=r+t$. We say an $r×d$ matrix A and a $t×d$ matrix B are perpendicular to each other if they satisfy

1) $\text{rank}\left(A\right)=r$, $\text{rank}\left(B\right)=t$ ;

2) $A{B}^{\text{T}}=B{A}^{\text{T}}$

Let $M=x\in {\Re }^{d}:Ax=0$ and $N=y\in {ℝ}^{d}:By=0$. Then we have M perpendicular to N and

${ℝ}^{d}=M\oplus N.$

Theorem 7: Let matrices A and B be perpendicular to each other. Suppose $\Phi :{ℝ}^{d}\to {ℝ}^{d}$ is strictly monotonic, then foe every $s\in {ℝ}^{d}$, the solution to the system of equations

$\left(\begin{array}{l}AQ=s\\ B\Phi \left(Q\right)=0\end{array}$ (13)

is unique.

Proof: Suppose both u and v are solutions of system (13), then

$\left(\begin{array}{l}A\left(u-v\right)=0\\ B\left(\Phi \left(u\right)-\varphi \left(v\right)\right)=0\end{array}$

Hence, $u-v\in M$, and $\Phi \left(u\right)-\varphi \left(v\right)\in N$. Thus $\left(u-v,\Phi \left(u\right)-\Phi \left(v\right)\right)=0$. Since $\Phi$ is strictly monotonic, the above equation implies that $u=v$. Hence, the solution is unique.

By applying the above theorem, and taking $A={A}_{f}$, $B={B}_{f}$,and $s={s}_{f}$, it is shown the solution of the system (11) is unique.

4.2. Existence of the Solution for a Gas Flow Network

In this section, we discuss the existence of solution to a general gas flow network that contains pipes and compressors. Let us consider a network of m pipes, t compressors, and n nodes. Assume the pressure is given, say at node 1, i.e. ${P}_{1}={S}_{p}$. We use the The Contraction Mapping Theorem to show existence of a solution for such gas flow network.

Definition: Let $\left(X,d\right)$ be a metric space. A function $T:X\to X$ is called a contraction mapping if there exists $k\in ℝ$ such that $0 and

$d\left(T\left(x\right),T\left(y\right)\right)\le kd\left(x,y\right)$

for any $x,y\in X$.

Theorem 8: (The Contraction Mapping Theorem): If $\left(X,d\right)$ is a complete metric space and $T:X\to X$ is a contraction mapping, then T has one and only one fixed point, (i.e., there exists exactly one $x\in X$ such that $T\left(x\right)=x$ ).

Sketch of the proof of the Contraction Mapping Theorem. The complete proof is found in . The proof proceeds in several steps:

1) starting from an arbitrary ${x}_{0}\in X$, construct the sequence ${x}_{n}\subset X$ by taking ${x}_{n}=T\left({x}_{n-1}\right)$ ;

2) prove that $d\left({x}_{n},{x}_{n+1}\right)\le {k}^{n}d\left({x}_{0},{x}_{1}\right)$, and therefore $d\left({x}_{m},{x}_{n}\right)=\frac{{k}^{n}-{k}^{m}}{1-k}d\left({x}_{0},{x}_{1}\right)$ for any natural numbers $n,m$ with $m>n$ ;

3) conclude that $\left\{{x}_{n}\right\}$ is a Cauchy sequence, and thus, since X is complete, converges;

4) prove that $x={\mathrm{lim}}_{n\to \infty }{x}_{n}$ is a fixed point of T;

5) prove that x is the only fixed point of T.

Now, we want to show that the equation

$T\left(x\right)=y$

has locally unique solution. The Contraction Mapping Theorem can be applied as follows: If we can find ${x}_{0}$ such that the jacobian determinant, J, of T at ${x}_{0}$ is different from zero, then we can introduce

$G\left(x\right)=x-{J}^{-1}\left({x}_{0}\right)\left(T\left(x\right)-y\right).$

The equation $x=G\left(x\right)$ is equivalent to $T\left(x\right)=y$. If G is a contraction mapping, then G has a unique fixed point $\stackrel{¯}{x}$ and thus $T\left(\stackrel{¯}{x}\right)=y$.

The jacobian $J\left(x\right)$ is expressed in the following form

$J\left(x\right)=\left(\begin{array}{cc}A\left(x\right)& D\left(x\right)\\ B\left(x\right)& C\left( x \right)\end{array}\right)$

where

$A\left(x\right)={\left({a}_{ij}\left(x\right)\right)}_{\left(m+t\right)×n}$

${a}_{ij}\left(x\right)=\left(\begin{array}{l}{a}_{iu}=\frac{\partial {F}_{i}\left(x\right)}{\partial {P}_{iu}}=-\left({P}_{iu}-\frac{C{Q}_{i}}{{P}_{iu}}\right),\text{ }\text{if}\text{\hspace{0.17em}}\text{Node}\text{\hspace{0.17em}}j\text{\hspace{0.17em}}\text{is}\text{\hspace{0.17em}}\text{the}\text{\hspace{0.17em}}\text{upstream}\text{\hspace{0.17em}}\text{node}\text{\hspace{0.17em}}\text{of}\text{ }\text{\hspace{0.17em}}\text{Pipe}\text{\hspace{0.17em}}\text{ }i\\ {a}_{iu}=\frac{\partial {F}_{m+c}\left(x\right)}{\partial {P}_{cu}}=-\left(\eta W\frac{\gamma -1}{\gamma }+b{Q}_{c}{P}_{cu}^{-\frac{1}{\gamma }}\right),\text{ }\text{if}\text{\hspace{0.17em}}\text{Node}\text{\hspace{0.17em}}j\text{\hspace{0.17em}}\text{is}\text{\hspace{0.17em}}\text{the}\text{\hspace{0.17em}}\text{upstream}\text{\hspace{0.17em}}\text{node}\text{\hspace{0.17em}}\text{of}\text{\hspace{0.17em}}\text{Compressor}\text{\hspace{0.17em}}c\\ {a}_{id}=\frac{\partial {F}_{i}\left(x\right)}{\partial {P}_{id}}=\left({P}_{id}-\frac{C{Q}_{i}}{{P}_{id}}\right),\text{ }\text{if}\text{\hspace{0.17em}}\text{Node}\text{\hspace{0.17em}}j\text{\hspace{0.17em}}\text{is}\text{\hspace{0.17em}}\text{the}\text{\hspace{0.17em}}\text{downstream}\text{\hspace{0.17em}}\text{node}\text{\hspace{0.17em}}\text{of}\text{\hspace{0.17em}}\text{Pipe}\text{\hspace{0.17em}}i\\ {a}_{id}=\frac{\partial {F}_{m+c}\left(x\right)}{\partial {P}_{cd}}=b{Q}_{c}{P}_{cd}^{-\frac{1}{\gamma }},\text{ }\text{if}\text{\hspace{0.17em}}\text{Node}\text{\hspace{0.17em}}j\text{\hspace{0.17em}}\text{is}\text{ }\text{\hspace{0.17em}}\text{the}\text{\hspace{0.17em}}\text{downstream}\text{\hspace{0.17em}}\text{node}\text{\hspace{0.17em}}\text{of}\text{\hspace{0.17em}}\text{Compressor}\text{\hspace{0.17em}}c\\ 0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{else}\end{array}$

$D\left(x\right)={\left({d}_{ij}\left(x\right)\right)}_{\left(m+t\right)×\left(m+t\right)}$

${d}_{ii}\left(x\right)=\left(\begin{array}{l}\frac{\partial {F}_{i}\left(x\right)}{\partial {Q}_{i}}=\left(C\mathrm{log}\frac{{P}_{iu}}{{P}_{id}}+\frac{CfL|{Q}_{i}|}{D}\right)\text{ }\text{if}\text{\hspace{0.17em}}i\le m\\ \frac{\partial {F}_{m+c}\left(x\right)}{\partial {Q}_{c}}=b\frac{\gamma }{\gamma -1}\left({\left(\frac{{P}_{cd}}{{P}_{cu}}\right)}^{\frac{\gamma -1}{\gamma }}-1\right)\text{ }\text{if}\text{\hspace{0.17em}}m

${d}_{ij}\left(x\right)=0\text{ }\text{\hspace{0.17em}}\text{ }\text{if}\text{ }\text{ }\text{ }i\ne j$

$B\left(x\right)={\left({b}_{ij}\left(x\right)\right)}_{n×n}$

${b}_{ij}\left(x\right)=\left(\begin{array}{l}1,\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\text{Node}\text{ }\text{\hspace{0.17em}}j\text{ }\text{\hspace{0.17em}}\text{is}\text{\hspace{0.17em}}\text{the}\text{\hspace{0.17em}}{i}^{\text{th}}\text{\hspace{0.17em}}\text{Pressure}\text{\hspace{0.17em}}\text{Node}\\ 0,\text{ }\text{ }\text{else}\end{array}$

$C\left(x\right)={\left({C}_{ij}\left(x\right)\right)}_{n×\left(m+t\right)}$

${C}_{ij}\left(x\right)=\left(\begin{array}{l}1,\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\text{flow}\text{\hspace{0.17em}}j\text{\hspace{0.17em}}\text{flows}\text{\hspace{0.17em}}\text{into}\text{ }\text{\hspace{0.17em}}\text{the}\text{\hspace{0.17em}}{i}^{\text{th}}\text{\hspace{0.17em}}\text{FlowNode}\\ -1,\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\text{flow}\text{\hspace{0.17em}}j\text{\hspace{0.17em}}\text{flows}\text{\hspace{0.17em}}\text{out}\text{\hspace{0.17em}}\text{of}\text{\hspace{0.17em}}\text{the}\text{\hspace{0.17em}}{i}^{\text{th}}\text{\hspace{0.17em}}\text{Flow}\text{\hspace{0.17em}}\text{Node}\\ 0,\text{\hspace{0.17em}}\text{else}\end{array}$

Choose ${x}_{0}$ such that ${d}_{ij}\left({x}_{0}\right)\ne 0$, ${a}_{iu}\left({x}_{0}\right)<0$ and ${a}_{id}\left({x}_{0}\right)>0$

Then $\mathrm{det}\left(J\left({X}_{0}\right)\right)\ne 0$

Example: Let us consider the following network with 3 pipes, 1 compressor, and 5 junctions.

The Jacobian of the matrix of the system of equations of the gas flow network dipicted in Figure 1 is given by

$J=\left(\begin{array}{ccccccccc}{a}_{1u}& {a}_{1d}& 0& 0& 0& {d}_{11}& 0& 0& 0\\ 0& 0& {a}_{2u}& {a}_{2d}& 0& 0& {d}_{22}& 0& 0\\ 0& 0& 0& {a}_{3u}& {a}_{3d}& 0& 0& {d}_{33}& 0\\ 0& {a}_{4u}& {a}_{4d}& 0& 0& 0& 0& 0& {d}_{44}\\ 1& 0& 0& 0& 0& 0& 0& 0& 0\\ 0& 0& 0& 0& 0& 1& 0& 0& -1\\ 0& 0& 0& 0& 0& 0& -1& 0& 1\\ 0& 0& 0& 0& 0& 0& 1& -1& 0\\ 0& 0& 0& 0& 0& 0& 0& 1& 0\end{array}\right)$

By performing elementary row and column operations on J

$J~\left(\begin{array}{cc}0& {I}_{\left(m+t\right)}\\ {B}^{\ast }& 0\end{array}\right)$

where ${B}^{*}={\left({b}_{ij}^{*}\right)}_{n×n}$

${b}_{ij}^{*}=\left(\begin{array}{l}{b}_{ij},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\text{node}\text{\hspace{0.17em}}i\text{\hspace{0.17em}}\text{is}\text{\hspace{0.17em}}\text{a}\text{\hspace{0.17em}}\text{pressure}\text{\hspace{0.17em}}\text{node}\\ \frac{-{a}_{ku}}{{d}_{kk}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\text{flow}\text{\hspace{0.17em}}k\text{\hspace{0.17em}}\text{flows}\text{\hspace{0.17em}}\text{from}\text{\hspace{0.17em}}\text{node}\text{\hspace{0.17em}}j\text{\hspace{0.17em}}\text{to}\text{\hspace{0.17em}}\text{node}\text{\hspace{0.17em}}i\\ \frac{{a}_{kd}}{{d}_{kk}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\text{flow}\text{\hspace{0.17em}}k\text{\hspace{0.17em}}\text{flows}\text{\hspace{0.17em}}\text{from}\text{\hspace{0.17em}}\text{node}\text{\hspace{0.17em}}i\text{\hspace{0.17em}}\text{to}\text{\hspace{0.17em}}\text{node}\text{\hspace{0.17em}}j\\ \underset{k\in ou{t}_{i}}{\sum }\frac{{a}_{ku}}{{d}_{kk}}-\underset{k\in I{n}_{i}}{\sum }\frac{{a}_{kd}}{{d}_{kk}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}i=j\text{ }\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{node}\text{\hspace{0.17em}}i\text{\hspace{0.17em}}\text{is}\text{\hspace{0.17em}}\text{a}\text{\hspace{0.17em}}\text{flow}\text{\hspace{0.17em}}\text{node}\\ 0,\text{ }\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{else}\end{array}$

$\overline{)det\left(J\right)\ne 0⇔det\left({B}^{*}\right)\ne 0}$

For our example given above

Figure 1. Example of a gas flow network.

${B}^{*}=\left(\begin{array}{ccccc}1& 0& 0& 0& 0\\ \frac{-{a}_{1u}}{{d}_{11}}& \frac{-{a}_{1d}}{{d}_{11}}+\frac{{a}_{4u}}{{d}_{44}}& \frac{{a}_{4d}}{{d}_{44}}& 0& 0\\ 0& \frac{-{a}_{4u}}{{d}_{44}}& \frac{-{a}_{4d}}{{d}_{44}}+\frac{{a}_{2u}}{{d}_{22}}& \frac{{a}_{2d}}{{d}_{22}}& 0\\ 0& 0& \frac{-{a}_{2u}}{{d}_{22}}& \frac{-{a}_{2d}}{{d}_{22}}+\frac{{a}_{3u}}{{d}_{33}}& \frac{{a}_{3d}}{{d}_{33}}\\ 0& 0& 0& \frac{-{a}_{3u}}{{d}_{33}}& \frac{-{a}_{3d}}{{d}_{33}}\end{array}\right)$

Let ${B}_{d}^{*}$ be the matrix obtained by deleting the rows and columns containing the Pressure Nodes, then

$\overline{)det\left({B}^{*}\right)\ne 0⇔det\left({B}_{d}^{*}\right)\ne 0}$

Properties of ${B}_{d}^{*}$ are

1) Diagonally dominant (column), strictly dominant at column j, where node j is connected to a Pressure Node( there is at least one Pressure Node).

2) Irreducible.

Hence, ${B}_{d}^{*}$ is invertible.

$⇒$ ${B}^{*}$ is invertible.

$⇒$ J is invertible.

Define $G\left(x\right)=x-{J}^{-1}\left({x}_{0}\right)\left(F\left(x\right)-y\right)$

$⇒$ ${G}^{\prime }\left(x\right)=I-{J}^{-1}\left({x}_{0}\right)J\left( x \right)$

Since ${G}^{\prime }\left(x\right)$ is continuous at ${x}_{0}$

$\forall ϵ>0$ $\exists \delta$ such that $‖x-{x}_{0}‖<\delta ⇒‖{G}^{\prime }\left(x\right)‖<ϵ$

Choose $ϵ$ such that $0<ϵ<1$

$\forall {x}_{1},{x}_{2}$ in $B\left({x}_{0},\delta \right)$

$‖G\left({x}_{1}\right)-G\left({x}_{2}\right)‖\le ‖{G}^{\prime }\left(\stackrel{¯}{x}\right)‖‖{x}_{1}-{x}_{2}‖<ϵ‖{x}_{1}-{x}_{2}‖|$

If $y\in B\left(F\left({x}_{0}\right),\frac{1-ϵ}{‖{J}^{-1}\left({x}_{0}\right)‖}\delta \right)$, then $G\left(B\left({x}_{0},\delta \right)\right)\subset B\left({x}_{0},\delta \right)$

$⇒$ G is a contraction in $B\left({x}_{0},\delta \right)$

$⇒$ G has a unique fixed point.

$⇒$ $F\left(x\right)=y$ has a unique solution in $B\left({x}_{0},\delta \right)$

5. Conclusion

In this paper, we discussed the existence and uniqueness of solution to the real gas flow network. A review of some basic concepts from graph theory is given as we discussed the flow network as a graph. The flow network is represented by a graph in which pipes and compressors are represented as edges and the junction points as nodes of the graph. The flow is assumed to be steady state and isothermal. The Contraction Mapping Theorem is applied to show the existence of the solution to the considered gas flow network, and its uniqueness is proved using monotonicity property of a mapping.

Nomenclature

$\rho$ Gas density

u Flow velocity

P Pressure

Pu Upstream pipe inlet pressure

Pd Downstream pipe outlet pressure

T Temperature

Q Volumetric flow rate in a pipe

x Position

D Pipe diameter

L Pipe length

$\eta$ Compressor efficiency

W Compressor power

Zcu Comprehensibility factor of a gas at compressor entry

R Gas constant

Qc Flow rate through the compressor

$\gamma$ Ratio of specific heats

Pcu Inlet compressor pressure

Pcd Outlet compressor pressure

Cite this paper: Atena, A. , Tekalign, W. and Muche, T. (2020) Steady State Gas Flow in Pipeline Networks: Existence and Uniqueness of Solution. Journal of Applied Mathematics and Physics, 8, 1155-1167. doi: 10.4236/jamp.2020.86087.
References

   Reigstad, G.A. (2015) Existence and Uniqueness of Solutions to the Generalized Riemann Problem for Isentropic Flow. SIAM Journal on Applied Mathematics, 75, 679-702.
https://doi.org/10.1137/140962759

   Deo, N. (2016) Graph Theory with Applications to Engineering and Computer Science, Mineola. Dover Publications, New York.

   Feistauer, M. (2003) Mathematical and Computational Methods for Compressible Flow. Oxford University Press, Oxford.

   Chorin, A.J. and Marsden, J.E. (2012) A Mathematical Introduction to Fluid Mechanics. 2nd Edition, Springer, New York.

   Brook, R.M. and Schmitt, K. (2009) The Contraction Mapping Principle and Some Applications. Electronic Journal of Differential Equations, 9.

Top