Manual Programa ALLIEVI
Manual Programa ALLIEVI
Manual Programa ALLIEVI
ITA
Universidad Politécnica de Valencia
Camino de Vera s/n - Edificio 5C 46022
Valencia (Spain)
Ph: (+34) 96 387 98 98
3 USERS MANUAL
System requirements
Intel ® Core Duo processor or higher (or equivalent processor in the case
of AMD);
Windows 7
At least 2 GB of RAM memory;
At least 100 MB available hard disk space;
1.280 x 1.024 resolution or higher, with 32 bits color (recommended)
Internet connection. Broadband recommended (1)
Adobe® Reader or compatible
.NET Framework 3.5
(1) This software requires internet connection. When executed for the first time, ALLIEVI will
have to be activated since it’s not possible to use ALLIEVI without activation. To do so, an
internet connection will be needed to validate your ALLIEVI copy with our servers and, in
addition, a periodic re-validation we’ll be required.
Registration
During the ALLIEVI installation process, the system will require some license
information. If this information is not provided, it won’t be possible to use the
program.
Activation
When executed for the first time, the program will have to be activated. To do so,
you’ll need an internet connection to validate your copy with our servers at ITA.
8 USERS MANUAL
1. Introduction
This ALLIEVI technical users guide has been formulated with the intention of
offering the user the necessary tools and concepts for the adequate usage of the
program, specifically for those who have basic hydraulic knowledge or have
experience in the usage of hydraulic transient (water hammer) simulation
computer programs. We advise you to follow the tutorial videos available at
www.allievi.net, as well as other useful network resources such as FAQs, also
available at this webpage.
To regulate the circulating flow through the conduits, the system may have
different types of valves installed either in line with these pipes, or with an exterior
discharge. Some of these valves can additionally function as protection systems
against hydraulic transients.
In hydraulic systems there are also other actions which can provoke transients,
such as the change in a tank´s level, the alterations in demand of a water supply
network and the breakage of a conduit at a given time. Of these three doings, the
first two are usually considered to be slow or very slow maneuvers, thus resulting
in minor transients. However, the breakage of a pipe may result in a rapid
transient, reaching its maximum when the breakage is significant and occurs in a
relatively short period of time.
9 USERS MANUAL
During the transient, and while the system is adapting from one steady state
regime to the other, there’s an exchange between pressure and kinetic energy,
and vice versa, with the corresponding oscillations in flow, velocity and pressure
inside the pipes. The oscillations in flow and velocity don’t have a direct
consequence on the system, but the fluctuations in pressure generate pressure
changes between maximum and minimum values. It’s the maximum pressure
values which can break the system’s pipes either if sufficient mechanical resistance
is not offered or if the installation is not appropriately protected.
On the other hand, minimum pressures tend not have any effect on the system, if
kept over atmospheric pressure. Although, if lowered beneath atmospheric
pressure, they can generate a depression inside the conduit which may lead to
the entrance of air or liquid through small cracks or fissures, with the related
danger of pathogen intrusion or conduit collapse if the ratio thickness/diameter
is sufficiently low.
In view of the harmful effects that may arise during the development of a
hydraulic transient, the solution is to, either provide the system with enough
mechanical resistance to bear the pressure oscillations or to protect the
installation in order to diminish the amplitude of such fluctuations and maintain
them under acceptable limits. To proceed in either way we need to make use of
a computer program capable of simulating transients and providing, among other
results, the oscillations in pressure originated after a particular maneuver. This is,
as we’ve already mentioned, the main objective of the ALLIEVI program.
10 USERS MANUAL
Consequently to what we’ve indicated, ALLIEVI is a software that can be used both
at design stages and at stages of verification of a hydraulic system’s functioning.
At design stages, the program allows to incorporate to the system several
protection devices (pressure vessels, surge tanks, relief valves, etc.) so as to
maintain the pressure oscillations caused by the different elements’ maneuvers
between previously determined limits. On the other hand, at verification stages,
the program enables the visualization of the transient effects when effectuating
the maneuvers of a previously designed installation.
H H a 2 V
V 0
t x g x
(1)
V V H VV
V g f 0
t x x 2D
Since this system doesn’t have an analytical solution, in order to integrate it, the
standard procedure is to assume that the solution will be obtained for specific
time instants (separated by a t) and in particular points of the pipe (separated
by a x), satisfying the following condition,
x
a (2)
t
where “a” is the celerity of the pressure wave inside the pipe. This celerity, in the
case of water, is calculated through the expression
9.900
a
D
47'6 C
e (3)
in which the C coefficient depends on the pipe’s material.
11 USERS MANUAL
f n x
H in 1 H in1
g
a n 1
Vi Vi n1 i 1 Vi n1 Vi n1 0
D
(4)
f x n n
n
a n 1
H in 1 H in1 Vi Vi n1 Vi 1 Vi 1 0
i 1
g D
from which we can calculate the hydraulic head H and the velocity V at a point i
of the pipe and in the calculation instant n+1, using the H and V values of the
points i-1 and i+1 at the n instant. Implementing the expressions at (4) is what is
known as the “METHOD OF CHARACTERISTICS”, and is the calculation basic for
the ALLIEVI program.
The system at (4) can be solved for all calculation points of all pipes, with the
exception of end points where we lack one of the two equations. In these
extremes, where supposedly the pipe is connected to a system’s element, the
missing equation is substituted by the equation, or system of equations, which
represent the behavior of said element, and are known as “boundary conditions”.
In this way, the ALLIEVI program enables us to simulate the transient functioning
of a pressurized hydraulic system including the following elements, connected to
each other through their corresponding pipe systems:
y y A V
V 0
t x T x (5)
V V y
V g g s0 s f 0
t x x
where the unknowns are the depth y = y(x, t) and the velocity V = V(x, t), both as
a function of the position along the conduit and time. In the expression at (5), A
is the occupying cross-section of the flow, T is the free surface width defined by
cross-section A, s0 is the bottom slope and sf the hydraulic slope, defined by the
expression
n 2V 2
s f 4/3 (6)
Rh
where n is the Manning coefficient and Rh the hydraulic radius of the fluid-
occupied cross-section. Assuming certain simplifications, system (5) can be solved
similarly to (1) through the method of characteristics, with an analogue
formulation to (4). In this case, the expression at (2) is substituted by
x
V c (7)
t
where c is the celerity of the conduit’s gravity wave, calculated with the expression
A
c g (8)
T
ALLIEVI allows the transient to be simulated in systems where there are both
pressurized and open-channel flows, applying the method of characteristics to
13 USERS MANUAL
both types of conduits, differentiated for each particular case. Since the
characteristic time for solving pressurized transients is in the tenths or hundredths
of a second, being considerably higher for open-channel flow solutions (seconds),
the time increment for the solving of open-channel flow transients must be a
multiple of the time increment for the solving of pressurized flow transients.
In this program, each open flow conduit can only be connected through their
ends to a tank or node linked to other free flow channels. In these conduits,
subcritical flow (or critical at the end point of a freefall conduit) is assumed, taking
rectangular, trapezoidal, circular or any type of cross-section.
In any hydraulic system, the transient regime initiates from a steady state regime
when a maneuver takes place in any of the system’s elements. In the program
ALLIEVI, the steady state regime is obtained by processing the transient regime
calculation using whichever operating conditions, with the system’s elements in
their initial position and without effectuating any maneuvers. When, in this
calculation, unaltered operating conditions through time are obtained, they’re
adopted as starting values to initiate the transient regime calculation.
NOTE: It’s of great importance to point out that in order to calculate any
transient, ALLIEVI calculates first the steady state solution, being available for
the user in table format for each of the project’s elements.
14 USERS MANUAL
Node
2. Tanks
15 USERS MANUAL
3. Pipes
4. Pumps
5. Protection structures
8. Francis Turbines
ALLIEVI’s graphic interface is intended so that the user can introduce, with just
one click, all the necessary elements to create from simple systems, e.g. a simple
impulse, to complex system like mesh networks.
Elements are grouped into five categories. On the one hand we have the basic
elements, which include pumps, pipes, nodes and open-channel flow conduits or
channels. On the other hand, we have the different types of tanks we can simulate.
Finally we have the group of valves, the group of protection elements and the
varied group, which contains flow and head rules, Francis turbines and injectors.
Go to the Image tab Basics Select between Node, Pipe, Pump, Channel.
Go to the Image tab Tanks Select the tank type and its location
16 USERS MANUAL
Go to the Image tab Valves Select the valve type from the possibilities
provided by ALLIEVI (Regulation, retention, automatic, Headloss valve,
among others)
Go to the Image tab Protection Select between Air vessel, Surge tank
or Reservoir.
3.2 Nodes
Once we’ve established each node, ALLIEVI will require the user to introduce the
values corresponding to:
At the nodes, ALLIEVI will show the type of element to which each of them
corresponds. In this way, we can corroborate the elevation values entered for the
respective elements.
NOTE: When two nodes belong to the same element, ALLIEVI will
automatically assign the same elevation value to both of them.
18 USERS MANUAL
3.3 Tanks
Tanks are elements of the system where water at atmospheric pressure is stored.
The types of tanks included on the program are:
Weir
In order to ease the graphical representation of the system, we have tank icons
with a single node (right or left) or with two nodes. In the case of the tank with
two nodes, they must have the same elevation. The elevation of a tank’s nodes
corresponds to its base elevation.
- Initial water level (Z0): It’s the water elevation level used to calculate the
steady state. Being a tank of great dimensions, this elevation will remain
constant throughout the entire simulation, independently of the established
inlet and outlet flows.
To declare a tank of great dimensions and input the initial water level on
ALLIEVI, the procedure is as follows:
The units for Initial water level can be modified by simply right clicking on the
cell once it’s in edit mode.
Right click on Z0 Change unit Units management: this assistant
will allow you to change the current unit to the desired one used in the
field you’re working on.
NOTE The way to change units is analogous for all entry data required for
calculation.
- Initial water level (Z0). It’s the water elevation level used to calculate the
steady state. This being a tank of small dimensions, this level will change
throughout the simulation, depending on the tank’s cross-section and the
established inflows and outflows.
IMPORTANT: In case we ignore the meaning of any entry field, it’s possible to
access an explanation by right clicking on the name of the field you’d like to
know about. With this action, ALLIEVI displays a label with all the necessary
information.
- Initial water level (Z0): It’s the water elevation level used to calculate the
steady state. This level might change along the transient simulation.
The way to introduce the data for a variable cross-section tank is as follows:
Go to Basic Data Tank cross-section S (m2) An assistant will be
displayed where you’ll have to input the elevations and their
corresponding accumulated volumes.
22 USERS MANUAL
3.3.4 Weirs
They’re tanks which have a rectangular transverse thin-walled weir that separates
the entry zone from the exit zone. Its characteristics are:
- Upstream section (Si). It’s the cross section at the entry zone.
- Downstream section (Sf). It’s the cross section at the exit zone.
- Initial and end section elevation (Z Si; Z Sf). It’s the elevation at the initial
and end sections of the tank, used to calculate the steady state solution.
This elevation may be equal or different for both of the tank’s zones, and
might change along the transient simulation.
- Weir length (L). It’s the length of the tank’s transverse weir or the maximum
level of spillway.
- Crest Elevation (Z Crest ). It’s the elevation of the weir’s crest, with regard
to an established reference level.
- Flow coefficient of weir (Cq). It’s the Cq coefficient of the rectangular thin-
walled tank’s weir, used to calculate the discharged flow through the
expression.
2 3
Q Cq L 2g H 2 (9)
3
23 USERS MANUAL
where L is the weir’s length and H the height of the water surface over the
Crest elevation(Z crest), before it initiates its rapid descent. Typical values
for CQ range between 0.60 and 0.80, and a maximun value of 0.85.
section
- Spillway length (L): It’s the length of the spillway. If this length is zero, the
corresponding tank doesn’t have one.
- Crest Elevation (Z Crest). It’s the elevation of the weir’s crest, with regard to
an established reference level.
24 USERS MANUAL
- Flow coefficient of the spillway (Cq). It’s the Cq coefficient of the rectangular
thin-walled or thick-walled tank’s spillway, used to calculate the discharged
flow through the expression at (9).
- Discharge tank: The discharged water through the spillway can be poured to
a receiving waterway or to a small or variable cross-section tank, installed in
the system. Through this variable, we’ll indicate the designation of this tank,
or, if left blank, specify that the pouring is done to a receiving waterway.
In case the tank includes a spillway, such spillway can be fixed or be constituted
by a gate with vertical movement and top pouring. In this last case, this gate might
do a single linear elevation or descent movement during the transient
development. The characteristics of this mobile spillway are the following ones:
- Crest (Fixed - Mobile). It indicates the crest characteristic, either mobile (Yes)
or fixed (No). In case it’s fixed, the following data are not required:
- Final crest elevation (Zf). It’s the elevation that the spillway crest will reach at
the end of its movement. The elevation for the start of the movement is the
previously indicated one through Z Crest.
weirs
In these tanks, as oppose in those of large dimensions, the water level will change
during the transient regime functioning. However, for the steady state calculation,
the water level in these tanks may adopt a predetermined value or vary until the
continuity equation between the outflows and inflows is verified.
In case that, for the steady state calculation, the water level in one of these tanks
is a predetermined value, it’ll result in a steady state regime which generally won’t
fulfill the continuity equation between outflows and inflows.
In most cases, this is done because is interesting to start the transient regime
calculations from certain levels in the system tanks.
If, on the contrary, we allow the water level of whichever of these tanks to fluctuate
during the steady state calculation, when said calculation is done the water level
reached in this tank will be such so as to fulfill the continuity equation between
its inflows and outflows.
- Variable water level for steady state calculation (Yes - No). In order to
calculate the steady state solution, the water inside the deposit might have
variable level (Yes) or fixed level (No), the latter being equal to the initial
level, previously indicated.
3.4 Pipes
- Pipe thickness (e): It refers to the pipe’s wall thickness, used to calculate
celerity.
Colebrook equation). The user must introduce the absolute roughness of the
pipe’s inner walls.
IMPORTANT
If the user introduces in this cell a negative value of roughness, for
instance (-1), the program effectuates calculations considering null
friction losses for this pipe.
Vt 2
hm k (10)
2g
- Speed (a): It’s the pressure wave speed. This datum is directly assigned to the
pipe or can be calculated from the inner diameter, the wall thickness and the
pipe’s material coefficient: If the given value is zero, the transient throughout
the pipe is treated as a mass oscillation.
This assistant is found right on the speed cell on the table of the pipe’s basic
data.
28 USERS MANUAL
- Material coefficient of pipe kmt. It’s the pipe’s material dependent coefficient
used to calculate speed through the expression(only valid for water)
9900
a (11)
Dt
47 '3 k mt
et
On the program, values for kmt are indicated for different pipe materials.
- Distance to origin. It’s the distance to which point i is located from the
pipe’s origin, measured along its axis. For the pipe’s initial point, the
program assigns automatically a zero distance, and for the end point
a distance equal to its length.
The distance and elevation columns, except for the extreme values that
are already assigned, can be copied and pasted from an Excel Sheet.
30 USERS MANUAL
Air valves are installed in order to avoid negative pressures inside pipes. They
admit air to its interior when the pressure at its location tends to be lower than
atmospheric pressure. When the pressure reaches values over the atmospheric
one, they release air. Nonetheless, since the entrance of air to a conduit may
cause, in some cases, more problems than solutions, the usage of air valves as
protection devices against hydraulic transients should only be done when no
other method to avoid negative pressures can be implemented.
ALLIEVI has an internal database with the characteristics of several air valves.
However, it also enables the user to define their own database with the desired
air valves and their characteristic curve.
Go to Pipes – Basic Data Profile Pipe settings Air valves. First, it has to be
selected the air valves file. It can be the air valves file of Allievi database, or a file
created by the user using Allievi.
The air valves file selected can be marked as default file (option “Mark file as
default”). This allows that while working with a project, it will be easier to access
to the desired air valves file. Therefore, it will not be necessary to select the folder
where the file is as the chosen database option will be marked as “Default” and
Allievi will present it automatically.
31 USERS MANUAL
Once the file has been selected, air valves can be chosen from the list of the air
valves database. For each air valve selected to be installed, the user has to specify
the distance in meters of the air valve to the pipe origin. These data will be stored
in a table. As air valves are introduced in the table, with their distance to the origin
pipe, they are represented in the profile of the pipe.
To delete an air valve, select the row where the air valve is and click the “Del” key
from the keyboard. If there is an air valve in the table without distance, Allievi will
not save it. It will be automatically deleted from the list of selected air valves.
In the simulation process, each air valve is moved to the nearest point of
calculation from its original situation. It is not possible to simulate more than one
air valve in the same point of calculation.
32 USERS MANUAL
IMPORTANT: It’s very important to bear in mind that, for calculation purposes,
if a pipe’s end is connected to an element, the installed vents must be located
at a minimum distance equal to:
∆𝑡 : Time interval
a: Speed
In the case of pipes whose end connects to other pipes, air valves can be
installed at the node associated to one of the pipes. For the others air valves,
a minimum distance to the node must be left:
For two consecutive air valves, there must be at least a minimum distance so
there can be a calculation point between them, this distance must be:
Allievi has his own air valves database file, but allows users to work with other an
air valve database file. The create it, Allievi allows the user to create it as an empty
file or as a copy from an existing file that can be modified afterwards. In this file,
the user can add new air valves with their admission and expulsion curves, or can
modify all the air valves that are in the file and their characteristics.
Go to the Settings BBDD: Manage the air valves files defined in ALLIEVI
33 USERS MANUAL
To define a new air valves database select “Create file” option. To define new air
valves in a file that already exists, to modify the air valves of a file, or to see the
air valves and its curves, select the “Setup” option.
If the option “Create empty file” is chosen, only the name of file and its path
is required.
34 USERS MANUAL
existing one, the “Duplicate” button has to be pressed. Air valves can be
removed by clicking on the “Delete” button.
It is also possible to import to the file an air valve and its characteristics if
it belongs to another air valves file introduced in Allievi, or from the Allievi’s
air valves file.
ALLIEVI enables you to simulate the effect of a pipe’s breakage on the system. To
do so, in the wizard where the pipe’s profile is introduced, select the “Burst” tab
and mark the option “Pipe with burst”. Then, the fields that define the burst are
enabled:
- Distance to origin. It refers to the distance from the pipe’s origin to the point
in which the burst will occur.
- Burst instant. It’s the instant of time in which the burst occurs.
- Kv final. It’s the burst’s final Kv coefficient that relates the discharged burst
flow Qrot with the inner pipe’s pressure at the point of breakage P tub,
according to the expression:
The final Kv of the burst could be assimilated to that of an open valve whose
passage section was that of the burst.
IMPORTANT: For the simulation, the burst is moved to the closest calculation
point to its location. However, the burst is not allowed to be moved to one
end of the pipe, where other element of the system will be connected.
Therefore, the minimum distance between the burst and any of the ends of
the pipe must meet
If a vent and a burst coincide in the same calculation point, the vent is
ignored and does not intervene in the simulation.
To calculate the steady state solution, the circulating pipe flow to initiate
calculations is obtained by assigning a supposed velocity in the order of 0.5 m/s
(or a different user-defined value). During the steady state calculation this velocity
will change with every iteration and for each pipe, until reaching a permanent
regime value.
39 USERS MANUAL
In case steady state flow in one or several pipes is zero (because of stopped pumps
and/or closed valves), the permanent regime conditions are obtained through a
great number iterations, which can consume important calculation time. In these
cases, for the steady state calculation you can directly assign a null flow to those
pipes which you know beforehand will have this regime flow value. For each of
the pipes, the necessary data are:
- Zero flow (Yes – No). It indicates whether we enforce a null initial flow.
- Number of parallel pumps. It refers to the number of parallel units, all equal,
present at the pumping station. These units will have the same exact
behavior during the transient and will have, in turn, the same regulation and
retention valves which will effectuate the same maneuvers.
- Pump characteristics curve shape. It indicates the way in which the pumps’
characteristics curve will be presented, Hp = Hp (Qp) and Pac = Pac (Qp), which
can be either Universal or Adjusted (By points).
- Regime rotation speed. It’s the rotation speed of the groups when they
operate at steady state regime conditions.
- Nominal rotation speed. It’s the pump’s rotating speed, corresponding to the
used characteristics curves.
- Inertia of each group. It’s the moment of inertia of the rotating masses of
each group, including the impeller, the axis and the driving engine’s rotor. If
42 USERS MANUAL
As we’ve previously mentioned, there are two ways of inputting the moment of
inertia. The direct way is to write it down in the corresponding box:
Rotation I (Kg.m2)
If you don’t know this datum, ALLIEVI gives you the option of calculating it, with
the optimal performance point of the pump, through an assistant, found in the
following way:
If pump has universal curve, the ‘Save’ option allows to copy the data used
to calculate the inertia to pump's grid. If pumps have characteristic curve
by points, the data from the optimal performance point to calculate inertia
are obtained from the characteristique curve previously introduced.
- Initial pump status. At initial conditions, the groups may be running (at
regime rotation speed) or stopped.
Pumping stations which are initially running may stop at a particular time and
then restart. In these cases, we have:
44 USERS MANUAL
- Stop (Yes - No). In pumping groups with initial running status, indicates if
they stop.
- Stop instant (Stop inst). Instant in which all pumping groups are
simultaneously stopped.
- New startup (Yes – No). It indicates whether the groups are going to be
restarted after their shut down.
- Startup instant. Instant in which the pumps begin the simultaneous startup
maneuver, if applicable.
- Duration. In case startup occurs, it’s the time taken until the nominal rotation
speed is achieved, considering a linear evolution.
Pumping stations which are initially stopped may startup at a given time and stop
yet again afterwards. In these cases we have:
- Startup (Yes – No). In pumping groups with initial stopped status, indicates if
they start run.
- Startup instant. Instant in which the pumps begin the simultaneous startup
maneuver.
- Duration. It’s the time taken until the nominal rotation speed is achieved,
considering a linear evolution.
45 USERS MANUAL
- New stop (Yes – No). It indicates whether the groups are going to be stopped
yet again after their startup.
- Stop instant. Instant in which the pumps begin the simultaneous stopping
maneuver, if applicable.
In this case, the pumps’ characteristics curves are obtained from the Marchal,
Flesch y Suter universal curves for a total of 14 different specific velocity values.
These curves enable the characterization of the pump’s behavior in whichever
functioning conditions, i.e. positive, negative or null flow, head or rotation speed
values.
- Power consumption Pcon0. Power consumption by the pump’s shaft for the
optimal performance point.
46 USERS MANUAL
Once these values are entered, the program shows for you to check the pump’s
efficiency at the optimal performance point.
The fitted characteristics curves are defined by points for nominal rotating speed,
obtained from the catalogue information provided by the pump’s manufacturer.
With this information the program is unable to characterize the pump’s behavior
for negative flow, head or rotating speed values. Hence, this manner of providing
the characteristics curve is only useful when the functioning of the valves is limited
to the first quadrant, with positive rotating speed.
In this case, the needed data to obtain the fitted characteristics curves are:
Once these three values are entered for each i operating point, the program
displays the pump’s efficiency for it to be verified.
47 USERS MANUAL
IMPORTANT: The pump’s operating points which define the fitted or point-
defined curves must cover the whole first quadrant, from Qp = 0 until Hp = 0,
with power consumptions that in no case should take negative or nil values.
Thus, the user must provide enough points so as to adequately define this
working interval, even if this means extrapolating the catalogue information
provided by the manufacturer. This is so because throughout the simulation
the pump may function in any point of its characteristics curve within the first
quadrant, and therefore the program must count on sufficient information
regarding the pump’s characteristics.
.
In order to input the characteristics curves’ defining data, ALLIEVI has an assistant
that can be reached in the following way:
Pumping groups with fitted curves must include, in all cases, a check valve. Groups
with universal curves may or may not incorporate this valve.
- Check valve’s minor head loss coefficient KVRt. It’s the check valve’s loss
coefficient when fully opened, with units in m/(m3/s)2, used to calculate
losses as a function of the pump’s boosted flow. This loss coefficient must
also include the losses concerning the pumping station’s remaining
accessories which have its same circulating flow, except those of the
regulation valve, if present.
The regulation valve is installed at the exit of each group, generally following the
check valve, if present. Its mission is to regulate the group’s boosted flow by
closing partially or, if fully closed, isolate the group from the discharge manifold.
For all these types, except the commercial valve, the program has a typical
characteristics curve which provides the dimensionless loss coefficient as a
function of the valve’s opening position.
For the case of the commercial valve, the characteristics valve curve will be
obtained from the valve database compiled using the catalogue information, and
can be uploaded through the program. In this case, the introduced data is:
- Opening position Ori. It’s the relative opening of the valve, in % with respect
to its maximum opening, at time instant ti.
This by-pass, if present, is installed between the inlet and outlet manifolds of the
pumping station, in parallel configuration with the other boosting groups. The
by-pass will include a check valve with the following data:
- Check valve’s minor head loss coefficient KVRbp. It refers to the minor head loss
coefficient of the check valve when fully opened, in m/(m3/s)2 used to
calculate the losses as a function of the circulating flow through the by-pass.
This coefficient must also include the losses concerning the rest of the
accessories that constitute the by-pass (Including the friction losses of the
corresponding pipe).
51 USERS MANUAL
2. Check valves
3. Relief valves
5. Injectors
7. Composite valves
8. Automatic valves
They are valves whose actioning allows the circulating flow through the
installation to be modified. In case they’re used as on/off valves, under ordinary
conditions they’re usually opened but might close to isolate a part of the system.
Their characteristics are:
- Elevation (Z). Elevation of the valve, equal for start and end nodes.
- Type of valve. It’s the type of valve to be installed, selecting from one of the
following types:
Butterfly
Spherical
Gate
Poppet
Diaphragm
Commercial
Poppet
Diaphragm
Commercial
For all these types, excluding commercial valves, the program possesses a typical
characteristics curve that provides the dimensionless headloss coefficient as a
function of the valve’s opening position.
For commercial valves, its characteristics curve will be obtained from the valve
database compiled using catalogue information. For these valves one must have
the dimensionless loss coefficient as a function of opening position, ranging from
0 % to 100 % with 5 % intervals of valve position.
When introducing a commercial valve, the user must decide which type of valve
he or she wants to introduce. The way to do it is as follows:
Regulation valves can effectuate opening and/or closing maneuvers during the
transient simulation. These maneuvers may be defined as:
Tabulated
Sinusoidal
In any of these maneuvers, the relative opening of the valve is limited to 0 - 100%
interval.
- Opening position Ori. It’s the relative opening of the valve, in % with respect
to its maximum opening, at time instant ti.
The user must introduce the valve’s maneuver in the following manner:
2
Apt Ap0 A sen ( t B) (151)
T
Where the relative opening in % with respect to the maximum one, depends on
the following parameters:
- Initial opening Ap0. It’s the initial valve opening in % with respect to the
maximum one.
- Phase B. It’s the phase angle of the valve’s opening position’s oscillation.
56 USERS MANUAL
These valves open when flow circulates with forward direction, and close when
backflow is established (inverse flow direction). They are considered as ideal
elements, instantaneously closing exactly when the valve’s outflow becomes
negative. Their characteristics include:
- Loss coefficient k for fully closed valve. It refers to the dimensionless headloss
coefficient for a fully closed valve. If, when the valve is closed, the backflow
is entirely prevented, this coefficient will take a value of 1·10252.
57 USERS MANUAL
- Loss coefficient k for open valve. It’s the dimensionless loss coefficient for a
fully opened valve.
They’re valves which open when pressure at their connection point exceeds a
given value, discharging to the exterior all or part of the circulating flow as long
as this pressure is kept over the aforesaid value. Once opened, their opening
position depends on inlet pressure until fully opened. On the contrary, when
pressures are low this valve will remain closed.
For each of the installed relief valves, the necessary data are:
- Type of valve. It’s the type of relief valve to be installed at the branch,
selecting from one of the following types:
Butterfly
Spherical
Gate
Poppet
Diaphragm
Commercial
Once we’ve introduced the number of branches, we then display the data table,
and introduce the following information for each of the branches:
For all these predefined types, excluding commercial valves, the program
possesses a typical characteristics curve that provides the dimensionless headloss
coefficient as a function of the valve’s opening position.
For commercial valves, its characteristics curve will be obtained from the valve
database compiled using catalogue information. For these valves one must have
the dimensionless loss coefficient as a function of opening position, ranging from
0 % to 100 % with 5 % intervals of valve position.
59 USERS MANUAL
- Valve maneuver: It relates valve opening with inlet node pressure. To define
this maneuver in a tabular manner the following data must be provided:
To eliminate the relief valve of one of the branches, select this by clicking on the
corresponding box on the left. Then the valve is eliminated by pressing the Delete
key.
- Type of valve. It refers to the type valve to be installed, selecting from one
of the following types:
Butterfly
Spherical
Gate
Poppet
Diaphragm
Commercial
For all these types, excluding commercial valves, the program possesses a typical
characteristics curve that provides the dimensionless headloss coefficient as a
function of the valve’s opening position.
For commercial valves, its characteristics curve will be obtained from the valve
database compiled using catalogue information. For these valves one must have
the dimensionless loss coefficient as a function of opening position, ranging from
0 % to 100 % with 5 % intervals of valve position.
- Instant ti. It’s the time instant corresponding to point i. This time is taken
from the instant in which the valve is triggered or when the closing maneuver
begins.
- Opening position Ori. It’s the relative opening of the valve, in % with respect
to its maximum opening, at time instant ti. For the initial instant, the valve’s
opening position will be 100 %.
When dealing with an anti-burst valve, its trigger may be conditioned by:
For valves which are triggered by minimum or maximum sped, the velocity which
triggers the opening maneuvers will be:
For those which are triggered by flow difference, initial and final flowmeters
whose difference will provoke the trigger have to be placed. To do so, the
necessary data are as follows:
- Flow difference. Flow difference reference value for the excess flow valve to
be triggered.
- Pipe where flowmeter is located. Designation of the pipes where we find both
the initial and final flowmeters. In general, both flowmeters are placed on
the same pipe, although we can locate them in different ones also.
- Q coefficient of the flow meter (+1, -1). This parameter indicates the flow’s
direction used to obtain the flow difference which in turn produces the
trigger. If this coefficient takes up the value +1, the flow is considered to be
positive in the direction of the pipe, and if on the contrary its value is -1 the
64 USERS MANUAL
NOTE: In case no pipes are available, the program shows the following
warning.
65 USERS MANUAL
For each calculation time instant, the program will calculate the difference in flows
obtained at the initial and final flowmeters’ locations, these being affected by the
corresponding flow coefficient. The absolute value of this difference is compared
to a reference one in order to decide whether the excess flow valve should be
triggered.
3.6.5 Injectors
They are valves, with discharge to the atmosphere, whose outlet diameter is less
than its nominal diameter. They admit a high output speed, which is why the
kinetic term corresponding to this speed is not negligible. Their characteristics are:
- Nozzle Diameter (D Boq). It is the outlet diameter of the injector. This diameter
must be smaller than the nominal diameter.
- Head loss coefficient (k) of the branch. It’s the dimensionless head loss
coefficient of the branch where the injector is connected. This head loss
coefficient is related to the speed referred to the ND of the injector.
The data on Type of Injector Valve and Maneuver Type are defined in the same
way as the corresponding data of the Regulation Valves, section 3.6.1.
- Head loss coefficient k.. It’s the dimensionless head loss coefficient
corresponding to the accessory.
- Elevation. Elevation of the composite valve, equal for start and end nodes.
- Number of check valves. It refers to the number of parallel check valves which
comprise the composite valve.
The data of each of the components which structure the composite valve are
analogue to those presented at subsection 3.6.1 (regulation valves), 3.6.2 (check
valves) and 3.6.6 (minor headlosses).
Display the Basic data table of the Composite valve( ) one must input
the data for each of the previously declared valves
To eliminate any of the components of the composite valve, select this one, by
clicking on the corresponding box on the left. The component is then deleted by
pressing the Delete key. It is the same procedure that is used to suppress branches
in the relief valve.
- Elevation. Elevation of automatic valve, equal for start and end nodes.
- Valve function. It’s the type of valve to be installed, according to the following
classification:
Pressure reducing
Pressure sustaining
Flow control
- Valve type. It’s the type of valve to be installed, selecting from one of the
following types:
70 USERS MANUAL
Butterfly
Spherical
Gate
Poppet
Diaphragm
Commercial
For each of these types, apart from commercial valves, the program
possesses a typical characteristics curve that provides the dimensionless
headloss coefficient as a function of the valve’s opening position.
For commercial valves, its characteristics curve will be obtained from the
valve database compiled using catalogue information. For these valves one
must have the dimensionless loss coefficient as a function of opening
position, ranging from 0 % to 100 % with 5 % intervals of valve position.
du t
K p et (16)
dt
Where,
71 USERS MANUAL
e = error between present parameter and setting value, divided by the setting
value for it to be dimensionless
u = Valve’s response to correct the error, in percentage of valve opening
In expression (16) we use a positive sign for pressure sustaining valves, whereas,
for pressure reducing and flow control valves, we use a negative one.
This option is only activated when pressure sustaining or pressure reducing valves
are declared.
- Flow setting. It’s the flow setting for flow control valves.
Again, this option is only enabled when flow control valves are declared.
Protection structures are used to reduce pressure oscillations once the transient
has been generated. They generally store a certain volume of water, and function
by supplying flow to the main conduit, to which they are connected, when the
circulating flow is reduced. Afterwards, when the flow direction changes, the latter
is gradually introduced into the protection structure, and in turn its breakage is
slowed down. In this manner, alterations in flow at the main pipe are slowed down,
reducing as a result the pressure oscillations.
Elements to be placed at the branch junction between the main pipe and the
protection structure should present in their whole a low headloss coefficient for
the structure’s outflow and a high one for its outflow. In such manner, the outflow
from the structure towards the main pipe is facilitated, which in turn reduces
depressions in said pipe, thus averting negative pressures. Moreover, when the
flow circulates from the main conduit to the protection structure important
headlosses are generated at the branch pipe which dissipate water’s energy,
slowing down the transient’s development.
2. Surge tank
3. Unidirectional tank
Regarding the necessary data, we’ll make a distinction between data of the
structure itself and that of the union branch. Since these last data are common to
all structures, specially the head loss coefficient for the inflow and outflow, their
regard will be done once the particular information is described for each structure.
- Elevation of the derivation node. It is the elevation of the node to which the
branch or union branch of the air vessel is connected.
- Number of air vessels. It’s the number of identical air vessels, installed in
parallel connection, which constitute the protection structure.
- Length. It refers to the air vessel’s length, supposedly with a cylindrical shape.
- Volume. It refers to the total volume of each air vessel, calculated by the
program once the vessel’s diameter and length are introduced.
- Position. The Air Vessel can be vertical (small Vessels) or horizontal (middle
and large Vessels).
- Initial inflation pressure. In the case of an air vessel with bladder, it’s the
pressure to which the bladder is inflated before the air vessel is installed. For
the bladderless air vessel, the inflation will be carried out closing a valve at
the vessel’s connection pipe (afterwards, this valve will open once the
installation is running and pumps are working at normal regime). To
calculate the inflation pressure, admitting that the gas suffers an isotherm
evolution between the inflation and the operation of the installation on
permanent regime condition, the expression used is the following one:
*
P0C = Gas pressure in permanent regime conditions, in absolute
value. In air vessels connected to the discharge manifold of a
pumping station, this pressure value is approximated as the
head provided by the pumps at permanent regime
conditions, deducting the difference in elevation between the
pumps and the average water level inside the air vessels, and
adding atmospheric pressure.
On ALLIEVI, the user can input the value for inflation pressure, in bar, or make use
of the assistant provided by ALLIEVI, which automatically calculates this parameter
once the values for initial volume of gas in cubic meters and initial gas pressure
in meters are given.
As indicated above the air vessel can include a bladder or not. Initially, the
behaviors of both types of air vessels are alike as long as they aren’t emptied
during the transient regime functioning. In case it is emptied, if the air vessel
has a bladder, when it’s emptied said vessel closes the connection orifice to
the branch union pipe. It’s at this time when the discharge of water from the
vessel to the main pipe is abruptly interrupted, thereby creating a depression
wave which can submit such pipe to negative pressures. Later, when in the
overpressure phase of the transient, water can flow back again to the air
vessel, repeating the cycle depending on the transient’s evolution.
Nonetheless, if the air vessel is bladder- free, in case it’s emptied, the gas (air
in this case) continues its expansion occupying the whole of the vessel and
seeping into the connection pipe, eventually reaching the main pipe. In this
case, the discharge to such pipe is not interrupted, thus preventing the
depression wave of the previous scenario, although part of the gas can travel
within the inside of the main pipe and not return back to the air vessel in the
subsequent overpressure phase. Besides, a certain volume of air travelling
uncontrolled within the main pipe’s interior may create a risk of breakage
during the following startup of the installation. Through the simulation, it’s
admitted that in case of emptying all the gas is confined within the air vessel
77 USERS MANUAL
and the branch connection pipe, without transferring to the main pipe. As a
consequence, the whole gas mass can go back to the air vessel’s interior
during the following overpressure phase of the transient.
The data of the union branch between the impulsion pipe and each of the Air
Vessels that constitute the protection structure are:
- Length. It’s the union branch length. This length can’t be less than the
difference between the Air Vessel base elevation and the derivation node
elevation.
- Absolute Roughness. It’s the absolute roughness of the inner walls of the
union branch.
- Head Loss Coefficient for the outlet flow of the Air Vessel. It is the loss
coefficient Ks of the union branch for the flow from the Vessel to the
impulse pipe, which provides these losses through the expression.
hs K s Qs2 (19)
- Coefficient of singular losses for the inlet flow of the Air Vessel. It is the loss
coefficient Ks of the union branch for the flow from the impulsion pipe to
the Vessel, which provides these losses through the expression.
he K e Qe2 (20)
If several air vessels have been defined in parallel, each one of them will have its
own union branch line with the characteristics that we have just indicated. In this
case, the singular losses are evaluated by expressions (19) and (20) with the output
or input flow rate to each air vessel.
air vessels
The head loss coefficients of the union branch in air vessels depend on the
characteristics of the said branch and the accessories installed in it. The value of
these coefficients can be introduced directly as indicated in the previous section,
or estimated using the assistant provided by the program.
The data necessary to estimate the head loss coefficients of the union branch are:
- Number of Air Vessels. It is the number of equal air vessels, in parallel, that
constitute the protection structure.
- Diameter of the impulsion pipe. It is the diameter of the pipe to which the
air vessels are connected by their union branches.
- Diameter of the union branch. The assistant estimates this diameter as half
the diameter of the discharge pipe divided by the square root of the
79 USERS MANUAL
number of air vessels. Its value can be modified by the user, according to
his own criteria, or to adapt it to the diameters of the commercial pipes
available.
- Diameter of the air vessels connection nozzle. It is accepted that the air
vessels has a shut-off valve with the same diameter as the connection pipe.
- Bypass diameter. This data is entered when union branch with bypass is
specified. Its value cannot exceed the diameter of the union branch.
This data is entered into ALLIEVI using one of the following screens:
80 USERS MANUAL
In the union branch, with or without by pass, it is possible to have or not an energy
dissipater for the inlet flow to the air vessel. This energy dissipater, if it exists, will
be characterized by a nozzle and will be located in the union branch when there
is no bypass available, or by a hole-plate or a nozzle that will be placed in the
bypass if it exists. The data to be introduced in this case are:
- Diameter of the nozzle, if it exists. This diameter must be smaller than the
diameter of the conduit where it is installed.
- Diameter of the plate orifice, if it exists. This diameter must be smaller than
the diameter of the conduit where the plate it is installed.
Once the indicated data have been entered, the following screen shows the
estimated head loss coefficients for the output and inlet flows to the air vessel.
With the option to Finish these coefficients can be copied to the data table of the
air vessel, along with the number of boilers and the diameter of the union branch
if they have changed with respect to the values previously introduced.
83 USERS MANUAL
From the last screen, and using the Advanced Assistant option, you can access the
data set introduced in the previous screens, as well as the head loss coefficients
calculated by the program for each of the elements that are part of the union
branch. These coefficients have been calculated from the existing information in
different publications of Fluid Mechanics.
The calculated coefficients are grouped into Output coefficients, Input coefficients
and Other coefficients. In this screen, you can modify the input data, which
modifies the value of the calculated coefficients. You can also directly modify the
head loss coefficients of the accessories, if more appropriate values obtained from
other sources or by catalog are available.
84 USERS MANUAL
With the option Accept the head loss coefficients can be copied to the data table
of the air vessel, together with the number of air vessels and the diameter of the
union branch if they have changed with respect to the values previously
introduced. On the other hand, with the Restore option you can restore the data,
or the loss coefficients of the accessories, to their original values without
modification.
ALLIEVI allows to carry out a pre-dimensioning of air vessels when the installation
consists of, or can be assimilated to, a pumping between tanks by means of simple
pipes. For this the Parmakian method is applied.
To carry out the pre-dimensioning of air vessels, previously the head loss
coefficients of the union branch have been estimated with the required
provisional data. This assistant is accessed through the pre-dimensioning option.
- Friction factor. It is the estimated friction factor for the impulse pipe. In
general, the friction factor takes values between 0.012 for large diameters,
and 0.020 for small diameters.
- Speed. It is the estimated value for the speed of the impulse pipe.
- Limit the maximum piezometric height at the exit of the pumping station.
- Limit the minimum piezometric height at the exit of the pumping station.
- Limit the minimum piezometric height at the midpoint of the impulsion pipe.
87 USERS MANUAL
If one of the cases of limiting the maximum piezometric height is selected, the
following data will be introduced:
If one of the cases of limiting the minimum piezometric height is selected, the
following data will be introduced:
Once the indicated data have been entered, the following screen shows the results
of the pre-dimensioning of the air vessel, such as length, diameter, and initial
volume of gas inside each air vessel. If necessary, the dimensions of the air vessels
can be modified by the user, to adapt them to the closest commercial values.
Similarly, the inflation pressure can be modified to adopt more convenient values.
With the option Finish the values of pre-dimensioning can be copied to the table
88 USERS MANUAL
of data of the air vessel, if they have changed with respect to other values
previously introduced.
- Elevation of the union branch node. It is the elevation of the node to which
the derivation or union branch of the surge tank is connected..
Constant Section
Variable Section
Differential
89 USERS MANUAL
In these surge tanks, it is admitted that the section varies stepwise, with segments
of constant section between the lower and upper levels of each step. A number
of sections Ntr of constant section, in which the surge tank is divided, is defined.
This results in Ntr + 1 points so that point 0 corresponds to the base of the tank,
point Ntr corresponds to the top of the tank, and the remaining Ntr- 1 points
correspond to the different section changes, ordered from bottom to top.
For each i point form i=0 until i= Ntr the data needed is:
- Accumulated Volume. It is the stored tank volume from the base to point i.
For i = 0, the program directly assigns a stored volume of zero and then
the user enters the volume for each of the elevations of the tank.
The volume of the surge tank of variable section is introduced in the following
way:
For differential surge tanks the necessary data are the following:
Main standpipe
- Diameter. It’s the inner diameter of the circular cross-section surge tank.
- Base elevation. It refers to the elevation of the surge tank’s base or bottom.
Inner tube
- Height. It’s the height of the surge tank’s inner tube, measured from the
bottom.
- Flow coefficient of the inner tube’s overflow section. It’s the CQ coefficient of
the weir formed by tube’s superior contour, used to calculate the discharged
flow according to expression at (9). To implement this expression, the
rectangular weir’s length L is substituted by the circumference’s length of the
tube’s top contour. Typical values for CQ thin wall weir range from 0.60 to
0.80, with a maximum value of 0.85. For thick-walled weirs, this coefficient is
0.577.
- Coefficient of singular head loss for outflow of the inner tube. It is the head
loss coefficient Ks for the flow from the inner tube to the impulsion pipe (or
forced pipe in case of hydraulic power plants), which provides these losses
by means of the expression (19) applied to the exit flow of the pipe.
- Coefficient of singular head loss for inflow of the inner tube. It is the head loss
coefficient Ks for the flow from the impulsion pipe to the inner tube, which
provides these losses by means of the expression (20) applied to the inflow
of the pipe.
K in (m/(m3/s)2)
The data of the union branch between the impulsion pipe and the surge tank,
whether it is of constant section, variable section or main standpipe in the
differential type are:
- Length. It’s the union branch length. This length can’t be less than the
difference between the surge tank base elevation and the derivation node
elevation.
- Absolute Roughness. It’s the absolute roughness of the inner walls of the
union branch.
- Head Loss Coefficient for the outlet flow of the surge tank. It is the loss
coefficient Ks of the union branch for the flow from the surge tank to the
impulse pipe, which provides these losses through the expression (19)
applied to the outflow of the surge tank.
- Coefficient of singular losses for the inlet flow of the Air Vessel. It is the loss
coefficient Ks of the union branch for the flow from the impulsion pipe to
the Vessel, which provides these losses through the expression (20) applied
to the inflow of the surge tank.
For the differential surge tank, it is possible to define one or several connecting
branches between the impulsion pipe and the main surge tank, so that the data
of the union branch indicated are for each one of those branches, all equal. For
the inner tube of this surge tank, the head loss coefficient for the output and inlet
flows is defined as indicated in the previous section, without specifying the data
of its own union branch.
The head loss coefficients of the union branch in surge tanks depend on the
characteristics of the said branch and the accessories installed in it. The value of
these coefficients can be entered directly as indicated in the previous section, or
estimated using the assistant provided by the program.
For the surge tanks of constant or variable section, the necessary data to estimate
the head loss coefficients of the union branch are:
- Diameter of the surge tank. It is the inside diameter of the surge tank, which
cannot be smaller than the diameter of the connection hole. This diameter
will be indicated in m, the previous diameters being in mm.
In the variable section surge tanks, as diameter for the calculation of the head loss
coefficients of the union branch, the one that corresponds to the lower section is
taken, assuming this as a circular section.
The next step indicates whether the connection between the main pipe and the
surge tank is direct or in derivation. In case of connection in derivation this one is
done by means of a union branch of a certain length, introducing the number of
elbows of 90º to dispose in the mentioned union branch.
Once the indicated data have been entered, the following screens show the
estimated head loss coefficients for the surge tanks outflows and inflows, in the
cases of direct connection or by means of a union branch. With the option Finish
these coefficients can be copied to the surge tank data table, together with its
diameter and the diameter of the union branch if they have changed with respect
to the values previously introduced.
97 USERS MANUAL
From the last screen, and using the Advanced Assistant option, you can access the
data set introduced in the previous screens, as well as the head loss coefficients
calculated by the program for each of the elements that are part of the union
branch. These coefficients have been calculated from the existing information in
different publications of Fluid Mechanics.
The calculated coefficients are grouped into Output coefficients, Input coefficients
and Other coefficients. In this screen, you can modify the input data, which
modifies the value of the calculated coefficients. You can also directly modify the
head loss coefficients of the accessories, if more appropriate values obtained from
other sources or by catalog are available.
98 USERS MANUAL
99 USERS MANUAL
With the option Accept the head loss coefficients can be copied to the data table
of the surge tank, together with the diameter of the surge tank and the diameter
of the union branch if they have changed with respect to the values previously
introduced. On the other hand, with the Restore option you can restore the data,
or the head loss coefficients of the accessories, to their original values without
modification.
The Allievi software has an assistant for the estimation of the head loss coefficients
in the union branches of the differential surge tanks.
The differential surge tanks can be connected to the impulsion pipe (or forced
pipe in hydraulic power plants), either directly or by derivation. The data necessary
to estimate the loss coefficients in both cases are:
- Number of 90º elbows in each union branch. It is the number of 90º elbows
of each connecting branch when the surge tank is connected in derivation.
- Diameter of the union branch. This diameter can be the same or different as
the diameter of the orifice in the main standpipe.
- Diameter of the inner tuve orifice. It is the diameter of the connection orifice
of the union branch with the inner tube.
Once the indicated data have been entered, the following screens show the
estimated head loss coefficients for the main standpipe and inner tube outflows
and inflows, in the cases of direct connection or by means of a union branch. With
the option Finish these coefficients can be copied to the surge tank data table,
together with its diameter, the diameter of the union branch and inner tube if
they have changed with respect to the values previously introduced.
103 USERS MANUAL
From the last screen, and using the Advanced Assistant option, you can access the
data set introduced in the previous screens, as well as the head loss coefficients
calculated by the program for each of the elements that are part of the union
branch. These coefficients have been calculated from the existing information in
different publications of Fluid Mechanics.
104 USERS MANUAL
105 USERS MANUAL
For the main standpipe and the inner tube, the calculated coefficients are grouped
into Output coefficients, Input coefficients and Other coefficients. In this screen,
you can modify the input data, which modifies the value of the calculated
coefficients. You can also directly modify the head loss coefficients of the
accessories, if more appropriate values obtained from other sources or by catalog
are available.
With the option Accept the head loss coefficients can be copied to the data table
of the surge tank, together with the diameter of the surge tank and the diameter
of the union branch if they have changed with respect to the values previously
introduced. On the other hand, with the Restore option you can restore the data,
or the head loss coefficients of the accessories, to their original values without
modification.
- Elevation of the derivation node. It is the elevation of the node that connects
the derivation or union branch of the unidirectional tank.
- Absolute roughness. It refers to the absolute roughness of the inner wall of the
tank.
- Maximum water level in unidirectional tank. It’s the maximum water level
allowed inside the tank, referenced to its bottom.
The data of the union branch between the impulsion pipe and the unidirectional
tank is:
- Length. It’s the union branch length. This length can’t be less than the
difference between the surge tank base elevation and the derivation node
elevation.
- Absolute Roughness. It’s the absolute roughness of the inner walls of the
union branch.
- Head Loss Coefficient for the outlet flow of the surge tank. It is the loss
coefficient Ks of the union branch for the flow from the surge tank to the
impulse pipe, which provides these losses through the expression (19)
applied to the outflow of the surge tank.
In this case, the singular head loss coefficient is not required for inflow to
the unidirectional tank, as the inlet flow to the tank is prevented by closure
of the check valve.
unidirectional Tanks
The head loss coefficients of the union branch for the outflow in unidirectional
tanks depend on the characteristics of the said branch and the accessories
installed in it. The value of these coefficients can be entered directly as indicated
in the previous section, or estimated using the assistant provided by the program.
In this case, the data necessary to estimate the head loss coefficient of the union
branch are:
- Diameter of the check valve. It is the nominal diameter of the check valve to
be installed in the union branch, which allows the unidirectional tank outlet
flow and prevents the inflow. This diameter cannot be greater than the
diameter of the union branch.
Once the indicated data have been entered, the following screens show the
estimated head loss coefficient for the outflow of the unidirectional tank. With the
option Finish these coefficients can be copied to the surge tank data table,
together with its diameter, the diameter of the union branch and inner tube if
they have changed with respect to the values previously introduced.
From the last screen, and using the Advanced Assistant option, you can access the
data set introduced in the previous screens, as well as the head loss coefficients
calculated by the program for each of the elements that are part of the union
branch. These coefficients have been calculated from the existing information in
different publications of Fluid Mechanics.
109 USERS MANUAL
The calculated coefficients are grouped into Output coefficients, Input coefficients
and Other coefficients. In this screen, you can modify the input data, which
modifies the value of the calculated coefficients. You can also directly modify the
head loss coefficients of the accessories, if more appropriate values obtained from
other sources or by catalog are available.
With the option Accept the head loss coefficients can be copied to the data table
of the unidirectional tank, together with the diameter of the unidirectional tank
and the diameter of the union branch if they have changed with respect to the
values previously introduced. On the other hand, with the Restore option you can
restore the data, or the head loss coefficients of the accessories, to their original
values without modification.
110 USERS MANUAL
Flow rules
Head rules
With this law a specific outflow is imposed to the corresponding node. The
following possibilities exist:
By Q law
By node pressure
By table
- Elevation. It’s the elevation to which the law is located, equal to the
connection node’s elevation.
The flow’s imposition by law forces the node’s outflow to follow the following lay
2
Qt A t 2 B t C D sen t E (21)
T
With t in seconds, E in degrees and Q in l/s.
When forcing the outflow, the program calculates in the node the hydraulic head
required in order to achieve this particular outflow. In this case, the required data
are:
- Constant term C. It refers to the constant term at the polynomial part of the
Q law.
- Coefficient D. It’s the coefficient for the sinusoidal term of the Q law.
- Phase angle E. Phase angle for the sinusoidal term of the Q law.
113 USERS MANUAL
- Limitation. It allows for the extreme flow values given by the Q law to be
limited. The different possibilities include:
No limitations
Lower limit
Upper limit
Upper and lower limit
In this case it’s assumed that there’s an element at the node which generates a
headloss so as to discharge a flow to the atmosphere following the law below
PN t KQ N t QN t (22)
where PN(t) is the node’s pressure, K the loss coefficient of the flow-controlling
element, and QN(t) the node’s discharged flow. The K coefficient will be constant
throughout the whole simulation, and it’s evaluated using the initial imposed flow
and the node’s outlet pressure in steady state conditions. The required datum is:
- Value of Q. It’s the discharged flow through the node at initial conditions.
A tabulated flow law is imposed. For each point i the data to enter are:
- Outflow Qi.. It’s the node’s discharged flow at time instant ti. If it’s an inflow,
its value will be negative.
Go to Rule – Basic data Q imp: through this button we can input the law’s
maneuver
116 USERS MANUAL
By means of this law, a certain hydraulic head is imposed only to a node where
several pipes converge. When imposing this piezometric head, the program
calculates the node’s discharged flow under these conditions.
To match up each head law to a particular tank, which the water is poured we
must:
NOTE The way to assign a law without tank discharge is by selecting the dotted
line.
.
117 USERS MANUAL
In case there’s no associated tank, the function which defines the node’s hydraulic
head imposition is
2
H t A t 2 B t C D sen t E (23)
T
With t in seconds, E in degrees and H in m.
- Constant term C. It refers to the constant term at the polynomial part of the
Q law.
- Coefficient D. It’s the coefficient for the sinusoidal term of the Q law.
- Phase angle E. Phase angle for the sinusoidal term of the Q law.
- Limitation. It allows for the extreme hydraulic head values provided by the
Q law to be limited. The different possibilities include:
No limitations
Lower limit
Upper limit
Upper and lower limit
These data are introduced in a likewise manner as shown for flow laws.
IMPORTANT: Once the head law is set, the program calculates the node’s
corresponding flow, which can be positive (exiting flow from the system) or
negative (inflow to the system). This should correspond to the action of a
subsystem, which isn’t accounted for on the simulation but imposes the
hydraulic head law represented at expression (23). This part of the system is
the one which would receive or provide the node’s outflow or inflow.
In case the head law is linked to a single pipe’s end node, such law can utilize an
associated tank to which the pipe’s flow is discharged, at atmospheric pressure.
The program then calculates the outflow keeping the node’s hydraulic head equal
to its elevation, and considering said flow as an inflow to the tank through its top
section.
If, during the transient simulation, the head law’s outflow becomes negative it
means that there’s a backflow through the law-connected pipe. From this point
on, the backflow generates an air column at the pipe’s end, which will increase or
decrease its length depending on the movement of water inside the conduit. At
all pipe points where there’s an air column null flow and atmospheric pressure will
be supposed, with air being able to enter or exit according to the head law.
Immediately upon the end of the air column’s expulsion, the outflow will be
reestablished following the head law. If the air column occupies at any time the
119 USERS MANUAL
whole of the pipe, the program stops the simulation process, since air movement
at elements which connect the pipes to their opposite law’s nodes can no longer
be evaluated
- Tank. It’s the law –linked tank. All system’s tank can be considered as
associated tanks. The program provides, for information purposes solely, the
elevation from which water is poured into the tank, which coincides with
law’s connection node.
- Area. If the tank associated to the law is a weir, one must indicate if the
discharge area is the initial or the final one.
A turbine must function keeping a constant rotating speed for any load case. In
this fashion, the turbine’s regulator will action the guide’s rotor blades, increasing
the turbined flow as the resistant turbine torque increases (load increments), or
decreasing the turbined flow when this torque is diminished (load declines). With
this, we’re looking to equalize the driving and the resistant torque for any load
combination, while keeping the rotation speed constant. Hence the frequency of
the generated alternate current will also be unaltered.
The expression that defines the functioning of the PID reagulator feedback is as
follows
det
ut K p et K i et dt K d K ra u t dt (24)
dt
d 2 et
K ra ut K i et K d
du(t ) de(t )
K p (25)
dt dt dt 2
With this formulation, any partial action of the regulator can be canceled making
its constant null.
For each turbine, the general data requested by the program are:
- Elevation Z. It’s the elevation to which the turbine rests, equal for start and
end nodes.
- Inertia. It’s the moment of inertia of the rotating masses in the turbine,
including the impeller, the shaft and the rotor of the electric generator.
- N regime or reference rotation velocity. It’s the rotation speed to which the
desired alternate current’s frequency value is produced. This rotation velocity
will have to be kept constant by means of the regulator’s actions.
The data necessary to characterize the turbine correspond to the optimal point of
operation (or point of maximum performance) at the nominal rotation speed, that
is:
- Height Ht0. It is the net height of the turbine at the optimal performance
point.
- Shaft Power Pejet0. It is the power in the shaft of the turbine at the optimal
performance point. .
Once these three values have been entered, the software displays the turbine
efficiency at the optimal performance point for verification.
- N nominal. It is the nominal rotation speed of the turbine, for which the
optimal performance point is obtained.
It represents the time evolution of the turbine’s imposed load, during which the
regulator will try to keep the reference rotation speed constant, regardless of this
evolution.
For each of the turbines, the time instants and their respective loads percentages
are entered in tabular form so as to define the load variations of the generator.
- Load Li with respect to normal regime. It’s the supplied power needed at time
instant i, or the resistant torque to overcome, in percentage with respect to
the normal regime’s power or resistive torque.
If, during the regulation maneuvers, the load on the turbine shaft decreases below
a certain value, the turbine generator will be disconnected from the network to
avoid operation at excessively low partial load. At that time, the distributor will
initiate a quick closing maneuver to avoid turbine runaway.
- Minimum disconnection load. It is the minimum load that will cause the
disconnection of the turbine, in percentage with respect to the regime
power.
The time instants and their respective guide openings for the closing process are
introduced in table form. The time reference is the turbine’s disconnection instant,
starting from the initial guide opening. These values are:
Instant ti. Time instant corresponding to point i for the guide’s closing
maneuver.
Guide opening Api. Guide opening at time instant i, during the process of
closure in order to avoid turbine runaway.
The closure of the guide vanes will follow this law, beginning the maneuver with
the guide’s opening at the instant of disconnection.
Time instants and guide opening during closure procedure are assigned in the
following manner:
For the feedback PDI controller in each turbine, the required data are:
-
Ki = Integral term constant of the regulator, with units in s-1
- Kra = Feedback constant of the regulator, with units in s-1. This constant can
take negative values.
NOTA: To set the value of the regulator constants it is advisable to use the
Ziegler and Nichols method. This method provides the values of the constants
of a PI or PID regulator type, with null feedback constant. In some cases, the
values provided by this method are adjusted manually to improve the
characteristics of the regulation.
Open channel flow conduits, or channels, can only be connected to tank nodes or
nodes which connect other channels. In case there’s a connection to a tank, this
specific node will solely link the tank and the channel. In case several channels are
connected together, the connection node will only join channels.
In order to define open- channel flow conduits, the necessary data are:
- Length. Length of the channel, supposed straight and with a uniform slope.
- Bottom elevation. Bottom elevations for start and extreme ends of the
channel. In case there’s a tank connection, the channel’s bottom elevation
will be higher or equal to the connection node’s elevation. In case its higher,
the channel’s bottom’s elevation will be higher than the tank’s bottom’s
elevation. In the case of channel connection, at this end all channels must
have the same bottom elevation, which will match the node’s elevation.
128 USERS MANUAL
The length of the channel, together with the bottom elevation of the initial
and final ends, will define the channel's slope according to the expression
Zsi Zsf
s0 % 100 (26)
L
Rectangular
Trapezoidal
Circular
Any
- Channel width. It’s the width of the channel, constant for all its cross-section.
- Maximum depth. It’s the maximum depth of the channel, from which
overflowing occurs and simulation is stopped.
- Channel width. It’s the width of the channel, constant for all its cross-section.
129 USERS MANUAL
- Maximum depth. It’s the maximum depth of the channel, from which
overflowing occurs and simulation is stopped.
- Lateral side slope, in degrees. It’s the incline angle of the lateral walls of the
channel, with respect to the horizontal plane, measured from the exterior
part of the channel.
- Channel diameter. It’s the inner diameter of the circular channel. The
maximum depth is indicated, which will coincide with this value, from which
said conduct becomes pressurized and simulation is stopped.
The data for the channel cross-section will be given in table form, with increasing
channel depth, these being for point i:
- Perimeter pi. Wetted perimeter for conduit cross-section, under depth yi.
- Section Ai. It refers to the area of the conduit’s cross section, under channel
depth yi.
As in pressurized conduits, it can occur that the steady state circulating flow
through the channel is zero. Thus, and also to avoid the need for a large number
of iterations in order to calculate the steady state, we can directly assign a null
initial flow to those channels to which this particular circumstance applies. So, for
each channel with null initial flow, the necessary information is:
131 USERS MANUAL
- Zero initial flow (Yes – No). It indicates whether we enforce a null initial flow
to the channel
IMPORTANT: We have to be very careful when inputting this data, since when
we indicate zero initial flow for a channel the program calculates the steady
state assigning this value to such channel, as well as the depth corresponding
to the water level one has previously introduced. The obtained steady state
solution will be incorrect if in fact these circumstances will not take place, as in
the case, for instance, of incompatible initial water levels between the
connected tanks and the channels with zero initial flow. In case we’re hesitant,
it’s better not to impose this type of data.
When declaring a new project, ALLIEVI will allow the user to define and save the
considered useful information, as in the case of the title, responsible Engineer,
date of creation, page size, amongst others. The objective is to enhance the
organization of the created projects.
132 USERS MANUAL
ALLIEVI allows the user to establish several program options. These options can
be changed for the present project as well as for all projects made with ALLIEVI.
Go to Project tab Project options ALLIEVI options Units: The user can
select the unit system to work with.
ALLIEVI counts with multiple tools which allow the user to adjust the work area to
their needs and/or preferences. We introduce below the way to access every one
of them.
In the first place, the default tool on ALLIEVI which allows us to introduce and edit
the elements, activate or deactivate tools or even introduce network data, is the
Selection tool.
- On the one hand, on the “View” tab, “Zoom”, we find the “Width” and
“Height” buttons.
- The screen adjustment buttons “Width” and “Height” can also be found on
the taskbar, located at the bottom of the work area.
- Magnifying glass: It’s the tool which allows us to zoom in or zoom out a
part of the drawn network. To access it we must:
- Slider: It’s the most direct way to zoom in and zoom out on the drawing.
It’s found on the taskbar and it indicates the zoom percentage.
- Drag: This tool acts as a complement of the zoom tool, since in many times,
the zoomed area doesn’t correspond to the one we are interested in.
Therefore, by activating the hand option we can adjust and locate the
desired screen without changing zoom or network drawing.
- Element lock: this tool enables the user to introduce the same element
several times on a row, being especially useful when we need to introduce
a large number of elements, for instance a network’s nodes.
- Lock drawing: It’s located on the taskbar and its function is to avoid
unwanted and unintentional modifications to the definitive drawing.
Schematic view allows the user to visualize the network in a much simpler manner,
which can be easily exported to either present as a report or to print it with the
desired format.
138 USERS MANUAL
On the Image tab, as we’ve previously mention, we can find all the available
elements to draw a network. Each set of elements is organized in a tab designed
to display on the left hand side the last user-selected element and on the right
hand side the miniature options of the remaining elements. Each of the elements
has a displayable menu with all the available elements.
139 USERS MANUAL
The entire toolbar is completely configurable, allowing the user to save the button
disposition which better suits his or her needs, having this combination available
at all times.
If we want to restore the default button disposition, push the restore button.
5.2 Labels
ALLIEVI grants every element with a default name, which is shown at the label. To
change the name of an element, we merely have to select it, either with the
selection tool or directly on the data window, and introduce the new name.
By selecting one or several elements and right clicking on one of them, the
program displays a menu with the option “Labels”, which allows the user to
change the label position according to his or her preferences. This action can also
be done through the computer’s numerical keyboard.
141 USERS MANUAL
ALLIEVI gives you the possibility of showing or hiding element labels effectuating
the following actions:
The quick drawing tool allows for the rapid introduction of connected elements,
saving time and clicks, since one can create an entire network without interrupting
the drawing process. The tool is found on the Image tab and in order to activate
it the listed steps must be followed:
Once the Quick drawing tool is activated, we have to select the first element to
be introduced to the network. Automatically, the cursor is set in pipe mode and
by right clicking we can change the selected element.
142 USERS MANUAL
The way to exit the quick drawing mode is either activating once more the
selection tool or from the actual displayable menu, when right clicking. We can
also use a shortcut, pressing the ESC key twice.
Once we’ve drawn an installation, it’s likely that modifications are required. For
instance, placing vertices to a pipe and converting it to a Polypipe. This is achieved
in the following manner:
Pressing the Control (Ctrl) key a pipe is selected, which displays a vertex edit
tool.
When clicking on any of the vertices, the cursor changes to and the
vertices can me moved to the desired place.
To eliminate vertices:
To change the direction of a pipeline, select it, and with the right mouse button
you can access the Change direction option.
144 USERS MANUAL
The change of direction of a pipe affects the definition of its profile, as well as the
position of vents and breakage that could be along its layout.
A pipe can be released from a node where it connects with other elements. To do
so, select the pipeline, and with the right mouse button you access the option
Release pipe. Next, the program requests from what extreme node it is desired to
release the pipe.
A pipe can be divided by a point at a certain distance from its origin. To do this,
select this pipeline, and with the right mouse button you access the Split pipe
option. Next, the program asks how far from the origin you want to divide the
pipeline.
145 USERS MANUAL
When dividing a pipe, both sections will have the same characteristics (diameter,
thickness, celerity and roughness), distributing the coefficient of minor head losses
of the original pipe proportionally to the length of each section. If there are vents
and breakage in the original pipe, its position will be defined in the corresponding
section after the division.
You can join two pipes if they are connected by a node in which there are no
other elements, are defined in the same direction, and have the same
characteristics (diameter, thickness, speed and roughness). To do this, the first
pipeline is selected and, by pressing the Ctrl button, the second pipeline is
selected. Then, with the right mouse button you access the Join pipes option.
The resulting pipe has the same characteristics as the original pipes, its same
direction, the sum of the lengths, and the sum of the minor head loss coefficients.
If there are vents in the original pipes and in one of them a breakage, all this is
relocated in the resulting pipe.
146 USERS MANUAL
ALLIEVI allows for the comparison of different alternatives to solve the same
problem to be done. In order to do so, it possesses a tool which enables the
creation of different operating scenarios from the main one.
For each scenario, we can change different protection elements, change data,
calculate and compare their results.
It’s important to highlight that the main scenario is the one which is shown when
we initialize the program. The steps needed to create a new one are the following
ones:
IMPORTANT: When a new scenario is created, ALLIEVI locks the drawing and
the user can’t change the typology of the network nor input any elements. For
such reason, we recommend that in the main scenario all necessary elements
are drawn, keeping in mind that one can activate or deactivate them, and also
change from one protection element to another using the keyboard shortcuts.
Each scenario differs from the other by their names and the background color
ALLIEVI applies when we create new scenarios, the latter being configurable by
the user in the following way:
6. Results Wizard
To explore these results, you need to select the data and by pressing Ctrl+C you
can copy them other programs like Microsoft Excel.
The table above shows piezometric heights and pressures in permanent regime
for the nodes of the scheme indicated below:
It should be noted that for the program, all elements of the system are connected
between two nodes, one initial and one final, although in the scheme some types
of deposits, protective devices and the laws of Q or H have only one node. These
hidden knots are included in the results table, each of them with its corresponding
150 USERS MANUAL
visible node. For the scheme, we are considering the initial and final knots of the
elements with a single visible node are those shown in the following table:
According to the values indicated in the results table, the pressure in the final
node of the air vessels(node N17_fin with pressure of 68'04 mca), is the pressure
of the gas inside the air vessel.
One the permanent regime is calculated, to access the results of the transient
regime calculation one must:
NOTE: Elements like channels or pipes have plots which are incompatible with
other elements. Hence, ALLIEVI doesn’t allow the selection of these elements
with a single result’s output, and shows an error message.
IMPORTANT: If the user considers that the selected series are the most used
ones for these elements, he or she can activate the option of setting graphs by
default. In this way, the results are saved and one can access them later, without
the need of going through the assistant yet again.
At the left hand side, the user selects the type of series to be represented.
By pressing the Finish button, ALLIEVI displays the graphical and numerical
result.
NOTE: In each of the wizard’s steps, the user can press the “Finish” button,
which exits the wizard and shows the plot without the need of undergoing all
the steps.
Once the graphs are generated, we can modify them without the need of going
back to the wizard:
155 USERS MANUAL
At the Series tab, those which are represented in the graph are shown:
By selecting More series options you can configure the colors of each series,
change the name and decide among other things if you want to draw the series
or not.
156 USERS MANUAL
At the Setting tab it’s possible to change the title of the graph, the axis, configure
the scale, define the font type, the size, etc. To do so, you must select the option
“More title and legend options”
IMPORTANT: Be aware that to change the default scale, you must make sure
that the Autoscale option is deactivated. Otherwise, the initial scale will remain
unchanged.
157 USERS MANUAL
Finally, at the Data tab the user finds the values of each of the plotted series. These
values can also be copied to external programs, like Microsoft Excel.
ALLIEVI has a tool which can quickly access to transient results by merely selecting
the elements that the user wants to analyze. To access them, the user has to right
click on the selected elements and choose the option on the
displayed menu.
Elements
At a point
Envelopes
Animations
158 USERS MANUAL
By default, the most representative plot is shown, or the one the user has
configured and saved on the results wizard.
6.2.3 SmartCharts
This tool enables the user to save a graph and have it always available, in such
way that if the user simulates the project once more, SmartCharts updates it and
shows it with the new obtained results.
Right click Add to “My SmartChart”: Input the name of the graph you
want to save.
159 USERS MANUAL
EPANET and AutoCad are the most widely used software in water engineering.
For this reason, ALLIEVI allows the user to import hydraulic networks from EPANET
(.inp) and AutoCAD (.dxf)
When opening a file with (.inp) extension, ALLIEVI asks for the network’s scale
factor
161 USERS MANUAL
IMPORTANT: Allievi not only imports the outline of the network, but also its
data. However, special care must be taken in those elements that are
incompatible between both programs, such as the case of vents, composite
valves or consumption in nodes, this last one represented in Alievi through
flow laws.
It must also be taken into account that in Epanet a long pipeline may have
been divided into relatively short sections to define its profile. In Allievi, the
short sections must be reduced to a minimum, as in the case of auxiliary pipes,
so there are no mistakes in the simulation of the transient regime. Therefore,
when importing an Epanet network, apart from reviewing and completing the
data, it is necessary to try to join those sections of pipe that, being contiguous,
have the same characteristics.
Once the conversion has been done, either from an EPANET or AutoCAD file, the
user can visualize the project on ALLIEVI with all the elements and their
respective data. In these cases it is always necessary to review the data and
imported elements, to detect possible incompatibilities and correct them
properly.
8. Keyboard shortcuts