An accurate time integration method for simplified overland flow models

An accurate time integration method for the diffusion-wave and kinematic-wave approximated models for the overland flow is proposed. The discretization of the first- and second-order spatial derivatives in the basic equation is obtained by using the second-order Lax–Wendroff and the three-point centred finite difference schemes, respectively. For the solution in time, the system of ordinary differential equations, obtained by the linearization of the celerity and of the hydraulic diffusivity by Taylor series expansions, is integrated analytically. The stability and the numerical dissipation and dispersion are investigated by the Fourier analysis. A proper Courant number, and the corresponding time step for the numerical simulations can be established. In addition, the proposed diffusion-wave and kinematic-wave models are straightforwardly extended to the two-dimensional flow. Test cases for both one- and two-dimensional problems, compare the solutions of the diffusion-wave and kinematic-wave models with analytical solutions, with experimental results and with numerical solutions obtained by the Saint–Venant equations. These simulations show that the proposed numerical–analytical models accurately predict the overland flow for several situations, in particular for unsteady rainfall rate and for spatial variations of the surface roughness.