SWMM5 Link Upstream Weighting

Purpose: The purpose of this note is to explain a significant dynamic wave routing difference between EPA SWMM 5.0.013 and EPA SWMM 5.0.011 and before. A few people have detected a difference. The previous solution(s) would use only the midpoint area (Amid) and hydraulic radius (Rmid) in the dynamic wave solution. The new solution will use a slider or linear combination of the midpoint area (Amid) and hydraulic radius (Rmid) and the upstream cross sectional area (A1) and hydraulic radius (R1). The slider is based on the Froude number in the link. The change involves the A and R link spacing in the two dominant terms of the St. Venant Equation:


The new method is a linear combination or slider that weights the value of A and R in the St. Venant Equation based on the value of rho (), or



where, Rho () is a function of the Froude number. The effect of this addition is that as the Froude number increases from 0.5 to 1.0 and beyond the area and hydraulic radius used as the pivot point in the St. Venant equation moves from the midpoint of the link to the upstream end of the link. When the Froude number is above 1.0 the St. Venant and Normal Flow equation both use the same cross sectional area and hydraulic radius which makes for a more stable model.

Just for reference, the equation for Qnorm or the Manning's Equation flow is



The equations for the calculation of Rho () as a function of the Froude Number (Fr) are:



If ALL of the follow conditions are true Rho ()is calculated:
  • the pipe is not full,
  • h1 >= h2, and
  • qLast > 0.

where,
h1 is the head at the upstream end of the link,
h2 is the head at the downstream end of the link and
qLast is the last flow value in the link.

If any of these conditions are true then rho = 1.0 and the value of A and R are the values Amid and Rmid, respectively.
The next graph shows the relationship between Rho and the Froude Number.



The value of Awtd and Rwtd move from the midpoint of the link to the upstream end of the link as the Froude number increases from 0.5 to 1.0.




Conclusion: This change should make the solution more stable because there is no longer an oscillation between the St. Venant Equation A and R and the Normal Flow Equation A and R.
Zemanta Pixie

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