## Abstract

Knowing the scarcity of water in the southeast of Spain and how the rain occurs, we considered the design of intake systems in ephemeral riverbeds in order to try to capture part of the runoff flow. The intake systems generally consist of a rack located in the bottom of a river channel, so that the water collected passes down the rack and leads to the side channel. This behaviour has been studied in the laboratory by several researchers. However, due to the many effects that occur on the bars, it is not possible to analyse the whole problem of characterization with traditional methodologies. For instance, the wetted rack length necessary to collect a required flow presents important differences depending on what each author has considered relevant. Computational fluid dynamics simulations have been done to improve the knowledge of the hydraulic phenomenon observed in different laboratory experiences, for which we have previously calibrated the numerical models using laboratory results. The ANSYS CFX code was selected. Several two-equation turbulence models have been considered. The results show differences smaller than 1% in the wetted rack length, and discharge coefficients also present good agreement.

## INTRODUCTION

Intake systems generally consist of a rack located in the bottom of the channel that allows water to pass through. These structures have been adopted in small mountain rivers with steep slopes and an irregular riverbed, intense sediment transport and flash flood. Their design is intended to satisfy two primary objectives: (1) to maximize water intake; and (2) to minimize sediment intake.

In designing intake systems, we need to consider geomorphologic, hydraulic, structural and economic aspects to avoid unnecessary maintenance and functionality problems during the lifetime of the project. The efficiency of racks depends on a number of factors, such as the number of bars, incoming flow conditions, longitudinal inclination, shape and spacing between bars.

The hydraulic behaviour of the racks is also influenced by the bars’ disposition. In longitudinal bars, the flow collected appears as a function of the local energy flow. However, in transversal bars or circular perforations, the flow collected is related to the local flow level (Mostkow 1957).

In the analysis of clear water flows it is assumed that the flux over the rack is one-dimensional, the flow decreases progressively, and the hydrostatic pressure distribution acts over the rack in the flow direction. Two broad approaches to dealing with the energy head over the rack are typically used and these are presented in Table 1.

Horizontal energy level . | Energy level parallel to the rack plane . |
---|---|

Kuntzmann & Bouvard (1954), Frank & Von Obering (1956), Free overfall (Henderson 1966), Vargas (1998), Drobir et al. (1999), Brunella et al. (2003) and Righetti & Lanzoni (2008) | Chaguinov (1937), Noseda (1956a, 1956b), Gherardelli (1956), Mostkow (1957), Dagan (1963) and Krochin (1978) |

Horizontal energy level . | Energy level parallel to the rack plane . |
---|---|

Kuntzmann & Bouvard (1954), Frank & Von Obering (1956), Free overfall (Henderson 1966), Vargas (1998), Drobir et al. (1999), Brunella et al. (2003) and Righetti & Lanzoni (2008) | Chaguinov (1937), Noseda (1956a, 1956b), Gherardelli (1956), Mostkow (1957), Dagan (1963) and Krochin (1978) |

The classical approach considers a two-dimensional perspective. However, when analyzing the flow near the solid edges, the flow becomes extremely three-dimensional, rendering the two-dimensional analysis tools less useful.

In numerical modelling, different approaches have been used. Some codes, such as HEC-RAS, MIKE, TELEMAC or IBER, consider one-dimensional and two-dimensional models to study the hydraulic behaviour of rivers and sediment transport phenomena (Castillo *et al.* 2015; USACE 2016). Nevertheless, flows occurring in hydraulic structures tend to be highly three-dimensional (Bayon *et al.* 2016; Castillo *et al.* 2016). For this reason, three-dimensional models were considered in this study.

The rest of this paper is organized as follows. The Purpose section indicates the challenges and objectives of the study. The Theoretical background section is a literature review to familiarize the reader with the empirical expressions obtained by several authors to calculate intake systems. The Methods section introduces the main characteristics of the laboratory equipment and of the numerical simulations carried out. The Results section compares the laboratory and simulated data in terms of mesh size sensitivity, turbulence models, flow profiles over the rack, wetted rack length, water collected and discharge coefficient. The Discussion section explains the main findings related to the purpose of the study. Finally, the Conclusions section briefly summarizes the overall conclusions of the work.

