SWMM 5 Controls in Matrix Form

Subject:  SWMM 5 Controls in Matrix Form

Comments

Post a Comment

Popular posts from this blog

SWMM 5 Error Messages

Engine Error Number Description ERROR 101: memory allocation error.                                                                                             ERROR 103: cannot solve KW equations for Link                                                    ...
Read more

More St Venant Equations in #SWMM5

Image
This blog shows the relationship between the terms dq1, dq2, dq3 and dq4 in the SWMM5 code and the St. Venant Partial Differential Equations. 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 Qnew = (Qold – dq2) / ( 1 + dq1) The dq4 term in dynamic.c uses the area upstream (a1) and area downstream (a2), the midpoint velocity, the sigma factor (a function of the link Froude number), the link length and the time step or 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...
Read more