Home      |     Contact Us      |     Users Site SEARCH:       


Evaluating Hydraulic Energy Losses and Total Hydraulic Head with
FLOW-3D
's Flux Baffles

This article was contributed by Laurent Bilodeau, ing. Conception des aménagements de production Hydro-Québec Équipement.

Introduction

This note reviews the computation of hydraulic energy flow and total hydraulic head provided by FLOW-3D, with a particular interest for the evaluation of hydraulic energy dissipation rates.

In FLOW-3D v9.3, total hydraulic head was included among the variables that can be visualized directly from model outputs. Figure 1 shows the total hydraulic head distribution across a cut through a river diversion tunnel. In version 10, FLOW-3D introduced hydraulic energy flow and total hydraulic head as integral values that can be computed with flux baffles, visualized as time series, and analyzed with external tools.

Simulation of the distribution of total hydraulic height

Figure 1. The distribution of total hydraulic height within a cross section through a river diversion's tunnel

Total Hydraulic Energy

Bernoulli’s Equation

Hydraulic energy, eG, is the sum of the potential and kinetic energies of particles of water in a flow. It is an energy density, expressed in J/m³, and is given by Bernoulli's equation (Eq. 1):

(1)Bernoulli’s equation

where the symbol meanings are

eG hydraulic energy density (J/m3)
p pressure (Pa ≡ N/m2 ≡ J/m3)
g gravity's acceleration (-9,81 m/s2)
ρ density (kg/m3)
u, v, w velocity in x, y and z (m/s)
z height (m) above some reference level, or altitude

The Hydraulic Energy Simplified Cascade

Energy in general is considered to be a conservative quantity that transforms itself but is never lost. For the purposes of representing water flow in civil engineering, it is often sufficient to conceive of energy transformations as a cascade starting with gravitational potential energy, transforming into kinetic energy and then into heat energy. It is also conventional to further limit the scope of the energy cascade by modeling explicitly the quantities of only the first two forms (potential and kinetic).

The hydraulic energy cascade is demonstrated well by the case of a flow plunging from a weir into a pool, as illustrated in Figure 2.

Partial trajectory of a water particle
Figure 2. Partial trajectory of a water particle going from an upper basin through a weir and to a stilling basin.

Following the track of a particle as indicated in the figure:

While travelling from A to B and C, viscous and turbulent processes usually have little influence on the flow. Total hydraulic energy eG can be treated as a conserved quantity, like mass would be, with small loss terms taken into account if need be. Downstream of C, this model of conservative hydraulic energy can be extended by taking into account energy loss terms of a larger magnitude and effect on the flow.

Mass and Energy Budgets

Control Volumes

Both the transport of eG and mass density ρ have to be monitored for reasons that will become clear later on; this is simply taken for granted now.

The transport of eG and mass density ρ by the flow are easy to analyze with the help of control volumes, denoted as CV below, and Gauss's law of divergence.

A spherical control volume defined so that it intersects a falling jet of water
Figure 3. A spherical control volume defined so that it intersects a falling jet of water

CVs are volumes enclosed by a defining surface of one choice as long as it follows these rules:

The volume of a CV is used to calculate integral, self-preserving quantities such as enclosed mass or energy.

The surface of a CV is used to define fluxes coming in and out, allowing one to make budgets of the enclosed quantities and to monitor their time history.

Figure 3 gives an example of a control volume that could be used to analyze the properties of a falling jet of water. The sole inflows and outflows into this control volume are the jet itself, coming in from the upper left and leaving at the lower right.

Control Volumes of Fixed Shapes in FLOW-3D

With FLOW-3D, CVs of fixed shape and position can be defined easily with the help of flux baffles of three basic shapes:

Figure 4 is an example drawn from an actual model, where one can see the three types of flux baffles, after they have been rendered into a computational mesh. They are rendered as opaque but they are completely permeable to the flow if they are defined as solely flux measuring baffles.

flux baffles of planar, cylindrical and spherical shape rendered within a sample mesh
Figure 4. Example of flux baffles of planar, cylindrical and spherical shape rendered within a sample mesh
Two control volumes defined by joining pairs of vertical cross-sections cutting through a tube or penstock

Figure 5. Two control volumes defined by joining pairs of vertical cross-sections cutting through a tube or penstock. The flow is coloured according to total hydraulic head.

Figure 5 shows how to use planar baffle surfaces to define two control volumes:

The vertical flux baffles in Figure 5 could be used to define several other CVs using permutations of the available vertical surfaces: DB, DA, CB, CA.

Energy Budget

