In this blog we show how the St Venant terms are used in SWMM5 as equations, table, graphs and units. We use a QA/QC version of SWMM 5 that lists many more link, node, system and Subcatchment variables than the default SWMM 5 GUI and engine. This also applies to #InfoSWMM and any software the uses the #SWMM5 engine.
SWMM5 is using is the most advanced equations as it takes into consideration the full dynamic (St. Venant) equations and not the more simplified kinematic wave / manning equations. The manning equation only considers the uniform flow conditions which represents a situation where the gravitational force on a column of water (due to the channel slope) balances out the frictional force. The full dynamic equations contains additional factors that affect the movement of water in a conduit or channel. These include the pressure force due to variation of depth along the length of the channel and the inertial (or convective acceleration) effect due to variation of flow area along the channel length. Because of these additional terms the flow/head relation you have in uniform flow conditions can be completely different according to the configuration of his network.
How are the St Venant Terms used in SWMM5?
dq2 = Time Step * Area wtd * (Head Downstream – Head Upstream) / Link Length or
dq2 = Time Step * Area wtd * (HGL) / Link Length Qnew = (Qold – dq2 + dq3 + dq4) / ( 1 + dq1) when the force main is full dq3 and dq4 are zero and
dq4 = Time Step * Velocity * Velocity * (a2 – a1) / Link Length * Sigma the dq3 term in dynamic.c uses the current midpoint area (a function of the midpoint depth), the sigma factor and the midpoint velocity
dq3 = 2 * Velocity * ( Amid(current iteration) – Amid (last time step) * Sigma
dq1 = Time Step * RoughFactor / Rwtd^1.333 * |Velocity| The weighted area (Awtd) is used in the dq2 term of the St. Venant equation:
dq2 = Time Step * Awtd * (Head Downstream – Head Upstream) / Link Length
You can also see the QA/QC report for SWMM 5 https://www.epa.gov/water-research/storm-water-management-model-swmm#downloads
How are the St Venant Units used in #SWMM5?The new flow (Q) calculated at during each iteration of time step as
(1) Q for the new iteration = (Q at the Old Time Step – DQ2 + DQ3 + DQ4 ) / ( 1.0 + DQ1 + DQ5)
In which DQ2, DQ3 and DQ4 all have units of flow (note internally SWMM 5 has units of CFS and the flows are converted to the user units in the output file, graphs and tables of SWMM 5).
The equations and units for DQ2, DQ3 and DQ4 are:
(2) Units of DQ2 = DT * GRAVITY * aWtd * ( H2 – H1) / Length = second * feet/second^2 * feet^2 * feet / feet = feet^3/second = CFS
(3) Units of DQ3 = 2 * Velocity * ( aMid – aOld) * Sigma = feet/second * feet^2 = feet^3/second = CFS
(4) Units of DQ4 = DT * Velocity * Velocity * ( aMid – aOld) * Sigma / Length = second * feet/second * feet/second * feet^2 / feet = feet^3/second = CFS
The equations and units for DQ1 and DQ5 are:
(5) Units of DQ1 = DT * GRAVITY * (n/PHI)^2 * Velocity / Hydraulic Radius^1.333 = second * feet/second^2 * second^2 * feet^1/3 * feet/second / feet^1.33 = Dimensionless
(6) Units of DQ5 = K * Q / Area / 2 / Length * DT = feet^3/second * 1/feet^2 * 1/feet * second = Dimensionless
|Figure 1. St Venant Terms in Table and Graphs for #SWMM5 for dq1, dq2, dq3, dq4, dq5, dq6|
|Figure 2. St Venant Equation in SWMM5|