# Numerical Simulation of the Unsteady 3D Flow in Vertical Slot Fishway-The Impact of Macro-Roughness

Gérard Pineau, Aurélien Ballu, Laurent David, Damien Calluaud

► To cite this version:

Gérard Pineau, Aurélien Ballu, Laurent David, Damien Calluaud. Numerical Simulation of the Unsteady 3D Flow in Vertical Slot Fishway-The Impact of Macro-Roughness. Water, 2025, 17, pp.1088. 10.3390/w17071088. [hal-05104168]

HAL Id: hal-05104168

https://hal.science/hal-05104168v1

Submitted on 9 Jun 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Numerical Simulation of the Unsteady 3D Flow in Vertical Slot Fishway—The Impact of Macro-Roughness

Gérard Pineau<fnref n="1" />*, Aurélien Ballu<fnref n="2" />, Laurent David<fnref n="1" /> and Damien Calluaud<fnref n="1" />

**Abstract:** Vertical slot fishways (VSFs) are crossing devices that are built on rivers or streams. They were initially designed to help salmons to complete their migratory cycle by crossing permanent obstructions. In order to favor the passage of smaller or benthic species, stones or concrete cylinders, called macro-roughnesses, are often inserted at the bottom of the fishway. To study the effects of macro-roughnesses on the flow inside a VSF, three-dimensional unsteady simulations were carried out using the volume of fluid method to model the free surface. In this paper, kinematic quantities obtained by CFD are used to detail the flow inside a VSF with and without macro-roughnesses. It can provide valuable information about the flow characteristics, especially in areas where the experimental measurements are difficult to implement.

**Keywords:** ecohydraulics; vertical slot fishway; numerical simulation; URANS; LES; macro-roughnesses

## 1. Introduction

To restore ecological continuity and allow fish species to cross structures that hinder their migration (dams, weirs), crossing devices are built on rivers and streams. Among these devices, the vertical slot fishway (VSF) is the most common kind of fishway in France. It is almost insensitive to changes in upstream and downstream water levels, allowing it to be functional for a wide range of flows. Originally designed for highly migratory species such as salmons, sea trouts or lampreys, this type of device must be adapted to small species with limited swimming capacities in order to restore biological connectivity. To this end, numerous studies have been conducted in recent years to characterize the flow in VSFs according to their main geometric features and to better understand the effect of mean and turbulent flow variables on fish behavior [1–5]. To modify the characteristics of the flow and adapt it to the swimming capacities of small species, the addition of obstacles in the flow has been proposed and studied [6–8]. An increasingly common solution to facilitate the passage of small or benthic species is to insert macro-roughnesses (blocks of cylindrical shape or rocks) at the bottom of the fishway [9–11]. In this type of configuration, experimental measurements are difficult to obtain because of the complex geometry. Numerical simulations are therefore a valuable alternative. Numerous authors [12] have studied the flows in VSFs through numerical simulation using RANS steady-state modeling, in 2D [13], 2D depth-averaged [11] [14–20] and in 3D [21–32]. However, the flow in VSFs is highly unsteady [5,33] and the presence of macro-roughness makes it strongly

<figure><img src="image_2.png" /></figure>

### check for updates

Academic Editor: Chin H Wu

Received: 28 February 2025

Revised: 28 March 2025

Accepted: 3 April 2025

Published: 5 April 2025

**Citation:** Pineau, G.; Ballu, A.; David, L.; Calluaud, D. Numerical Simulation of the Unsteady 3D Flow in Vertical Slot Fishway—The Impact of Macro-Roughness. *Water* 2025, 17, 1088. https://doi.org/10.3390/w17071088

**Copyright:** © 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Pprime Institute, CNRS-Université de Poitiers-Isae Ensma, 86000 Poitiers, France; damiencalluaud@univ-poitiers.fr (D.C.); laurent.david@univ-poitiers.fr (L.D.)</footnote>

<footnote n="2">² Institute of Earth Surface Dynamics, University of Lausanne, 1015 Lausanne, Switzerland; aurelien.ballu@unil.ch</footnote>

<footnote>* Correspondence: gerard.pineau@univ-poitiers.fr</footnote>
three-dimensional [10]. Consequently, steady 2D or 2D depth-averaged models are there-
fore insufficient. Moreover, it is essential to account for non-uniform and unsteady flows
when adapting fishway hydrodynamics to the needs of target species. Therefore, the use of
3D unsteady models such as URANS, LES [24,34–36] or DES then becomes necessary to
correctly model these flows and, in particular, to access the flow fluctuations in statistical
turbulent quantities such as Reynolds stresses or turbulent kinetic energy.

This paper aims to complement the work of Fuentes-Pérez [31] on unsteady numerical flow modeling. For unsteady and 3D flows inside a VSF with macro-roughnesses, URANS and LES models are compared with experimental data to assess their ability to accurately predict the complex flow pattern. In addition, this numerical study highlights the influence of macro-roughness uniformly distributed on the bottom of a vertical slot fishway on the characteristics of the 3D flow.

**2. Materials and Methods**

2.1. *VSF Geometry*

