==================== Example Applications ==================== Example 1 - End-valve closure ----------------------------- This example shows how to simulate the closure of a valve located at the boundary of a network. The first example network is shown below in :numref:`tnet1`, adopted from [[STWY67],WOLB05]_. It comprises 9 pipes, 8 junctions, one reservoir, 3 closed loops, and one valve located at the downstream end of the system. There are five steps that the user needs to take to run the transient simulation using the TSNet package: .. _tnet1: .. figure:: figures/Tnet1.PNG :scale: 100 % :alt: tnet1 Tnet1 network graphics 1. Import TSNet package, read the EPANET INP file, and create transient model object. .. literalinclude:: ../examples/Tnet1_valve_closure.py :lines: 1-4 2. Set the wave speed for all pipes to :math:`1200m/s`, time step to :math:`0.1s`, and simulation period to :math:`60s`. .. literalinclude:: ../examples/Tnet1_valve_closure.py :lines: 6-11 3. Set valve operating rules, including how long it takes to close the valve (:math:`tc`), when to start close the valve (:math:`ts`), the opening percentage when the closure is completed (:math:`se`), and the shape of the closure operating curve (:math:`m`, :math:`1` stands for linear closure, :math:`2` stands for quadratic closure). .. literalinclude:: ../examples/Tnet1_valve_closure.py :lines: 13-19 4. Compute steady state results to establish the initial condition for transient simulation. .. literalinclude:: ../examples/Tnet1_valve_closure.py :lines: 21-24 5. Run transient simulation and specify the name of the results file. .. literalinclude:: ../examples/Tnet1_valve_closure.py :lines: 26-28 After the transient simulation, the results at nodes and links will be returned and stored in the transient model (tm) instance. The time history of flow rate on the start node of pipe P2 throughout the simulation can be retrieved by: >>> print(tm.links['P2'].start_node_flowrate) To plot the head results at N3: .. literalinclude:: ../examples/Tnet1_valve_closure.py :lines: 31-42 yields :numref:`tnet1_node`: .. _tnet1_node: .. figure:: figures/tnet1_node.png :width: 600 :alt: tnet1_node Tnet1 - Head at node N3. Similarly, to plot the flow rate results in pipe P2: .. literalinclude:: ../examples/Tnet1_valve_closure.py :lines: 45-58 yields :numref:`tnet1_pipe`: .. _tnet1_pipe: .. figure:: figures/tnet1_pipe.png :width: 600 :alt: tnet1_pipe Tnet1 - Flow rate at the start and end node of pipe P2. Example 2 - Pump operations --------------------------- This example illustrates how the package models a transient event resulting from a controlled pump shut-off , i.e., the pump speed is ramped down. This example network, Tnet2, is shown below in :numref:`tnet2`. Tnet2 comprises 113 pipes, 91 junctions, 2 pumps, 2 reservoir, 3 tanks, and one valve located in the middle of the network. A transient simulation of 50 seconds is generated by shutting off PUMP2. There are five steps user needs to take: .. _tnet2: .. figure:: figures/Tnet2.PNG :width: 600 :alt: tnet2 Tnet2 network graphics 1. Import TSNet package, read the EPANET INP file, and create transient model object. .. literalinclude:: ../examples/Tnet2_pump_shutdown.py :lines: 1-4 2. Set the wave speed for all pipes to be :math:`1200m/s` and simulation period to be :math:`50s`. Use suggested time step. .. literalinclude:: ../examples/Tnet2_pump_shutdown.py :lines: 6-10 3. Set pump operating rules, including how long it takes to shutdown the pump (:math:`tc`), when to the shut-off starts (:math:`ts`), the pump speed multiplier value when the shut-off is completed (:math:`se`), and the shape of the shut-off operation curve (:math:`m`, :math:`1` stands for linear closure, :math:`2` stands for quadratic closure). .. literalinclude:: ../examples/Tnet2_pump_shutdown.py :lines: 12-18 4. Compute steady state results to establish the initial condition for transient simulation. .. literalinclude:: ../examples/Tnet2_pump_shutdown.py :lines: 20-23 5. Run transient simulation and specify the name of the results file. .. literalinclude:: ../examples/Tnet2_pump_shutdown.py :lines: 25-27 After the transient simulation, the results at nodes and links will be returned to the transient model (tm) instance, which is then stored in **Tnet2.obj**. The actual demand discharge at JUNCTION-105 throughout the simulation can be retrieved by: >>> print(tm.nodes['JUNCTION-105'].demand_discharge) To plot the head results at JUNCTION-105: .. literalinclude:: ../examples/Tnet2_pump_shutdown.py :lines: 30-41 yields :numref:`tnet2_node`: .. _tnet2_node: .. figure:: figures/tnet2_node.png :width: 600 :alt: tnet2_node Tnet2 - Head at node JUNCTION-105. Similarly, to plot the velocity results in PIPE-109: .. literalinclude:: ../examples/Tnet2_pump_shutdown.py :lines: 44-55 yields :numref:`tnet2_pipe`: .. _tnet2_pipe: .. figure:: figures/tnet2_pipe.png :width: 600 :alt: tnet2_pipe Tnet2 - Velocity at the start and end node of PIPE-109. Example 3 - Burst and leak --------------------------- This example reveals how TSNet simulates pipe bursts and leaks. This example network, adapted from [OSBH08]_, is shown below in :numref:`tnet3`. Tnet3 comprises 168 pipes, 126 junctions, 8 valve, 2 pumps, one reservoir, and two tanks. The transient event is generated by a burst and a background leak. There are five steps that the user would need to take: .. _tnet3: .. figure:: figures/Tnet3.PNG :width: 600 :alt: tnet3 Tnet3 network graphics 1. Import TSNet package, read the EPANET INP file, and create transient model object. .. literalinclude:: ../examples/Tnet3_burst_leak.py :lines: 1-4 2. The user can import custom wave speeds for each pipe. To demonstrate how to assign different wave speed, we assume that the wave speed for the pipes is normally distributed with mean of :math:`1200 m/s` and standard deviation of :math: `100m/s`. Then, assign the randomly generated wave speed to each pipe in the network according to the order the pipes defined in the INP file. Subsequently, set the simulation period as :math:`20s`, and use suggested time step. .. literalinclude:: ../examples/Tnet3_burst_leak.py :lines: 6-12 3. Define background leak location, JUNCTION-22, and specify the emitter coefficient. The leak will be included in the initial condition calculation. See WNTR documentation [WNTRSi]_ for more info about leak simulation. .. literalinclude:: ../examples/Tnet3_burst_leak.py :lines: 14-16 4. Compute steady state results to establish the initial condition for transient simulation. .. literalinclude:: ../examples/Tnet3_burst_leak.py :lines: 24-27 5. Set up burst event, including burst location, JUNCTION-20, burst start time (:math:`ts`), time for burst to fully develop (:math:`tc`), and the final emitter coefficient (final_burst_coeff). .. literalinclude:: ../examples/Tnet3_burst_leak.py :lines: 18-22 6. Run transient simulation and specify the name of the results file. .. literalinclude:: ../examples/Tnet3_burst_leak.py :lines: 29-31 After the transient simulation, the results at nodes and links will be returned to the transient model (tm) instance, which is subsequently stored in **Tnet3.obj**. To understand how much water has been lost through the leakage at JUNCTION-22, we can plot the leak discharge results at JUNCTION-22: .. literalinclude:: ../examples/Tnet3_burst_leak.py :lines: 34-45 yields :numref:`tnet3_leak`: .. _tnet3_leak: .. figure:: figures/tnet3_leak.png :width: 600 :alt: tnet3_leak Tnet3 - Leak discharge at node JUNCTION-22. Similarly, to reveal how much water has been wasted through the burst event at JUNCTION-20, we can plot the burst discharge results at JUNCTION-20: .. literalinclude:: ../examples/Tnet3_burst_leak.py :lines: 48-58 yields :numref:`tnet3_burst`: .. _tnet3_burst: .. figure:: figures/tnet3_burst.png :width: 600 :alt: tnet3_burst Tnet3 - Burst discharge at node JUNCTION-20. Additionally, to plot the velocity results in LINK-40: .. literalinclude:: ../examples/Tnet3_burst_leak.py :lines: 61-72 yields :numref:`tnet3_pipe`: .. _tnet3_pipe: .. figure:: figures/tnet3_pipe.png :width: 600 :alt: tnet3_pipe Tnet3 - Velocity at the start and end node of LINK-40. Moreover, we can plot head results at some further nodes, such as JUNCTION-8, JUNCTION-16, JUNCTION-45, JUNCTION-90, by: .. literalinclude:: ../examples/Tnet3_burst_leak.py :lines: 75-90 The results are demonstrated in :numref:`tnet3_multi`. It can be noticed that the amplitude of the pressure transient at JUNCTION-8 and JUNCTION-16 is greater than that at other two junctions which are further away from JUNCTION-20, where the burst occurred. .. _tnet3_multi: .. figure:: figures/tnet3_multi.png :width: 600 :alt: tnet3_multi Tnet3 - Head at multiple junctions.