Accuracy of kinematic wave and diffusion wave approximations for space-independent flows with lateral inflow neglected in the momentum equation

Error equations for the kinematic wave and diffusion wave approximations with lateral inflow neglected in the momentum equation are derived under simplified conditions for space-independent flows. These equations specify error as a function of time in the flow hydrograph. The kinematic wave, diffusion wave and dynamic wave solutions are parameterized through a dimensionless parameter γ which is dependent on the initial conditions. This parameter reflects the effect of initial flow depth, channel-bed slope, lateral inflow and channel roughness when the initial condition is non-vanishing; and it reflects the effect of bed slope, channel roughness and acceleration due to gravity when the initial condition is vanishing. The error equations are found to be the Riccati equation. The structure of the error equations in the case when the momentum equation neglects lateral inflow is different from that when the lateral inflow is included.     

Accuracy of kinematic wave and diffusion wave approximations for space-independent flows on infiltrating surfaces

Error equations for the kinematic wave and diffusion wave approximations were derived under simplified conditions for space-independent flows occurring on infiltrating planes or channels. These equations specify error as a function of time in the flow hydrograph. The kinematic wave, diffusion wave and dynamic wave solutions were parameterized through a dimensionless parameter γ which is dependent on the initial conditions. This parameter reflects the effect of initial flow depth, channel bed slope, lateral inflow and channel roughness when the initial condition is non-vanishing; it reflects the effect of bed slope, channel roughness, lateral inflow and infiltration when the initial condition is vanishing. The error equations were found to be the Riccati equation.     

ERRORS OF KINEMATIC-WAVE AND DIFFUSION-WAVE APPROXIMATIONS FOR SPACE-INDEPENDENT FLOWS ON INFILTRATING SURFACES

Error equations for the kinematic-wave and diffusion-wave approximations were derived under simplified conditions for space-independent flows occurring on infiltrating planes or channels. These equations specify error as a function of time in the flow hydrograph. The kinematic-wave, diffusion wave and dynamic-wave solutions were parameterized through a dimensionless parameter γ which is dependent on the initial conditions. This parameter reflects the effect of initial flow depth, channel-bed slope, lateral inflow and channel roughness when the initial condition is non-vanishing; and it reflects the effect of bed slope, channel roughness, lateral inflow and infiltration when the initial condition is vanishing. The error equations were found to be the Riccati equation.     

A practical method for analysis of river waves and for kinematic wave routing in natural channel networks

The behaviour of river waves is described using a simplified dimensionless form of the momentum equation in conjunction with the continuity equation. Three dimensionless parameters were derived based on a quantitative linear analysis. These parameters, which depend on the Froude number of the steady uniform flow and the geometric characteristics of the river, permit quantification of the influence of inertia and pressure in the momentum equation. It was found that dynamic and diffusion waves occur mainly on gentle channel slopes and the transition between them is characterized by the Froude number. On the other hand, the kinematic wave has a wide range of applications. If the channel slope is greater than 1%, the kinematic wave is particularly suitable for describing the hydraulics of flow. Since slopes in natural channel networks are often greater than 1%, an analytical solution of the linearized kinematic wave equation with lateral inflow uniformly distributed along the channel is desirable and was therefore derived. The analytical solution was then implemented in a channel routing module of an existing simple rainfall–runoff model. The results obtained using the analytical solution compared well with those obtained from a non‐linear kinematic wave model. Copyright © 2000 John Wiley & Sons, Ltd.     

ACCURACY OF KINEMATIC WAVE AND DIFFUSION WAVE APPROXIMATIONS FOR TIME-INDEPENDENT FLOW WITH MOMENTUM EXCHANGE INCLUDED

Errors in the kinematic wave and diffusion wave approximation for time-independent (or steady-state) cases of channel flow with momentum exchange included were derived for three types of boundary conditions: zero flow at the upstream end, and critical flow depth and zero depth gradient at the downstream end. The diffusion wave approximation was found to be in excellent agreement with the dynamic wave approximation, with errors of less than 1% for KF²₀≥7·5 and up to 12% for KF²₀≤0·75 for the upstream boundary condition of zero discharge and finite depth, where K is the kinematic wave number and F₀ is the Froude number. The kinematic wave approximation was reasonably accurate except at the channel boundaries and for small values of KF²₀ (≤1). The accuracy of these approximations was significantly influenced by the downstream boundary condition both in terms of the error magnitude and the segment of the channel reach for which these approximations would be applicable. © 1997 by John Wiley & Sons, Ltd.     

Errors of kinematic wave and diffusion wave approximations for time‐independent flows with infiltration and momentum exchange included