This paper investigates a VSF design with a pool geometry similar to the design
studied by Calluaud et al. [7] (Figure 1). The length of the pools is L = 0.75 m, and the
width of the vertical slots is b = 0.075 m, giving a ratio of L/b equal to 10. The width of the
pool is B = 0.675 m, i.e., B/b = 9.

<figure><img src="image_3.png" /><figcaption>Figure 1. Geometry of the pool of the VSF.</figcaption></figure>

The slope of the VSF was set to s = 7.5%. The average flow depth measured at the center of the pools is equal to 0.30 m. The discharge is fixed to Q = 0.023 m³/s. The macro-roughnesses arranged on the bottom of the VSF model are equally spaced cylindrical studs with a diameter of 0.035 m and a height h<sub>r</sub> = 0.05 m.

The density d<sub>r</sub> is defined as the ratio of the elementary surface covered by three ele-
ments of macro-roughnesses (S<sub>r</sub>) to the total surface of the patch of cylinders (S<sub>t</sub>) (Figure 2).
In the studied configuration, the density is set to d<sub>r</sub> = S<sub>r</sub>/S<sub>t</sub> = 15%.

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Density of macro-roughnesses.</figcaption></figure>
## 2.2. *Numerical Simulations*

The VSF was simulated using Star CCM+ for two methods: URANS and LES. The URANS method (Unsteady Reynolds-Averaged Navier-Stokes) is based on the same principle as Reynolds decomposition, which divides the instantaneous quantities into a sum of a mean and a fluctuating value. Since the averaging process leads to a loss of information, closure models are required to complete the system of equations. Here, the low Reynolds number $k-\epsilon$ model is used. The URANS method takes the unsteady nature of the flow into account, while the RANS method does not.

The LES method (Large Eddy Simulation), however, takes a different approach. It uses a spatial filter to resolve large eddies (resulting from the domain geometry) while modeling only the smaller ones. The Smagorinsky model was used in this simulation. In both numerical simulations, the volume of fluid (VOF) method was implemented to emulate the free surface, which assigns volume fractions of both air and water to each computational cell. For the unsteady simulations, an implicit temporal discretization scheme is used for both methods. A system of two nested loops is employed. The first corresponds to the physical simulation time and captures the unsteady behavior of the flow. The second is nested within the first and aims to reach a provisional steady state. To exit the inner loop and proceed to the next time step, velocity convergence criteria at several monitoring points are applied.

Boundary conditions are schematized in Figure 3. The water levels are imposed at the inlet and outlet. Hydrostatic pressure and the volume fractions of each phase (water and air) are specified at the inlet and outlet of the computational domain. For all solid walls, no-slip wall boundary conditions are applied. The computational domain is extended above the fishway to move the upper boundary away from the region of the free surface. A symmetry condition is imposed on the boundaries of this domain.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Boundary conditions used for LES and URANS simulations.</figcaption></figure>

The size of the cells ($T^*$) in the various parts of the mesh was defined based on the width of the slot ($T^*/b$). A mesh sensitivity analysis was performed for three different mesh sizes, from the coarsest to the finest: $T^*/b = 1/2, 1/4$, and $1/8$ [10]. To estimate the errors due to spatial discretization, the Grid Convergence Index (GCI) method was used [37]. For the URANS method, the $T^*/b$ ratio was set to $1/4$ (Figure 4), which provides both a good geometric definition and a reduced spatial discretization error. The mean position of the free surface in each of the basins was estimated from the experimentally determined water level. The mesh was refined ($T^*/b = 1/8$) in the region around this surface position (+/- 20%). The part of the domain containing only air (above the free surface) was meshed with $T^*/b = 1$. To simulate the VSF flow near macro-roughnesses, a refinement ($T^*/b = 1/8$) was applied in a region bounded by the height of the cylinders. For the LES simulation, the
mesh size used for the URANS simulations was retained for all pools except the third one, which served as the reference pool. In this pool, the core mesh size is T*/b = 1/8 (Figure 4).

<figure><img src="image_6.png" /><figcaption>Figure 4. Mesh generated for simulation of the flow in a VSF with macro-roughnesses.</figcaption></figure>

The k-ε turbulence model used for URANS simulations is valid down to the wall (low Reynolds turbulence model). Thus, the first inner node of the boundary layer mesh should be located within the viscous sub-layer, at $y^+ = 1$, where $y^+$ is the non-dimensional wall distance [38], defined as follows:

$$y^+ = \frac{y_w \cdot u_\tau}{\nu} \qquad (1)$$

where $y^+$ is the non-dimensional wall distance, $y_w$ is the normal distance from the wall (m), $u_\tau$ is the wall friction velocity (m/s) and $\nu$ is the kinematic viscosity (m²/s).

In URANS simulation, $u_\tau$ is estimated from the friction coefficient given by the Schultz-Grunow expression:

$$C_f = \frac{0.37}{\log\left(\frac{U_\infty x}{\nu}\right)^{2.584}} \qquad (2)$$

where $C_f$ is the friction coefficient (-), $U_\infty$ is the velocity of the flow (m/s), $x$ is the distance in the flow direction (m) and $\nu$ is the kinematic viscosity (m²/s).

The same condition must be met in LES to correctly resolve the near-wall region. However, the wall shear velocity was approximated using the URANS results.

