Next Article in Journal
Flash Flood Prediction by Coupling KINEROS2 and HEC-RAS Models for Tropical Regions of Northern Vietnam
Previous Article in Journal
Non-Stationary Flood Frequency Analysis in the Ouémé River Basin, Benin Republic
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Migrating Storms and Optimal Control of Urban Sewer Networks

Department of Civil Engineering, Faculty of Engineering, Sri Lanka Institute of Information Technology, Malabe 10115, Sri Lanka
Hydrology 2015, 2(4), 230-241; https://doi.org/10.3390/hydrology2040230
Submission received: 1 September 2015 / Revised: 13 October 2015 / Accepted: 1 November 2015 / Published: 4 November 2015

Abstract

:
Uniform storms are generally applied in most of the research on sewer systems. This is for modeling simplicity. However, in the real world, these conditions may not be applicable. It is very important to consider the migration behavior of storms not only in the design of combined sewers, but also in controlling them. Therefore, this research was carried out to improve Rathnayake and Tanyimboh’s optimal control algorithm for migrating storms. Promising results were found from the model improvement. Feasible solutions were obtained from the multi-objective optimization and, in addition, the role of on-line storage tanks was well placed.

1. Introduction

Uniform storms are generally applied in most of the research on sewer systems [1]. This is for modeling simplicity. However, these conditions may not be applicable in the real world. The spatial distribution of rainfall events (such as land use) leads to non-uniform hydrographs. These non-uniform storms affect the runoff hydrograph and the magnitude of the peak discharge of the runoff hydrograph. The overall peak discharge of a downstream migrating storm exceeds that of an upstream migrating storm [1,2,3,4]. The safety level is low compared to a uniform rainfall event, with equal directions for both wind and flow [1], and can easily cause combined sewer overflows (CSOs). This notable difference has been excluded when considering the uniform rainfall conditions for the simulations for simplicity.
CSOs are a major environmental concern for many cities in the world. There are many approaches in reducing these CSOs. Sustainable urban drainage systems (SUDS) are one way of addressing the issue. SUDS try to reduce the storm water flow to the combined sewer networks during the stormy periods. This is done by introducing natural infiltration techniques to infiltrate the majority of the storm water. In addition, providing more storage facilities to store wastewater during stormy periods is another alternative. Controlling the existing combined sewer networks is another option and it has taken the interest of the world.
Optimal control of the sewer system is not a new topic; however, it is still a challenging topic in introducing the integration of optimal control of urban sewer systems to real-time control with the consideration of water quality aspects. Therefore, a holistic approach is still to be tabled for discussions. Rathnayake and Tanyimboh [5,6,7,8] have presented an optimal control model for urban sewer systems considering the water quality effects in receiving water due to CSOs. However, that was for uniform storms. In other words, runoff storms have reached the corresponding interceptor sewers at the same time. This paper fills that identified gap in research by introducing the migrating behavior of storms. Rathnayake and Tanyimboh’s optimal control model [5,6,8] was improved with two different migrating storms, including migrating downstream and migrating upstream storms.

2. Mathematical Formulations

