Non Linear Term in the Saint Venant Equation of SWMM 5
The flow equation has six components that have to be in balance at each time step:
1. The unsteady flow term or dQ/dt
2. The friction loss term (normally based on Manning's equation except for full force mains),
3. The bed slope term or dz/dx
4. The water surface slope term or dy/dx,
5. The non linear term or d(Q^2/A)/dx and
6. The entrance, exit and other loss terms.
All of these terms have to add up to zero at each time step. If the water surface slope becomes zero or negative then the only way the equation can be balanced is for the flow to decrease. If the spike is due to a change in the downstream head versus the upstream head then typically you will a dip in the flow graph as the water surface slope term becomes flat or negative, followed by a rise in the flow as the upstream head increases versus the downstream head.
You get more than the normal flow based on the head difference because in addition to the head difference you also get a push from the non linear terms or dq3 and dq4 in this graph.
Comments
The QA/QC report from the EPA Web Sitehas further information on the SWMM5 solution and of course the open source SWMM 5 code can be viewed for free
map of the code is shown here http://swmm2000.com/page/swmm5-code-mindmaps
//-----------------------------------------------------------------------------
// dynwave.c
//
// Project: EPA SWMM5
// Version: 5.0
// Date: 3/11/08 (5.0.013)
// 1/21/09 (5.0.014)
// 4/10/09 (5.0.015)
// 6/22/09 (5.0.016)
// 10/7/09 (5.0.017)
// 11/18/09 (5.0.018)
// 07/30/10 (5.0.019)
// 04/20/11 (5.0.022)
// Author: L. Rossman
// R. Dickinson
//
// Dynamic wave flow routing functions.
//
// This module solves the dynamic wave flow routing equations using
// Picard Iterations (i.e., a method of successive approximations)
// to solve the explicit form of the continuity and momentum equations
// for conduits.
Downloads
thank you very muh Robert! what i mean to know is how to solve the St. Venant Terms,not the explanation of the four St. Venant Terms. I want to know the method of Picard Algorithm in swmm5 and how to difference the St. Venant Terms. thank you!