The boundary layer mesh thickness $\delta(x)$ was estimated to 0.03 m, based on the equation describing the evolution of the turbulent boundary layer thickness over a flat plate [39]:

$$\frac{\delta(x)}{x} \approx \frac{0.37}{R_e(x)^{1/5}} \qquad (3)$$

where $\delta(x)$ is the boundary layer thickness (m), $x$ is the distance in the flow direction (m) and $R_e(x)$ is the Reynolds number (-).

In contrast to URANS simulations, the anisotropy of the near-wall mesh must be minimized in LES. To resolve the inner-layer vortices, the streamwise and spanwise grid sizes in wall unit $\Delta x^+ \approx 100$ and $\Delta z^+ \approx 20$, respectively, were used [40].

In URANS, a starting water level is imposed at the output of the simulation domain (Figure 5). This results in a residual water level in the VSF. This initial state corresponds to a configuration in which the fishway entrance is closed and the exit remains open. In LES, the initial conditions are derived from the URANS calculations. These include pressure, velocity, and fluid volume fraction.
<figure><img src="image_7.png" /><figcaption>Figure 5. Initial conditions of the volume fraction of water (in blue) for URANS and LES methods.</figcaption></figure>

To meet precision requirements, the computational time step was set to $\Delta t = 10$ ms, based on the cell size of the computational domain and the flow kinematics. In the case of the URANS simulation, equilibrium of the water levels in each pool is reached after 20 s. The calculation of statistical turbulent quantities is performed only from this equilibrium state. The overall simulation time was set to 300 s corresponding to the acquisition time of the experimental measurements. As the flow is already fully developed at the beginning of the LES simulations, the total simulation time is based on the convergence of the mean and fluctuating component of velocity in the center of each pool.

## 2.3. Validation

Experiments were carried out to provide a basis of comparison to the numerical study. They were undertaken in a laboratory VSF, located at the Prime Institute of the University of Poitiers, France. The physical and numerical models of the VSF share identical characteristics (5 pools, $B/b = 9$, $Q = 0.023 \text{ m}^3/\text{s}$ and $s = 7.5\%$). A detailed description of the experimental apparatus is provided in [33]. In the configuration with macro-roughnesses, flow velocity measurements were performed using an Acoustic Doppler Velocimeter (ADV). The sampling rate of the ADV was 50 Hz over a 300 s acquisition time. A phase/space filter was applied [41]. A transverse profile (T) was defined in a plane parallel to the floor at a distance of $Z/b = 2$ and $X/b = 4.67$. For this profile, the first and last measurement points were located 50 mm from the walls. A vertical profile (V) was also measured along the Z-axis at $Y/b = 2.95$ from $X/b = 4.5$ at the first point, located 10 mm above the macro-roughnesses, to $X/b = 4.8$ at the free surface. The longitudinal profile $L$ ($Y/b = 2.95$) was defined along the X-axis at $Y/b = 2.95$. The sampling interval for all profiles is 15 mm. The location of the profiles is shown in Figure 6.

A comparative analysis of the flow characteristics was conducted to determine which simulation method—URANS or LES—is more appropriate to model the flow. The velocity magnitude and the turbulent kinetic energy were normalized using the maximum velocity defined by

$$V_d = \sqrt{2 \cdot g \cdot s \cdot L} \qquad (4)$$

where $V_d$ is the maximum velocity (m/s), $g$ is the acceleration due to gravity (m/s²), $s$ is the slope of the VSF and $L$ is the length of the pools (m).
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Velocity and turbulent kinetic energy profiles obtained using URANS, LES and ADV experimental measurements in the presence of macro-roughnesses. Profile T: transverse profile defined in a plane parallel to the floor (Z/b = 2 and X/b = 4.67); profile V: vertical profile measured along the Z-axis (Y/b = 2.95 and X/b = 4.5); profile L: longitudinal profile defined in a plane parallel to the floor (Z/b = 2 and Y/b = 2.95).</figcaption></figure>

The 3D non-dimensional magnitude of the flow velocity and turbulent kinetic energy profiles obtained from experimental measurements and from both URANS and LES simulations are plotted in Figure 6. As mentioned above, the estimated errors due to spatial discretization in the URANS simulations were assessed using the Grid Convergence Index (GCI) method [37]. For the experimental measurements, uncertainties were estimated on the specifications provided by the ADV manufacturer. The method used to assess these uncertainties follows the Guide to the Expression of Uncertainty in Measurement [42]. The total relative uncertainty associated with each mean velocity component was calculated using the general uncertainty propagation equation, with a coverage factor of k = 2 (corresponding to a 95% confidence interval). Based on the uncertainties in both the mean and
instantaneous velocity, the uncertainties in the mean velocity magnitudes (4%) and in the turbulent kinetic energy (12%) were determined [43]. These uncertainties are represented as error bars in Figure 6.

To compare the experimental and numerical results, the correlation coefficient $r^2$ and the relative standard deviation $\sigma_m$ between the curves were calculated. The correlation coefficient was used to quantify the degree of agreement between the experimental and numerical data. The values of $\sigma_m$ and $r^2$ are reported in Table 1.

