Non Linear Term in the Saint Venant Equation of SWMM 5

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.

If you use the Keep Inertial term option then you use all 5 terms at all times, if you use the Dampen option you use all five terms until the Froude number of the link greater than 1 and then you do not use terms dq3 and dq4, if you use the Ignore option then the dq3 and dq4 terms are not used at all for the calculation of the flow in the link.

Comments

Post a Comment

Popular posts from this blog

Orifice and Weir flow calculations

Note: Orifice and Weir Flow Computations The orifice flow calculation proceeds as follows: 1. Initially and whenever the setting (i.e., the fraction opened) changes, flow coefficients for both orifice and weir behavior are computed as follows: a. For side orifices: Define Hcrit = h/2 where h is the opening height. b. For bottom orifices: i. For a circular orifice, compute area over length (i.e., circumference) as AL = h /4. ii. For a rectangular orifice compute AL = h*w/(2*(h+w)) where w is the opening width. iii. Compute Hcrit = Cd*AL/0.414 where Cd is the orifice discharge coefficient. At step 1b, the critical head for the bottom orifice, where orifice flow turns into weir flow, is found by equating the result of the orifice equation to that of the weir equation Cd*Area*sqrt(2g)*sqrt(Hcrit) = Cw*Length*sqrt(Hcrit)*Hcrit or Hcrit = (Cd * Area) / (Cw/sqrt(2g) * Length) The value of Cw/sqrt(2g) for a sharp crested weir is 0.414.
Read more

SWMM 5 Error Messages

Engine Error Number Description ERROR 101: memory allocation error.                                                                                             ERROR 103: cannot solve KW equations for Link                                                                                   ERROR 105: cannot open ODE solver.                                                                                              ERROR 107: cannot compute a valid time step.                                                                                     ERROR 108: ambiguous outlet ID name for Subcatchment                                                                            ERROR 109: invalid parameter values for Aquifer                                                                                  ERROR 110: ground elevation is below water table for Subcatchment                                                         
Read more

EMC Washoff in SWMM5

Image
There are four steps to using EMC concentrations in your network: 1. Define your pollutant by adding a pollutant using the Data=>Quality=>Pollutant command : 2. Define the Land Use by using the Data=>Land Uses command or the Land Use Editor: 3. Define Buildup to be None by clicking on the None Tab: 4. Define the EMC Washoff concentration by clicking on the Washoff Tab: This is where you would add a GW concentration More: http://www.swmm2000.com/
Read more