## PURPOSE

As a result of the lack of water resources in several parts of the world, bottom intake systems are a point of great interest. Research to date has been developed in the laboratory. Thanks to the computational fluid dynamics (CFD) programs, researchers and designers can evaluate different effects with a smaller cost than that incurred by building scale models. However, references to numerical simulations of bottom intake systems are scarce.

As a result of the lack of numerical simulations in this field, the present paper is focused on filling the gap between laboratory results and CFD simulations in T-shape bottom racks. ANSYS CFX software (version 16.2) is used to analyse the flow profiles over the rack, and the collected water along the intake system. Due to the variety of turbulence models, and in order to compare their effect over different parameters, three distinct two-equation turbulence models are examined.

Knowing the parameters analysed, designers will be able to design the bottom intake system with greater certainty.

## THEORETICAL BACKGROUND

### Wetted rack length

Several researchers have estimated the theoretical wetted rack length *L* necessary to derive a defined flow rate *q*_{1}. The different required lengths are due to the variation of experimental conditions used to adjust the discharge coefficient, such as the shape of the bars, their separation and width, the void ratio, the approximation flow conditions, the initial flow depth *h*_{1} or the longitudinal rack slope *θ*.

Figure 1 differentiates the wetted rack length in two magnitudes: the distance along the rack where the nappe enters directly through the racks (measured between the bars), called *L*_{1}, and the maximum distance where the bars are wet, called *L*_{2}.

The theoretical longitudinal rack slope, *θ*, has been considered in several ways. Some authors consider the influence of the slope in the required wetted rack length (Garot 1939; Orth *et al.* 1954; White *et al.* 1972; Righetti & Lanzoni 2008). Other researchers experimentally found that there is no further influence for slopes greater than 19° (34.43%) (Brunella *et al.* 2003).

Castillo & Lima (2010) and García (2016) compared the wetted rack length proposed by several authors, for the case of critical conditions at the beginning of the rack, void ratio *m* = 0.60, and longitudinal rack slope of 20% (Figure 2). In each case, the horizontal energy level remains constant along the rack. Differences may be double in some cases.

In Figure 2, the length proposed by Righetti & Lanzoni (2008) approaches the free overfall (Henderson 1966). These authors introduced the approximation level conditions in their formula and critical approximation conditions were adopted. These differ from other experimental results. There are also some field measurements in a prototype made with circular bars (Drobir *et al.* 1999).

To avoid rack occlusion, several authors proposed design recommendations from prototype observations in mountain rivers (e.g., Orth *et al.* 1954; Ract-Madoux *et al.* 1955; Krochin 1978; Drobir 1981; Bouvard 1992; Raudkivi 1993). Castillo *et al.* (2016), from experimental measurements in a flow with gravel-sized materials, determined that a longitudinal rack slope around 30% minimizes the occlusion effect.

### Discharge coefficient

*q*

_{d}*=*

*dq/dx*(m

^{3}/s/m), may be calculated using the orifice equation. The collected flow through the rack plane is influenced by the velocity distribution close to it. The deviation from the uniform distribution is included in the velocity coefficient,

*C*. In a similar way, there is a change in the available section, which involves the flow contraction. In this case, a contraction coefficient,

_{v}*C*, may be considered. Both coefficients require experimental measurement and depend on the shape of the bars and the spacing between them. The orifice equation may be written as: where

_{c}*m*represents the relation void area and the total area (void ratio),

*H*the total energy available refers to the plane of the rack,

*x*the longitudinal coordinate along the rack, and

*C*the discharge coefficient depending on the energy height (

_{qH}*C*=

_{qH}*C*).

_{c}C_{v}Table 2 shows the discharge equation, considered by several authors.

*Note:**h(x)* is the normal water surface height depending on the *x* coordinate, *C _{qh}(x)* the discharge coefficient depending on

