Publications

Different methods of two-dimensional and three-dimensional numerical resolution models have been used to predict the air–water interaction in pipe systems in the early twenty-first century, where reliable and adequate results have been obtained when compared with experimental results. However, the study of the drainage process in pressurized systems with air admitted through openings has not been studied using this type of model due to the complexity that this represents. In this research, a two-dimensional numerical model is developed in the open-source software OpenFOAM; this model represents the drainage of an irregular pipe with air admitted by an air valve, defined by a structured mesh. A validation of the numerical model related to the air admitted by the variation of the air valve diameter is also performed.

Entrapped air pockets in water pipelines play a significant role in influencing transient over-pressures during filling procedures. Several research is focused on highlighting the attenuation of pressure peaks in pipes with single air pockets. This research studies the air-water interaction during rapid water filling processes in an irregular pipeline and air pockets in different branches, and how the trapped air can attenuate the over-pressure peaks. A three-dimensional computational fluid dynamics (CFD) model was developed, and numerical results of the model were validated through experimental measurements. For a given initial air pocket condition upstream of the high point, the maximum air pocket over-pressure was 11% to 32% lower when the descending pipe segment initially contains air compared to when it contains water. In sum, it was found that entrapped air pockets at high points of water pipelines can help mitigate transient over-pressures considering specific initial hydraulic conditions prior to filling operations.

Some countries grapple with data scarcity for calibration purposes when establishing current hydrodynamic models, which often require many parameters. In this context, this research presents a practical simulation methodology for hydrodynamic modelling suitable for application in bay and estuarine systems based on mass and momentum equations and requiring only one parameter for calibration—bed friction. The proposed simulation methodology is applied to a linear open channel measuring 200,000 m long. A sensitivity analysis of the bed friction is conducted to assess the proposed methodology’s response to the maximum water levels achieved. The results are compared to linear theory, indicating that the proposed simulation methodology effectively represents the water phase. In all simulations, the maximum root mean square error is less than 2.1% when neglecting bed friction and 4.69% when a bed friction of 0.005 is considered. The proposed simulation methodology can be a practical tool for hydrodynamic modelling in shallow waters.

Filling and emptying processes are common maneuvers while operating, controlling and managing water pipeline systems. Currently, these operations are executed following recommendations from technical manuals and pipe manufacturers; however, these recommendations have a lack of understanding about the behavior of these processes. The application of mathematical models considering transient flows with entrapped air pockets is necessary because a rapid filling operation can cause pressure surges due to air pocket compressions, while an uncontrolled emptying operation can generate troughs of sub-atmospheric pressure caused by air pocket expansion. Depending on pipe and installation conditions, either situation can produce a rupture of pipe systems. Recently, reliable mathematical models have been developed by different researchers. This paper reviews and compares various mathematical models to simulate these processes. Water columns can be analyzed using a rigid water column model, an elastic water model, or 2D/3D CFD models; air–water interfaces using a piston-flow model or more complex models; air pockets through a polytropic model; and air valves using an isentropic nozzle flow or similar approaches. This work can be used as a starting point for planning filling and emptying operations in pressurized pipelines. Uncertainties of mathematical models of two-phases flow concerning to a non-variable friction factor, a polytropic coefficient, an air pocket sizes and an air valve behavior are identified.

Accurately estimating storm hydrographs for various return periods is crucial for planning and designing hydrological infrastructure, such as dams and drainage systems. A key aspect of this estimation is the separation of baseflow from storm runoff. This study proposes a method for deriving storm hydrographs for different return periods based on hydrological station records. The proposed approach uses three baseflow separation methods: constant, linear, and master recession curve. A significant advantage of the proposed method over traditional rainfall–runoff approaches is its minimal parameter requirements during calibration. The methodology is tested on records from the Lengupá River watershed in Colombia, using data from the Páez hydrological station, which has a drainage area of 1090 km2. The results indicate that the linear method yields the most accurate hydrograph estimates, as demonstrated by its lower root mean square error (RMSE) of 0.35%, compared to the other baseflow separation techniques, the values of which range from 2.92 to 3.02%. A frequency analysis of hydrological data was conducted using Pearson Type III and Generalized Extreme Value distributions to identify the most suitable statistical models for estimating extreme events regarding peak flow and maximum storm hydrograph volume. The findings demonstrate that the proposed methods effectively reproduce storm hydrographs for return periods ranging from 5 to 200 years, providing valuable insights for hydrological design, which can be employed using the data from stream gauging stations in rivers.

