An Asymptotic-Preserving Scheme for Isentropic Flow in Pipe Networks

Abstract

We consider the simulation of isentropic flow in pipelines and pipe networks. Standard operating conditions in pipe networks suggest an emphasis to simulate low Mach and high friction regimes -- however, the system is stiff in these regimes and conventional explicit approximation techniques prove quite costly and often impractical. To combat these inefficiencies, we develop a novel asymptotic-preserving scheme that is uniformly consistent and stable for all Mach regimes. The proposed method for a single pipeline follows the flux splitting suggested in [Haack et al., Commun. Comput. Phys., 12 (2012), pp. 955--980], in which the flux is separated into stiff and non-stiff portions then discretized in time using an implicit-explicit approach. The non-stiff part is advanced in time by an explicit hyperbolic solver; we opt for the second-order central-upwind finite volume scheme. The stiff portion is advanced in time implicitly using an approach based on Rosenbrock-type Runge-Kutta methods, which ultimately reduces this implicit stage to a discretization of a linear elliptic equation. To extend to full pipe networks, the scheme on a single pipeline is paired with coupling conditions defined at pipe-to-pipe intersections to ensure a mathematically well-posed problem. We show that the coupling conditions remain well-posed in the low Mach/high friction limit -- which, when used to define the ghost cells of each pipeline, results in a method that is accurate across these intersections in all regimes. The proposed method is tested on several numerical examples and produces accurate, non-oscillatory results with run times independent of the Mach number.

References in corpus (13)

• All speed scheme for the low mach number limit of the Isentropic Euler equation
• An Asymptotic-Preserving all-speed scheme for the Euler and Navier-Stokes equations
• A structure-preserving staggered semi-implicit finite volume scheme for continuum mechanics
• Nonlinear waves in networks: a simple approach using the sine-Gordon equation
• Second order Implicit-Explicit Total Variation Diminishing schemes for the Euler system in the low Mach regime
• Simulation of non-Newtonian viscoplastic flows with a unified first order hyperbolic model and a structure-preserving semi-implicit scheme
• A staggered semi-implicit hybrid finite volume / finite element scheme for the shallow water equations at all Froude numbers
• High order well-balanced asymptotic preserving finite difference WENO schemes for the shallow water equations in all Froude numbers
• Wave dynamics on networks: method and application to the sine-Gordon equation
• Stability and asymptotic analysis for instationary gas transport via relative energy estimates
• Coupling Conditions for Water Waves at Forks
• An asymptotic-preserving discretization scheme for gas transport in pipe networks
• An asymptotic preserving hybrid-dG method for convection-diffusion equations on pipe networks