Extended Period Dynamic Simulation (Unsteady Flow) in InfoSewer SA
InfoSewer SA tracks the movement of wastewater flowing through the network over an extended period of time under varying wastewater loading and operating conditions. The extended period simulation (EPS) model implemented in InfoSewer SA is a fully dynamic (unsteady) model and is predicated on solving a simplified form of the full 1D Saint-Venant equations neglecting local acceleration.
The Saint-Venant equations or full dynamic wave equations for open channel flow routing consist of the conservation of momentum equation and the equation of continuity. The momentum equation is:
The continuity (mass conservation) equation is:
x = distance along the pipe (longitudinal direction of sewer)
A = flow cross sectional area normal to x
y = coordinate direction normal to x on a vertical plane
d = depth of flow of the cross section, measured along y direction
Q = discharge through A
V = cross sectional average velocity along x direction
S0 = pipe slope, equal to sin θ
θ = angle between sewer bottom and horizontal plane
Sf = friction slope
g = gravitational acceleration
t = time
β = Boussinesqmomentum flux correction coefficient for velocity distribution
These complete unsteady flow equations (momentum together with continuity) along with appropriate initial and boundary conditions are rather tedious and computationally expensive to solve, especially for large sewer collections systems. As a result, acceptable simplifications and improved solution methods have been proposed including non-inertial, kinematic wave and dynamic wave simplifications. Hydraulically, the dynamic wave approach is the most accurate model among the approximations. The Muskingum-Cunge explicit diffusion wave dynamic flow routing model, obtained by neglecting local acceleration term in the momentum equation, is the most commonly used dynamic wave model.
In InfoSewer SAH20Map SewerPro, unsteady open channel (free surface) flow is simulated using Muskingum-Cunge technique whereas pressurized flow in any pipe is modeled assuming the pipe is flowing full and the energy equation is applied to the entire pipe section.
Here c is the dynamic wave celerity and B is the top width at normal depth for discharge Q. This highly efficient and accurate flow routing algorithm is used by InfoSewer SA to track the spatial and temporal variation of flows throughout the collection system.
In this method (a.k.a., one sweep explicit solution method), the network flow dynamic equations are formulated by using an explicit finite difference scheme such that the flow depth, discharge, or velocity at a given location and the current time can be solved explicitly from the known information at the previous locations at the same time level, as well as known information at the previous time level. Thus, the solution is obtained segment by segment, pipe by pipe, over a given time interval for the entire sewer network before progressing to the next interval for another sweep of individual solutions of the network flow equations for the entire network. A variable time step approach (based on the Courant number is used to minimize numerical dispersion and ensure robustness and stability of the numerical scheme. Complex flow attenuation calculations can be explicitly carried out to more accurately simulate the movement and transformation of sanitary sewer flows in the collection system.
An excellent review and comparison between simulated and observed hydrographs of the various numerical methods for solving unsteady flow in simple and compound channels was presented by Chatila (Chatila 2003). In terms of overall performance, the Muskingum unsteady solution scheme compared favorably and proved to be a simple and reliable method avoiding complicated mathematical and numerical computations for the cases considered.
Flooding at manholes and wet-wells in InfoSewer is not modeled during an extended period dynamic simulation. Instead, the flows at the flooded structures are conserved and are not lost by the occurrence of flooding at the manholes. In actual flooding situations, flows may be diverted away from the flooded structures and out of the sewer collection system. However, a surcharged pipe or manhole is generally an indication of poor hydraulic performance of the sewer system. InfoSewer SA Pro assumes that the downstream pipes of flooded manholes are flowing full.
Sanitary sewer systems are typically designed to flow less than full and have an upper pressure limit of 4 to 6 psi. Sewer systems operating under pressurized flow condition may run the risk of violating local, state, and federal health codes. The USEPA regulations would also be in violation if raw sewage were discharged into the ground, potentially affecting groundwater. For these reasons, pressurized flows in sanitary sewers not designed to sustain pressures can be dangerous and in some cases can present an unlawful activity.
Sewer pipes can flow full with water under pressure, which is often known as surcharge flow. Surcharge flow occurs in under-designed pipes (or under extreme flows) when the flow rate Q exceeds the full pipe capacity Qf .
Flow conditions are unstable at the transition between open-channel (free surface) flow and full pipe flow. A wave or surge can induce full flow in the pipe in the unstable range. Surcharge in sewer pipes is modeled in InfoSewer SA using energy and continuity principles. The energy equation between sections 1 and 2 in a pipe can be written as:
Here z denotes the invert elevation; d represents the water depth; and HL designates the head loss between sections 1 and 2. The energy equation is used to determine the difference in hydraulic grade line elevation (which is added at the upstream manholes) needed to pass downstream flows under the surcharge condition.
The procedure for analyzing surcharge in sewer pipes is illustrated using the figure below as a reference.
Assuming that pipe 4 (between manholes 4 and 5) is under-designed, Q4 will exceed its full flow capacity and the hydraulic grade line at manhole 4 will increase based on energy consideration to allow Q4 to pass through pipe 4 (note that water always flows from higher to lower energy) as continuity must be satisfied. This forces the hydraulic grade line at manhole 3 to increase in order for Q3 to pass through pipe 3. The procedure continues upstream until the slope of the energy grade line needed to transport the flow allows open-channel flow condition to occur in the pipe. The projected hydraulic grade line will then intersect the uniform water surface flow to complete the backwater curve.
The energy equation is also used to model the flow in siphons, which can occur in adverse pipes. InfoSewer assumes that the siphon flows full, with a continuous liquid column throughout it. This image shows how an adverse sloped link InfoSewer looks
When a flow hydrograph is injected and propagates downstream in sewer pipes the bulk of the water will normally travel slower than its induced disturbance or wave. That is, if the water is injected with a tracer then the tracer lags behind the disturbance. The speed of the disturbance depends on parameters such as depth, width and flow velocity. This disturbance will tend to flatten, or spread out, the peak flow in the downstream direction along the sewer pipes.
Flow attenuation in a sewer system is defined as the process of reducing the peak flow rate by redistributing the same volume of flow over a longer period of time as a result of friction (resistance), internal storage and diffusion along the sewer pipes. InfoSewer SA uses the distributed Muskingum-Cunge flow routing method based on diffusion analogy, which is capable of accurately predicting hydrograph attenuation or peak flow damping effects (peak subsidence). The method is attractive since the routing parameters can be directly calculated as a function of pipe and flow properties, is applicable for a wide range of flow conditions, and does not require calibration or any iterative scheme. The Muskingum coefficients are derived from the pipe diameter, length, discharge, dynamic wave celerity, and slope of the flow. The magnitude of attenuation depends on parameters such as the peak discharge, the curvature of the hydrograph, and the width of flow. An example of flow attenuation process as a hydrograph is routed through a sewer system is illustrated in the figure below.
Hydrograph Aggregation/Flow Accumulation
Proper aggregation of multiple hydrographs with distinct time steps is essential in a sewer collection system as the flows are routed in both time and space. Aggregation normally occurs when laterals are merging around manholes and wet-wells. This can create offset of time-steps, which can affect accurate determination of flow peaks and volumes. InfoSewer SA utilizes a highly accurate dynamic hydrograph aggregation method that allows preservation of both flow peaks and flow volumes when multiple hydrographs with different time steps are added. The method is Lagrangian in nature and tracks the hydrograph ordinates as they are transported along the sewer pipes and mix together at manholes and wet-wells. A variable time step is used to minimize numerical dispersion, enhance stability, and maximize computational efficiency. See the User Guide for more information on Extended Period Simulations.