# Experimental flows through an array of emerged or slightly submerged square cylinders over a rough bed

Marina Oukacine, Sébastien Proust, Frédérique Larrarte, Nicole Goutal

► To cite this version:

Marina Oukacine, Sébastien Proust, Frédérique Larrarte, Nicole Goutal. Experimental flows through an array of emerged or slightly submerged square cylinders over a rough bed. Scientific Data, 2021, 8 (1), 10.1038/s41597-020-00791-w. hal-03130105

HAL Id: hal-03130105

https://hal.science/hal-03130105

Submitted on 3 Feb 2021

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
SCIENTIFIC DATA 110110
0111101
11011110
011101101

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

Check for updates

OPEN

DATA DESCRIPTOR

# Experimental flows through an array of emerged or slightly submerged square cylinders over a rough bed

Marina Oukacine ID<fnref n="1" />✉, Sébastien Proust ID², Frédérique Larrarte ID<fnref n="1" />,³ & Nicole Goutal<fnref n="1" />,⁴

The experimental dataset presented was collected in an 18 m long and 1 m wide laboratory flume. Low to high flood flows through an urbanized floodplain were modelled. The floodplain bed is rough, modelled with dense artificial grass. A square cylinder array, representing house models, was set on the rough bed. The cylinder immersion rate was varied: cylinders are emerged for three flow cases and slightly submerged for one case. The experimental dataset comprises water levels, measured using an ultrasonic transit time probe, velocities across the channel measured using an Acoustic Doppler Velocimetry with a side looking probe, and velocities in longitudinal-vertical planes measured using Particle Image Velocimetry. These data could help understanding the physical processes associated with high flood flows through urbanized floodplains, with a focus on the transition from emerged to submerged obstacles. They could also be used as benchmark data to assess the ability of numerical models from one to three-dimensions to estimate the flood hazard (water depth, velocity) over a wide range of flood event magnitudes.

## Background & Summary

As a result of climate change, extreme floods will become more frequent and more intense. People and properties, such as housing and industrial facilities, must therefore be protected against these floods for which data are very scarce or even non-existent<fnref n="1" />. When moving from low to extreme flooding, the flood extent over the floodplain greatly varies. This study aimed at investigating if there is a change in the physical processes when the houses initially emerged become weakly submerged and at providing well documented data for validation of simulations. To our knowledge, very few experimental studies have focused on this issue by considering a dense urban area, as depicted in Fig. 1a². Experimental study exists³⁻¹⁰ and used blocks up to 25. Most of them provide only water level measurements and are not interested in the study of the vertical confinement. Another complementary area of research concerns the development of simulation tools for extreme flows in a congested area. The field data are very scarce so the experimental data are essential to validate the different ways of modellisation. One-dimensional models are reasonably convenient to simulate flood propagation in straight streets, except near-street intersections where the flow is typically highly perturbed, and strongly 2-D or even 3-D¹¹. 2-D modelling was improved by verifying the accuracy of the roughness approach against full buildings incorporation in flood simulation¹² or by representing the urban area topography¹³ or by introducing porosity in the model¹⁴⁻¹⁷ in order to take into account the building blockage effect⁷,¹⁸. Recently, more information can be found in a review of the data available¹⁹.

As extreme flood data are very scarce or even non-existent in the field, the present dataset could help understanding the physical processes associated with extreme flood flows in a large simplified urban district area. Moreover, these data could be used as benchmark data to assess the ability of numerical models from one-dimension (1D) to three-dimensions (3D) to estimate the flood hazard over a wide range of flood event magnitudes²,²⁰.

<footnote n="1">¹Saint-Venant Laboratory for Hydraulics (LHSV), Chatou, France. ²Riverly Research Unit, INRAE, Villeurbanne, France. ³Univ. Gustave Eiffel, Marne la Vallée, France. ⁴EDF R&D, National Laboratory for Hydraulics and Environment (LNHE), Chatou, France. ✉e-mail: marina.oukacine@yahoo.fr</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** (a) Open-channel flume at INRAE Lyon-Villeurbanne (view looking upstream); (b) sketch of a cross-section (view looking upstream) with a zoom of the artificial grass; (c) sketched top view of the flume without the square cylinders: 1. Inlet tank, 2. Vertical linear ramp with a 15% slope rising the water until the synthetic grass bed level, the water level is nearly flat in the upstream part of the flume, 3. Working channel length, 4. Downstream tail weir; (d) sketch of the lateral view of the entire flume.</figcaption></figure>