<table><caption>**Table 1.** $\sigma_m$ and $r^2$ values calculated between the velocity and turbulent kinetic energy profiles obtained from numerical simulations and those derived from the experimental measurements.</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="2">Profile T</th><th colspan="2">Profile V</th><th colspan="2">Profile L</th></tr><tr><th>σ<sub>m</sub> (%)</th><th>r<sup>2</sup> (-)</th><th>σ<sub>m</sub> (%)</th><th>r<sup>2</sup> (-)</th><th>σ<sub>m</sub> (%)</th><th>r<sup>2</sup> (-)</th></tr></thead><tbody><tr><td rowspan="2">||V||<sub>3D</sub>/V<sub>d</sub></td><td>URANS</td><td>28.8</td><td>0.87</td><td>15.7</td><td>0.56</td><td>22.6</td><td>0.92</td></tr><tr><td>LES</td><td>28.5</td><td>0.89</td><td>14.2</td><td>0.91</td><td>11.3</td><td>0.95</td></tr><tr><td rowspan="2">k<sub>3D</sub>/V<sub>d</sub><sup>2</sup></td><td>URANS</td><td>66.0</td><td>0.95</td><td>55.4</td><td>0.55</td><td>65.0</td><td>0.26</td></tr><tr><td>LES</td><td>28.8</td><td>0.98</td><td>27.0</td><td>0.74</td><td>38.0</td><td>0.64</td></tr></tbody></table>

The mean velocity curves $||V||_{3D}/V_d$ of the transverse profiles T obtained from both URANS and LES simulations are similar to those obtained experimentally. The jet flow induces a sharp increase in both velocity and turbulent kinetic energy at $Y/b = 3$. On either side of the jet, two lower-intensity zones with local velocity minima reveal the presence of two recirculation zones indicating the existence of the first flow pattern [1]. Both URANS and LES simulations clearly capture this flow topology.

This visual observation is confirmed by high correlation coefficients, for both the URANS and LES simulations (Table 1). The errors on the mean velocity are nearly equivalent in URANS and LES and remain relatively low for all profiles T, V and L ($\sigma_m = 18.7\%$ for the LES and $\sigma_m = 19.3\%$ for the URANS). The turbulent kinetic energy $k_{3D}/V_d^2$ is underestimated by both methods. However, LES models turbulent kinetic energy more accurately than URANS for all profiles with $\sigma_m = 27\%$ compared to $\sigma_m = 58.5\%$ for URANS. It is worth noting that the correlation coefficients between the LES simulation and the experimental profiles are higher for profiles T, V and L than those obtained using the URANS method. This indicates that LES provides a more accurate distribution of turbulent kinetic energy in the flow within the VSF.

Regarding computation time, it is intrinsically dependent on the methods and the respective mesh resolution. For the VSF configuration studied, 320 s of unsteady flow was simulated. For the URANS method, the computational time per iteration is 1.9 s on a machine with 128 cores. For LES, this time increases to 5.2 s per iteration. This corresponds to a total computation time of 17 h for the URANS method, compared with 46 h for LES. Turbulent kinetic energy plays a key role in fish behavior within the basins and in their ability to pass through them. The LES simulation provides more accurate results for TKE than URANS, justifying its use despite the higher computational cost, which remains reasonable (under two days).

## 3. Results

During the validation of the numerical simulations, the velocity profiles discussed in the previous section highlighted the three-dimensional nature of the flow. This is also clearly illustrated in Figures 7 and 8, which show the influence of macro-roughnesses on the three-dimensionality of the flow. These figures display the mean flow streamlines in a micro-pool basin with and without macro-roughnesses, simulated using LES. In the smooth configuration, the two-dimensional nature of the flow is confirmed, with vortex structures
whose centers remain aligned along a vertical axis. In configuration with $dr = 15\%$, the lower vortex also retains a vertical axis. However, the axis of the upper vortex becomes steeply inclined, turning horizontal near the macro-roughnesses. This development highlights a significant modification of the mean flow, which becomes three-dimensional well beyond the immediate influence of the canopy. The jet remains, on average, two-dimensional. It feeds the upper recirculation zone by plunging towards the central deflector, primarily due to the strong velocity gradient induced by the presence of macro-roughness.

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Evolution of the position of the vortex centers (in red) of the main recirculation zones in configuration without macro-roughnesses. Streamlines from the mean velocity field are shown in green to illustrate the spatial evolution of the upper recirculation zone, in red for the jet, in dark blue for the lower recirculation zone and in light blue for the small recirculation zone located behind the small baffle (LES).</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 8.** Evolution of the position of the vortex centers (in red) of the main recirculation zones in configuration with macro-roughnesses. Streamlines from the mean velocity field are shown in green to illustrate the spatial evolution of the upper recirculation zone, in red for the jet, in dark blue for the lower recirculation zone and in light blue for the small recirculation zone located behind the small baffle (LES).</figcaption></figure>

The LES model was used to study the flow generated both through and over the macro-roughnesses. The flow is primarily described using two normalized heights (Z/hr) and (Z/b) relative to the roughness height (hr) and the slot width (b) in order to facilitate comparisons with results from the literature.

