Numerical Solution of Saint-Venant Equations

The general, one-dimensional Saint-Venant equations are presented for a rigid open channel of arbitrary form, not necessarily prismatic, containing a flow that may be spatially varied. The theoretical basis for the method of characteristics is reviewed and used to show that, in the general case, the speed of long-wave disturbances is given by the slope of the characteristic curves. Finite-difference schemes on a rectangular net in the x - t plane and based on the characteristic forms of the Saint-Venant equations, as well as on the direct forms, are given and examined for their stability. The von Neumann technique for stability analysis is presented in detail. Explicit numerical schemes, which are simple, but require small steps in time because of stability problems, are contrasted with implicit schemes that permit numerical solution over large time steps but require the solution of large sets of simultaneous algebraic equations at each step. The double-sweep or progonka method, an exact time- and space-saving technique for solving these (locally linearized) equations, is also given in detail.