Introducing Excel Spreadsheet Calculations and Numerical Simulations With Professional Software Into An Undergraduate Hydraulic Engineering Course
Introducing Excel Spreadsheet Calculations and Numerical Simulations With Professional Software Into An Undergraduate Hydraulic Engineering Course
Introducing Excel Spreadsheet Calculations and Numerical Simulations With Professional Software Into An Undergraduate Hydraulic Engineering Course
DOI: 10.1002/cae.22185
RESEARCH ARTICLE
1
Environmental Fluid Dynamics Group,
Andalusian Institute for Earth System Abstract
Research, University of Granada, This paper presents a comprehensive set of computer tools for implementing
Granada, Spain
an active teaching–learning strategy on an undergraduate hydraulic engineer-
2
Department of Civil Engineering,
ing course, within the subject blocks of pressure flow and free surface flow.
Environmental and Water Engineering
Group, University of A Coruña, A The tools comprise a number of simple Excel‐based applications, to enhance
Coruña, Spain students’ understanding of basic concepts, plus a series of practical exercises
Correspondence
using the free professional models EPANET and IBER, aimed at familiarizing
María Bermúdez, Andalusian Institute for students with the tools used in real practice. The application of the new
Earth System Research, Av. Del methodology to a hydraulics course at the Civil Engineering School of the
Mediterráneo s/n, 18006 Granada, Spain.
Email: mariabermudez@ugr.es; University of A Coruña (Spain) resulted in increased student satisfaction and
mbermudez@udc.es engagement, as well as greater interaction with instructors and classmates.
The applications and practical exercises are easily adaptable and freely
Funding information
H2020 Marie Sklodowska‐Curie Actions, available to the community to facilitate their application to other courses and/
Grant/Award Number: 754446 or degrees.
KEYWORDS
engineering education, EPANET model, hydraulic engineering, IBER model, Microsoft‐Excel
spreadsheet
laboratory practical work. On previous courses teachers programming, but this is not what we are looking for
have been able to verify the usefulness of these as a here.
means of understanding and reinforcing complex con- With the professional models, simple cases are chosen
cepts learned from a theoretical perspective [4]. (d) Use of for a first approximation to their use in real practice.
computer applications, be they developed by instructors Working with relatively simple examples also allows
or professional ones, adapted to the knowledge and skills students to critically evaluate their results and thus
of the students. These final activities, sometimes linked to avoids these results becoming a “black box.”
practical laboratory work, are the main subject of the
present article. 3.1 | Applications associated with
pressure flow
FIGURE 1 Calculation of head loss along a pipe system with three segments: (a) Discharge calculation, (b) pressure calculation, and
(c) diagram of the system
defined with three points, on which a second degree analysis of pipe segments is completed with a spreadsheet
function is interpolated (this being the usual practice). for the coupling of pumps in series or in parallel.
The inlet and outlet pressures of the pump are also The next and final step in this set of simple
analyzed for different locations, this to illustrate the applications is the analysis of networks, which is
concept of net positive suction head and to detect illustrated with the problem of three tanks (Figure 4).
whether cavitation problems may occur. Using Excel’s optimization procedures, the equation of
We work with a spreadsheet in which the pump energy conservation in the node is imposed, maintaining
characteristic curves are synthetic (they do not corre- a consistent distribution of flows. In this application,
spond to real pumps), to understand the concepts, and which is the preamble to the next block, where a
subsequently with another spreadsheet in which a pump professional model will be used, the effect of time is also
must be chosen from those available in a commercial introduced. The circulating flows fill and empty the
catalog (Figure 3), to understand that the options are tanks, and this modifies the distribution of flows, in a
limited by the market (although they are wide). The similar way to how the EPANET model works.
FIGURE 2 Calculation of a pipe network with pump: (a) Calculation of the system curve, (b) pump characteristic curve considered
(synthetic) and resulting system curve, (c) operating point obtained, and (d) diagram of the system
BERMÚDEZ ET AL. | 197
FIGURE 3 Pump selection for a pipe system using a pump catalog: (a) Selection of a typology, and (b) choice of a specific pump
considering the set of commercial curves for that typology
FIGURE 4 Branched pipe system with three reservoirs: (a) Calculation of flows for each branch, and (b) diagram of the system
These simple models are used in parallel with the Once students have assimilated and managed the
theoretical training and the problems that are solved simple concepts, it is important that they manage a
manually, but also with laboratory practices on the pipe professional program for the calculation of pressure
network panel in the UDC School of Civil Engineering’s pipelines. The current standard is EPANET (https://
laboratory, which has pipes of a certain magnitude www.epa.gov/water‐research/epanet), which is very well
(Figure 5) and equipment to evaluate the performance of known in the engineering community. Students receive
pumps. The pumps included in the computer application interactive training on this program, analyzing basic
correspond to the same series as those available in the problems, such as some of those mentioned in the
laboratory. previous paragraphs, and subsequently they are asked to
FIGURE 5 Diagram of the small‐scale pipe network panel in the University of A Coruña School of Civil Engineering’s laboratory
198 | BERMÚDEZ ET AL.
FIGURE 6 Photograph of the full‐size pipe network in the Center for Technological Innovation in Construction and Civil Engineering
(left). Graphical user interface of the remote control application
model a system of pipelines with real characteristics. For calculated with Allievi’s formula) and slow closure (in
this, the existing research infrastructure of the CITEEC is which the pressure peaks are attenuated).
used, which includes a cast iron pipeline with similar These two spreadsheets are “playable,” in the sense
characteristics to that of the water supply of the city of A that one can modify all the parameters and observe the
Coruña (Figure 6). distribution of pressures in very different configurations.
This installation is used by the GEAMA research As a final point to this block, a laboratory practice is
group and by the water supply company of the city carried out in the hydraulic laboratory of the School of
(EMALCSA) for development work and for the calibra- Civil Engineering where a rapid closure is generated in a
tion of equipment, and is also used by students for PVC pipeline by means of an electrovalve. The pressure
working on pipes with a real appearance. The pipe distribution is analyzed, it being possible here to estimate
network is amply instrumented with pressure sensors parameters such as the velocity of the pressure wave. The
and flow meters. The purpose of the laboratory practice is entire process is tele‐controlled and the pressure field is
to create the pipe network in EPANET (the pipeline’s exported to a text file where students can analyze the
dimensions are indicated in the infrastructure itself), to results (Figure 9). In addition, the validity of the test can
establish flow conditions (which can be modified), and to be verified by means of the graphic interface presented in
check that the flow regime and the distribution of the figure, where the pressure peaks and their attenua-
pressures modeled coincide with reality. The experimen- tion over time can be seen.
tal facility has a computer platform for remote control
and operation (Figure 6).
The final topic related to the pressure flow is the 3.2 | Applications related to free surface
analysis of hydraulic transients. This is very complex and flow
requires all the students’ attention. In addition to
studying it theoretically and through class discussions This introductory subject explains in some detail perma-
among students, an analysis is made through some nent and uniform movement, as well as gradually varied
simple applications of the compression and decompres- permanent flow (backwater curves). For this purpose, the
sion wave fronts in a simple pipe, in several locations. In theoretical materials are accompanied by small Excel
the first of the applications, we consider the effect of the applications, which allow for the reproduction of results
total closure of a valve both instantaneously and as the observed in the laboratory channel of the School of Civil
sum of partial closures, so that the cumulative effect of Engineering (Figure 10). This channel has a length of
these partial closures can be perceived, which in some 15 m, a width of 50 cm, and a maximum depth of 50 cm
cases have an additive effect and in others a canceling and both the flow and the slope can be varied. In addition,
effect (Figure 7). it has a variable level reservoir at its downstream end that
The second application (Figure 8), in which linear allows various boundary conditions to be imposed to
valve closures are considered, these being essentially simulate different types of profile on the free surface.
staggered closures with a very small time increment, The first application presents the basic characteristics
seeks to illustrate the concepts of rapid closure (in which of a channel with different geometries, such as its normal
the pressure increase reaches the maximum that can be depth, its critical depth, and variables derived from these,
BERMÚDEZ ET AL. | 199
F I G U R E 7 Water hammer calculation considering a stepwise closure: (a) Input data and main calculations, (b) diagram of the system,
and (c) pressure fluctuations caused by the closure
such as its velocity or the conjugate depth of the normal The optimal shape of a channel section is the one that
depth, the meaning of which will be understood during transports a maximum flow for a given cross‐sectional
the final classes of the course. Rectangular, triangular, area. This means that the wetted perimeter is minimal. In
trapezoidal, and circular sections are proposed. addition to the exploration of theoretical issues, students
This spreadsheet allows for interaction with the are provided with a small Excel application that allows
laboratory channel, imposing on it a constant depth by them to optimize sections for the four geometries
manipulating the boundary condition, and verifying that indicated (Figure 12), and to verify that the optimum
this depth effectively fulfills the equation of the normal shapes are the half of a square, in the rectangular and
depth for different flows and slopes. Because the section triangular cases, the half of a hexagon, in the trapezoidal
of the laboratory channel is rectangular, the tab case, and the semicircular one in the circular case [12].
corresponding to this geometry is used, and is reproduced To do this, Excel optimization tools are used (Goal Seek
in Figure 11. and Solver).
FIGURE 8 Water hammer calculation considering a linear valve closure: (a) Input data and main calculations, (b) diagram of the
system, and (c) pressure changes caused by the closure
200 | BERMÚDEZ ET AL.
FIGURE 9 Laboratory experiment used to generate and visualize water hammer in pipes
The last simple application is the one that allows drawing regime), and the integration interval is also indicated. The
backwater curves, and is used to reproduce the curves that application integrates the equation of the backwater curves
are observed in the laboratory channel (Figure 13). This by means of a fourth‐order Runge–Kutta algorithm, which
application therefore makes it possible to “test” (in is explained in the spreadsheet (the calculation columns
simulation form) before, during, and after the practical are hidden in Figure 13 but can be displayed) and which
work in the laboratory. Students set different slopes, flows has been previously analyzed in the theoretical classes.
and initial conditions, calculate the normal and critical The result is a graph showing the bottom of the
depths, and draw on paper the curves observed in the channel, the normal depth, the critical depth, the free
laboratory, extracting from them a maximum of parameters surface profile, and each of the conjugate depths. The
(the critical depth if there is a free fall condition, or the application makes it possible to include the measured
normal depths if a constant depth is reached, for example). data in the channel and to check the average quadratic
At the same time, the situation is reproduced on a error between the measured data and the calculated data,
screen in the laboratory itself using the application, in and to calculate the Manning coefficient that optimizes
which the flow parameters are set, the direction of the adjustment, using the solver tool.
integration is indicated (from downstream to upstream in The determination of the conjugate depths in a
subcritical flow regime, and vice versa in supercritical flow hydraulic jump, which is measured in the laboratory,
can also be analyzed by means of this application, in applications, a step is taken towards professional models
which the supercritical and subcritical reaches are through the presentation and use of the two‐dimensional
reproduced, as well as their respective backwater (2D) hydraulic model Iber (http://www.iberaula.es/). The
curves and their conjugate depths, plus the point of Iber model allows the analysis of rivers reaches or
compatibility. estuaries and has thousands of registered users. It is
In a parallel way to that proposed in the case of developed by the UDC and the Polytechnic University of
pressure flow, once the basic concepts have been Catalonia, together with the Center for Studies and
assimilated in the laboratory and by means of simple Experimentation on Public Works (CEDEX) and the
FIGURE 12 Best hydraulic sections for four cross‐section shapes: rectangle, triangle, trapezoid, and semicircle
202 | BERMÚDEZ ET AL.
F I G U R E 1 3 Calculation of backwater curve observed in the laboratory experiment: (a) Input data, (b) plotted results, and (c) numerical
calculation of the water surface profile (partial view of the table with calculations)
FIGURE 14 Calculation of backwater curves in a rectangular channel with four sections of different slopes
These positive results have encouraged the authors to available to the community at Zenodo (https://doi.org/10.
share their experiences and extend them to other courses 5281/zenodo.3522670). The two professional software
(e.g., within the Civil Engineering Degree at the packages used in the examples are distributed free of
University of Granada). Introductory courses in fluid charge and can be obtained from their respective
mechanics are common to various university degrees websites. The aim here is to eliminate the external
within engineering (e.g., mechanical, civil, or chemical). barriers usually encountered during the integration of
However, it is undeniable that the introduction of new technology in the classroom, such as time constraints and
active methodologies in the university classroom requires lack of resources [16]. Implementing the methodology
an effort on the part of teaching staff, who must produce does not entail any economic cost and the tools provided
new high‐quality content adapted to the new teaching can be used “as is.”
strategy, often with limited institutional support. To It should be noted that the material provided
facilitate its application or adaptation to other courses constitutes a comprehensive set of computer tools
and/or degrees, the tools developed have been made that covers the basics of pressure flow and free surface
204 | BERMÚDEZ ET AL.
FIGURE 15 Calculation of backwater curves in a rectangular channel with three sections of different slope and width, and with a step
flow, starting from very simple applications based on these tools into an undergraduate hydraulic engineer-
Excel up to numerical simulations with professional ing course. Following the proposed methodology, it is
software. All applications are easily adaptable, and also possible to incorporate other problems or tools as
geometries can be easily varied according to the a means of exploring in greater depth those topics
laboratory facilities available, or new cases can be that may be of more interest, according to the focus of
generated from them. It is also worth noting that this the subject (e.g., the study of branched systems
study goes beyond the development of a set of of pipes [13,26,27,34], or the calculation of transients
computer tools by designing and implementing an in pipes with the professional tool Allievi [www.
innovative teaching–learning strategy that introduces allievi.net]).
BERMÚDEZ ET AL. | 205
5 | C ON C LU S I O N S RE FER E NCES
1. M. Álvarez et al., Two‐dimensional dam‐break flood analysis in
This paper presents an active teaching–learning strategy of
data‐scarce regions: The Case study of Chipembe Dam, Mozambi-
hydraulics for engineering students within an under- que, Water 9 (2017), no. 6, 432. https://doi.org/10.3390/w9060432
graduate curriculum, and a comprehensive set of computa- 2. L. W. Anderson, and D. R. Krathwohl, A Taxonomy for Learning,
tional tools for its application. The tools comprise a number Teaching, and Assessing: A Revision of Bloom’s Taxonomy of
of simple Excel‐based applications, plus a series of practical Educational Objectives, Longman, New York, NY, 2001.
exercises using the free professional models EPANET and 3. M. Bermúdez et al., Quantifying local rainfall dynamics and
IBER. The Excel applications cover the following content uncertain boundary conditions into a nested regional‐local flood
modeling system, Water Resour. Res. 53 (2017), no. 4,
within the study of pressure flow and free surface flow: (a)
2770–2785. https://doi.org/10.1002/2016WR019903
Evaluation of friction and minor losses in a pipe segment, 4. M. Bermúdez et al., Laboratorio virtual para la adquisición de
(b) selection of pumps (and coupling of pumps), (c) analysis conceptos umbrales en ingeniería hidráulica, Contextos Uni-
of pipe networks, (d) analysis of hydraulic transients, (e) versitarios Tranformadores: Retos e Ideas Innovadoras. II
calculation of normal and critical depths, (f) determination Xornadas de Innovación Docente (E. de la Torre Fernández,
of the optimum hydraulic cross‐section of a channel, and ed.), Servizo de publicacións, 2018, pp. 391–393. Universidade
(g) calculation of free surface profiles in a channel. In da Coruña https://doi.org/10.17979/spudc.9788497496780
5. D. Berrett, How “Flipping” the Classroom Can Improve the
practical exercises with EPANET a system of pipes with real
Traditional Lecture, Educ. Dig. Essent. Read. Condens Quick
characteristics is modeled, while in exercises with IBER the
Rev. 78 (2012), no. 1, 36–41. https://eric.ed.gov/?id=EJ987290.
effect of changes of slope and cross‐section, and of the Accessed April 30, 2019.
presence of steps and other singularities on the flow in 6. E. Bladé et al., Iber: Herramienta de simulación numérica del
rectangular channels is reproduced. flujo en ríos, Revista Internacional de Métodos Numéricos para
The results of this experience show a greater Cálculo y Diseño en Ingeniería 30 (2014), no. 1, 1–10. https://
participation of students, with more interaction with doi.org/10.1016/j.rimni.2012.07.004
instructors, and this of a higher quality, as well as higher 7. B. S. Bloom, Taxonomy of Educational Objectives: The
Classification of Educational Goals. Handbook 1: Cognitive
motivation for learning. The perception of a more
Domain, Longmans, Green & Co., London, 1956.
pleasant experience compared with the traditional class
8. J. M. Bodoque et al., Improvement of resilience of urban areas by
model was also found. integrating social perception in flash‐flood risk management, J.
Hydrol. 541 (2016), 665–676. https://doi.org/10.1016/J.
JHYDROL.2016.02.005
ACKNOWLEDGMEN TS
9. D. Brkic, Solution of the implicit colebrook equation for flow
The authors would like to thank Dr Félix Sánchez‐ friction using Excel, Spreadsheets Educ. 10 (2017), no. 2, 4663.
Tembleque Díaz‐Pache, and Dr Luís Pena for their help in https://sie.scholasticahq.com/article/4663‐solution‐of‐the‐
implicit‐colebrook‐equation‐for‐flow‐friction‐using‐excel.
the laboratory practices. M. B. gratefully acknowledges
Accessed September 4, 2019.
funding from the European Union’s Horizon 2020 10. L. Cea, M. Bermúdez Pita, and B. Sobral Areán, Cálculo de
research and innovation program under the Marie Curvas de Remanso y Fenómenos Locales Con Iber, University
Skłodowska‐Curie grant agreement No. 754446 and UGR of Coruña, Coruña, Spain, 2018, https://ruc.udc.es/dspace/
Research and Knowledge Transfer Found—Athenea3i. handle/2183/20804?locale‐attribute=es. Accessed April 26, 2019.
11. L. Cea et al., IberWQ: New simulation tool for 2D water quality
modelling in rivers and shallow estuaries, J Hydroinformatics. 18
DATA A VA I L A B I L I T Y ST A T EM E N T (2016), no. 5, 816–830. https://doi.org/10.2166/hydro.2016.235
12. V. Te Chow, Open‐Channel Hydraulics, McGraw‐Hil, New
The MS Excel spreadsheets developed in this study are
York, NY, 1959.
openly available in Zenodo at https://doi.org/10.5281/ 13. S. Demir, N. Manav Demir, and A. Karadeniz, An MS Excel tool
zenodo.3522670. Detailed instructions on how to use them for water distribution network design in environmental engineer-
are also provided. The two professional software packages ing education, Comput. Appl. Eng. Educ. 26 (2018), no. 2,
used in the examples are distributed free of charge and can 203–214. https://doi.org/10.1002/cae.21870
be obtained from their respective websites. 14. S. Demir et al., An MS excel add‐in for teaching hydraulics of
pipe flow in engineering curricula, Comput. Appl. Eng. Educ. 26
(2018), no. 3, 449–459. https://doi.org/10.1002/cae.21898
ORCID 15. T. F. Driscoll, Flipped learning and democratic education. Grad-
uate thesis Teachers College, Columbia (2012).
María Bermúdez http://orcid.org/0000-0003-3189-4791 16. P. A. Ertmer, Addressing first‐ and second‐order barriers to
Jerónimo Puertas http://orcid.org/0000-0001-6502-0799 change: Strategies for technology integration, Educ. Tech. Res.
Luis Cea http://orcid.org/0000-0002-3920-0478 Dev. 47 (1999), no. 4, 47–61. https://doi.org/10.1007/
BF02299597
206 | BERMÚDEZ ET AL.