*h*,

*C*the discharge coefficient depending on the energy and the

_{qH}(x)*x*coordinate,

*H*

_{0}the energy height at the beginning of the rack,

*C*the discharge coefficient measured under static conditions (with negligible approaching velocity), and

_{q0}*θ*the angle of the rack with the horizontal.

*et al.*(2000) considered that it is possible to calculate the differential flow of the water collected as: where

*m*is the void ratio,

*dx*the increment longitudinal in the flow direction,

*H*

_{0}the specific energy at the beginning of the rack, Δ

*z*the vertical difference between the initial rack section and the analysed section, and

*C*is the discharge coefficient.

_{qH}The same authors proposed that *C _{qH}* ≈ sin

*θ*, with

*θ*being the angle between the velocity vector of water derived and the plane of the rack (Figure 3).

### Flow profile over the rack

The flow profile over the rack has also been analysed by several authors (Table 3).

Author . | Flow profile over the rack (dh/dx)
. |
---|---|

Garot (1939) | |

De Marchi (1947) | |

Noseda (1956a, 1956b) | |

Frank & Von Obering (1956) | ; |

Dagan (1963) | |

García (2016) | |

García (2016) |

Author . | Flow profile over the rack (dh/dx)
. |
---|---|

Garot (1939) | |

De Marchi (1947) | |

Noseda (1956a, 1956b) | |

Frank & Von Obering (1956) | ; |

Dagan (1963) | |

García (2016) | |

García (2016) |

With velocity measurements in the free surface, Brunella *et al.* (2003) found that the dissipation effects are insignificant. However, in the final part of the racks these effects cannot be neglected since the local effects generate friction effects. Differences between measured and calculated depth profiles at the beginning of the rack are due to the consideration of hydrostatic pressure distribution.

## METHODS

### Physical device

An intake system based on the physical model analysed by Noseda (1956a, 1956b) was built in the Hydraulic Laboratory of the Universidad Politécnica de Cartagena (Figure 4). It consists of a 5.00 m long and 0.50 m wide approximation channel, a rack with different slopes (from horizontal to 33%), the discharge channel, and the collected discharge water channel. The model was made with methacrylate walls allowing for direct observation of the flow.

The racks were located at the bottom of the channel. They were built with aluminium bars with T profiles (T 30/25/2 mm) with the spacing parallel to the direction of the flow. The bars used had the same width, but the longitudinal layout was modified to consider different spacing between them.

The rack was 0.90 m in length. The spacing between bars was 11.70 mm. With those considerations, the void fraction analysed was *m* = 0.28. The inlet flow rates were measured by an electromagnetic flowmeter Endress Häuser Promag 53 W of 125 mm with an error of 0.5% of the full scale. The rejected flow rates were measured by a 90 V-notch weir. The collected or derived flows were obtained as the difference between them. Once steady conditions were reached, a vertical point gauge with an accuracy of 0.5 mm was used to measure the free surface flow profile. Further details on the experimental setup can be found in Castillo *et al.* (2017b).

The inflow was subcritical in all the cases at the beginning of the inlet channel. The flow reached supercritical conditions at the beginning of the rack. In each experiment, the inlet and exit flows and the longitudinal flow profile were measured. *q*_{1} is the entrained specific flow, while *q*_{2} is the specific discharge flow, and *q _{d}* is the specific discharge flow collected in the intake system.

### Numerical modelling

CFD programs allow for the simulation of the interaction among different fluids as a two-phase air-water flow, or flows with different concentrations in the case of sediment transport. The programs solve the fluid mechanic problem in various geometric configurations, providing a great deal of data, with greater flexibility and speed than that obtained with experimental procedures. However, mathematical models still present accuracy issues when modelling some hydraulic phenomena (Chanson & Gualtieri 2008). For correct use, it is necessary to contrast and to calibrate with data obtained in prototypes or physical models.