Error equations for kinematic wave and diffusion wave approximations were derived for time‐independent flows on infiltrating planes and channels under one upstream boundary and two downstream boundary conditions: zero flow at the upstream boundary, and critical flow depth and zero depth gradient at the downstream boundary. These equations specify error in the flow hydrograph as a function of space. The diffusion wave approximation was found to be in excellent agreement with the dynamic wave approximation, with errors below 2% for values of KF (e.g. KF ≥ 7·5), where K is the kinematic wave number and F is the Froude number. Even for small values of KF (e.g. KF = 2·5), the errors were typically less than 3%. The accuracy of the diffusive approximation was greatly influenced by the downstream boundary condition. For critical flow depth downstream boundary condition, the error of the kinematic wave approximation was found to be less than 10% for KF ≥ 7·5 and greater than 20% for smaller values of KF. This error increased with strong downstream boundary control. The analytical solution of the diffusion wave approximation is adequate only for small values of K. Copyright © 2005 John Wiley & Sons, Ltd.     

Accuracy of kinematic wave and diffusion wave approximations for time-independent flows

Errors in the kinematic wave and diffusion wave approximations for time-independent (or steady-state) cases of channel flow were derived for three types of boundary conditions: zero flow at the upstream end, and critical flow depth and zero depth gradient at the downstream end. The diffusion wave approximation was found to be in excellent agreement with the dynamic wave approximation, with errors in the range 1–2% for values of KF (? 7.5), where K is the kinematic wave number and F₀ is the Froude number. Even for small values of KF (e.g. KF²₀ = 0.75), the errors were typically less than 15%. The accuracy of the diffusion wave approximation was greatly influenced by the downstream boundary condition. The error of the kinematic wave approximation was found to be less than 13% in the region 0.1 ? x ? 0.95 for KF = 7.5 and was greater than 30% for smaller values of KF (? 0.75). This error increased with strong downstream boundary control.     

Accuracy of the kinematic wave and diffusion wave approximations for time-independent flows with infiltration included

Errors in the kinematic wave and diffusion wave approximation for time-independent (or steady-state) cases of channel flow with infiltration were derived for three types of boundary conditions: zero flow at the upstream end, and critical flow depth and zero depth gradient at the downstream end. The diffusion wave approximation was found to be in excellent agreement with the dynamic wave approximation, with errors of less than 1·4% for KF²₀≥7·5, and up to 14% for KF²₀≤0·75 for the upstream boundary condition of zero discharge and finite depth, where K is the kinematic wave number and F₀ is the Froude number. The kinematic wave approximation was reasonably accurate except at the channel boundaries and for small values of KF²₀ (≤1). The accuracy of these approximations was significantly influenced by the downstream boundary, both in terms of the magnitude of the error and the segment of the channel reach for which these approximations would be applicable. © 1998 John Wiley & Sons, Ltd.     

Flood routing based on diffusion wave equation using mixing cell method

A one-dimensional non-linear diffusion wave equation is derived from the Saint Venant equations with neglect of the inertia terms. This non-linear equation has no general analytical solution. Numerical schemes are therefore employed to discretize the space and time axes and convert the differential equation to difference form. In this study, the mixing cell method is used to convert the diffusion wave equation to difference form, in which the difference term can be eliminated by selecting an optimal space step size Δx when time step size Δt is given. When the time step size Δt→0, the space step size Δx=Q/(2S₀BC]_k) where Q is discharge, S₀ is bed slope, B is channel width and C_k is kinematic wave celerity, which is the same as the characteristic length proposed by Kalinin and Milyukov. The results of application to two cases show that the mixing cell and linear channel flow routing methods produce hydrographs that are in agreement with the observed flood hydrographs. © 1997 John Wiley & Sons, Ltd.     

Stability and accuracy of two-dimensional kinematic wave overland flow modeling

A two-dimensional finite element based overland flow model was developed and used to study the accuracy and stability of three numerical schemes and watershed parameter aggregation error. The conventional consistent finite element scheme results in oscillations for certain time step ranges. The lumped and the upwind finite element schemes are tested as alternatives to the consistent scheme. The upwind scheme did not improve on the stability or the accuracy of the solution, while the lumped scheme provided stable and accurate solutions for time steps twice the size of time steps needed for the consistent scheme. A new accuracy based dynamic time step estimate for the two-dimensional overland flow kinematic wave solution is developed for the lumped scheme. The newly developed dynamic time step estimates are functions of the mesh size, and time of concentration of the watershed hydrograph. Due to lack of analytical solutions, the time step was developed by comparing numerical solutions of various levels of discretization to a reference solution using a very fine mesh and a very small time step. The time step criteria were tested on a different set of problems and proved to be adequate for accurate and stable solutions. A sensitivity analysis for the watershed slope, Manning’s roughness coefficient and excess rainfall rate was conducted in order to test the effect of parameter aggregation on the stability and accuracy of the solution. The results of this analysis show that aggregation of the slope data resulted in the highest error. The roughness coefficient had a smaller effect on the solution while the rainfall intensity did not show any significant effect on the flow rate solution for the range of rainfall intensity used. This work pioneers the challenge of providing guidelines for accurate and stable numerical solutions of the two-dimensional kinematic wave equations for overland flow.     

