Numerical issues in gas flow dynamics with hydraulic shocks using high order finite volume WENO schemes

Ferdinand Uilhoorn

Abstract

Accurate and efficient simulation of the hydraulic shock phenomenon in pipeline systems is of paramount importance. Even though the conservation-law formulation of the governing equations is here strongly advocated, the nonconservative form is still frequently used. This also concerns its mathematical conservative form. We investigated the numerical consequences of using the compressible gas flow model in the latter form while simulating a hydraulic shock. In this context, we also solved two Riemann problems. For the investigation, we used the third-, fifth- and seventh-order accurate weighted essentially non-oscillatory (WENO) scheme along with the Lax–Friedrichs solver at the cell interfaces. Both the classical finite volume WENO scheme and its modification WENO–Z have been implemented. A procedure based on the method of manufactured solutions has been developed to verify whether the numerical code solved correctly the hyperbolic set of equations. We demonstrated that the solutions of the conservative and nonconservative formulations are similar if we have smooth variations in the solution domain. The convective inertia term in the momentum equation should not be ignored. In the presence of shocks, differences in oscillating behavior and slope steepness near the discontinuities were observed. For the hydraulic shock problem, spurious oscillations appeared while using the nonconservative formulation in combination with the WENO–Z reconstruction.
Author Ferdinand Uilhoorn (FEE / DPEGHS)
Ferdinand Uilhoorn,,
- Department of Power Engineering and Gas Heating Systems
Journal seriesJournal of Computational Physics, ISSN 0021-9991, e-ISSN 1090-2716
Issue year2020
Vol404
No1
Pages1-26
Publication size in sheets1.85
Article number109137
Keywords in Englishhydraulic shock Riemann problems weighted essentially non-oscillatory method of manufactured solutions nonconservative scheme
ASJC Classification1706 Computer Science Applications; 3101 Physics and Astronomy (miscellaneous)
Abstract in original languageAccurate and efficient simulation of the hydraulic shock phenomenon in pipeline systems is of paramount importance. Even though the conservation-law formulation of the governing equations is here strongly advocated, the nonconservative form is still frequently used. This also concerns its mathematical conservative form. We investigated the numerical consequences of using the compressible gas flow model in the latter form while simulating a hydraulic shock. In this context, we also solved two Riemann problems. For the investigation, we used the third-, fifth- and seventh-order accurate weighted essentially non-oscillatory (WENO) scheme along with the Lax–Friedrichs solver at the cell interfaces. Both the classical finite volume WENO scheme and its modification WENO–Z have been implemented. A procedure based on the method of manufactured solutions has been developed to verify whether the numerical code solved correctly the hyperbolic set of equations. We demonstrated that the solutions of the conservative and nonconservative formulations are similar if we have smooth variations in the solution domain. The convective inertia term in the momentum equation should not be ignored. In the presence of shocks, differences in oscillating behavior and slope steepness near the discontinuities were observed. For the hydraulic shock problem, spurious oscillations appeared while using the nonconservative formulation in combination with the WENO–Z reconstruction.
DOIDOI:10.1016/j.jcp.2019.109137
URL https://www.sciencedirect.com/science/article/pii/S0021999119308423?via%3Dihub
Languageen angielski
Score (nominal)140
Score sourcejournalList
ScoreMinisterial score = 140.0, 10-06-2020, ArticleFromJournal
Publication indicators WoS Citations = 0; Scopus Citations = 0; Scopus SNIP (Source Normalised Impact per Paper): 2018 = 1.737; WoS Impact Factor: 2018 = 2.845 (2) - 2018=3.321 (5)
Citation count*
Cite
Share Share

Get link to the record


* presented citation count is obtained through Internet information analysis and it is close to the number calculated by the Publish or Perish system.
Back
Confirmation
Are you sure?