The hydraulic energy balance is leaky by definition because it explicitly labels viscous heat production as a loss. Ideally, the hydraulic energy cascade should lose energy to no other cause, and should not gain any. Here again, practice with numerical models paints a slightly different picture. Any numerical model will have some artificial source or sink of hydraulic energy.

For instance, there is a flow interference with the computational mesh when the cell size is not much smaller than energy carrying flow features. When the cell size is not small enough, velocity contrasts are spread over larger spatial extents than in the natural flow. The spread makes the kinetic energy slightly smaller and acts as an energy dissipation that is due to grid effects rather than to natural phenomena.

Monitoring the energy budget may yield clues about the reliability of the model and can be used to compare runs that use different parameter values or grid cell sizes. When artificial gains and losses are under control, the rate of hydraulic energy dissipation is often one of the important results taken from numerical models and important for discriminating for between design variations.

Total Hydraulic Head

Total Hydraulic Head as an Energy Density

Total hydraulic head, symbolized as hG below, is the total hydraulic energy eG of Eq. 1 recast into an altitude by simply dividing it by (-g ρ ):

(2) Equation for total hydraulic head

 

(3)  Final equation for total hydraulic head

 

where the symbols have all been introduced already except for:

hG, which is total hydraulic head (m)

Total hydraulic head is qualified as total because it is the sum of:

The piezometric head measured in the flow is considered to provide a good measurement of the water's local free surface elevation.

In reservoirs and tranquil reaches of rivers, the flow's velocity may be low enough that the kinetic energy head is considered negligible so that hG is sometimes taken to be equal to the piezometric head.

Total hydraulic head hG is sometimes called the stagnation height: given a particle of fluid within a flow, it is the ultimate height the particle would reach if all its velocity was suddenly directed upwards and the surrounding fluid was not an obstacle.

Total hydraulic head hG is an ubiquitous variable in hydraulic design, such as at bridge piers and abutments. It also stands in for the flow’s hydraulic energy whenever that energy has to be conveyed or dissipated in a managed fashion, such as along in channels and penstocks. hG is a variable of choice because it can appear in engineering drawings as annotations using the same altitude scales as other engineering work’s key heights.

Hydraulic Energy Dissipation from Integral Values of Total Hydraulic Head

The integral approach to energy dissipation occurring between two flow cross-sections A and B is calculated from the decrease in HG in the downward direction of the flow going.

To calculate the energy dissipation between A and B with the help of HG, the value of HG at each cross section is first multiplied by -ρg to cast it into an energy density flow, and then multiplied by Q to cast it into a total hydraulic energy flow.
Taking the difference in energy flows of two cross-sections, in a situation where the volume flow rate Q is the same at both cross-sections, yields the hydraulic energy loss incurred by the flow in moving from the upstream cross-section to the downstream one, as expressed below:

(4)  Energy dissipation

where the newly introduced symbols have the following meanings:

P (W) Loss of hydraulic energy undergone by the flow while travelling between cross-sections A and B
Q (m³/s) Flow of water volume going from cross-section A to cross-section B
HGA, HGB (m) Total hydraulic head at cross-sections A and B

Example 1 - A River Diversion

Figures 6 and 7 illustrate a river diversion through a tunnel. Only the water is shown. A river diversion is used on dam construction sites so that the dam’s foundations can be easily accessed and constructed.

In the example, water flows from left to right. The river is in a state of flood and it is flowing into the upper reach at a high rate. At times when the rate of incoming water is higher than the rate of outflow through the diversion, water accumulates in the upstream reach, the water level rises. The rise in water level increases the pressure on water in the tunnel, which increases the rate of flow in the tunnel, which decreases the rate of accumulation in the upstream reach. This negative feedback eventually equalizes and the rates of inflow and outflow stabilizes. The water level that is reached at peak flood is of much concern for managing the construction site.

A river diversion through a tunnel modeled with FLOW-3D
Figure 6. A river diversion through a tunnel, flowing from left to right, colored according to the flow velocity.
River diversion simulation, showing only the water and the flux measuring cross-sections.

Figure 7. The same river diversion, showing only the water and the flux measuring cross-sections.

The upstream reach presents an example of a volume balance where the water level fluctuates according to the difference between the incoming and outgoing flows. At the same time, the rate of flow through the tunnel results from a pressure balance between the upstream and downstream ends in which friction on the tunnel walls and flow energy losses at discrete structures and transitions along the passages play a large part. Figure 10 shows numerous flow measuring cross sections, known as flux baffles in FLOW-3D. Their uses are many:

Example 2 – High Velocity Free Run Over a Natural Rock Surface

Figure 8 shows an example of a spillway’s tailrace and free run over natural rock.