Scale invariance and self‐similarity in kinematic wave overland flow in space and time

Fractals are famous for their self‐similar nature at different spatial scales. Similar to fractals, solutions of scale invariant processes are self‐similar at different space–time scales. This unique property of scale‐invariant processes can be utilized to translate the solution of the processes at a much larger or smaller space–time scale (domain) based on the solution calculated on the original space–time scale. This study investigates scale invariance conditions of kinematic wave overland flow process in one‐parameter Lie group of point transformations framework. Scaling (stretching) transformation is one of the one‐parameter Lie group of point transformations and it has a unique importance among the other transformations, as it leads to the scale invariance or scale dependence of a process. Scale invariance of a process yields a self‐similar solution at different space–time scales. However, the conditions for the process to be scale invariant usually dictate various relationships between the scaling coefficients of the dependent and independent variables of the process. Therefore, the scale invariance of a process does not assure a self‐similar solution at any arbitrary space and time scale. The kinematic wave overland flow process is modelled mathematically as initial‐boundary value problem. The conditions to be satisfied by the system of governing equations as well as the initial and boundary conditions of the kinematic wave overland flow process are established in order for the process to be scale invariant. Also, self‐similarity of the solution of the kinematic wave overland flow under the established invariance conditions is demonstrated by various numerical example problems. Copyright © 2011 John Wiley & Sons, Ltd.     

A MATHEMATICAL QUASI-2D OVERLAND FLOW MODEL

1 INTRODUCTIONThe prediction of future impacts on terrestrial ecosystems by atmospheric, climatic and land-usechanges is the aim of watershed management. Meeting these requirements scientists, managers and policymakers try to achieve the sustainable management of the vitally important resources of watersheds due toan integrated ecosystem approach at the catchment scale. As composite landscapes often have a highdegree of contingency between its elements, the transport over these landscape s…     

A primitive pseudo wave equation formulation for solving the harmonic shallow water equations

A finite element method formulation for solving the harmonic shallow water equations in their primitive or unmodified form is developed and analysed. The scheme, referred to as the Primitive Pseudo Wave Equation Formulation (PPWE), involves developing a weak weighted residual form of the continuity equation and furthermore forming a pseudo wave equation by substituting the discretized form of the momentum equation into the discretized form of the continuity equation. The final set of equations to be solved, the pseudo wave equation and the primitive momentum equations, deceptively resemble the discretized equations of the wave equation formulation of Lynch and Gray. Despite this resemblance, Fourier analysis indicates that the PPWE scheme is still fundamentally primitive.However, application of the PPWE scheme to a set of stringent test problems results in very good solutions with well controlled nodal oscillations. It is shown that this low degree of spurious oscillations is due to the treatment of the boundary conditions such that elevation is taken as an essential condition and normal flux is taken as a natural condition. This particular boundary condition treatment is suggested by the formation of the pseudo wave equation. Furthermore, even though the equation re-arrangement does not in itself effect the solutions, it does make the scheme much more efficient.     

Implication of the flow resistance formulae on the prediction of flood wave propagation

ABSTRACT

This study examines the difference in the predictions of flood wave propagation in open channels depending on the flow resistance formulae, such as the Chézy and Manning’s equation. The celerity and diffusion coefficient are functions of the channel geometry, slope, roughness as well as the resistance formulae. The results suggest that substituting the Chézy equation with Manning’s equation results in different characteristics of flood propagation, which are consistent regardless of the cross-sectional geometry except for a circular cross-section: increasing celerity and decreasing diffusion coefficient. The celerity is more sensitive to the selection of resistance formulae than the diffusion coefficient. Geometry has a greater effect on the celerity and diffusion coefficient, and consequently on the resulting hydrographs. Manning’s equation results in a larger difference in celerity and diffusion coefficient compared to Chézy equation regardless of the water depth. Overall, this study shows that the selection of resistance formulae is important in terms of the resulting hydrographs and peak flow.

EDITOR Z.W. Kundzewicz ASSOCIATE EDITOR not assigned     

ANALYTICAL HAYAMI SOLUTION FOR THE DIFFUSIVE WAVE FLOOD ROUTING PROBLEM WITH LATERAL INFLOW

