Search

Accurate flood modeling is crucial for effective analysis and forecasting. Full momentum hydrodynamic models often require extensive computational time, sometimes exceeding the forecast horizon. In contrast, low-complexity models, like local-inertial approximations, provide accurate results in subcritical flows but may have limited skillfulness in supercritical conditions. This paper explores two main aspects: (i) the impact of urban infrastructure on 2D hydrodynamic modeling without detailed sewer and drainage data, and (ii) the accuracy of 2D local-inertial modeling using three numerical schemes (original formulation, s-centered, and s-upwind) in a dam-break scenario on complex, flat terrain. The HydroPol2D model is benchmarked against HEC-RAS 2D full momentum solver. We present one numerical case study and three real-world scenarios in S\~ao Paulo, Brazil: a detention pond with a $1$ in $100$-year inflow, a highly urbanized catchment with a $1$ in $50$-year hyetograph, and a dam-break scenario threatening a coastal city of nearly 200,000 residents. Results show that the model accurately simulates internal boundary conditions, achieving peak errors under 5\% compared to HEC-RAS 2D. However, neglecting urban infrastructure can lead to a 17.5\% difference in peak discharges at the outlet and significant mismatches in hydrographs, with computational times nearly doubling. The dam-break scenario demonstrates good predictive performance for maximum flood depths (CSI = $0.95$ for the original model, $0.92$ for s-centered, and $0.89$ for s-upwind), though the model's lack of convective inertia results in faster flood wave propagation than the full momentum solver. Notably, HydroPol2D is 23 times faster than HEC-RAS 2D, making it well-suited for simulating dam collapses in forecasting systems and capable of modeling urban drainage infrastructure such as orifices, weirs, and pumps., Comment: 1 table and 10 figures, under review

Cost-effective water distribution network (WDN) design with acceptable pressure performance is crucial for the management of drinking water in cities. This paper presents a Microsoft Excel tool to model, simulate, and optimize WDNs with looped pipelines under steady-state incompressible flow simulations. Typically, the hardy-cross method is applied using spreadsheet calculations to estimate discharges. This method requires mass-conservative initial estimates and requires successive iterations to converge. In this paper, however, we develop an alternative method that uses the built-in solver capabilities of Excel, does not require initial mass-conservative estimation, and is free of flow corrections. The main objective of this paper is to develop an open-source accessible tool for simulating hydraulic networks also adapted for teaching and learning purposes. The governing equations and the mathematical basis for the hydraulic modeling of the system are mathematically described, considering the topology of the network, mass and energy conservation, cost of tank material, foundation, and cost of pumping energy to fill the tank. The use of this tool is encouraged at the undergraduate and graduate engineering levels, as it offers the opportunity to address complex concepts in a comprehensive way using a spreadsheet that does not require coding expertise. Hence, users can debug all cells and understand all equations used in the hydraulic model, as well as modify them. To demonstrate the model capabilities, three practical examples are presented, with the first one solved step by step, and the results are compared with the EPANET and with the results reported in the literature. Using the optimization method presented in this paper, it was possible to achieve a cost reduction of 151,790 USD (9.8% of the total cost) in a network that supplies a 44,416 population.

Detention ponds can mitigate flooding and improve water quality by allowing the settlement of pollutants. Typically, they are operated with fully open orifices and weirs (i.e., passive control). Active controls can improve the performance of these systems: orifices can be retrofitted with controlled valves and spillways can have controllable gates. The real-time optimal operation of its hydraulic devices can be achieved with techniques such as Model Predictive Control (MPC). A distributed quasi-2D hydrologic-hydrodynamic coupled with a reservoir flood routing model is developed and integrated with an MPC algorithm to estimate the operation of valves and movable gates. The control optimization problem is adapted to switch from a flood-related algorithm focusing on mitigating floods to a heuristic objective function that aims to increase the detention time when no inflow hydrographs are predicted. The case studies show the potential of applying the methods developed in a catchment in Sao Paulo, Brazil. The performance of MPC compared to alternatives with either fully or partially open valves and gates are tested. Comparisons with HEC-RAS 2D indicate volume and peak flow errors of approximately 1.4% and 0.91% for the watershed module. Simulating two consecutive 10-year storms shows that the MPC strategy can achieve peak flow reductions of 79%. In contrast, passive control has nearly half of the performance. A 1-year continuous simulation results show that the passive scenario with 25% of the valves opened can treat 12% more runoff compared to the developed MPC approach, with an average detention time of approximately 6 hours. For the MPC approach, the average detention time is nearly 14 hours indicating that both control techniques can treat similar volumes; however, the proxy water quality for the MPC approach is enhanced due to the longer detention times.