To test the hydraulic behaviour of the intake system, laboratory measurements were used to model and calibrate CFD simulations. The Finite Volume Scheme program ANSYS CFX (version 16.2) has been used. This program was previously used for solving intake systems with accurate results (Castillo *et al.* 2014b).

*i*and

*j*are indices,

*x*represents the coordinates directions (

_{i}*i*= 1 to 3 for

*x*,

*y*,

*z*directions, respectively),

*ρ*the flow density,

*t*the time,

*U*the velocity vector,

*p*the pressure, presents the turbulent velocity in each direction (

*i*= 1 to 3 for

*x*,

*y*,

*z*directions, respectively),

*μ*is the molecular viscosity, is the mean strain-rate tensor and is the Reynolds stress.

Although the RANS equations can be applied to variable-density flows, in this case Navier–Stokes equations are considered in their incompressible form.

### Turbulence models

To reach closure of the Navier–Stoke equations, one-equation to direct simulation turbulence models can be used.

*μ*being the eddy viscosity or turbulent viscosity, the turbulent kinetic energy and

_{t}*δ*the Kronecker delta function.

The choice of the turbulence model may have great importance in some studies (Castillo & Carrillo 2017). In this work, some of the most usual RANS turbulence models have been tested.

#### Standard k-ɛ model

*C*is an empirical coefficient and

_{μ}*ɛ*is the dissipation rate of turbulent kinetic energy.

#### Re-normalization group k-ɛ model

*k-ɛ*model is more accurate than the standard

*k-ɛ*model. The RNG

*k-ɛ*turbulence model is based on a RNG analysis of the Navier–Stokes equations. The transport equations for turbulence generation and dissipation are the same as those for the standard model, although the model constants differ and the constant is replaced by the function (Yakhot & Smith 1992): with this being: where = 0.085 is a closure coefficient, a constant with a value of 0.012,

*P*the turbulence production due to viscous forces (ANSYS Inc. 2015),

_{k}*ρ*the flow density and

*ε*the dissipation rate of turbulent kinetic energy.

#### k-ω based shear-stress transport model

*k-ω*turbulence models assume that the turbulence viscosity is linked to the turbulence kinetic energy,

*k*, and the turbulent frequency,

*ω*, as: The shear-stress transport (SST) model takes into account the accuracy of the

*k-ω*model in the near wall region and the free stream independence of the

*k-ɛ*model in the outer part of the boundary layer. To do this, the original

*k-ω*model (Wilcox 2006) is multiplied by a blending function

*F*

_{1}, while the

*k-ɛ*model (Launder & Sharma 1972) is transformed to a

*k-ω*formulation and multiplied by a function 1−

*F*

_{1}(Menter 1994).

*F*

_{1}is designed to be 1 inside the boundary layer and decreases to a value of 0 away from the surface. If Φ

_{1}represents any constant in the original

*k-ω*model and Φ

_{2}represents any constant in the transformed

*k-ɛ*model, then the corresponding constant in the model Φ

_{3}may be written as: From the three two-equation turbulence models considered, this model has obtained the best results in the study of free surface falling jets (see Castillo

*et al.*2014a, 2017a).

### Convergence criteria

In judging the convergence of a solution in a finite volume scheme, a widely used method entails monitoring the residuals. Residuals are defined as the imbalance in each conservation equation following each iteration. The solution is said to have converged if the scaled residuals are smaller than prefixed values ranging between 10^{−3} and 10^{−6}. In this work, the residual values were set to 10^{−4} for all the variables.

### Free surface modelling

To solve the air-water two-phase flow, the Eulerian–Eulerian multiphase flow homogeneous model was selected. In each control volume, the sum of volume fraction of all phases (*r _{α}*) is the unit. Three conditions are possible for each cell:

*r*= 0. The cell is empty of the_{α}*α*phase.*r*= 1. The cell is full of the_{α}*α*phase.0 <

*r*< 1. The cell contains the interface between the_{α}*α*phase and one or more other phases.

It may be assumed that the free surface is on the 0.5 air volume fraction.

### Boundary conditions

