Node Overflowing, Surface Flooding and Overland/Street Flow In 1D XPSWMM Modeling

Node Overflowing, Surface Flooding and Overland/Street Flow In 1D XPSWMM Modeling

Even if 2D modeling by XPSWMM, HEC-RAS, and InfoWorks has become popular these days, 1D modeling still holds its ground on storm water management and flooding study due to several advantages over 2D:

  1. Mainstream storm water modeling approach for the last half century and having a large user base in the engineering community
  2. Less amount of data required to establish a 1D model and smaller model file sizes
  3. Faster model run which means more time for model debugging, validation and calibration, and alternative study under a limited budget and a tight schedule
  4. Capable of generating relatively reliable results when applying proper modeling skills and reasonable assumptions
  5. Methodology and results widely acceptable by reviewing and regulating agencies

SWMM including EPA SWMM and the proprietary software XPSWMM was originally developed to model 1D flow in a storm sewer network or an open channel system and extra considerations must be given when a 1D model is to deal with node overflowing or surface flooding. Different methodologies compute HGL in an overflowing node differently which in turn would impact 1D flow in a pipe or a channel. When a storm sewer network is overwhelmed, any storm water exceeding its capacity will rise up in nodes until the water level reaches the ground surface (spill crest) to cause node overflowing and at this point XPSWMM has different ways to handle the situation:

1. Set up Spill Crest as the actual ground or top of node elevation while the Ponding is “None” (Figure 1). Any excessive water above node spill crest elevation will get lost from the system through nodes (no water loss through links in XPSWMM). This option probably should be avoided at all costs since the lost water will not be accounted for by the model. The maximum HGL at a node is its spill crest under this setting.

Figure 1

2. Set up Spill Crest to an arbitrarily high value to prevent the water from getting lost while the Ponding is “None” (Figure 2). The extremely high Spill Crest will contain the water in the system (HGL can go up as high as Spill Crest), but it should only be used sparsely in a screening level study since it may result in unrealistically high HGL and flow rates in the storm sewer network. Only a few field situations can be modeled properly by this method such as the nodes at each end of a forcemain or a bolted or sealed manhole.

Figure 2

3. Similar to #2 above, sealing a node will also prevent water loss (Figure 3) in which HGL can go up as high as needed. A bolted manhole is an ideal candidate for applying the “Sealed” ponding option.

Figure 3

In a real world, if a storm sewer system is overwhelmed, the excessive water will spill over from the inlet or manhole openings and then ponds on the ground or flows downstream along an overland path such as street gutters. In 1D XPSWMM, this is modeled by allowing ponding, a storage node with depth measured from spill crest, a dual drainage system, or their combination.

4. Setting the ponding option as “Allowed” will allow excessive water to be kept temporarily in a ponding area above the overflowing node and only return to the subsurface storm sewer system at a later time (Figure 4). This is a quick way to account for surface flooding without causing an unrealistically high HGL and can be used for preliminary analysis. According to XPSWMM manual, the ponding area is calculated as A * exp (* Ponding Depth), which can be further edited in Hydraulic Job Control/Junction Defaults.

Figure 4

5. For a more detailed study, use a storage node with its ponding option set as “Allowed” and its storage depth being measured from Spill Crest (Figure 5). This option works in a similar way to #4, with only the ponding area being user-defined to better describe the field situation. In the stepwise linear storage method window (Figure 6), the maximum defined depth probably should be greater than the expected ponding depth; otherwise, XPSWMM will re-use the last surface area by assuming a vertical wall surrounding it to calculate ponding depths above the maximum defined value.

Figure 5
Figure 6

The default stepwise linear storage method requires a paired data set of depth and surface area, however, if desired, a modeler can override the program setting using “VS” configuration parameter to enter depth and volume values (Figure 7).

Figure 7

Using a storage node to model a ponding area above node spill crest can be applied to anywhere there is a node overflowing issue and the spilled water needs to be stored temporarily according to the local topography, although the most ideal locations may be at depressed areas such as a sag inlet or a low area with an inlet in the middle. To facilitate inter-storage area flow transfer, a weir can be set up at the high ground boundary of two neighboring storage areas (Figure 8).

Figure 8

6. The Multi Link option of XPSWMM can be used to model a dual drainage system (minor system – subsurface storm sewer network; major system – open channel flow in street/overland/ditch). Usually the node spill crest should be raised to match the highest elevation that the open channel cross section is representing. For example, if the major system is represented by a link using a natural section from road ROW to ROW, the node spill crest should be raised from its normal position by D as shown in Figure 9. Under some rare situations, the spill crest may need to be raised by more than D to contain higher than expected HGL and in this scenario two vertical glass walls will be assumed at each end of the open channel cross section to match spill crest (assuming VERT_WALLS=ON, refer to this post for more details).

Figure 9

A single link can be converted to a multi-link manually by right clicking the selected single link and choosing Multi Link as shown in Figure 10.

Figure 10

XPSWMM also has a tool to help convert a single link (pipe) to a dual drainage system if the major system is to be defined by a natural section shape. First select as many as single links that need to be converted, and then go to Tools —>Dual Drainage Batch Converter… to make the conversion (Figure 11).

Figure 11

In Figure 11, the max depth is usually left as 0 to let XPSWMM calculate the natural section shape depth as default (Figure 12). After conversion, the major drainage system conduit inverts are to be set at the original node spill crests, while the new node spill crests will be automatically increased by the calculated natural section shape maximum depths.

Figure 12

In a dual drainage system, the node’s inlet capacity feature can be turned on (Figure 13, this example uses Maximum Capacity – the simplest inlet capacity setting for demonstration purpose) to model an inlet’s restriction to incoming flow: only the flow not exceeding the inlet capacity will drop into the inlet chamber to be carried away by the subsurface storm sewer, while the remaining will flow down the major system – a street or an overland open channel.

Figure 13

Enabling a node’s inlet capacity will split the original node into two by XPSWMM internally: a “top” node (connecting major system) and a “bottom” subsurface node (connection minor system) which is denoted by “$I” (Figure 14 and Figure 15).

Figure 14
Figure 15

The two internal nodes are connected by a rating curve (Node name + “$R“) defined by the inlet capacity setting (Figure 16). The flow captured by the inlet capacity rating curve can be viewed in the result window (Figure 17).

Figure 16
Figure 17

When the inlet capacity feature is enabled in a XPSWMM model, extra attention should be paid to model stability. In this example, the maximum inlet capacity had to be lowered from 8 cfs to 4 cfs to get a stable captured flow result as shown in Figure 17.

If the inlet capacity is not enabled in a dual drainage system, an internal node with “$I” will not be created. In this case, a surface open channel will only receive inflows from the nodes it is connected to when the water levels in the nodes rise above the open channel invert.

Leave a Reply

Your email address will not be published. Required fields are marked *