Figure 1 illustrates the schematic diagram of an interceptor sewer system. These schematics were used to develop the multi-objective optimization approach in controlling the urban sewer systems.
Figure 1. Schematic diagram of interceptor sewer system.
Figure 1. Schematic diagram of interceptor sewer system.
Hydrology 02 00230 g001
Equation (1) gives the formulated first objective function (F1) for the multi-objective optimization approach. This was developed to minimize the pollution load from CSOs. The effluent quality index (EQI) was used to formulate this pollution load. Concentrations of total suspended solids (TSS), chemical oxygen demand (COD), five-day biochemical oxygen demand (BOD), total Kjeldahl nitrogen (TKN), and nitrates/nitrites (NOX) are integrated to develop this EQI. A detailed explanation of this EQI can be found in Rathnayake and Tanyimboh [5,6,7].
M i n i m i z e    F 1 = i = 1 n P i
where n and Pi are the number of interceptor nodes or CSO chamber points and the pollution load to the receiving water from the ith CSO chamber, respectively. The second objective function (F2), given in Equation (2), was formulated to minimize the wastewater treatment cost at downstream wastewater treatment plants.
M i n i m i z e    F 2 = C T
where CT (€/year) is the treatment cost at the treatment plant. The corresponding control settings of the sewer networks are the decision variables of these two objective functions. The CT is expressed as a function of the wastewater volume flow rate (qT) to the wastewater treatment plant and given in the following equation (Equation (3)).
C T = { 916.862 × ( 86400 × q T ) 0.659 q T 3 × D W F 916.862 × ( 3 × D W F ) 0.659 + 2 3 ( 1.69 ( q T 3 × D W F ) + 11376 ) , 6 × D W F q T 3 × D W F 916.862 × ( 3 × D W F ) 0.659 + 2 3 ( ( 1.69 × 3 × D W F ) + 11376 ) , q T > 6 × D W F
where qT (m3/s) is the treated wastewater volume flow rate. A detailed explanation on the derivation of this generic cost function is given in Rathnayake and Tanyimboh [5,6,7].
Referring to Figure 1, the following continuity equations can be listed. Equation (4) is for the continuity in the interceptor node, whereas Equations (5) and (6) are the conditional continuity equations for the sewer chamber. These two conditional continuity equations are based on the capacity of the ith sewer chamber.
Q i + q i 1 q i = 0
A C Δ h C Δ t = I i Q i ;   h C < h S
A C Δ h C Δ t = I i Q i O i ;   h C > h s
where AC is the surface area of the CSO chamber. Flows inside the interceptor sewers are constrained. These constraints are inserted in the multi-objective optimization approach. The mathematical formulations for the constraints are given in the following equation (Equation (7)).
0 q i q max , i
where qmax,i is the maximum flow rate at the ith conduit. Often, combined sewer networks are constructed with storage tanks. They are to provide additional capacities to the sewer network. Formulations for on-line storage tanks were carried out in this research. More information on these formulations can be found in Rathnayake and Tanyimboh [8]. Furthermore, the above-stated continuity Equations (4)–(6) and constraints (Equation (7)) can handle the spatial and temporal behavior of the inflow. Therefore, the implementation of migration storms can be managed without any disruptions.

3. Solution Technique

U.S. EPA SWMM 5.0 [9] is a powerful hydraulic and water quality simulation model. This software is capable of simulating stormwater runoff and routing processes, including water quality routing. SWMM has been successfully used in many practical design and control problems in both combined and separate sewer systems [10,11]. Gradually varied, unsteady flows inside sewers are routed using the mass conservation and momentum equations, i.e., one-dimensional St. Venant equations. In addition, the sewers are assumed to behave as continuously stirred tank reactors for water quality routing. Mass conservation is used to calculate the concentrations of each water quality parameter leaving the conduits at the end of a time-step. The same principle is applied at the storage nodes.
SWMM 5.0 was linked to Non-dominated Sorting Genetic Algorithm (NSGA II) [12] using the C programming language. NSGA II is a multi-objective optimization module and it has already been successfully applied to many practical optimization problems in various disciplines.
NSGA II is a computationally fast and elitist multi-objective evolutionary algorithm. Briefly, an initial parent population of size N is generated randomly. An offspring population also of size N is then generated from the parent population. This involves binary tournament selection and crossover operations that create two children from two parents. On completion of the crossover, a mutation operator that alters a selection of individuals in the offspring population is applied. Then, the parent and offspring populations are merged to form a population of size 2N. The merged population is then sorted into various levels or fronts of non-domination based on the relevant fitness criteria. Finally, a new population of size N is created by taking the best solutions first, according to the rank order of the respective non-domination levels. If any non-domination level has more solutions than required, then preference is given to solutions in the parts of the solution space that are under-represented. This cycle is repeated until the algorithm convergences and/or the specified termination criteria are reached. At the end of the optimization run, the solutions to the objective functions can be plotted as “Pareto optimal fronts”.
A rectangular orifice was placed at the bottom of each CSO chamber. This orifice was used to control the wastewater flow to the interceptor sewer. The orifice openings are, therefore, the control settings or decision variables of the developed multi-objective optimization approach. These orifice openings were initially generated randomly. Then, a full hydraulic simulation, including water quality routing, was carried out using SWMM 5.0. In order to carry out the hydraulic simulations, runoff hydrographs from migrating storms and their corresponding pollutographs were fed to the respective CSO chambers. Generic pollutographs for total suspended solids (TSS), chemical oxygen demand (COD), biochemical oxygen demand (BOD), nitrates and nitrites (NOX), and total Kjeldahl nitrogen (TKN) were generated based on the previous literature [13,14,15,16,17,18,19]. Hydraulic and water quality results from the simulations were used to calculate the two objective functions (F1 and F2).
Continuity equations given in Equations (4)–(6) are automatically satisfied by the hydraulic model, SWMM 5.0. However, the constraints given in Equation (7) are satisfied inside the multi-objective optimization module, NSGA II. Constraints in NSGA II are satisfied by a tournament constraint handling approach [12].

4. Case Study

The interceptor sewer system used for this analysis can be seen in Figure 2. The interceptor sewer system found in Thomas [20] was modified for this study.
Figure 2. Interceptor sewer system.
Figure 2. Interceptor sewer system.
Hydrology 02 00230 g002
Two on-line storage tanks (T8 and T9) were introduced to T2 and T5 CSO chambers. These two on-line storage tanks were controlled inside the hydraulic model. The control rules in SWMM 5.0 were used to control the on-line storage tank. When the wastewater level in the CSO chamber reaches its spill level, the wastewater is transferred to the attached on-line storage tank. However, when the on-line storage tank researches its maximum capacity, this diversion of wastewater from the CSO chamber to the storage tank stops. This avoids having CSOs from any on-line storage tanks. However, when there is less stress on the combined sewer network, the stored wastewater is released back to combined sewer system.
The lengths of the C1C7 conduits are 895, 740, 460, 15, 719, 19, 710, 350, and 196 m, respectively. More information on the geometric data of the interceptor sewer system is found in Thomas [20] and Rathnayake and Tanyimboh [6]. Flow rates inside the conduits were constrained to 3.26 m3/s in C1C3 conduits and 7.72 m3/s in C4C7 conduits.
Average dry weather flows were fed to T1T7 CSO chambers. Runoff hydrographs for migrating downstream and migrating upstream runoff hydrographs were separately fed to the T1T7 CSO chambers for two analyses. These migrating storms are illustrated in Figure 3a,b. For simplicity, only a set of generic hydrographs were fed to CSO chambers. The hydrographs presented in Thomas [20] and Rathnayake [6] were modified to represent the migration behavior of the storms. Therefore, Figure 3a presents the migration downstream behavior of the storms whereas Figure 3b shows the migration upstream behavior. However, there is a flexibility in applying the true rainfall data to generate migrating storms in the hydraulic model. Therefore, the model can be applied for any variation of runoff hydrographs.
Figure 3. (a) Stormwater runoff hydrograph for migrating downstream storms; (b) Stormwater runoff hydrograph for migrating upstream storms.
Figure 3. (a) Stormwater runoff hydrograph for migrating downstream storms; (b) Stormwater runoff hydrograph for migrating upstream storms.
Hydrology 02 00230 g003
The peak duration of the hydrographs for the same CSO chamber variation can be clearly seen from Figure 3a,b. The peak flow of the Rimrose catchment to the T1 CSO chamber is around the 30th minute for the migrating downstream storms. However, that of the migrating upstream storms is around the 90th minute. In addition, the peak of the Northern is around the 90th minute in the migrating downstream storms (Figure 3a). Therefore, the two peak flows have a 45 min difference. These features show the migration features of the storms. In addition to the runoff hydrographs, corresponding pollutographs for five different water quality constituents (TSS, COD, BOD, NOX, and TKN) for different land uses (residential, commercial, industrial, agricultural, and mid-urban) were fed to each CSO chamber for the two cases.
A few pollutographs generated based on previous literature [13,14,15,16,17,18,19] for two different migrating storms can be seen in Figure 4. Figure 4a,b show the BOD pollutographs and the corresponding stormwater runoffs for the Millers Bridge and Sandhills Lane catchments. The concentration levels are different in the two figures due to the different land uses. However, the migrating effects can be clearly visualized in the two figures. The peak of the BOD pollutograph of Figure 4a is roughly at the 30th min of the storm. However, the peak of the BOD pollutograph in Figure 4b has migrated roughly to the 1 h time duration. This is according to the migration of the corresponding stormwater runoff hydrograph. In contrast, Figure 4c,d show the COD pollutographs and the corresponding stormwater runoffs for the Rimrose and Bankhall catchments. The concentration levels are different in the two catchments due to the different land uses. However, the migrating upstream effects can be clearly visualized in the two figures. The peak of the COD pollutograph of the Rimrose catchment is roughly at the 1 h mark of the storm. However, the peak of the COD pollutograph at the Bankhall catchment has migrated roughly to the 30th min of the time duration. This is again according to the migration of the corresponding stormwater runoff hydrograph. It is noted herein that the Bankhall catchment is at the downstream of the combined sewer network.
Figure 4. Pollutographs for migrating storms. (a) BOD pollutograph for Millers Bridge; (b) BOD pollutograph for Sandhills Lane; (c) COD pollutograph for Rimrose; (d) COD pollutograph for Bankhall.
Figure 4. Pollutographs for migrating storms. (a) BOD pollutograph for Millers Bridge; (b) BOD pollutograph for Sandhills Lane; (c) COD pollutograph for Rimrose; (d) COD pollutograph for Bankhall.
Hydrology 02 00230 g004
Combined sewer network modeling was carried out using the above-stated information. The obtained results from the hydraulic and water quality modeling were input to the NSGA II optimization module. A real-coded NSGA II program was used in this study. The optimization process was done with a population of 100,100 generations for both migrating storm conditions. Many optimization runs with different random seeds were conducted. The routing time-step in SWMM 5.0 was kept at 30 s, and the results were obtained at each 15 min interval for the total duration of the storms. Then, the NSGA II optimization module was run using the obtained results. Each Genetic Algorithm (GA) run took about 42 to 47 min on an Intel® Core™ i3 desktop personal computer with a 3.40 GHz processor and 4 GB of RAM.

5. Results and Discussion

Figure 5a,b show the best Pareto optimal fronts for both storm conditions after 3 h. For the calibration purposes many optimization runs were carried out with different mutation rates. It was found out that the best Pareto optimal fronts can be achieved when the mutation rate is at 0.4 for both migrating downstream and upstream storms. These two Pareto optimal fronts look similar in shape. However, they have different maximums and minimums. For example, the minimum pollution load for the migrating downstream storm condition is 22.67 kt/day whereas that for the migrating upstream storm condition is 21.98 kt/day. The corresponding treatment costs for those mentioned above are 49.87 M.€/year and 49.96 M.€/year, respectively. Therefore, these numerical values clearly show the effects of migrating storms.
Figure 5. Pareto optimial fronts after 3 h runoff. (a) For migrating downstream; (b) For migrating upstream.
Figure 5. Pareto optimial fronts after 3 h runoff. (a) For migrating downstream; (b) For migrating upstream.
Hydrology 02 00230 g005
Five optimal solutions (A to E and P to T) for each case were selected for further analysis. Solutions A to E are for the migrating downstream storm conditions and solutions P to T are for the migrating upstream storm conditions. Solutions A and P correspond to the minimum pollution loads in each case and solutions E and T give minimum treatment costs in each case. Control settings (orifice openings) were obtained for these 10 optimal solutions. Figure 5a,b presents the orifice openings for four optimal solutions.
Figure 6a shows the orifice openings for minimum pollution load solutions, whereas Figure 6b gives the orifice openings for minimum treatment cost solutions. Even though the Pareto optimal fronts (shown in Figure 5a,b) show similar patterns, the corresponding orifice openings are different. For example, the orifice opening for the O1 orifice in the migrating downstream storm condition is about 21 cm; however, that of the migrating upstream storm condition is about 6 cm. This observation can be seen for all other orifices in Figure 6a,b. Therefore, the developed model is capable of handling migrating storms and then producing sound control settings. These obtained orifice openings were used to carry out the hydraulic simulations for each case.
Figure 6. Orifice openings for optimal solutions. (a) Minimum pollution load case; (b) Minimum treatment cost case.
Figure 6. Orifice openings for optimal solutions. (a) Minimum pollution load case; (b) Minimum treatment cost case.
Hydrology 02 00230 g006
Figure 7a,b present the flow rates inside the interceptor sewers for solution A. As it was stated in the case study section, the flow rates in the conduits were constrained. The maximum allowable flow rate in C1 to C3 conduits was 3.26 m3/s and in C4 to C7 was 7.72 m3/s. Red-color dashed lines in Figure 7a,b graphically represent these maximum allowable flow rates. It can be clearly seen herein that the flow rates in the conduits are lower than the maximum allowed. This shows the developed multi-objective optimization approach produces feasible solutions. In addition, flow rates for the rest of the solutions were checked and the same conclusion can be determined.
Figure 7. Flow rates in interceptor sewers. (a) C1, C2 and C3; (b) C4, C5, C6 and C7.
Figure 7. Flow rates in interceptor sewers. (a) C1, C2 and C3; (b) C4, C5, C6 and C7.
Hydrology 02 00230 g007
Table 1 presents the wastewater heights in the CSO chambers and storage tanks for solution A. The geometric depths of these tanks from T1 to T9 are 5.42, 6.91, 7.95, 8.04, 8.18, 8.47, 9.26, 6.91, and 8.18 m, respectively. The grey color values shown in the following table are higher than the geometric heights of the corresponding CSO chambers. They represent the CSOs. However, the bold values for the T8 and T9 storage tanks are exactly the same as the geometric depths. They are kept at the maximum possible storage limits and that shows that the role of the storage tanks was well developed for the multi-objective optimization approach. In addition, as it is expected, there are no CSOs for these two storage tanks. This further suggests the robustness of the developed control algorithm.
Table 1. Water heights in CSO chambers and storage tanks.
Table 1. Water heights in CSO chambers and storage tanks.
Time (h:min:s)Water Height (m)
T1T2T3T4T5T6T7T8T9
0:15:005.765.555.647.27.714.643.0601.67
0:30:005.997.168.218.228.498.594.966.918.18
0:45:006.147.268.498.58.748.646.296.918.18
1:00:006.17.248.478.899.068.769.426.918.18
1:15:005.877.148.359.049.188.869.566.918.18
1:30:005.647.068.219.049.188.839.546.918.18
1:45:005.617.048.018.889.148.729.466.918.18
2:00:005.617.047.578.698.928.629.346.918.18
2:15:005.617.047.468.358.628.69.26.918.18
2:30:005.617.047.448.228.498.69.096.918.18
2:45:005.617.047.438.228.498.696.918.18
3:00:005.617.037.438.218.488.598.946.918.18
Table 2 presents the total treated wastewater volumes and total CSO volumes for the four extreme solutions (A, E, P and T). Solution A corresponds to the minimum pollution load approach in migrating downstream storms and shows much higher treated wastewater volumes compared to that of solution E, which corresponds to the minimum treatment cost. A similar observation can be visualized in the migrating upstream storm case. In contrast, total CSO volumes in solutions A and P are lower than those in solutions E and T. Therefore, the optimization algorithm is robust even in volumetric measures.
Table 2. Total treated wastewater and CSOs.
Table 2. Total treated wastewater and CSOs.
AEPT
Total treated wastewater volume (m3)5432.112.35394.511.5
Total CSO volume (m3)6876.27737.36682.87392.7

6. Conclusions

Rathnayake and Tanyimboh’s optimal control algorithm was successfully applied for migrating storm conditions. Hydraulic analysis shows that the developed multi-objective optimization control algorithm produces feasible results. In addition, the role of the storage tanks was well placed for the control algorithm. Further research should be carried out in order to apply this optimal control algorithm in real-time controls.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vaes, G.; Clemens, F.; Willems, P.; Berlamont, J. Design rainfall for combined sewer system calculations. In Proceedings of the 9th International Conference on Urban Drainage, Portland, OR, USA, 8–13 September 2002.
  2. Nelen, F.; Mooijman, A.; Jacobsen, P. The importance of rainfall distribution in urban drainage operation. Nord. Hydrol. 1992, 23, 121–136. [Google Scholar]
  3. Hochedlinger, M.; Sprung, W.; Kainz, H.; König, K. CSO Modelling considering moving storms and tipping bucket gauge failures. In Proceedings of the 2nd International IWA Conference on Sewer Operation and Maintenance, Vienna, Austria, 26–28 October 2006.
  4. Zawilski, M.; Brzezińska, A. Spatial rainfall intensity distribution over an urban area and its effect on a combined sewerage system. In Proceedings of the 12th International Conference on Urban Drainage, Porto Alegre, Brazil, 11–16 September 2011.
  5. Rathnayake, U.S.; Tanyimboh, T.T. Optimal control of combined sewer systems using SWMM5.0. WIT Trans. Built Environ. 2012, 122, 87–96. [Google Scholar]
  6. Rathnayake, U.S.; Tanyimboh, T.T. Integrated optimal control of urban wastewater systems. In Proceedings of the IWA World Congress on Water, Climate and Energy, Dublin, Ireland, 13–18 May 2012.
  7. Rathnayake, U.S.; Tanyimboh, T.T. Multi-objective optimization of urban wastewater systems. In Proceedings of the 10th International Conference on Hydroinfomatics (HIC2012), Hamburg, Germany, 14–18 July 2012.
  8. Rathnayake, U.; Tanyimboh, T.T. Multi-objective optimization of combined sewer systems using SWMM 5.0. In Proceedings of the 13th International Conference on Urban Drainage (13ICUD), Sarawak, Malaysia, 7–12 September 2014.
  9. Rossman, L.A. US EPA SWMM 5.0 User’s Manual, Manual EPA/600/R-05/040, Water Supply and Water Resources Division; National Risk Management Research Laboratory: Cincinnati, OH, USA, 2009. [Google Scholar]
  10. Vojinovic, Z.; Sanchez, A.T.; Barreto, W.J. Optimising sewer systems rehabilitation strategies between flooding overflow emissions and investment costs. In Proceedings of the 11th International conference on Urban Drainage, Edinburgh, UK, 31 August–5 September 2008.
  11. Shinma, T.A.; Reis, L.F.R. Analysis of different objective functions sets applied to automatic calibration of the storm water management model (SWMM). In Proceedings of the 11th International Conference Computing and Control for the WATER Industry, Exeter, UK, 5–7 September 2011.
  12. Deb, K.; Pratap, A.; Agrawal, S.; Meriyarivan, T. A fast elitist multi-objective genetic algorithm: NSGA II. IEEE Trans. Evol. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef]
  13. Morris, G.M.; Agg, A.R.; Clifforde, I.T.; Ridout, K.J.; Singleton, P.; Tyson, J.M.; Wilson, A. Urban Pollution Management Manual, 2nd ed.; Foundation for Water Research; Allen House: Buckinghamshire, UK, 1998. [Google Scholar]
  14. Duncan, P. Urban Stormawater Quality: A Statistical Overview; Cooperative Research Centre for Catchment Hydrology: Victoria, Australia, 1999. [Google Scholar]
  15. Yusop, Z.; Tan, L.W.; Ujang, Z.; Mohamed, M.; Nasir, K.A. Runoff quality and pollution loading from a tropical urban catchment. Water Sci. Technol. 2005, 52, 125–132. [Google Scholar] [PubMed]
  16. Li, L.; Yin, C.; He, Q.; Kong, L. First flush of storm runoff pollution from an urban catchment in China. J. Environ. Sci. 2007, 19, 295–299. [Google Scholar] [CrossRef]
  17. Kim, G.; Yur, J.; Kim, J. Diffuse pollution loading from urban stormwater runoff in Daejeon city, Korea. J. Environ. Manag. 2007, 85, 9–16. [Google Scholar] [CrossRef] [PubMed]
  18. Nazahiyah, R.; Yusop, Z.; Abustan, I. Stormwater quality and pollution loading from an urban residential catchment in Johor, Malaysia. Water Sci. Technol. 2007, 56, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Qin, H.; Khu, S.; Yu, X. Spatial variations of storm runoff pollution and their correlation with land-use in a rapidly urbanizing catchment in China. Sci. Total Environ. 2010, 408, 4613–4623. [Google Scholar] [CrossRef] [PubMed]
  20. Thomas, N.S. Optimal Pollution Control Models for Interceptor Sewer Systems. Ph.D. Thesis, Department of Civil Engineering, University of Liverpool, Liverpool, UK, 2000. [Google Scholar]

Share and Cite

MDPI and ACS Style

Rathnayake, U. Migrating Storms and Optimal Control of Urban Sewer Networks. Hydrology 2015, 2, 230-241. https://doi.org/10.3390/hydrology2040230

AMA Style

Rathnayake U. Migrating Storms and Optimal Control of Urban Sewer Networks. Hydrology. 2015; 2(4):230-241. https://doi.org/10.3390/hydrology2040230

Chicago/Turabian Style

Rathnayake, Upaka. 2015. "Migrating Storms and Optimal Control of Urban Sewer Networks" Hydrology 2, no. 4: 230-241. https://doi.org/10.3390/hydrology2040230

APA Style

Rathnayake, U. (2015). Migrating Storms and Optimal Control of Urban Sewer Networks. Hydrology, 2(4), 230-241. https://doi.org/10.3390/hydrology2040230

Article Metrics

Back to TopTop