The model boundary conditions corresponded to the flow, the turbulence at the inlet condition obtained with an acoustic Doppler velocimeter (located 0.50 m upstream of the front edge of the rack), the upstream and downstream water levels, and their hydrostatic pressure distributions. In the bottom of the exit channel of water collected, opening boundary conditions were used. For simplicity, it was considered that all the longitudinal bars work in the same mode in the intake system. For this reason, the domain fluid considers three bars and two spaces between bars. Symmetry conditions were used in the central plane of the extreme bars (Figure 5).

For the walls and the rack, no slip wall conditions and smooth walls were considered. The atmospheric condition was simulated as an opening condition with a relative pressure of 0 Pa, air volume fraction of 1 and water volume fraction of 0.

## RESULTS

Before studying all the range of flows and slopes in the intake system, mesh sensitivity analysis and turbulence model comparisons were performed. Once the best options were obtained, the different settings of the physical model were analysed. Laboratory and simulated data were compared at the same geometric scale after the steady state was reached.

### Mesh size independence

In the study of intake systems, there is flow separation and high turbulence that need high quality mesh elements to solve the problem with the highest accuracy. For these reasons, hexahedral mesh elements were used.

To determine the accuracy of the numerical simulations data, in the first place, the flow profiles along the rack were compared by using three different mesh sizes (0.002, 0.004 and 0.008 m). Table 4 shows the number of elements and the mean time required to solve the simulations using an 8-core Intel(R) Xeon(R) CPU at 2.40 GHz and with 12 GB of RAM.

Mesh size (m) . | Number of elements . | Mean required time . |
---|---|---|

0.002 | 2,501,741 | 22 h |

0.004 | 313,342 | 4 h |

0.008 | 51,393 | 2 h |

Mesh size (m) . | Number of elements . | Mean required time . |
---|---|---|

0.002 | 2,501,741 | 22 h |

0.004 | 313,342 | 4 h |

0.008 | 51,393 | 2 h |

Figure 6 compares the flow profiles measured in the laboratory over the centre of the bars with the simulated data obtained for the three mesh sizes considered. Table 5 shows the comparison of water depths for the *q*_{1} = 155.4 l/s/m. In all the cases, the water profiles obtained with the CDF methodology were similar to the laboratory measurements (differences smaller than 0.70% of the characteristic length over bars, *L*_{2}, for the 0.004 and 0.002 m mesh sizes and smaller than 1.30% for the 0.008 m mesh size). The water depths obtained with the 0.008 m mesh size tend to be slightly smaller than the measurements and the other mesh sizes. This seems to be related to the way in which CFX obtains the interface capture approach of the free surface (ANSYS Inc. 2015). When the mesh sizes increases, the precision of the interface tends to be smaller. There are no outstanding differences between the results obtained with the 0.004 and 0.002 m mesh sizes.

Rack section (m) . | Water depth (m) . | |||
---|---|---|---|---|

Laboratory . | Mesh = 0.002 m . | Mesh = 0.004 m . | Mesh = 0.008 m . | |

−0.50 | 0.142 | 0.142 | 0.142 | 0.142 |

0.00 | 0.097 | 0.095 | 0.095 | 0.094 |

0.15 | 0.070 | 0.066 | 0.066 | 0.060 |

0.30 | 0.045 | 0.039 | 0.039 | 0.035 |

0.45 | 0.029 | 0.024 | 0.023 | 0.017 |

0.60 | 0.018 | 0.014 | 0.013 | 0.005 |

0.75 | 0.009 | 0.007 | 0.007 | 0.002 |

0.90 | 0.002 | 0.002 | 0.002 | 0.000 |

Rack section (m) . | Water depth (m) . | |||
---|---|---|---|---|

Laboratory . | Mesh = 0.002 m . | Mesh = 0.004 m . | Mesh = 0.008 m . | |

−0.50 | 0.142 | 0.142 | 0.142 | 0.142 |

0.00 | 0.097 | 0.095 | 0.095 | 0.094 |