The diffusive wave equation is generally used in flood routing in rivers. The two parameters of the equation, celerity and diffusivity, are usually taken as functions of the discharge. If these two parameters can be assumed to be constant without lateral inflow, the diffusive wave equation may have an analytical solution: the Hayami model. A general analytical method, based on ‘Hayami’s hypothesis, is developed here which resolves the diffusive wave flood routing equation with lateral inflow or outflow uniformly distributed over a channel reach. Flood routing parameters are then identified using observed inflow and outflow and the Hayami model used to simulate outflow. Two examples are discussed. Firstly, the prediction of the hydrograph at a downstream section on the basis of a knowledge of the hydrograph at an upstream section and the lateral inflow. The second example concerns lateral inflow identification between an upstream and a downstream section on the basis of a knowledge of hydrographs at the upstream and downstream sections. The new general Hayami model was applied to flood routing simulation and for lateral inflow identification of the River Allier in France. The major advantages of the method relate to computer simulation, real-time forecasting and control applications in examples where numerical instabilities, in the solution of the partial differential equations must be avoided.     

MÉTHODE D'ÉTUDE DES NAPPES LIBRES DANS LES AQUIFÈRES ANISOTROPES

Abstract

The necessary and sufficient conditions for non-zero phase shift and non-zero attenuation in linear flood routing can be derived from the continuity equation alone and are found to depend on the existence of an imaginary part in the expression for frequency or in the expression for wave number. It is shown that in linear flood routing the phase lag between flow rate and area of flow is directly related to the attenuation per unit wave length. The effects of using various forms of the momentum equation, in addition to the continuity equation, are exemplified by deriving analytical expressions in terms of the frequency, both for attenuation per unit channel length and for phase shift, for the kinematic wave, the general diffusion analogy, and the complete St. Venant equation.     

A fuzzy dynamic wave routing model

In this paper a fuzzy dynamic wave routing model (FDWRM) for unsteady flow simulation in open channels is presented. The continuity equation of the dynamic wave routing model is preserved in its original form while the momentum equation is replaced by a fuzzy rule based model which is developed on the principle that during unsteady flow the disturbances in the form of discontinuities in the gradient of the physical parameters will propagate along the characteristics with a velocity equal to that of velocity of the shallow water wave. The model gets rid off the assumptions associated with the momentum equation by replacing it with the fuzzy rule based model. It overcomes the necessity of calculating friction slope (S_f) in flow routing and hence the associated uncertainties are eliminated. The robustness of the fuzzy rule based model enables the FDWRM to march the solution even in regions where the aforementioned assumptions are violated. Also the model can be used for flow routing in curved channels. When the model is applied to hypothetical flood routing problems in a river it is observed that the results are comparable to those of an implicit numerical model (INM) which solves the dynamic wave equations using an implicit numerical scheme. The model is also applied to a real case of flow routing in a field canal. The results match well with the measured data and the model performs better than the INM. Copyright © 2007 John Wiley & Sons, Ltd.     

Algorithms for solving the diffusive wave flood routing equation

Generally, the diffusive wave equation, obtained by neglecting the acceleration terms in the Saint-Venant equations, is used in flood routing in rivers. Methods based on the finite-difference discretization techniques are often used to calculate discharges at each time step. A modified form of the diffusive wave equation has been developed and new resolution algorithms proposed which are better adapted to flood routing along a complex river network. The two parameters of the equation, celerity and diffusivity, can then be taken as functions of the discharge. The resolution algorithm allows the use of any distribution of lateral inflow in space and time. The accuracy of the new algorithms were compared with a traditional algorithm by numerical experimentation. Special attention was given to the instability caused by the inflow signal which constitutes the upstream boundary condition. For the fully diffusive wave flood routing problem, all three algorithms tested gave good results. The results also indicate that the efficiency of the new algorithms could be significantly improved if the position of the x-axis is modified by rotation. The new algorithms were applied to flood routing simulation over the Gardon d'Anduze catchment (542 km²) in southern France.     

Diffusive wave solutions for open channel flows with uniform and concentrated lateral inflow

The diffusive wave equation with inhomogeneous terms representing hydraulics with uniform or concentrated lateral inflow into a river is theoretically investigated in the current paper. All the solutions have been systematically expressed in a unified form in terms of response function or so called K-function. The integration of K-function obtained by using Laplace transform becomes S-function, which is examined in detail to improve the understanding of flood routing characters. The backwater effects usually resulting in the discharge reductions and water surface elevations upstream due to both the downstream boundary and lateral inflow are analyzed. With a pulse discharge in upstream boundary inflow, downstream boundary outflow and lateral inflow respectively, hydrographs of a channel are routed by using the S-functions. Moreover, the comparisons of hydrographs in infinite, semi-infinite and finite channels are pursued to exhibit the different backwater effects due to a concentrated lateral inflow for various channel types.