Note: The Keep and Dampen options and their effect on the four main terms of the St Venant equation.
Note: The Keep and Dampen opt
The four terms are are used in the new flow for a time step of Qnew:
Qnew = (Qold – dq2 + dq3 + dq4) / ( 1 + dq1)
when the force main or gravity 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
where Sigma is a function of the Froude Number and the Keep, Dampen and Ignore Inertial Term Options. Keep sets Sigma to 1 always and Dampen set Sigma based on the Froude number, Ignore sets Sigma to 0 all of the time during the simulation
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 or
dq2 = Time Step * Awtd * (Head Downstream – Head Upstream) / Link Length
Normally, dq1 (Friction Loss / Maroon in the Graph) balances dq2 (Water Surface Slope Term or Green in the Graph) but often for links with a large difference between upstream and downstream depths dq4 (Red in the Graph) can have a significant value. If dq4 or dq3 are important then the depth of water to increases to pass the same flow using the Keep option over the Ignore. If you have a link with a Froude number near or over 1.0 (Supercritical) then using Keep or Dampen for the Options may result in depth differences. The effect of K
Post a Comment