0.15 | 0.070 | 0.066 | 0.066 | 0.060 |

0.30 | 0.045 | 0.039 | 0.039 | 0.035 |

0.45 | 0.029 | 0.024 | 0.023 | 0.017 |

0.60 | 0.018 | 0.014 | 0.013 | 0.005 |

0.75 | 0.009 | 0.007 | 0.007 | 0.002 |

0.90 | 0.002 | 0.002 | 0.002 | 0.000 |

As there are no remarkable differences between the results obtained with the intermediate and the smaller mesh size, but the required time is around five times longer with the 0.002 m mesh size, the 0.004 m mesh size was used to analyse the different specific flows and rack slopes.

### Turbulence model

The influence of the turbulence model has also been analysed. Three different turbulence models have been considered: the standard *k-ɛ*, the RNG *k-ɛ*, and the SST models.

Figure 7 compares the flow profiles over the centre of the bars, measured and simulated. For the cases considered, there are no remarkable differences between the three turbulence models.

Although there are no outstanding differences, the SST turbulence model was selected to analyse the different specific flows and rack slopes, based on previous studies comparing the three turbulence models in free falling flows (Castillo *et al.* 2014a, 2017a).

### Flow profiles over the rack

Once the mesh size and the turbulence model had been tested, simulations were carried out with five specific flows and five slopes.

In order to know the accuracy of the numerical simulations data, the longitudinal flow profiles calculated with CFD were compared with the experimental measurements in the laboratory. Figures 8 and 9 show the longitudinal flow profiles over the rack for three specific flows (77.0, 114.6 and 155.4 l/s/m), measured and simulated over the centre of the bars. In both cases, the values show a good agreement.

### Wetted rack length

Knowing the required length of the rack necessary to derive a determinate flow is essential for the designers. Experimental measurements and CFD simulated values of *L*_{1} and *L*_{2} have been used to compare the required rack length for several specific flows and slopes (Tables 6 and 7). Differences between measured and calculated rack lengths are around 1% of the laboratory value for all the cases considered. Those differences tend to increase as the input discharge decreases. This seems to be related to the associated reduction of the water depth, more accurate results being obtained with higher water depth situations.

Rack slope = 10% . | ||||||
---|---|---|---|---|---|---|

q_{1} (l/s/m)
. | L_{1_lab} (m)
. | L_{2_lab} (m)
. | L_{1_CFD} (m)
. | L_{2_CFD} (m)
. | Error L_{1} (%)
. | Error L_{2} (%)
. |

155.4 | 0.812 | >0.90 | 0.808 | >0.90 | 0.485 | – |

138.88 | 0.744 | >0.90 | 0.740 | >0.90 | 0.532 | – |

114.6 | 0.636 | >0.90 | 0.633 | >0.90 | 0.609 | – |

77.00 | 0.457 | 0.839 | 0.453 | 0.835 | 0.849 | 0.473 |

53.8 | 0.334 | 0.664 | 0.330 | 0.660 | 1.183 | 0.598 |

Rack slope = 10% . | ||||||
---|---|---|---|---|---|---|

q_{1} (l/s/m)
. | L_{1_lab} (m)
. | L_{2_lab} (m)
. | L_{1_CFD} (m)
. | L_{2_CFD} (m)
. | Error L_{1} (%)
. | Error L_{2} (%)
. |

155.4 | 0.812 | >0.90 | 0.808 | >0.90 | 0.485 | – |

138.88 | 0.744 | >0.90 | 0.740 | >0.90 | 0.532 | – |

114.6 | 0.636 | >0.90 | 0.633 | >0.90 | 0.609 | – |

77.00 | 0.457 | 0.839 | 0.453 | 0.835 | 0.849 | 0.473 |

53.8 | 0.334 | 0.664 | 0.330 | 0.660 | 1.183 | 0.598 |

Rack slope = 30% . | ||||||
---|---|---|---|---|---|---|