The model was aimed at evaluating the rate of hydraulic energy dissipation per unit area of ground surface. This rate, in W/m², was an input for evaluations of the erosion potential of the rock surface along the free run.

A spillway’s tailrace and free run over rock
Figure 8. An example of how flux baffles can be used is given by a spillway's tailrace over natural rock. The aim is to evaluate the rate of hydraulic energy dissipation per unit area of ground surface.
Instrumentation installed in the spillway’s tailrace and free run for assessing rates of hydraulic energy dissipation per unit surface

Figure 9. This simulation indicates where planar and cylindrical flow baffles were located to provide estimates of energy dissipation.

Figure 9 presents the metering apparatus that was used to ascertain the tools themselves and to perform the evaluations of the dissipation rate. The mesh blocks are also outlined.

An assortment of planar and cylindrical flux baffles were deployed for measuring the desired rate of dissipation and, just as importantly, for assessing the quality of the flow and measuring tools.

The planar flux baffles allow one to construct control volumes and use them to monitor the steadiness of the volume flow and the dissipation of energy within the CV. In the tailrace, the flow is well contained by sidewalls and fairly uniform across the cross-sections. The energy dissipation rates were of the order of 25 to 50 kW/m².

A cylindrical flow baffle is located at the foot of the discharge gates. The mean volume flow rate through the cylinder was variable through time because of normal flow fluctuations, but it did tend to zero when a suitable averaging interval was taken. When the same average was taken for the net hydraulic energy flow through the baffle, it yielded as expected a negative value which, divided by the area, was close to 30 kW/m². Another cylindrical flow baffle, extended into an elliptical horizontal cross-section, was also located near the end of the tail race. Similar verifications made there also showed a similar agreement.

It was concluded that the cylindrical flow baffles worked as expected for measuring hydraulic energy dissipation. Attention was then turned to the cylindrical baffles placed in the downward slope, where the dissipation was expected to be the highest.

Figure 10 shows the time series of the net volume and energy flows through cylindrical baffle number 3, located over the free run, near the transition from the tailrace to the natural rock surface. The figure shows that high amplitude fluctuations of both flows are present and any tendency is rather well hidden by them.

Net volume flow and net hydraulic energy flow
Figure 10. Net volume flow (m³/s) and net hydraulic energy flow (W) through cylindrical flux baffle number 3 of the second example, as a function of time.

Figure 11 presents the time integration of the net volume and energy fluxes of Figure 10. The time integration yields values of volume (m³) and energy (J). The volume time series allows one to select a time interval where the regime is stationary, and to choose integration time boundaries so that the net volume change is close to zero. The energy time series shows a regular downwards trend, as expected for energy dissipation. The slope of the trend, in W/s, gives an estimate of the rate of dissipation. One can then divide it by the surface area of the base of the cylindrical baffle enclosure to get the desired dissipation rate per unit of area. The radii of the cylinder were chosen so that the area would be close to 100 m². In this case, the rate of dissipation found was 286 kW/m².

Accumulated water volume and hydraulic energy
Figure 11. Volume and total energy within the volume of the cylindrical baffles as a function of time, obtained by time integration, offset by an arbitrary value so that the zero value is simply the initial quantity present, and unknown.

These results were then used in the discussions with engineers from other disciplines. It was made clear that the margin of uncertainty was large since, among other factors, the flow computation was not aerated ant the rock surface resulted from preliminary estimations. It was also pointed out by the modeller that hydraulic energy dissipates in the water, not in the rock, and that the location of maximum dissipation is not necessarily that of maximum action on the rock. The minutia of the analysis presented above was a significant burden but was deemed necessary as the baffles were being used in a way that was novel for the modeller. The cross-disciplinary discussions and the order of magnitude of the values were the most useful results of the model, rather than the numerical values themselves.

Conclusion

FLOW-3D's flux baffles offer precise evaluations of volume and hydraulic energy net flows through them. Their computation algorithm is finely adjusted to the underlying numerical scheme of FLOW-3D used with the control volume approach, they are well designed for numerous assays including verifying the performance of FLOW-3D itself with regards to mass conservation in situations demanding a high level of consistency.

The computation of total hydraulic head can be made in numerous ways, a situation that is not surprising in view of the very high usefulness of the quantity to civil and hydraulics engineers. One of the ways offered by FLOW-3D is the calculation over the flow area of a baffle of the flux averaged total hydraulic head. It is shown here that hydraulic energy rates of loss in a flow measured with the differences between the values at two flux baffles crossing a given flow tube will be exactly those that can be calculated with a control volume approach linked to primitive flow variables by the Gauss divergence theorem.

^ back to top