Numerical methods for hydraulic transients in visco-elastic pipes

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.

References in corpus (1)

• Well balancing of the SWE schemes for moving-water steady flows

Cited by in corpus (8)

• Hyperbolic models for the spread of epidemics on networks: kinetic description and numerical methods
• Modeling blood flow in viscoelastic vessels: the 1D augmented fluid-structure interaction system
• Modeling blood flow in networks of viscoelastic vessels with the 1-D augmented fluid-structure interaction system
• Computational hemodynamics in arteries with the one-dimensional augmented fluid-structure interaction system: viscoelastic parameters estimation and comparison with in-vivo data
• Uncertainty quantification of viscoelastic parameters in arterial hemodynamics with the a-FSI blood flow model
• Dam break in rectangular channels with different upstream-downstream widths
• Dynamic Cavitation Inception by Wave Propagation Across Solid-Fluid Interface with Varying Solid Surface Wettability
• Multiscale constitutive framework of 1D blood flow modeling: Asymptotic limits and numerical methods