Numerical methods for hydraulic transients in visco-elastic pipes

TL;DR

This paper compares three numerical methods for modeling hydraulic transients in visco-elastic pipes, incorporating viscoelastic wall behavior, flow unsteadiness, and cavitation, validated against experimental data and analytical solutions.

Contribution

It introduces original extensions for simulating viscoelastic pipe behavior within established numerical schemes and compares their accuracy and efficiency in complex transient scenarios.

Findings

Method of Characteristics shows high efficiency in standard cases.

Viscoelastic damping dominates over unsteady friction damping.

All methods agree well with experimental data in tested scenarios.

Abstract

A wide and critical comparison of the capability of Method of Characteristics, Explicit Path-Conservative Finite Volume Method and Semi-Implicit Staggered Finite Volume Method is presented and discussed, in terms of accuracy and efficiency. The viscoelastic behaviour of the pipe wall, the effects of the unsteadiness of the flow on the friction losses, cavitation and cross-sectional changes are considered. The analyses are performed comparing numerical solutions obtained using the three models against experimental data and analytical solutions. Water hammer studies in high density polyethylene pipes, for which laboratory data have been provided, are used as test cases. Considering the viscoelastic mechanical behaviour of plastic materials, a 3-parameter and a multi-parameter linear viscoelastic rheological model are adopted and implemented in each numerical scheme. Original extensions of existing techniques for the numerical treatment of such viscoelastic models are introduced in this work for the first time. After a focused calibration of the viscoelastic parameters, the different performance of the numerical models is investigated. A comparison of the results is presented considering the unsteady wall-shear stress, with a new approach proposed for turbulent flows, or simply considering a quasi-steady friction model. A predominance of the damping effect due to viscoelasticity with respect to the damping effect related to the unsteady friction is confirmed in these contexts. All the numerical methods show a good agreement with the experimental data and a high efficiency of the Method of Characteristics in standard configuration is observed. Three Riemann Problems are chosen and run to stress the numerical methods, considering cross-sectional changes, more flexible materials and cavitation cases. In these demanding scenarios, the weak spots of the Method of Characteristics are depicted.