<table><caption>**Table 1.** Flow parameters of the test cases: *H* is water depth, *k* is square cylinder height, *Q* is flow rate, *U*<sub>Q</sub> is bulk velocity, *H*<sub>weir</sub> is the height of the downstream weir, *λ*<sub>f</sub> is frontal density per unit horizontal surface, *F*<sub>RUQ</sub>(*H*) is the Froude number and *Re*<sub>UQ</sub>(*H*) is the Reynolds number.</caption><thead><tr><th>N°</th><th>H/k</th><th>H</th><th>Q</th><th>U<sub>Q</sub></th><th>H<sub>weir</sub></th><th>λ<sub>f</sub></th><th>F<sub>RUQ</sub>(H)</th><th>Re<sub>UQ</sub>(H)</th></tr><tr><th>[-]</th><th>[-]</th><th>[m]</th><th>[m³·s⁻¹]</th><th>[m·s⁻¹]</th><th>[m]</th><th>[-]</th><th>[-]</th><th>[-]</th></tr></thead><tbody><tr><td>1</td><td>0.42</td><td>0.025</td><td>0.0016</td><td>0.064</td><td>0.013</td><td>0.08</td><td>0.130</td><td>1590</td></tr><tr><td>2</td><td>0.93</td><td>0.055</td><td>0.0033</td><td>0.060</td><td>0.037</td><td>0.17</td><td>0.082</td><td>3280</td></tr><tr><td>3</td><td>0.98</td><td>0.058</td><td>0.0038</td><td>0.065</td><td>0.038</td><td>0.18</td><td>0.086</td><td>3770</td></tr><tr><td>4</td><td>1.48</td><td>0.088</td><td>0.0103</td><td>0.117</td><td>0.059</td><td>0.185</td><td>0.120</td><td>101000</td></tr></tbody></table>

The experiments were carried out in an 18 m long and 1 m wide open-channel flume (Fig. 1a), which is located in the Hydraulic and Hydro-morphology Laboratory (HHLab) at INRAE Lyon-Villeurbanne, France. The channel bed slope in the longitudinal direction, *S*<sub>0</sub>, is equal to 1.05 mm/m, and the working length is 17.25 m (Fig. 1c). The cross-section is rectangular (Fig. 1b) with the right-hand sidewall made of glass, and the left-hand sidewall (removable) made of plexiglass. The channel bed is covered with dense artificial grass (5 mm high rigid blades with a density of 256 blades per square centimeter), see Fig. 1a,b. Over the artificial grass, 833 cylinders were placed with an in-line distribution (Fig. 1a). Each cylinder has a square section (side length *ℓ* = 6.4 cm) and has a height *k* = 5.92 cm when measured from the average top of the grass blades. The distance between two adjacent square cylinders is constant, the same in the transverse and longitudinal directions, and equals to 7.9 cm (Fig. 1a). The distance between the sidewall and the closest square cylinder is equal to 3.9 cm. The most upstream transverse square cylinder row is located at longitudinal position *x* = 0.17 m (see the origin of *x* in Fig. 1a). The array comprises 119 transverse rows of square cylinders, with 7 square cylinders in each row. At the downstream end of the flume (*x* = 17.25 m), an adjustable vertical weir enables the water surface to be controlled.

This physical model is at scale 1/100. At real scale, the spacing between two adjacent houses, *L* − *ℓ*, (where *L* is the distance between two square cylinder centres) corresponds to a two-way motorway including sidewalks, and the height of the house, *k*, corresponds to that of a single-storey house (approximately 6 m high). The aligned configuration was chosen to simplified a suburban area. Note that, because of the presence of the square cylinders, the classical bed friction laws are not applicable to these highly confined flows with small immersion rates.

Four streamwise uniform flows on average were investigated, varying the immersion rate *H*/*k*, where *H* is water depth and *k* is square cylinder height, with *H*/*k* = 0.42; 0.93; 0.98; 1.48. The hydraulic parameters are reported in Table 1, in which *F*<sub>RUQ</sub>(*H*) = *U*<sub>Q</sub>/√g*H* is the Froude number, *Re*<sub>UQ</sub>(*H*) = *U*<sub>Q</sub>*H*/ν is the Reynolds number, *U*<sub>Q</sub> = *Q*/*A* is bulk velocity, *A* = *BH* being the wet area and *B* the channel width. Note that, the data with immersion rate *H*/*k* = 0.93 and 1.48 are more exhaustive than the other two flow cases because the velocity measurement is more detailed (see section Methods subsection ADV and PIV velocity measurements positions).
<figure><img src="image_4.png" /><figcaption>Fig. 2 Average height of the square cylinders in each transverse row measured from the top of the grass blades at various longitudinal positions. The solid line indicates the spatial average in the horizontal plane $<k(x, y)>_{i,j} = 0.0592$ m and the dotted lines represent $\pm 2 \sigma(x, y)$, where $\sigma(x, y)$ is the standard deviation of the local average height $k(x, y)$.</figcaption></figure>