Regarding the macro-roughness configuration, the flow topology on the plane furthest from the floor ($Z/hr = 4.5$) is consistent from one basin to the next. The jet follows a curved
trajectory, generating two counter-rotating recirculation zones on either side. Compared to the smooth floor configuration, the upper vortex is deformed and its center is shifted downstream. At Z/hr = 1.2, the flow remains nearly identical in each basin, although it appears more disturbed, particularly in the upper vortex, which exhibits a less well-defined structure. This flow pattern, referred to as type 1, was clearly identified by Wang [5].

The Z/hr = 0.75 plane provides a useful view of the average flow within the canopy's macro-roughness elements. When encountering the first macro-roughness, the jet splits into two asymmetrical parts. At the bottom of the various basins, most of the flow is directed to the right of the first macro-roughness, flowing directly to the second element located along the same alignment. It then divides into two fairly symmetrical streams and continues toward the downstream slot, following the path formed by the aligned macro-roughnesses. The second part of the flow is deflected to the left of the first macro-roughness towards the opposite wall of the slot. This portion then encounters other macro-roughnesses aligned along a diagonal centered on the first cylinder and oriented at an angle of 60° to the longitudinal axis. As a result, the jet is distributed over a large portion of the available canopy surface, generating numerous wake zones and promoting the dissipation of turbulent kinetic energy. In Figure 9, the iso-contours of the mean velocity field in the vertical planes passing through the center of the slots show that the jet experiences greater velocity variations across the water column in this configuration than in the configuration without macro-roughnesses.

<figure><img src="image_11.png" /><figcaption>**Figure 9.** Mean velocity fields. **Left:** vertical cross-section located at 45° in middle of slots 2 and 3. **Right:** horizontal sectional planes at Z/b = 0.5 (Z/hr = 0.75), Z/b = 0.8 (Z/hr = 1.2) and Z/b = 3 (Z/hr = 4.5).</figcaption></figure>

Figure 10 shows the spatial distributions of the mean velocity magnitude and the turbulent kinetic energy on a plane close to the bottom with and without macro-roughnesses. It is clearly observed that the presence of roughness reduces the regions of high intensity, both in terms of velocity and turbulent kinetic energy. The roughness element near the
bottom diffuses the jet, causing it to widen and fragment across the entire basin surface. The vertical extent of the flow influenced by the roughness is also significant, as illustrated in Figure 11. These figures demonstrate that the peaks in velocity or TKE magnitude are dampened by the presence of roughness over a substantial layer near the bottom of the basins. As a result, in this near-bottom region, the local flow characteristics become more favorable for the movement and development of fish with more limited swimming capabilities, particularly benthic species.

<figure><img src="image_12.png" /><figcaption>**Figure 10.** LES results: (a,b) spatial distribution of the amplitude of the mean velocity field and (c,d) spatial distribution of the turbulent kinetic energy in the Z/hr = 0.5 plane (mid-height of the macro-roughness) for configurations without macro-roughnesses (**left**) and with macro-roughnesses (**right**). The dashed line indicates the plane shown in Figure 11.</figcaption></figure>

This three-dimensionality of the flow also influences the flow dynamics as described in [7]. Unsteady phenomena, such as jet oscillation (or “jet beating”), do not necessarily exhibit the same frequency and amplitude throughout the entire water column. To verify this (Figure 12), instantaneous velocity fields of the three phases of the flow were plotted on a plane at Z/hr = 3 and simultaneously on a plane close to the roughness element at Z/hr = 1.2. Between these two planes (Z/hr = 1.2 and Z/hr = 3), the unsteady flow structures differ significantly. In particular, during phase (b), the velocity burst feeding the upper recirculation zone is not observed in the plane closest to the bottom. This observation aligns with the previous analysis of the mean flow field, which showed that the upper vortex was strongly inclined, suggesting it is primarily fed “from above”.

In the configuration with dr = 15%, the flow is generally less intense in the upper recirculation zone, regardless of the phase observed. This reduction in kinetic energy in the main flow is a consequence of the increased dissipation near the bottom within the rough
underlayer. The high rate of turbulent kinetic energy dissipation within the canopy and in the rough underlayer is illustrated in Figure 13.

<figure><img src="image_13.png" /><figcaption>**Figure 11.** LES results: (a,b) spatial distribution of the amplitude of the mean velocity field and (c,d) spatial distribution of the turbulence kinetic energy in the vertical plane for configurations without macro-roughnesses (**left**) and with macro-roughnesses (**right**). The dotted line indicates the plane in Figure 10.</figcaption></figure>

