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 value of dq4 increases when there is a significant difference in the cross sectional  area of the downstream end of the link and the upstream end of the link.  In this  example, the downstream storage node causes a backflow in the link (Figure 1).   The flow may look unstable in the link  flow time series but the change in flow is simply due to the water sloshing back and forth.  There is no continuity error as the term dq4 keeps the water in the link  in balance (Figure 2).