The computation of hydrographs in large watersheds necessitates utilizing channel routing, which calculates the movement of hydrographs along channel branches. Routing methods rely on an implicit scheme to facilitate numerical resolution, which requires more computational time than the explicit scheme. This study presents an explicit scheme channel routing model that offers a versatile approach to open channel flow analysis. The model is based on mass conservation principles and Manning equations, and it can accommodate varying bed slopes, making it highly adaptable to diverse hydraulic scenarios. In addition, the proposed model considers backwater effects, which enhances its applicability in practical scenarios. The model was tested in a practical application on a rectangular channel with a width of 7 m, and the results showed that it can accurately predict outflow hydrographs and handle different flow conditions. Comparative analyses with existing models revealed that the proposed model’s performance in generating water flow oscillations was competitive. Moreover, sensitivity analyses were performed, which showed that the model is highly responsive to parameter variations, such as Manning’s coefficient, bed slope, and channel width. The comparison of peak flows and peak times between the proposed model and existing methods further emphasized the model’s reliability and efficiency in simulating channel routing processes. This research introduces a valuable addition to the field of hydrology by proposing a practical and effective channel routing model that integrates essential hydraulic principles and parameters. The results of the proposed model (lumped routing) are comparable with the solution provided by the Muskingum–Cunge method (distributed routing). It is of utmost importance to note that the proposed model applies to channel branches with bed slopes below 6°.

Air pockets in water distribution networks can cause various operational issues, as their expansion during drainage operations leads to sub-atmospheric conditions that may result in pipeline collapse depending on soil conditions and pipe stiffness. This study presents an analytical solution for calculating air pocket pressure, water column length, and water velocity during drainage operations in a pipeline with an entrapped air pocket and a closed upstream end. The existing system of three differential equations is reduced to two first-order nonlinear differential equations, enabling a rigorous analysis of the existence and uniqueness of solutions. The system is then further reduced to a single secondorder nonlinear ordinary differential equation (ODE), providing an intuitive framework for examining the physical behaviour of the hydraulic and thermodynamic variables. Furthermore, through a change of variables, the second-order ODE is transformed into a first-order linear ODE, facilitating the derivation of an analytical solution. The analytical solution is validated by comparing it with a numerical solution. Additionally, a practical application demonstrates the effectiveness of the developed tool in predicting the extreme pressure values in the air pocket during the water drainage process in a pipe, within a controlled environment.

Entrapped air pockets can cause failure in water distribution systems if air valves have not been appropriately designed for expelling air during filling manoeuvres performed by water utilities. One-dimensional mathematical models recently developed for studying this phenomenon do not consider the effect of blocking columns inside water pipelines. This research presents the development of a mathematical model for analysing the filling process in a pipeline with an undulating profile with various air valves, including blocking columns during starting-up water installations. The results show how different air pocket pressure peaks can be produced over transient events, which need to be analysed to ensure a successful procedure that guarantees pipeline safety during the pressure surge occurrence. In this study, an experimental set-up is analysed to observe the behaviour of two blocking columns during filling by comparing the air pocket pressure pulses.

The occurrence of sub-atmospheric pressure in the drainage of pipelines containing an air pocket has been known as a major cause of several serious problems. Accordingly, some system malfunction and pipe buckling events have been reported in the literature. This case has been studied experimentally and numerically in the current research considering objectives for a better understanding of: (i) the emptying process, (ii) the main parameters influencing the drainage, and (iii) the air-water interface deformation. Also, this research demonstrates the ability of a computational fluid dynamic (CFD) model in the simulation of this event. The effects of the air pocket size, the percentage and the time of valve opening on the pressure variation have been studied. Results show the pipeline drainage mostly occurs due to backflow air intrusion. The worst case scenario is associated with a fast valve opening when a tiny air pocket exists in the pipeline.

The time of concentration is the time it takes a drop of water in a basin to travel from the most distant point to the outlet, and is one of the most important parameters, along with the morphometric characteristics, for determining the design flow rate in rainfall-runoff models. This study aims to determine the sensitivity of the parameters included in different equations for the calculation of the time of concentration. A case study was conducted on small, urbanized watersheds in the city of Montería, Colombia. The study uses information obtained through field work using GPS equipment and electronic total station, supplemented by geographic information contained in the city drawings of the local sewage company, which includes data on elevations above sea level with sub-metric precision. The time of concentration determined by the 12 empirical equations was compared to the results obtained from the equation proposed by the Natural Resources Conservation Service (NRCS), which was considered as a baseline formulation for the intricacy of calculation. Based on this comparison, it was found that the Carter equation is the one that best fits the results obtained from the NRCS equation because it displayed highly significant goodness of fit values. Even though the equations by Kirpich, Ventura, California Culvert Practice, Simas-Hawkins and TxDOT provide a relatively good fit compared to other empirical equations, they tend to over-estimate time of concentration values, which could lead to the under-estimation of the design flow rates. For this reason, sensitivity analysis of the parameters of these equations represents an alternative for improving the calculation of the time of concentration. The current research analyses deepen the influence of some parameters in the estimation of time of concentration. The research can also be used by designers and engineers in the city of Montería, Colombia, as an important reference to compute time of concentrations in urbanized watersheds.