Beyond the zones of high turbulent kinetic energy (TKE) dissipation located near the walls, it is essential to evaluate the capacity of the numerical simulation to capture the full range of turbulent scales. To this end, the temporal evolution of the velocity magnitude was compared at a representative point (X/b = 2.4 and Y/b = 3.2, Z/hr = 3) where TKE levels are high (Figure 14). A significant difference is observed between the results depending on the simulation methods. The amplitude of velocity fluctuations is considerably reduced in the URANS model compared to the LES model. A Fourier analysis of the velocity data highlights the intrinsic characteristics of each method. The URANS approach explicitly resolves only the largest scales of motion and models the smaller-scale fluctuations representative of turbulent agitation. It thus behaves as a spatial or temporal filter. The cut-off frequency of this filter is defined by the integral scale (length or time) of the turbulence model, which, for the chosen mesh and time step, corresponds to 1 Hz. In contrast to the RANS method—which models all scales of turbulence—the LES (Large Eddy Simulation) method resolves the large-scale structures, which are strongly influenced by the geometry, and only models the smaller more universal scales. A spatial filter is implicitly defined by the mesh size, determining which scales are resolved and which are modeled. In our LES simulation, the FFT spectrum of the velocity magnitude shows
significant amplitudes at high frequencies, indicating that the mesh accurately resolves the energy cascade down to small-scale eddies, as discussed by Quaresma et al. (2025) [36].

<figure><img src="image_14.png" /><figcaption>**Figure 12.** Instantaneous velocity fields from LES simulations in $dr = 15\%$ configuration illustrating the three different phases (a–c) of the flow. Left hand column: $Z/h_r = 1.2$; right hand column: $Z/h_r = 3$.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Figure 13.** Isocontours of the turbulent kinetic energy dissipation rate.</figcaption></figure>
<figure><img src="image_16.png" /><figcaption>**Figure 14.** (a) One point temporal evolution of the velocity magnitude; (b) FFT from URANS (black) and LES (gray).</figcaption></figure>

**4. Conclusions**

The objective of this study was to validate unsteady 3D flow simulations inside a vertical slot fishway with and without roughness elements placed on the bottom of the pools. LES and URANS simulations were performed in the presence of a free surface and compared to ADV measurements conducted in the same geometry. Overall, the LES results showed better agreement with the experimental data than those obtained with URANS. Therefore, the analysis of flow behavior—with or without macro-roughnesses—should primarily rely on LES results. Moreover, LES provides a more accurate description of the unsteady nature of the flow. While the average velocity tends to be slightly overestimated and TKE is underestimated by approximately 30%, the overall flow topology and the spatial distribution of the kinematic quantities closely match experimental observations.

The hydraulic analysis revealed that adding macro-roughnesses inside the VSF alters the flow structure. To quantify this influence, velocity and turbulent kinetic energy profiles were compared between the smooth configuration and the roughened configuration with $dr = 15\%$ using LES simulations. For the studied basin width ($B/b = 9$), the presence of macro-roughnesses enhances the three-dimensionality of the mean flow by tilting the axis vortex of the upper recirculation zone. The analysis of unsteady flow demonstrated that the three characteristic flow phases observed in a conventional VSF are still present when macro-roughnesses are introduced. However, in the roughened sub-layer, the dynamics and intensity of these unsteady phases are reduced due to increased energy dissipation. Since the overall dynamics in the water column are not negatively affected, macro-roughnesses do not impair the hydraulic efficiency of the VSF.

On the contrary, the flow dynamics near the bottom are significantly modified by the roughness elements, suggesting an increase in passage efficiency for benthic species. The size and spacing of the macro-roughnesses will need to be optimized based on the target fish species present in the river system. To confirm the conclusions of this study, field trials on existing vertical slot fishways will be conducted by the Office Français de la Biodiversité (OFB) over the next two years.

**Author Contributions:** Methodology, A.B., L.D., D.C. and G.P.; validation, A.B., L.D., D.C. and G.P.; investigation, A.B., L.D., D.C. and G.P.; writing—original draft, A.B. and G.P.; writing—review and editing, A.B., L.D., D.C. and G.P. All authors have read and agreed to the published version of the manuscript.

**Funding:** This project received funding from the "Office Français de la Biodiversité" (OFB), "Voies Navigables de France" (VNF) and the European Union's Horizon 2020 research and innovation program FITHydro (www.fithydro.eu), under grant agreement [No 727830].
**Data Availability Statement:** The data presented in this study are available on request from the
corresponding author.

**Acknowledgments:** The authors would like to acknowledge the Nouvelle Aquitaine CPER and the Environmental Hydrodynamic Platform of the University of Poitiers and CNRS for their facilities and equipment.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## Abbreviations

The following abbreviations are used in this manuscript:

<table><tr><td>VSF</td><td>vertical slot fishway</td></tr><tr><td>CFD</td><td>Computational Fluid Dynamics</td></tr><tr><td>URANS</td><td>Unsteady Reynolds-Averaged Navier-Stokes</td></tr><tr><td>LES</td><td>Large Eddy Simulation</td></tr><tr><td>VOF</td><td>volume of fluid</td></tr><tr><td>ADV</td><td>Acoustic Doppler Velocimeter</td></tr></table>

## References

1. Wu, S.; Rajaratnam, N.; Katopodis, C. Structure of Flow in Vertical Slot Fishway. *J. Hydraul. Eng.* **1999**, *125*, 351–360. [CrossRef]

2. Puertas, J.; Pena, L.; Teijeiro, T. Experimental Approach to the Hydraulics of Vertical Slot Fishways. *J. Hydraul. Eng.* **2004**, *130*, 10–23. [CrossRef]

3. Liu, M.; Rajaratnam, N.; Zhu, D.Z. Mean Flow and Turbulence Structure in Vertical Slot Fishways. *J. Hydraul. Eng.* **2006**, *132*, 765–777. [CrossRef]

