An inexact Petrov-Galerkin approximation for gas transport in pipeline networks

Herbert Egger,Thomas Kugler,Vsevolod Shashkov

DOI: https://doi.org/10.48550/arXiv.1811.05215

2018-11-13

Abstract:This paper studies the discretization of gas transport in pipeline networks by an inexact Petrov-Galerkin method. A full convergence analysis is presented for single pipes under the assumption of a linear friction law and the possible extension to pipe networks is discussed. The generalization to nonlinear gas transport models and the efficient implementation by hybridization is investigated numerically.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the discretization problem of gas transmission in pipeline networks. Specifically, the authors studied an inexact Petrov - Galerkin method for simulating gas transmission in a single pipeline and pipeline networks, and carried out a comprehensive convergence analysis.

### Main problem description:

1. **Discretization of the gas transmission equation**:
- The paper mainly focuses on how to use the inexact Petrov - Galerkin method to discretize the gas transmission in pipelines. For a single pipeline, assuming that the linear friction law holds, the following basic equations are derived:
\[
A\frac{\partial \rho}{\partial t}+\frac{\partial m}{\partial x} = 0
\]
\[
\frac{\partial m}{\partial t}+\frac{\partial}{\partial x}\left(\frac{m^{2}}{A\rho}+Ap\right)=-\frac{\lambda|m|}{2D}\frac{m}{A\rho}
\]
where \(A\) and \(D\) are the cross - sectional area and diameter of the pipeline respectively, \(\lambda\) is a dimensionless friction parameter, and \(\rho\), \(p\) and \(m\) represent the density, pressure and mass flow rate of the gas respectively.

2. **Simplified model under isothermal flow conditions**:
- Under isothermal flow conditions, \(p = c^{2}\rho\), where \(c\) is the speed of sound. After ignoring the nonlinear terms, a simplified evolution problem can be obtained:
\[
a\frac{\partial p}{\partial t}+\frac{\partial m}{\partial x}=0
\]
\[
b\frac{\partial m}{\partial t}+\frac{\partial p}{\partial x}=-dm
\]
where \(a\) and \(b\) are positive constants, and \(d\) is a state - dependent friction coefficient.

3. **Pipeline network model**:
- For pipeline networks, the flow equations of a single pipeline are coupled through algebraic conditions. For example, at the pipeline junctions, the flow conservation and pressure continuity conditions are satisfied.

4. **Numerical implementation and efficient solution**:
- The paper also explores how to simplify the implementation through the hybridization strategy and allows the handling of non - standard coupling conditions. In addition, this method can also be extended to nonlinear gas transmission models.

### Solutions and contributions:

- **Stability and convergence analysis**: The authors proved that the proposed discretization method is stable and can obtain the optimal convergence rate.
- **High - order approximation**: This method has potential high - order accuracy and is suitable for practically relevant scales.
- **Extension to nonlinear models**: Through numerical experiments, it is verified that this method is still effective under the nonlinear friction model and maintains second - order convergence.

In summary, the main purpose of this paper is to develop an efficient numerical method to simulate gas transmission in pipeline networks and verify its effectiveness and accuracy through strict mathematical analysis and numerical experiments.