Finite volume methods for multi-component Euler equations with source terms. (12th October 2017)

• Record Type: Journal Article
• Authors: Bermúdez, Alfredo; López, Xián; Vázquez-Cendón, M. Elena
• Is Part Of: Computers & fluids. Volume 156(2017)
• Page Start: 113
• Page End: 134
• Publication Date: 2017-10-12
• Subjects: Multi-component gas flow -- Non-linear hyperbolic systems with sources -- Segregated scheme -- Finite volume method -- Well-balanced scheme; Fluid dynamics -- Data processing -- Periodicals; 532.050285
• DOI: 10.1016/j.compfluid.2017.07.004
• Languages: English
• ISSNs: 0045-7930

Abstract

Highlights: Euler equations coupled with species equations are solved with a segregated scheme. Flux vectors of both systems depend on the conservative variables, time and space. New discretizations are made so as to ensure that the full scheme is wellbalanced. The positivity of the new method for the species concentrations is proved. Numerical results show the performance of the new method also in a real gas network. Abstract: A first-order well-balanced finite volume scheme for the solution of a multi-component gas flow model in a pipe on non-flat topography is introduced. The mathematical model consists of Euler equations with source terms which arise from heat exchange, and gravity and viscosity forces, coupled with the mass conservation equations of species. We propose a segregated scheme in which the Euler and species equations are solved separately. This methodology leads to a flux vector in the Euler equations which depends not only on the conservative variables but also on time and space variables through the gas composition. This fact makes necessary to add some artificial viscosity to the usual numerical flux which is done by introducing an additional source term. Besides, in order to preserve the positivity of the species concentrations, we discretize the flux in the mass conservation equations for species, in accordance with the upwind discretization of the total mass conservation equation in the Euler system. Moreover, as proposed in a previous reference by the authors, [5], the discretizations of the flux and source terms are made so as to ensure that the full scheme is well-balanced. Numerical tests including both academic and real gas network problems are solved, showing the performance of the proposed methodology.