q_{1} (l/s/m)
. | L_{1_lab} (m)
. | L_{2_lab} (m)
. | L_{1_CFD} (m)
. | L_{2_CFD} (m)
. | Error L_{1} (%)
. | Error L_{2} (%)
. |

155.4 | 0.711 | 0.781 | 0.707 | 0.777 | 0.520 | 0.485 |

138.88 | 0.653 | 0.731 | 0.650 | 0.727 | 0.577 | 0.521 |

114.6 | 0.569 | 0.646 | 0.565 | 0.643 | 0.664 | 0.589 |

77.00 | 0.408 | 0.512 | 0.404 | 0.508 | 0.918 | 0.741 |

53.8 | 0.292 | 0.397 | 0.289 | 0.393 | 1.065 | 0.954 |

Rack slope = 30% . | ||||||
---|---|---|---|---|---|---|

q_{1} (l/s/m)
. | L_{1_lab} (m)
. | L_{2_lab} (m)
. | L_{1_CFD} (m)
. | L_{2_CFD} (m)
. | Error L_{1} (%)
. | Error L_{2} (%)
. |

155.4 | 0.711 | 0.781 | 0.707 | 0.777 | 0.520 | 0.485 |

138.88 | 0.653 | 0.731 | 0.650 | 0.727 | 0.577 | 0.521 |

114.6 | 0.569 | 0.646 | 0.565 | 0.643 | 0.664 | 0.589 |

77.00 | 0.408 | 0.512 | 0.404 | 0.508 | 0.918 | 0.741 |

53.8 | 0.292 | 0.397 | 0.289 | 0.393 | 1.065 | 0.954 |

### Water collected

Figures 10 and 11 compare the water collected along the rack length for several specific flows and two different rack slopes. In both cases, the ratio flow entrained–flow collected is almost the same. In the laboratory, it was only possible to measure in the first 50 cm of the rack length and the total water collected flow. After this distance, the water flows were very small and it was therefore complicated to split the increments in the collected and the rejected flows. As a result of the shape of the T bars, the surface tension phenomena tend to cause high values of *L*_{2}, even when more than 95% of *q*_{1} is derived in the vicinity of *L*_{1}.

### Discharge coefficient

*θ*, has been obtained in the numerical simulations as an estimator of the discharge coefficient along the rack. Figures 12 and 13 compare the results obtained for different specific flows and rack slopes. The numerical results have been compared with the empirical expression proposed by García (2016) for T shape bars:

Equation (12) was adjusted by using experimental measurements of the derived flow in Equation (3). The hypothesis of *H _{0}* =

*H*for the incoming flow was adopted for this adjustment. The discharge coefficient proposed in Equation (12) avoids the requirement of the water depth along the rack to calculate the derived flow.

_{min}Table 8 shows the differences between the energy head measurements at the beginning of the rack, *H*_{0}, and the minimum energy, *H _{min}*. To obtain the empirical discharge coefficient related to the real energy head, Equation (12) has been multiplied by the square root of the ratio

*H*

_{min}/H_{0}.

tan θ
. | H_{min}/H_{0}
. |
---|---|

0.1 | 0.902 |

0.3 | 0.762 |

tan θ
. | H_{min}/H_{0}
. |
---|---|

0.1 | 0.902 |

0.3 | 0.762 |

Despite the fact that we have used different bars types, and flow conditions from Righetti & Lanzoni (2008), similar results are obtained in terms of verifying *C _{qH}* ≈ sin

*θ*. At the beginning of the rack, the behaviour is similar to a free overfall, requiring some distance to change the streamlines, and to reach the maximum values of the

*C*. In each rack slope, the maximum values seem to obtain similar results with the different inlet flows tested. The rack slope tends to modify the maximum value. Hence, bigger rack slopes tend to obtain smaller maximum

_{qH}*C*values. After the maximum, the discharge coefficient tends to reduce with the decreasing of the water depth over the rack. Differences to the empirical formulae proposal by García (2016), corrected by the factor shown in Table 8, are due to the author considering that the energy remains constant along the rack.

_{qH}## DISCUSSION