For the design of stormwater pumping stations, there is often uncertainty regarding the selection of an appropriate rainfall event to determine the required pumping capacity and temporary storage volume for managing extreme events of a given magnitude. To account for the risk of system failure, the return period is considered, as recommended based on the size of the catchment’s drainage area or other considerations, depending on the local regulations of a country. This study focused on analysing the direct runoff volume from the catchment, the storage volume required for the operation of the pumping system, and the order of magnitude of the design flow rate. The results indicate that a rainfall event with a duration of at least twice the time of concentration should be used. The design flow rate should range between 50% and 70% of the peak discharge, and designing for flow rates near the peak is not advisable, as it can lead to intermittent pump operation and result in an oversized installed capacity. The methodology developed in this research was applied to the Coastal Protection Project located in the city of Cartagena, Colombia, which includes a 2045.6-m-long box culvert with a cross-sectional area of 2 × 2 m, and three pumping stations, each equipped with three pumps rated at 0.75 m3/s, for a total installed capacity of 6.75 m3/s.

The conventional practice in the design of storm drainage systems is based on statistically stationary load and resistance conditions that remain invariant over time. However, uncertainties in the variables affect the design accuracy and the satisfactory performance of these hydrosystems during their operation and service. To overcome this limitation, a design methodology for a storm drainage channel was proposed using a probabilistic framework that incorporates uncertainty analysis of random variables and estimates the system’s probability of failure in terms of design depth and maximum allowable velocity. This methodology employs the Monte Carlo simulation technique and offers an alternative design and analysis approach to strengthen the conventional sizing method for drainage channels in urban watersheds. Based on uncertainty criteria associated with hydraulic design, operation, and prospective changes in the watershed and the channel, appropriate dimensions were estimated regarding design depth and freeboard. The results of this study demonstrate that the annual probability of failure of a channel, when considering uncertainty, is significantly higher than the yearly exceedance probability associated with the hydrological design return period event. Therefore, the proposed methodology is appropriate for estimating the system’s capacity and potential failure risk. This methodology may also be applied to sizing other stormwater drainage structures.

Pipeline filling and emptying are critical hydraulic procedures involving transient two-phase air–water interactions, which can cause pressure surges and structural risks. Traditional Digital Twin models rely on one-dimensional (1D) approaches, which cannot capture air–water interactions. This study integrates Computational Fluid Dynamics (CFD) models into a Digital Twin framework for improved predictive analysis. A CFD-based Digital Twin is developed and validated using real-time pressure measurements, incorporating 2D and 3D CFD models, mesh sensitivity analysis, and calibration procedures. Key contributions include a CFD-driven Digital Twin for real-time monitoring and machine learning (ML) techniques to optimise pressure surges. ML models trained with experimental and CFD data reduce reliance on computationally expensive CFD simulations. Among the 31 algorithms tested, decision trees, efficient linear models, and ensemble classifiers achieved 100% accuracy for filling processes, while k-Nearest Neighbours (KNN) provided 97.2% accuracy for emptying processes. These models effectively predict hazardous pressure peaks and vacuum conditions, confirming their reliability in optimising pipeline operations while significantly reducing computational time.

Tidal inlet structures are engineering projects with associated benefits related to flood control, water quality enhancement, and coastal protection. This study analyzes the performance of hydraulic gates on a stabilized inlet in estuarine systems by developing a simplified hydraulic model that considers inlet and outlet water levels. The proposed model was applied to the stabilized tidal inlet structure in Cartagena de Indias, Colombia. This model offers a practical tool for engineers and designers operating estuarine systems. The analysis focuses on the coastal lagoon of Ciénaga de la Virgen. The proposed model was successfully calibrated using two water sensors, with extreme input and outlet flow rates of approximately 260 m3/s and 110 m3/s, respectively. The average daily output volume in the system is 3,361,000 m3, while the average daily input volume is 3,200,000 m3. Consequently, the manipulation of the opening gates results in a decrease in the estuarine water level, potentially by as much as 25 cm, which local authorities can use to make decisions to reduce extreme water levels during flooding events.