4. Tarrade, L.; Texier, A.; David, L.; Larinier, M. Topologies and measurements of turbulent flow in vertical slot fishways. *Hydrobiologia* **2008**, *609*, 177–188. [CrossRef]

5. Wang, R.; David, L.; Larinier, M. Contribution of experimental fluid mechanics to the design of vertical slot fish passes. *Knowl. Manag. Aquat. Ecosyst.* **2010**, *396*, 2. [CrossRef]

6. Mallen-Cooper, M.; Zampatti, B.; Stuart, I.; Baumgartner, L. *Innovative Fishways—Manipulating Turbulence in the Vertical Slot Design to Improve Performance and Reduce Cost*; Technicol Report, Fishway Consult; Murray Darling Basin Authority: Sydney, NSW, Australia, 2008.

7. Callaud, D.; Pineau, G.; Texier, A.; David, L. Modification of vertical slot fishway flow with a supplementary cylinder. *J. Hydraul. Res.* **2014**, *52*, 614–629. [CrossRef]

8. Zhao, H.; Xu, Y.; Lu, Y.; Lu, S.; Dai, J.; Meng, D. Numerical Study of Vertical Slot Fishway Flow with Supplementary Cylinders. *Water* **2022**, *14*, 1772. [CrossRef]

9. Ballu, A.; Pineau, G.; Callaud, D.; David, L. Experimental-Based Methodology to Improve the Design of Vertical Slot Fishways. *J. Hydraul. Eng.* **2019**, *145*, 04019031.

10. Ballu, A. Étude Numérique et Expérimentale de L'écoulement Turbulent au Sein des Passes à Poissons à Fentes Verticales. Analyse de L'écoulement Tridimensionnel et Instationnaire. Ph.D. Thesis, Université de Poitiers, Poitiers, France, 2017.

11. Bombač, M.; Rak, G.; Novak, G. Mathematical modeling of flow in vertical slot fishways. *Acta Hydrotech.* **2020**, *33*, 97–112. [CrossRef]

12. Tarena, F.; Nyqvist, D.; Katopodis, C.; Comoglio, C. Computational fluid dynamics in fishway research—A systematic review on upstream fish passage solutions. *J. Ecolhydraulics* **2024**, *10*, 107–126. [CrossRef]

13. Tarrade, L. Etude des Écoulements Turbulents dans les Passes à Poissons à Fentes Verticales: Adaptation aux Petites Espèces. Ph.D. Thesis, Université de Poitiers, Poitiers, France, 2007.

14. Cea, L.; Pena, L.; Puertas, J.; Vázquez-Cendón, M.E.; Peña, E. Application of several depth-averaged turbulence models to simulate flow in vertical slot fishways. *J. Hydraul. Eng.* **2007**, *133*, 160–172. [CrossRef]

15. Chorda, J.; Maubourguet, M.M.; Roux, H.; Larinier, M.; Tarrade, L.; David, L. Two-dimensional free surface flow numerical model for vertical slot fishways. *J. Hydraul. Res.* **2010**, *48*, 141–151.

16. Bermúdez, M.; Puertas, J.; Cea, L.; Pena, L.; Balairón, L. Influence of pool geometry on the biological efficiency of vertical slot fishways. *Ecol. Eng.* **2010**, *36*, 1355–1364. [CrossRef]

17. Bombač, M.; Novak, G.; Rodič, P.; Četina, M. Numerical and physical model study of a vertical slot fishway. *J. Hydrol. Hydromechanics* **2014**, *62*, 150–159. [CrossRef]
18. Bombač, M.; Četina, M.; Novak, G. Study on flow characteristics in vertical slot fishways regarding slot layout optimization. *Ecol. Eng.* **2017**, *107*, 126–136. [CrossRef]

19. Bousmar, D.; Zorzan, G.; Gillet, A.; Delhoulle, J.; Baugnée, A.; De Greef, J.C. Refurbishing an old fish pass: Physical and numerical modeling. In *Proceedings of the River Flow, Lausanne, Switzerland, 3–5 September 2014*; pp. 2427–2435.

20. Zhang, D.; Qu, Y.; Shi, X.; Liu, Y.; Jiang, C. Design of a novel multislot and pool-weir combined fishway based on hydraulic properties analysis and fish-passage experiments. *J. Hydraul. Eng.* **2024**, *150*, 04024004. [CrossRef]

21. Heimerl, S.; Hagmeyer, M.; Echteler, C. Numerical flow simulation of pool-type fishways: New ways with well-known tools. *Hydrobiologia* **2008**, *609*, 189–196. [CrossRef]

22. Chang, W.; Lien, H.; Tsai, W.; Lai, J.; Guo, W.; Hu, T. 3D Flow Simulations In A Pool-weir-type Fishway In The Tonghou River. In *Proceedings of the 33rd IAHR World Congress, Vancouver, BC, Canada, 9–14 August 2009*.

23. Barton, A.F.; Keller, R.J.; Katopodis, C. Verification of a numerical model for the prediction of low slope vertical slot fishway hydraulics. *Australas. J. Water Resour.* **2009**, *13*, 53–60. [CrossRef]