Intake systems are used in stepped streams to obtain water, minimizing the sediment collection. The basic idea may be extrapolated to semiarid regions, in which the lack of water is a problem to deal with.

Boundary layer separation and high turbulence appear in bottom intake systems. In addition, the water depths in some sections of the rack and the spacing between bars make it difficult to use intrusive instrumentation.

Numerical simulations can help in the choice of the right intake system design but prior to this they should be compared with experimental and field data.

In this paper, the accuracy of ANSYS CFX to solve an intake system has been tested. Simulations have been compared with laboratory data, obtaining good agreement in several parameters. This study complements previous laboratory works (Noseda 1956a, 1956b; Brunella *et al.* 2003; Righetti & Lanzoni 2008; Castillo *et al.* 2014b, 2016, 2017b; García 2016) in bottom intake systems with T-shaped bars disposed longitudinally.

From the two classical assumptions shown in Table 1, numerical results show that the energy level is not constant along the rack. The energy level is somewhere between both simplifications, closer to being parallel to the rack plane than to the horizontal energy level (see Table 8).

According to different authors (see Figure 2), the rack length required to collect a desirable flow may vary up to two-fold. This range of possible values generates uncertainty in designers. In these situations, numerical models previously validated may reduce the lack of confidence during the design of intake systems and complement the laboratory studies.

Numerical simulations require a mesh sensitivity analysis. With mesh sizes of 0.004 m (around 1/3 of the spacing between bars in this case), the results are independent of the mesh size. Smaller mesh sizes would require more computational effort without remarkable differences in the results. Larger mesh sizes tend to obtain flow profiles below the laboratory measurements.

Regarding the choice of turbulence models, three different turbulence models have been considered. The results obtained are almost the same, with no significant differences between them. Based on previous studies in free-falling jets, the SST turbulence model was selected.

Once the sensitivity analysis was done, different comparisons were been considered. Regarding the flow profiles over the rack, the results are in agreement with the laboratory measurements for all the cases analysed. In the same way, the maximum wetted rack lengths over the bars and over the spacing diverge around 1% from the laboratory measurements. The results obtained for the water collected along the rack are also quite similar.

Finally, the sinus of the angle of the velocity vector of water collected with the rack plane has been considered as a discharge coefficient along the rack. The results are in agreement with previous studies published by Righetti & Lanzoni (2008) and by García (2016).

## CONCLUSIONS

The results obtained with ANSYS CFX offered good agreement with the laboratory measurements and the empirical formulae (differences smaller than 1%). In this way, simulations can be used to complement laboratory and empirical data, allowing a better design to be obtained.

Numerical results are sensitive to the mesh size considered. In this way, the mesh size needs to be reduced until no significant differences are obtained between the results.

The extensive typology of bar shapes and settings of bottom intake systems do not allow to directly extrapolate results to other rack types. To reach accurate numerical results, several laboratory data are required, such as the water depth at the inlet channel or the flow profile along the rack.

With the aim of improving the design of bottom intake systems, it would be necessary to provide advances in the knowledge and characteristics of the flow through the rack. More experimental studies, both in physical models and prototypes, are necessary in characterizing the flow, combining measurements of pressures, velocities and flow profiles in T-shape and other types of bars.

## ACKNOWLEDGEMENTS

The authors are grateful for the financial support received from the Seneca Foundation of Región de Murcia (Spain) through the project ‘Optimización de los sistemas de captación de fondo para zonas semiáridas y caudales con alto contenido de sedimentos. Definición de los parámetros de diseño’, reference: 19490/PI/14.

## REFERENCES

*.*

*Prise d'eau du type tyrolien (Tyrolean Water Intake)*

*Estudio experimental y numérico de los sistemas de captación de fondo (Experimental and Numerical Study of the Bottom Intake Systems)*

*PhD Thesis*

*Energia Elettrica*, 1347 (in Italian)

*,*2nd edn.

*.*

*HEC-RAS River Analysis System. Hydraulic Reference Manual*. Version 5.0