A simplified flood routing model: variable prameter muskingum (VPM)

Summary Flood routing methods are numerical methods for estimating the movement of a flood wave along a channel reach, on the basis of the knowledge of the discharge hydrograph at the upstream end and of the hydraulic characteristics of the reach and, usually, in the hypothesis that no perturbation is coming from downstream (free boundary condition). The flood routing method wich is proposed is similar to the Muskingum one, but with variable and hydraulic parameters; it is able to estimate water levels too; is effective even if kinetic terms are not completely negligible; take advantage of the insignificance of the downstream condition and make it possible to obtain results starting upstream and proceeding downstream; for simplicity's sake, take advantage of the fact that the discharge loop of normal flood waves is quite small. Obtained results are much better that those obtainable from constant parameters methods and indeed, if the flood loop is less that 10%, very similar to those obtainable from more complex and time consuming models.

---

Sommario I metodi di flood routing studiano la propagazione di un'onda di piena lungo un tratto di un corso d'acqua, assegnato l'andamento temporale della portata nella sezione di monte e le caratteristiche dell'alveo, e usualmente nell'ipotesi di assenza di perturbazioni provenienti da valle (condizione di valle passiva). Viene qui proposto un procedimento di flood routing, formalmente simile ad un Muskingum ma con i parametri variabili e calcolati per via idraulica; idoneo a stimare anche i livelli idrici; valido anche se i termini cinetici non sono del tutto trascurabili; che sfrutta l'irrilevanza della condizione di valle procedendo a cascata da monte a valle; che sfrutta, a vantaggio della semplicità, il fatto che per le normali onde di piena dei corsi d'acqua il cappio di portata è di dimensioni modeste. I risultati ottenuti sono molto migliori di quelli ottenibili con metodi a parametri costanti e, almeno per i casi in cui il cappio relativo è inferiore al 10%, paragonabili a quelli ottenuti con metodi molto più complessi ed onerosi.

---

List of symbols x, t channel distance, starting upstream; time - z water surface height above datum - Q volumetric rate of discharge - P(x, z) steady rating curve - q=Q–P flood loop - g acceleration of gravity - A, B cross section wetted area and free surface width - I, S water surface slope and friction slope - c kinematic wave velocity - F Froude number - L,L ₀,L ₁,L ₂,L ₃ characteristic lengths of the channel - T=L/c characteristic time of the channel - D diffusion - p, l time and space steps - K, X Muskingum parameters - C ₁,C ₂,C ₃,C ₄ Muskingum coefficients - f _x=f/x,f _t=f/t etc. for the partial derivatives Paper presented at the First Italian Meeting of Computational Mechanics held in Milan, June 24–26, 1986.