24. Musall, M.; Schmitz, C.; Oberle, P.; Nestmann, F.; Henning, M.; Weichert, R. Analysis of flow patterns in vertical slot fishways. In *Proceedings of the River Flow, Lausanne, Switzerland, 3–5 September 2014*; pp. 2421–2426.

25. An, R.; Li, J.; Liang, R.; Tuo, Y. Three-dimensional simulation and experimental study for optimising a vertical slot fishway. *J. Hydro-Environ. Res.* **2016**, *12*, 119–129. [CrossRef]

26. Duguay, J.; Lacey, R.; Gaucher, J. A case study of a pool and weir fishway modeled with OpenFOAM and FLOW-3D. *Ecol. Eng.* **2017**, *103*, 31–42. [CrossRef]

27. Umeda, Y.C.L.; de Lima, G.; Janzen, J.G.; Salla, M.R. One- and three-dimensional modelling of a vertical-slot fishway. *J. Urban Environ. Eng.* **2017**, *11*, 99–107.

28. Stamou, A.I.; Mitsopoulos, G.; Rutschmann, P.; Bui, M.D. Verification of a 3D CFD model for vertical slot fish-passes. *Environ. Fluid Mech.* **2018**, *18*, 1435–1461. [CrossRef]

29. Li, Y.; Wang, X.; Xuan, G.; Liang, D. Effect of parameters of pool geometry on flow characteristics in low slope vertical slot fishways. *J. Hydraul. Res.* **2019**, *58*, 395–407. [CrossRef]

30. Quaranta, E.; Katopodis, C.; Comoglio, C. Effects of bed slope on the flow field of vertical slot fishways. *River Res. Appl.* **2019**, *35*, 656–668. [CrossRef]

31. Fuentes-Pérez, J.F.; Silva, A.T.; Tuhtan, J.A.; García-Vega, A.; Carbonell-Baeza, R.; Musall, M.; Kruusmaa, M. 3D modelling of non-uniform and turbulent flow in vertical slot fishways. *Environ. Model. Softw.* **2018**, *99*, 156–169. [CrossRef]

32. Ahmadi, M.; Kuriqi, A.; Nezhad, H.M.; Ghaderi, A.; Mohammadi, M. Innovative configuration of vertical slot fishway to enhance fish swimming conditions. *J. Hydrodyn.* **2022**, *34*, 917–933. [CrossRef]

33. Tarrade, L.; Pineau, G.; Callaud, D.; Texier, A.; David, L.; Larinier, M. Detailed experimental study of hydrodynamic turbulent flows generated in vertical slot fishways. *Environ. Fluid Mech.* **2011**, *11*, 1–21. [CrossRef]

34. Fuentes-Pérez, J.F.; Quaresma, A.L.; Pinheiro, A.; Sanz-Ronda, F.J. OpenFOAM vs FLOW-3D: A comparative study of vertical slot fishway modelling. *Ecol. Eng.* **2022**, *174*, 1064–1069. [CrossRef]

35. Quaresma, A.L.; Pinheiro, A.N. Modelling of Pool-Type Fishways Flows: Efficiency and Scale Effects Assessment. *Water* **2021**, *13*, 851. [CrossRef]

36. Quaresma, A.L.; Romão, F.; Pinheiro, A.N. A Comparative Assessment of Reynolds Averaged Navier-Stokes and Large-Eddy Simulation Models: Choosing the Best for Pool-Type Fishway Flow Simulations. *Water* **2025**, *17*, 686. [CrossRef]

37. Celik, I.B.; Ghia, U.; Roache, P.J.; Freitas, C.J. Procedure for estimation and reporting of uncertainty due to discretization in CFD applications. *J. Fluids Eng. Trans. ASME* **2008**, *130*, 078001.

38. Schlichting, H.; Gersten, K. *Boundary-Layer Theory*, 9th ed.; Springer: Berlin/Heidelberg, Germany, 2017.

39. Chassais, P. *Turbulence en Mécanique des Fluides: Analyse du Phénomène en vue de sa Modélisation à L'usage de L'ingénieur*, Cépaduès-Editions: Toulouse, France, 2000.

40. Piomelli, U.; Balaras, E. Wall-layer models for large-eddy simulations. *Annu. Rev. Fluid Mech.* **2002**, *34*, 349–374. [CrossRef]

41. Goring, D.G.; Nikora, V.I. Despiking Acoustic Doppler Velocimeter Data. *J. Hydraul. Eng.* **2002**, *128*, 117–126. [CrossRef]

42. Joint Committee for Guides in Metrology—*Guide to the Expression of Uncertainty in Measurement* (JCGM-2023). Available online: https://www.bipm.org/documents/20126/2071204/JCGM_GUM-1.pdf/74e7aa56-2403-7037-f975-cd6b555b80e6 (accessed on 27 February 2025).

43. Beaulieu, C.; Pineau, G.; Ballu, A.; David, L. Estimation process of uncertainties measurement in a research laboratory: Contribution and perspectives. Example of a research laboratory in hydrology for aquatic environments. In *Proceedings of the 17th International Congress of Metrology*, Paris, France, 21–24 September 2015.

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.