## Methods

A right-handed Cartesian coordinate system was used (Fig. 1), in which the x-axis is aligned with the longitudinal direction, parallel to the flume bottom, the y-axis is aligned with the spanwise direction, and the vertical z-axis is aligned with the vertical direction (normal to channel bed). The origin $x=0$ is located at the beginning of working length (Fig. 1c), $y=0$ at the right-hand sidewall, and $z=0$ at the top of the grass blades (Fig. 1b).

**Determining the flow rate.** The values of immersion rate $H/k$ and water depth $H$ were chosen in collaboration with the other partners of the FlowRes project²¹ (https://flowres.inrae.fr/en/), with experiments in four laboratories. $H/k$-values are thus input data while flow rate ($Q$) values are data to be determined.

The inlet flow rate is controlled by a control valve Samson 3310 with servomotor PSQ, and is monitored with an electromagnetic flow meter (Krohne Waterux 3000 IFC 100). The discharge measurement uncertainty is 0.3% of the maximum range set either here $\pm 0.16$ L.s$^{-1}$, calculated with manufacturer data for a maximum flow rate of 50 L.s$^{-1}$. For the four flows, the flow rate was found after successive iterations of the couple 'discharge/height' of the downstream weir. More details can be found in estimating the flow rate with presence of obstacle in²².

**Determining the square cylinder height.** At the flume bottom, the artificial grass was glued to 1 cm thick PVC sheets, which in turn were glued to the glass bottom of the flume. The square cylinders were then screwed to the PVC sheets. Each cylinder partly crushed the grass at its location. Considering that the flow within the dense artificial grass is negligible (256 blades per square centimetre with 1 mm wide blades), the vertical distance from the cylinder top to the top of blades was to be determined. Measurements of this distance were made at 6 positions along the y axis for 230 positions along the x axis, i.e. a total of 1380 positions. 115 positions along the x axis correspond to the centres of cylinder rows (except for 4 rows that are not reachable by the measuring devices set on the traversing system). 115 other positions along the x axis correspond to the centres of rows without cylinders. The difference between the air depth measurement at the cylinder top and that taken at the bed level (dense meadow) is the cylinder height. We obtained a spatially averaged height $<k(x, y)>_{x,y} = 0.0592$ m as shown in Fig. 2, with $i$ the index along the streamwise direction and $j$ along the spanwise direction. This value is obtained using the following equation:

$$<k(x, y)>_{i,j} = \frac{1}{115} \sum_{i=1}^{115} \frac{1}{6} \sum_{j=1}^{6} k(x, y) \qquad (1)$$

where $k(x, y)$ is the average height at position $(x, y)$

**Determining the water depth.** The average water depth is defined as the difference of air depths with and without water (both measured using the ultrasonic probe). The average water depth $H(x, y)$ was measured at 24 positions along the x axis and at 3 lateral positions $y/(L/2) = 4,8$ and 12 as shown in Fig. 3.

The spatial average of the water depth $<H(x, y)>_{x,y}$ will be merely noted as $H$, with $i$ the index along the streamwise direction and $j$ along the spanwise direction, calculated as follows:

$$H = \frac{1}{24} \sum_{i=1}^{24} \frac{1}{3} \sum_{y=1}^{3} H(x, y) \qquad (2)$$

The convergence time of the time-averaged water depth was estimated at two different locations, in a central vein, and behind an obstacle model, where the vertical boundary layer was assumed to be established at 12 m. The convergence time is the time needed for the value to be stationary and was determined when the relative deviation tends toward zero. For each flow case, water depth is measured independently for a duration of $t=10, 20, 30, 40, 60, 90$ and 100 s. The convergence time found for each flow case is equal to 40 s. Once the convergence time has been established, the measurement is then carried out. It is ensured that all measurements are within the
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Positions (red markers) of the flow depth measurements using the ultrasonic probe. Top: global top view. Bottom: zoomed section indicating the positions along the y axis of the three longitudinal transects.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 4** Positions of the PIV measurement $y_1$, $y_2$ and $y_3$, and ADV measurements across the channel and central pattern; with $L = 14.3$ cm, $l = 6.4$ cm and $L - l = 7.9$ cm for the 4 flows.</figcaption></figure>

<table><caption>**Table 2.** Summary of the duration *t* for the convergence of time-averaged velocity and turbulence statistics for flow cases with *H/k* = 93% and 148% Measurements for elementary pattern located in the central vein at *y*/(L/2) = 6, at the edge of the square cylinder at *y*/(L/2) = 6.55 and behind the middle of the obstacle *y*/(L/2) = 7.</caption><thead><tr><th rowspan="3">H/k [%]</th><th colspan="3">Convergence time t [s]</th></tr><tr><th colspan="3">Pattern: central vein</th></tr><tr><th>y/(L/2) = 6</th><th>y/(L/2) = 6.55</th><th>y/(L/2) = 7</th></tr></thead><tbody><tr><td>93</td><td>400</td><td>400</td><td>300</td></tr><tr><td>148</td><td>300</td><td>300</td><td>300</td></tr></tbody></table>

interval defined as $2\sigma$ with $\sigma$ the standard deviation of the mean value $H(x, y)$. Obtaining a constant water level flow flushing the obstacle tops is difficult because of the variation in the local cylinder height $k(x, y)$ (Fig. 2). The average water depth at $x = 12$ m (study area) for case n°3 is $H = 58.3$ mm.

**Determining the velocity.** Velocity measurements were performed in the central vein at $y/(L/2) = 6$ and at positions $x = 4, 8, 12$ and $15$ m to investigate the longitudinal flow development. The emerged obstacles flows $H/k = 0.42, 0.93, 0.98$ are developed from $12$ m and the submerged obstacle flow $H/k = 1.48$ is developed from $8$ m. This is determined when the vertical profiles of average velocity overlap.

PIV measurements were performed in the vertical-longitudinal plane while ADV measurements were performed across the channel (vertical-lateral plane). They thus complement each other to explore the physical processes involved and can be compared. The convergence time of the time-averaged velocities and turbulence statistics was estimated for the ADV data at positions $y/(L/2) = 6, 6.55, 7$ (Table 2), and at 3 elevations $z$: near the channel bed, at $z/H = 0.40$ and near the free surface (Table 3). For the emerged obstacles flows $H/k = 0.42, 0.93, 0.98$, the convergence time is $t = 400$ s in the central vein and near the obstacle edge and $t = 300$ s behind the obstacle. For the submerged obstacle flow $H/k = 1.48$, the convergence time $t = 300$ s everywhere. In these experiments an elementary pattern extends from the middle of the prism $P_3$ to the middle of the prism $P_4$ (Fig. 4) and is located between $5 \le y/(L/2) \le 7$. This elementary pattern will be termed central pattern.

In the same open-channel flume, with only artificial grass over the bed (no house models) and higher flow rates, the experimental work of<sup>23</sup> showed duration *t* for the convergence of time-averaged velocity and turbulence statistics around 120 s.
<table><caption>**Table 3.** Elevations of the ADV velocity measurements across the channel for the four flows studied with $H/k = 42\%$, 93\%, 98\% and 148\%.</caption><thead><tr><th colspan="3">ADV Vertical measurement elevations in the transverse plan (Oy, Oz) at x = 12 m</th></tr></thead><tbody><tr><td colspan="3">42% immersion</td></tr><tr><td>z [mm]</td><td>z/H [-]</td><td>z/k [-]</td></tr><tr><td>11</td><td>0.44</td><td>0.186</td></tr><tr><td colspan="3">93% immersion</td></tr><tr><td>z [mm]</td><td>z/H [-]</td><td>z/k [-]</td></tr><tr><td>11</td><td>0.20</td><td>0.186</td></tr><tr><td>24.2</td><td>0.44</td><td>0.41</td></tr><tr><td>35.2</td><td>0.64</td><td>0.60</td></tr><tr><td colspan="3">98% immersion</td></tr><tr><td>z [mm]</td><td>z/H [-]</td><td>z/k [-]</td></tr><tr><td>11</td><td>0.19</td><td>0.186</td></tr><tr><td>25.6</td><td>0.44</td><td>0.43</td></tr><tr><td colspan="3">148% immersion</td></tr><tr><td>z [mm]</td><td>z/H [-]</td><td>z/k [-]</td></tr><tr><td>11</td><td>0.125</td><td>0.186</td></tr><tr><td>36.96</td><td>0.42</td><td>0.62</td></tr><tr><td>58.96</td><td>0.67</td><td>1.00</td></tr><tr><td>79.33</td><td>0.90</td><td>1.34</td></tr></tbody></table>

**Measurement of the water level fluctuations.** Free surface oscillations measurements were carried out along and across the flume with the ultrasonic transit time probe. Two types of measurements were carried out, the first one is a measure across the flume of $t = 200$ s for each position with a total of 35 positions and the same measure was made 5 hours after the first measured position, the second one is a measure of 12 hours at the position $x = 12$ m and $y/(L/2) = 13$. These oscillations are mainly in the transverse direction and are caused by the vortex shedding behind the obstacles, and are termed seiching phenomenon. Seiching is the most important for the flow with $H/k = 0.42$ with a maximum standard deviation of the level fluctuations equal to 1.5% in the centre of the channel for measurements taken at a time $t = 0$ h and after 5 hours. All flows have the same oscillation mode but the higher the water height, the lower the oscillation. The water level fluctuations stabilize after various durations depending on the flow cases: after 4 hours for $H/k = 0.42$, 1 hour for $H/k = 0.93$ and 0.98 and after 3 hours for $H/k = 1.48$. For more details about seiching phenomenon due to vortex shedding in the flow through cylinder arrays see$^{24–27}$ in this open-channel flume in particular, see the experiments of$^{28}$.

**ADV and PIV velocity measurements positions.** The positions of the velocity measurements using PIV and ADV are summarized in Fig. 4. The PIV measuring planes for all flow cases are located at lateral positions $y_1/(L/2) = 6$, $y_2/(L/2) = 6.55$ and $y_3/(L/2) = 7$ and at $x = 12.143$ m.

ADV measurements were carried out at $x = 12$ m, with 1 cm spacing the spanwise direction between two measurements across almost the entire width of the channel (Table 3), and at various elevations.

The central pattern, from $y/(L/2) = 5$ to 7, is indeed a large pattern (Fig. 4), as it corresponds to two patterns in the literature. Usually, in the literature a pattern ranges from $y/(L/2) = 5$ to 6 or 6 to 7 corresponding to half a vein and half an obstacle. An elementary pattern in the literature includes half a prism and half a central vein. Our elementary pattern is twice wider to verify the symmetry usually assumed in a central vein. In addition, the fact that a longitudinal row of square cylinders is present in the middle of the channel, i.e. in its privileged central vein (for a channel without obstacles), prevents certainty on the symmetry hypothesis.

At the position, $x = 12$ m, an ADV punctual velocity profiles are measured with a spanwise spacing of 1 cm between two profiles at the positions $5 \le y/(L/2) \le 7$, 15 positions along the y axis between the centre of the prism $P_3$ and the centre of the prism $P_4$ with 8 positions along the z axis for $H/k = 93\%$ and 15 positions along the z axis for $H/k = 148\%$ (Fig. 5). The measuring volume of the ADV has a cylindrical shape of 6 mm diameter and 7 mm length.

**Data Records**

The experimental data have been uploaded on the Zenodo website$^{29}$, under the DOI number https://doi.org/10.5281/zenodo.3871773. Four Excel files termed according to the value of immersion rate $H/k$ can be found. In each files, several sheets present the water depth, the ADV time-averaged velocities and turbulence statistics across the channel at $x = 12$ m (described in Table 3), the ADV velocities at $y = 6$ for 4 different positions $x = 4,8,12,16$ m. The last sheet describes the central pattern ADV velocities and turbulence statistics for the two flows with $H/k = 93\%$ and 148% flows.
<figure><img src="image_7.png" /><figcaption>Fig. 5 ADV measuring mesh realized within for the central pattern ($5 \le y/(L/2) \le 7$) for cases (left) $H/k=93\%$ and (right) 148%, the grey areas indicate the blind zones for measurements.</figcaption></figure>

<table><caption>Table 4. Rotation angle $\theta_z$ around the vertical axis of the ADV probe as a function of the position along the x axis.</caption><thead><tr><td>Position x [m]</td><td>Angle θ<sub>z</sub> [°C]</td></tr></thead><tbody><tr><td>4</td><td>0.073</td></tr><tr><td>8</td><td>-0.078</td></tr><tr><td>12</td><td>-0.0627</td></tr><tr><td>12</td><td>0.230</td></tr><tr><td>12</td><td>-0.220</td></tr><tr><td>15</td><td>-0.715</td></tr></tbody></table>

## Technical Validation

**Flow sensor and displacement system.** The inlet flow rate is controlled by a control valve Samson 3310 with servomotor PSQ, and is monitored with an electromagnetic flow meter (Krohne Waterux 3000 IFC 100). The discharge measurement uncertainty is 0.3% of the maximum range set either here $\pm 0.16 \text{ L.s}^{-1}$, calculated with manufacturer data for a maximum flow rate of $50 \text{ L.s}^{-1}$. The flow sensor was used throughout the measurement campaign showing that the flow rate used for a given water level was replicated and reproducible within the device's margin of error.

The channel is equipped with a motorized displacement system (Siemens S7-1200), carrying the measuring devices. The accuracy along the three space directions is of 0.1 mm. These measuring devices are an ultrasonic probe with transit time and an Acoustic Doppler Velocimetry (ADV) with a 3D side looking probe.

Note that the laboratory ambient air temperature is maintained by a thermostat between 19°C and 23°C.

**Water level.** The ultrasonic probe with transit time measures the distance between the probe and the first reflecting surface. It is used between the transmitter and the free surface and also between the transmitter and the bottom of the channel without water. This difference gives the water depth. This apparatus is manufactured by Baumer (UNDK 20I6903/S35A). Accuracy is better than 0.3 mm and the reproducibility is better than 0.5 mm according to the manufacturer. The chosen acquisition rate is 50 Hz (maximum value permitted). The height of water determined for a flow was tested during the entire measurement campaign. Replicability was therefore verified.

**Velocity measurements.** The ADV velocity probe chosen acquisition rate is 100 Hz for a given acquisition duration. The measurement uncertainty on velocity is $\pm 0.5\%$ of the measured value $\pm 1\% \text{ mm.s}^{-1}$ (manufacturer data). The centre of the measuring volume is located at 5 cm away from the probe transmitter.

For all the flows studied, the ADV measurement volumes are independent, with vertical intervals of 6 mm. It is essential to ensure that the ADV probe is orthogonal to the planes of space. The probe verticality is obtained by using a plumb line. The axis $Oz$ perpendicular to the bottom of the channel is therefore confused with the vertical. Since the exact alignment of the ADV probe with the transverse axis is difficult to obtain³⁰, a first test consists in placing the measurement volume near the wall and rotating the probe along the vertical axis in order to cancel the transverse mean velocity component. The second step consists of post-processing the velocity measurements assuming that the vertical average of the transverse component of the time-averaged velocity is nil at the wall. Hence, a rotation angle is deduced to ensure that the depth-averaged transverse velocity over the water level is zero. This rotation angle should be lower than 1 degree according to³⁰.

This rotation angle $\theta_z$ is applied to all lateral and longitudinal time-averaged velocity values as a function of the longitudinal position $x$, because this angle also depends on the verticality of the wall. The Table 4 summarizes the variation of the angle $\theta_z$ according to the position $x$. The longitudinal position $x=12 \text{ m}$ has 3 different angles because the probe has been taken off twice.
<table><caption>**Table 5.** Frequencies and time between two image burst used for the PIV measurements as a function of the lateral position y for each flow case.</caption><thead><tr><th rowspan="2">H/k [%]</th><th colspan="3">Frequencies [Hz]</th></tr><tr><th>y/(L/2) = 6</th><th>y/(L/2) = 6.55</th><th>y/(L/2) = 7</th></tr></thead><tbody><tr><td>42</td><td>4</td><td>3</td><td>3.5</td></tr><tr><td>93</td><td>5</td><td>3</td><td>3.5</td></tr><tr><td>98</td><td>5</td><td>3</td><td>3.5</td></tr><tr><td>148</td><td>2.5</td><td>3</td><td>4</td></tr><tr><th colspan="4">Time between two image burst [µs]</th></tr><tr><th>H/k [%]</th><th>y/(L/2) = 6</th><th>y/(L/2) = 6.55</th><th>y/(L/2) = 7</th></tr><tr><td>42</td><td>4000</td><td>8000</td><td>—</td></tr><tr><td>93</td><td>1000</td><td>20000</td><td>16500</td></tr><tr><td>98</td><td>1000</td><td>8000</td><td>16500</td></tr><tr><td>148</td><td>1000</td><td>6000</td><td>9000</td></tr></tbody></table>

<figure><img src="image_8.png" /><figcaption>**Fig. 6** Vertical profile of the longitudinal time-averaged velocity $\bar{u}$ at $x = 12$ m and $y/(L/2) = 6$, towards the glass sidewall located at $y = 0$, or towards the plexiglass sidewall positioned at $y = 1$; after taking into account the correction alignment angle of the probe (Table 4) for the flow $H/k = 0.93$.</figcaption></figure>

The particles used were polyamide powder of 50 µm median diameter (Evonik Vestosint 1164 white). They were used to seeding the flow during the ADV measurements. It enables to increase the signal to noise ratio and the correlations inside the measuring volume. For some ADV measurements, hollow glass particles (Hologlass sphere) of 10 µm diameter were also used. The inlet tank has not been cleaned up because the mixture of particles is not problematic for ADV measurements whereas it has been for PIV measurements. Indeed, for PIV measurements, two different particle sizes imply that the illumination will not be homogeneous. Only hollow glass particles were used for the PIV measurements.

The Fig. 6 shows that with a different rotation angle correction $\theta_z$ (Table 4) and less than 1°$^{30}$, it is possible at the same point to obtain at the measurement uncertainty almost the same velocity profile whatever the direction of the probe (towards $y = 1$ or $y = 0$).

As regards the PIV measurements, the artificial grass was cut so that the laser could go through it, so some perturbations are expected, and also we made sure that the laser sheet was straight and the calibration with the camera at a very low angle. The Table 5 presents the frequency and the time between two image bursts used for each flow case.

**Data usage caution.** It should be noted that aligning 833 square cylinders was not an easy task. We estimated a maximum shift of 5 mm in both longitudinal and transverse directions which represents 8% of the square cylinder width.

## Usage Notes

This dataset is intended to bring insight into the transition from emergence to slight submergence of house models focusing on their impact on the velocity field in a dense urban area. As extreme flood data are very scarce or non-existent in the field, this dataset could be used as a benchmark data for other experimental studies and numerical modelling of flows through a simplified urban area. This is relevant in the flood risk management.

Another use of this dataset is to compare the data with other configurations and to conclude upon the physical processes involved.

## Code availability

WinADV software was used to process the ADV velocity data. Data were filtered using the despiking technique of$^{31}$. The data measured with the PIV system were post-processed with Davis software from LaVision company.
Received: 16 June 2020; Accepted: 1 December 2020;

Published online: 11 January 2021

## References

1. European Union. Directive 2007/60/ce du parlement européen et du conseil, du 23 octobre 2007, relative à l'évaluation et à la gestion des risques d'inondation. *Official journal of the European Union* 27-34 (2007).

2. Oukacine, M. Etude expérimentale et numérique des écoulements à surface libre en présence d'obstacles émergés et faiblement submergés. *PhD Thesis, Université Paris-Est* (2019).

3. Herbich, J. B. & Shulits, S. Large-scale roughness in open-channel flow. *Journal of the Hydraulics Division* 90(6), 203-230 (1964).

4. Testa, G., Zuccalà, D., Alcrudo, F., Mulet, J. & Soares-Frazão, S. Flash flood flow experiment in a simplified urban district. *J. Hydraul. Eng.* 45, 37-44 (2007).

5. LaRocque, L. A., Elkholy, M., Chaudhry, M. H. & Imran, J. Experiments on urban flooding caused by a levee breach. *J. Hydraul. Eng.* 139(9), 960-973 (2013).

6. Huang, C.-J., Hsu, M.-H., Teng, W.-H. & Wang, Y.-H. The impact of building coverage in the metropolitan area on the flow calculation. *Water* 6(8), 2449-2466 (2014).

7. Zhou, Q., Yu, W., Chen, A. S., Jiang, C. & Fu, G. Experimental assessment of building blockage effects in a simplified urban district. *Procedia Engineering* 154, 844-852 (2016).

8. Velickovic, M., Zech, Y. & Soares-Frazão, S. Steady-flow experiments in urban areas and anisotropic porosity model. *Journal of Hydraulic Research* 55(1), 85-100 (2017).

9. Finaud-Guyot, P., Garambois, P. A., Dellinger, G., Lawniczak, F. & François, P. Experimental characterization of various scale hydraulic signatures in a flooded branched street network. *Urban Water Journal* 16(9), 609-624 (2019).

10. Li, X. et al. Laboratory modelling of urban flooding: strengths and challenges of distorted scale models. *Hydrol. Earth Syst. Sci.* 23, 1567-1580 (2019).

11. Neary, V. S., Sotiropoulos, F. & Odgaard, A. J. Three-dimensional numerical model of lateral-intake inflows. *J. Hydraul. Eng.* 125(2), 126-140 (1999).

12. Beretta, R., Ravazzani, G., Maiorano, C. & Mancini, M. Simulating the influence of buildings on flood inundation in urban areas. *Geosciences* 8(2), 77 (2018).

13. El Kadi Abderrezzak, K., Paquier, A. & Mignot, E. Modelling flash flood propagation in urban areas using a two-dimensional numerical model. *Nat. Hazards* 50, 433-460 (2009).

14. Guinot, V. & Soares-Frazão, S. Flux and source term discretization in two dimensional shallow water models with porosity on unstructured grids. *Int. J. Numer. Methods Fluids* 50, 309-345 (2006).

15. Guinot, V. Multiple porosity shallow water models for macroscopic modelling of urban floods. *Adv. Water Resour.* 37, 40-72 (2012).

16. Viero, D. P. Modelling urban floods using a finite element staggered scheme with an anisotropic dual porosity model. *J. Hydrol.* 568, 247-259 (2019).

17. Ferrari, A. & Viero, D. P. Floodwater pathways in urban areas: A method to compute porosity fields for anisotropic subgrid models in differential form. *J. Hydrol.* 589, 125-193 (2020).

18. Turki, S., Abbassi, H. & Nasrallah, S. Effect of the blockage ratio on the flow in a channel with a built-in square cylinder. *Computational Mechanics* 33, 22-29 (2003).

19. Mignot, E., Li, X. & Dewals, B. Experimental modelling of urban flooding: A review. *J. Hydrol.* 568, 334-342 (2019).

20. Oukacine, M. et al. Large Eddy Simulation for flows through emerged or slightly submerged square cylinders. In *Proc. 10th RiverFlow Conference on Fluvial Hydraulics* (Ed. Taylor & Francis.) 8 (Delft, 2020).

21. Proust, S. et al. Predicting the flow on the floodplains with evolving land occupations during extreme flood events (FlowRes ANR project). In *Proc. 3rd European Conference on Flood Risk Management FLOODrisk* (Ed. EDP Sciences) 04004 (2016).

22. Guillen-Ludena, S., Lopez, D., Mignot, E. & Rivière, N. Flow resistance for varying density of obstacles on smooth and rough bed. *Journal of Hydraulic Engineering* 146(2), 04019059 (2020).

23. Dupuis, V. Experimental investigation of flows subjected to a longitudinal transition in hydraulic roughness in single and compound channels. *PhD thesis Université Claude Bernard Lyon* 1 (2016).

24. Zima, L. & Ackermann, N. L. Wave generation in open channels by vortex shedding from channel obstructions. *Journal of Hydraulic Engineering* 128(6), 596-603 (2002).

25. Defina, A. & Pradella, I. Vortex-induced cross-flow seiching in cylinder arrays. *Advances in Water Resources* 71, 140-148 (2014).

26. Zhao, K., Cheng, N.-S. & Huang, Z. Experimental study of free-surface fluctuations in open-channel flow in the presence of periodic cylinder arrays. *Journal of Hydraulic Research* 52(4), 465-475 (2014).

27. Viero, D. P., Pradella, I. & Defina, A. Free surface waves induced by vortex shedding in cylinder arrays. *J. Hydraul. Eng.* 55, 16-26 (2017).

28. Chetibi, M., Proust, S. & Benamar, S. Transverse surface waves in steady uniform and non-uniform flows through an array of emergent and slightly submerged square cylinders. *Journal of Hydraulic Research* 58(4), 1-13 (2020).

29. Oukacine, M., Proust, S., Larrarte, F. & Goutal, N. Experimental flows through an array of emerged or slightly submerged square cylinders over a rough bed. *Zenodo* https://doi.org/10.5281/zenodo.3871773 (2020).

30. Peltier, Y. et al. Estimation of the error on the mean velocity and on the Reynolds stress due to a misoriented ADV probe in the horizontal plane: Case of experiments in a compound open-channel. *Flow Measurement and Instrumentation* 34, 34-41 (2013).

31. Goring, D.-G. & Nikora, V.-I. Despiking acoustic doppler velocimetry data. *Journal of Hydraulic Engineering* 128(1), 117-126 (2002).

## Acknowledgements

The authors acknowledge the support of the FlowRes project²¹ (https://flowres.inrae.fr/en/) by the French National Research Agency (ANR), under grant N° ANR-14-CE 03-0010.

## Author contributions

S.P. designed the configuration. M.O. conducted all experiments and processed all data. S.P. and N.G. provided support for conducting the tests and processing the corresponding data. F.L. provided scientific support. M.O. wrote the first draft of the data descriptor. All authors revised successive versions of the data descriptor.

## Competing interests

The authors declare no competing interests.
**Additional information**

**Correspondence** and requests for materials should be addressed to M.O.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver http://creativecommons.org/publicdomain/zero/1.0/ applies to the metadata files associated with this article.

© The Author(s) 2021