Explicit Iteration Hydraulic Computation and Implicit Time Step Hydraulic Computations in SWMM 5

Explicit Iteration Hydraulic Computation and Implicit Time Step Hydraulic Computations in SWMM 5
The dynamic wave solution in SWMM 5 and InfoSWMM uses an interlocking explicit iteration in an overall time step implicit calculation of the node depths and link flows. The node depths and links flows are computed based on the last iteration values for 2 to 4 iterations. It explicitly calculates the node depths and link flows at each iteration based on the last iteration value but uses an underrelaxation or Gauss-Seidel method averaging to compute the final new iteration values. The underrelaxation method averages the last iteration value and the current iteration estimate to calculate the final iteration node depth and links flow.
At each time step this is the sequence of computations:
1. The time step is based on the Courant–Friedrichs–Lewy (CFL) condition for the most restrictive link (the CFL condition is based on the link length, current depth and current velocity),
2.Link flows and node depths are calculated at the 1st iteration based on the last time steps flows and depths,
3. The new iterations link flow and node depths are calculated by averaging the node depths and link flows found in step 2 with the original iteration node depths and link flows,
4. The iteration process continues for at least 2 iterations for all nodes and links,
5. The iteration process may stop after 2 iterations for those links in which the upstream and downstream node depths have converged,
6. The upstream or downstream node is considered converged when the absolute depth difference between iterations is less than 0.005 feet,
7. The maximum number of iterations for all nodes and links is set to 4.
8. Once the number of iterations has reached 4 or all nodes and links are converged the time step calculations are considered finished and the program moves on to a new time step.
9. The image shown below shows how the node continuity equation is solved for each iteration but the link combined continuity and momentum equation is calculated for either 2, 3 or 4 iterations depending on the upstream and downstream node convergence of 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