Hydrodynamic models with rain-on-the-grid capabilities are usually computationally expensive. This makes the use of automatic calibration algorithms hard to apply due to the large number of model runs. However, with the recent advances in parallel processing, computational resources, and increasing high-resolution climatologic and GIS data, high-resolution hydrodynamic models can be used for optimization-based calibration. This paper presents a global optimization-based algorithm to calibrate a fully distributed hydrologic-hydrodynamic and water quality model (HydroPol2D) using observed data (i.e., discharge, or pollutant concentration) as input. The algorithm can find a near-optimal set of parameters to explain observed gauged data. The modeling framework presented here, although applied in a poorly-gauged catchment, can be adapted for catchments with more detailed observations. We applied the algorithm in different cases of the V-Tilted Catchment, the Wooden-Board catchment, and in an existing urban catchment with heterogeneous data. The results of automatic calibration indicate $\mathrm{NSE} = 0.99$ for the V-Tilted catchment, $\mathrm{RMSE} = 830~\mathrm{mgL^{-1}}$ for salt concentration pollutographs (i.e., 8.3% of the event mean concentration), and $\mathrm{NSE} = 0.89$ for the urban catchment case study. This paper also explores the issue of equifinality in modeling calibration (EqMC). Equifinality is defined as the set of different parameter combinations that can provide equally good or accepted results, within the physical parameter ranges. EqMC decreases with the number of events and increases with the choice of partially or nonproducing runoff ones. Furthermore, results indicate that providing more accurate parameter ranges based on a priori knowledge of the catchment is fundamental to reduce the chances of finding a set of parameters with equifinality., Comment: Preprint submitted to Journal of Hydrology

The increasing access to non-expensive sensors, computing power, and more accurate forecasting of storm events provides unique opportunities to shift flood management practices from static approaches to an optimization-based real-time control (RTC) of urban drainage systems. Recent studies have addressed a plethora of strategies for flood control in stormwater reservoirs; however, advanced control theoretic techniques are not yet fully investigated and applied to these systems. In addition, there is an absence of a coupled integrated control model for systems composed of watersheds, reservoirs, and channels for flood mitigation. To this end, we develop a novel state-space model of hydrologic and hydrodynamic processes in reservoirs and one-dimensional channels. The model is tested under different types of reservoir control strategies based on real-time measurements (reactive control), and based on predictions of the future behavior of the system (predictive control) using rainfall forecastings. We apply the modeling approach in a system composed by a single watershed, reservoir, and a channel connected in series, respectively, for the San Antonio observed rainfall data. Results indicate that for flood mitigation, the predictive control strategy outperforms the reactive controls not only when applied for synthetic design storm events, but also for a continuous simulation. Moreover, the predictive control strategy requires smaller valve operations, while still guaranteeing efficient hydrological performance. From the results, we recommend the use of the model predictive control strategy to control stormwater systems due to the ability to handle different objective functions, which can be altered according to rainfall forecasting and shift the reservoir operation from flood-based control to strategies focused on increasing detention times, depending on the forecasting.

Chlorine is a widely used disinfectant and proxy for water quality (WQ) monitoring in water distribution networks (WDN). Chlorine-based WQ regulation and control aims to maintain pathogen-free water. Chlorine residual evolution within WDN is commonly modeled using the typical single-species decay and reaction dynamics that account for network-wide, spatiotemporal chlorine concentrations only. Prior studies have proposed more advanced and accurate descriptions via multi-species dynamics. This paper presents a host of novel state-space, control-theoretic representations of multi-species water quality dynamics. These representations describe decay, reaction, and transport of chlorine and a fictitious reactive substance to reflect realistic complex scenarios in WDN. Such dynamics are simulated over space- and time-discretized grids of the transport partial differential equation and the nonlinear reaction ordinary differential equation. To that end, this paper (i) provides a full description on how to formulate a high fidelity model-driven state-space representation of the multi-species water quality dynamics and (ii) investigates the applicability and performance of different Eulerian-based schemes (Lax-Wendroff, backward Euler, and Crank- Nicolson) and Lagrangian-based schemes (method of characteristics) in contrast with EPANET and its EPANET-MSX extension. Numerical case studies reveal that the Lax-Wendroff scheme and method of characteristics outperform other schemes with reliable results under reasonable assumptions and limitations.

Real-time water quality (WQ) sensors in water distribution networks (WDN) have the potential to enable network-wide observability of water quality indicators, contamination event detection, and closed-loop feedback control of WQ dynamics. To that end, prior research has investigated a wide range of methods that guide the geographic placement of WQ sensors. These methods assign a metric for fixed sensor placement (SP) followed by \textit{metric-optimization} to obtain optimal SP. These metrics include minimizing intrusion detection time, minimizing the expected population and amount of contaminated water affected by an intrusion event. In contrast to the literature, the objective of this paper is to provide a computational method that considers the overlooked metric of state estimation and network-wide observability of the WQ dynamics. This metric finds the optimal WQ sensor placement that minimizes the state estimation error via the Kalman filter for noisy WQ dynamics -- a metric that quantifies WDN observability. To that end, the state-space dynamics of WQ states for an entire WDN are given and the observability-driven sensor placement algorithm is presented. The algorithm takes into account the time-varying nature of WQ dynamics due to changes in the hydraulic profile -- a collection of hydraulic states including heads (pressures) at nodes and flow rates in links which are caused by a demand profile over a certain period of time. Thorough case studies are given, highlighting key findings, observations, and recommendations for WDN operators. Github codes are included for reproducibility.

State estimation (SE) of water distribution networks (WDNs) is difficult to solve due to nonlinearity/nonconvexity of water flow models, uncertainties from parameters and demands, lack of redundancy of measurements, and inaccurate flow and pressure measurements. This paper proposes a new, scalable successive linear approximation to solve the SE problem in WDNs. The approach amounts to solving either a sequence of linear or quadratic programs---depending on the operators' objectives. The proposed successive linear approximation offers a seamless way of dealing with valve/pump model nonconvexities, is different than a first order Taylor series linearization, and can incorporate with robust uncertainty modeling. Two simple testcases are adopted to illustrate the effectiveness of proposed approach using head measurements at select nodes., Comment: Published in IEEE 58th Conference on Decision and Control, Nice, France, December 11--13, 2019