Numerical methods for hydraulic transients in visco-elastic pipes

Abstract

In technical applications involving transient fluid flows in pipes the convective terms of the corresponding governing equations are generally negligible. Typically, under this condition, these governing equations are efficiently discretised by the Method of Characteristics (MOC). Only in the last years the availability of very efficient and robust numerical schemes for the complete system of equations, such as recent Finite Volume Methods (FVM), has encouraged the simulation of transient fluid flows with numerical schemes different from the MOC, allowing a better representation of the physics of the phenomena. In this work, a wide and critical comparison of the capability of Method of Characteristics, Explicit Path-Conservative Finite Volume Method (DOT solver) and Semi-Implicit (SI) Staggered Finite Volume Method is presented and discussed, in terms of accuracy and efficiency. To perform the analysis in a framework that properly represents real-world engineering applications, 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 taken into account. The analyses are performed comparing numerical solutions obtained using the three models against experimental data and analytical solutions. In particular, water hammer studies in high density polyethylene (HDPE) pipes, for which laboratory data have been provided, are used as test cases. Considering the visco-elastic mechanical behaviour of plastic materials, a 3-parameter and a multiparameter linear visco-elastic rheological models are adopted and implemented in each numerical scheme. Original extensions of existing techniques for the numerical treatment of such visco-elastic models are introduced in this work for the first time. After a focused calibration of the visco-elastic parameters, the different performance of the numerical models is investigated. A comparison of the results is presented taking into account 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 visco-elasticity with respect to the damping effect related to the unsteady friction is confirmed in these contexts. Moreover, all the numerical methods show a good agreement with the experimental data and a high efficiency of the MOC in standard configuration is observed. Finally, three Riemann Problems (RP) are chosen and run to stress the numerical methods, taking into account cross-sectional changes, more flexible materials and cavitation cases. In these demanding scenarios, the weak spots of the Method of Characteristics are depicted.

Modeling blood flow in viscoelastic vessels: the 1D augmented fluid–structure interaction system

“…Recently developed Runge-Kutta schemes overcome this issue. Thus, a formally implicit finite volume discretization is adopted, applying a second-order L-stable diagonally implicit Runge-Kutta method (DIRK) to the stiff part and a second-order explicit strong-stability-preserving (SSP) method to the non-stiff terms, with the addition of the path-conservative Dumbser-Osher-Toro (DOT) Riemann solver, as applied in [7] for compressible flows in polymer tubes.…”

Transient Liquid Flow in Plastic Pipes

“…For quantitative analysis, it is essential to define a minimum of two coefficients, whose role is to mathematically determine the compliance of the simulated pressure histories with respect to the experimental tests. To the authors' knowledge, only few studies in the literature [24], [28], [36] and [37] in the field of water hammer flows are concerned on such quantitative analysis. Urbanowicz's analysis carried out in [28] was based on pressure histories in dimensional form.…”

Dam break in rectangular channels with different upstream-downstream widths

“…The intermediate time step is defined as t n+1/2 = t n + ∆t/2. A second-order FVM in space and time is chosen, as described in [27,28]. The Dumbser-Osher-Toro (DOT) approximate Riemann solver [25] is adopted to evaluate fluctuations at the cell boundaries related to the nonconservative part of the system (2).…”