Operator splitting method for simulation of dynamic flows in natural gas pipeline networks
Abstract
Here, we develop an operator splitting method to simulate flows of isothermal compressible natural gas over transmission pipelines. The method solves a system of nonlinear hyperbolic partial differential equations (PDEs) of hydrodynamic type for mass flow and pressure on a metric graph, where turbulent losses of momentum are modeled by phenomenological DarcyWeisbach friction. Mass flow balance is maintained through the boundary conditions at the network nodes, where natural gas is injected or withdrawn from the system. Gas flow through the network is controlled by compressors boosting pressure at the inlet of the adjoint pipe. Our operator splitting numerical scheme is unconditionally stable and it is second order accurate in space and time. The scheme is explicit, and it is formulated to work with general networks with loops. We test the scheme over range of regimes and network configurations, also comparing its performance with performance of two other state of the art implicit schemes.
 Authors:
 Univ. of Illinois, UrbanaChampaign, IL (United States). Dept. of Mathematics
 Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
 Univ. of New Mexico, Albuquerque, NM (United States). Dept. of Mathematics & Statistics; L. D. Landau Inst. for Theoretical Physics, Moscow (Russia)
 Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Energy System Center, Skoltech, Moscow (Russia)
 Publication Date:
 Research Org.:
 Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
 Sponsoring Org.:
 USDOE Office of Electricity Delivery and Energy Reliability (OE); USDOE Advanced Research Projects Agency  Energy (ARPAE)
 OSTI Identifier:
 1394973
 Report Number(s):
 LAUR1629073
Journal ID: ISSN 01672789
 Grant/Contract Number:
 AC5206NA25396; AR0000673
 Resource Type:
 Journal Article: Accepted Manuscript
 Journal Name:
 Physica. D, Nonlinear Phenomena
 Additional Journal Information:
 Journal Volume: 361; Journal ID: ISSN 01672789
 Publisher:
 Elsevier
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 97 MATHEMATICS AND COMPUTING; Computer Science; Energy Sciences; Mathematics; Natural Gas; Computational Methods
Citation Formats
Dyachenko, Sergey A., Zlotnik, Anatoly, Korotkevich, Alexander O., and Chertkov, Michael. Operator splitting method for simulation of dynamic flows in natural gas pipeline networks. United States: N. p., 2017.
Web. doi:10.1016/j.physd.2017.09.002.
Dyachenko, Sergey A., Zlotnik, Anatoly, Korotkevich, Alexander O., & Chertkov, Michael. Operator splitting method for simulation of dynamic flows in natural gas pipeline networks. United States. doi:10.1016/j.physd.2017.09.002.
Dyachenko, Sergey A., Zlotnik, Anatoly, Korotkevich, Alexander O., and Chertkov, Michael. 2017.
"Operator splitting method for simulation of dynamic flows in natural gas pipeline networks". United States.
doi:10.1016/j.physd.2017.09.002.
@article{osti_1394973,
title = {Operator splitting method for simulation of dynamic flows in natural gas pipeline networks},
author = {Dyachenko, Sergey A. and Zlotnik, Anatoly and Korotkevich, Alexander O. and Chertkov, Michael},
abstractNote = {Here, we develop an operator splitting method to simulate flows of isothermal compressible natural gas over transmission pipelines. The method solves a system of nonlinear hyperbolic partial differential equations (PDEs) of hydrodynamic type for mass flow and pressure on a metric graph, where turbulent losses of momentum are modeled by phenomenological DarcyWeisbach friction. Mass flow balance is maintained through the boundary conditions at the network nodes, where natural gas is injected or withdrawn from the system. Gas flow through the network is controlled by compressors boosting pressure at the inlet of the adjoint pipe. Our operator splitting numerical scheme is unconditionally stable and it is second order accurate in space and time. The scheme is explicit, and it is formulated to work with general networks with loops. We test the scheme over range of regimes and network configurations, also comparing its performance with performance of two other state of the art implicit schemes.},
doi = {10.1016/j.physd.2017.09.002},
journal = {Physica. D, Nonlinear Phenomena},
number = ,
volume = 361,
place = {United States},
year = 2017,
month = 9
}

A model for simulating the dynamic behavior of gas distribution pipeline networks is developed. An integral part of this work is a presentation of a new implicit integration method and a special ''crossingbranch method.'' The first method can accomodate systems exhibiting both slow and fast responses and prevents numerical oscillations in the solution. In the crossingbranch method, the solution of very large systems describe by sparse matrices is provided. It is based upon special properties of the matrix, which in comparison with methods assuming general matrices, makes the numerical solution more effective. Based on these methods a universal model formore »

A model for simulating adaptive, dynamic flows on networks: Application to petroleum infrastructure
Simulation models can greatly improve decisions meant to control the consequences of disruptions to critical infrastructures. We describe a dynamic flow model on networks purposed to inform analyses by those concerned about consequences of disruptions to infrastructures and to help policy makers design robust mitigations. We conceptualize the adaptive responses of infrastructure networks to perturbations as market transactions and business decisions of operators. We approximate commodity flows in these networks by a diffusion equation, with nonlinearities introduced to model capacity limits. To illustrate the behavior and scalability of the model, we show its application first on two simple networks, thenmore » 
Spectral implementation of a new operator splitting method for solving partial differential equations
We extend and review a method of operator splitting for the numerical solution of a class of partial differential equations introduced by us recently, which we call the {ital odd}{ital even} {ital splitting}. This method extends the applicability of higher order composition methods that have appeared recently in the literature. These composition methods can be considerably more efficient than conventional methods and additionally require at least a factor of 2 less storage capacity. In addition, in some cases, it is possible to use an explicit treatment of certain terms that need to be treated implicitly in other time stepping schemes.more » 
On the use of operatorsplitting method for the equations of combustion
A model problem is analyzed to illustrate the operatorsplitting method for combustion problems with low Machnumber flows. In the limit of high Damkohler number, it is shown that the solution of the combustion equations can be separated at each time step into two parts, a spatially homogeneous explosion calculation over the initial part of the time step and a transport calculation over the remainder of the time step that redistributes the chemical energy released in the first part. The limit of high activation energy allows approximate analytical solution of the explosion equations. The method is applicable primarily for problems withmore » 
Advanced operator splittingbased semiimplicit spectral method to solve the binary phasefield crystal equations with variable coefficients
We present an efficient method to solve numerically the equations of dissipative dynamics of the binary phasefield crystal model proposed by Elder et al. [K.R. Elder, M. Katakowski, M. Haataja, M. Grant, Phys. Rev. B 75 (2007) 064107] characterized by variable coefficients. Using the operator splitting method, the problem has been decomposed into subproblems that can be solved more efficiently. A combination of nontrivial splitting with spectral semiimplicit solution leads to sets of algebraic equations of diagonal matrix form. Extensive testing of the method has been carried out to find the optimum balance among errors associated with time integration, spatialmore »