The Effects of Hydraulic Jumps Instability on a Natural River Conﬂuence: The Case Study of the Chiaravagna River (Italy)

: The occurrence and the effects of hydraulic jump instabilities on a natural river conﬂuence in a small river basin in Liguria (Italy) is here investigated. Hydraulic jump instability has been extensively studied in controlled and simpliﬁed laboratory rectangular ﬂumes. In the present study, a scaled physical model of the Chiaravagna River and Ruscarolo Creek conﬂuence has been used, retaining the realistic geometry of the reaches. This reach has been subject to frequent ﬂoods in the last twenty years and the entire area of the conﬂuence has been redesigned to decrease the ﬂood risk. A series of experiments has been performed varying the discharge on the two reaches and the geometrical conﬁgurations. Free surface levels and two dimensional horizontal velocities have been measured in several positions along the physical model. The analysis of the water levels and velocities reveals that oscillations characterised by large amplitude and low frequency occur under particular hydraulic conditions. These oscillations have been found to be triggered by the hydraulic jump toe instability of the smallest reach of the conﬂuence. Aiming at reducing the amplitude of the oscillations, which can be of the order of the ﬂow depth, possible constructive solutions have been tested to control or damp the oscillations. Indeed, the insertion of a longitudinal dyke at the conﬂuence has proven to be an effective solution to limit the amplitude of the transversal oscillations.


Introduction
Hydraulic jumps are common phenomena in free surface flows in natural rivers, artificial canals and industrial applications. They occur whenever a supercritical high velocity flow impacts a subcritical low-velocity flow and the abrupt transition leads to a strong energy dissipation. The role of the jumps in hydraulics [1] and environmental engineering [2][3][4] is well known and this motivated a great number of scientific contributions. Starting from the seminal studies of Bakhmeteff [5] and Rouse and Siao [6], a considerable effort has been devoted to understand the complex turbulent flow involved, the prediction of the free surface profile and its stability. Applying the integral form of the conservation of momentum, it is found that the momentum of the supercritical flow balances the momentum of the subcritical one downstream in correspondence of the hydraulic jump. In the case of a prismatic rectangular channel, neglecting the flow resistance and velocity distribution effects, the balance provides the classical Bélanger relationship between the flow depth upstream (d 1 ) and design configuration in terms of flood risk. The main cause of the free surface oscillations is rooted in an instability mechanism similar to the hydraulic jump toe instability discussed in the above mentioned contributions. The high flow coming from the Ruscarolo Creek generated an oblique hydraulic jump, which was found to be unstable under certain flow discharge conditions and the resulting toe oscillation frequency was close to the transversal free mode of oscillation of the free surface downstream the confluence. The understanding of this mechanism suggested some modification of the originally designed configuration in order to limit the effect of the hydraulic jump instability.
The paper proceeds with Section 2.1, regarding the description of the study area, then the physical model, measuring techniques and data analysis are described in the rest of Section 2. Section 3 describes the results of the experiments performed on the different geometrical and hydrological conditions, followed by a discussion. Eventually, the conclusions are drawn in Section 5.

Study Area
The basin of the Chiaravagna River is placed on the Tyrrhenian side of the Ligurian Apennines, close to Genova Sestri Ponente, as shown in Figure 1. It is characterised by a narrow and elongated shape and it extends over an area of about 11 km 2 . The main reach is only 3.3 km long and it is characterised by limited ramifications with the exception of the Ruscarolo Creek; downstream the confluence it flows towards the airport of Genova, on the west side of the Polcevera River. Along its entire flow, the Chiaravagna River presents various obstacles: mining activities, the largest landfill in the Region and intense and chaotic urban development. After an intense flooding event that occurred in October 2010, the local Administration decided to start works on the riverbed in order to increase the hydraulic capacity and, thus, to decrease the flooding risk.

Physical Model Description
In order to test the designed configuration, a physical scaled model of the confluence of the Chiaravagna River and the Ruscarolo Creek has been prepared in the hydraulic Laboratory of the Department of Civil, Chemical and Environmental Engineering of the University of Genova. The geometric reduction scale, defined as the ratio of the lengths on the prototype (indicated with the P subscript) compared to the lengths on the model (indicated with the M subscript), has been chosen as λ = L P /L M = 35. Scale reproduction of free surface flows is typically achieved through a partial similarity of Froude, defined as Fr = V/ gd where V and d are the average velocity and depth of the current, g the acceleration of gravity, thus accepting that the Reynolds number is not preserved in the model [27]. Once the geometrical reduction scale λ has been set, Froude's similarity determines the scales for other quantities. In particular, fluid velocity, time and volumetric flow rate assume the following reduction scales: In this case, the correct representation of the flow resistance requires special attention, owing to the reduced value of the Reynolds number on the model. In fact, the model Reynolds number is reduced by about two hundred times compared to the one in the prototype. This implies, for example, that the flow in the physical model might be in the hydrodynamically rough regime, as it typically occurs in real river flows. On the contrary, a smooth flow regime is observed on the model, leading to a conductance coefficient C that depends on the Reynolds number alone and not on the dimensional roughness. In the present case, owing to the geometrical scale adopted and the dimensional roughness of the material used for the model construction, the flow falls in the smooth case and, using the Blasius resistance law valid in this condition, it turns out that the mean non dimensional smooth conductance coefficient of the model C M = 17 corresponds on the prototype (in the rough regime) to a Gauckler-Strickler coefficient of 55 m 1 3 /s, which is adequate for a concrete pavement. Moreover, when dealing with free surface flows, Chanson and Gualtieri (2008) noted that possible air entrapment in hydraulic jumps could not be correctly reproduced under a Froude similarity, owing to the different processes that regulate aeration, i.e., surface tension. Even though this is probably not the main mechanism driving the hydraulic jump instability.
A schematic view of the physical model is shown in Figure 2, where it is clearly visible the part of the model representing the confluence and the hydraulic circuit used to feed the model. In particular, the model represents the following reaches (see Figure 3, for section references): The Chiaravagna and Ruscarolo branches of the model have been fed with two separate piping starting from a common tank, placed on the ground floor of the laboratory. The flowing discharges are regulated through two butterfly valves manually operated. Moreover, in Figure 3, some other features are reported: in red and white are shown the underpinning walls (see the Figure S1 in the Supplementary Materials), in brown the sewage box (Figures S1-S3) and in filled black the ELSAG building' pillars (painted in yellow and blue in Figures S4 and S5) above which Via Manara, the building itself and a pedestrian bridge are placed. All the superstructures are represented with dashed lines: proceeding left to right, Via Giotto insists both on the Chiaravagna and Ruscarolo reaches without any obstacles on the flow (i.e., pillars), then, on the Ruscarolo Creek, an old pedestrian arch bridge is present (see Figure S6) and then, as already mentioned, Via Manara, ELSAG building and the pedestrian bridge. The black transversal element ( Figure 4 and Figure S7 in the Supplementary Materials), are the supporting beams of the ELSAG Building. Weirs can be found between Section 14.85 and 14.80, and Section 14 and 13.5 (see Figure S8), Section 30.1 and 29 and Section 12 and 11.9, respectively shown in Figues S9 and S5 (see Figure 3 for section references). Between Section 11.2 and 11, a car park has been built on the left bank of the reach. It leans on eight pillars placed directly on the riverbed (see Figure S10). In the bottom panel of Figure 3, there are three different position of a dyke: in Figure S11 position B is shown.

Measuring Techniques and Data Analysis
The physical model has been equipped with several measuring instruments. In particular, the flow discharge has been monitored using two electromagnetic flowmeters PROMAG P, one on each feeding pipeline. During each experiment, the water level and the velocities in different position of the model were recorded. In order to perform the water level measurements, an ultrasound gauge (Honeywell model 946-A4V-2D-2C0-380E, with 30 cm range and an accuracy of 0.2% of the full scale) has been used. The velocities have been recorded through a micropropeller (velocity range 0-1.2 m/s) and an acoustic velocimeter (Nortek Vectrino ADV probe with a velocity range ±0.01-4 m/s and an accuracy ±0.5% of measured value ±1 mm/s). All the instruments have been installed on a movable carrier: this allowed to place the instruments on several and repeatable positions. In Figure 4 the measuring point of the water level (panel c) and the velocity (panel a) are shown. Note that in the latter, all the velocity measuring points are represented. Indeed, a complete map of the velocity has been obtained only in few runs, owed to the complexity of these measurements in an unsteady situation, as the one here examined. In the following, most of the discussion of the velocity spectra is focused on signals measured in few stations along a single cross section located in the confluence area. Moreover, not all features of the model are here illustrated, e.g., just few ELSAG Building's pillars are shown in black. The water level measuring points have been used to produce longitudinal free surface profiles along several paths indicated in Figure 4b. All instruments have been connected to a DAQ unit (National Instrument model 6250). The ultrasound gauge and the micropropeller have been used with an acquisition frequency of 100 Hz, whereas the acquisition frequency of the acoustic velocimeter was set at 200 Hz. The acquired signals have been analysed in both time and frequency domains. In particular, probability density function and its moments have been calculated for each signal, providing estimates of its mean and standard deviation values. In the frequency domain, a spectral analysis has been performed, by computing the power spectral density (PSD) estimate of the longitudinal velocity v x (S v x ( f )) and, when appropriate, of the transversal velocity v y (S v y ( f )), using Welch's overlapped segment averaging estimator. The signals are divided into the longest possible segments to obtain as close to but not exceeding 8 segments with 50% overlap. Each segment is windowed with a Hamming window. The modified periodograms are averaged to obtain the PSD estimate.

Experimental Conditions
In the present study, a total of 28 runs have been performed, varying the volumetric discharge and the geometrical configuration of the model. In Table 1, the controlling parameters, i.e., the flow rates, the water depth of the Chiaravagna and Ruscarolo Creek and the corresponding Froude numbers are reported. In the final column, it is indicated whether free surface oscillations were observed during the runs.
To investigate the response of the model in its original design (see Figure 3 top panel), a first set of experiments were performed in different hydraulic scenarios, labelled 1C and 1R. In fact, the discharge downstream the confluence was set to the value corresponding to a return period of 200 years. Two combinations of discharges on the upstream reaches producing the same output discharge were tested, one assuming a 200 return time event on the Chiaravagna river and setting the discharge on the Ruscarolo by difference(1C), the other assuming a bicentennial event on the Ruscarolo and setting the discharge on the Chiaravagna by difference (1R). The experiment 2R has been performed with the same hydraulic condition of experiment 1R, but with the insertion of a longitudinal dyke. The experiments labelled A, B and C have been performed with the hydraulic condition of experiment 2R with the aim to test the dyke orientation that better damp the free surface oscillations associated with the toe jump instability. In particular, the longitudinal dyke, tested in the physical model, consists of wall of height equal to 6 cm, with a rounded tip triangular base of 2 cm width and 32.5 cm length, positioned upstream of the main beam of the ELSAG building (see Figure 3 bottom panel).
The series of experimental runs, labelled from 1 to 23, has been performed with the specific aim to find the critical conditions for the onset of instability. The procedure required to progressively increase the flow discharges starting from very low values. The experiments have been performed keeping the discharge ratio Q ru /Q ch constant and equal to 26%, i.e., the same hydraulic condition of experiment 1C, in which the maximum discharge is set on the Chiaravagna River. For each run the following quantities have been measured: the discharges in the two branches, the water levels upstream the confluence, both in the Ruscarolo Creek and in the Chiaravagna River, and the relative Froude numbers. Moreover, the velocity at three different locations has been recorded. The locations have been selected along the cross section of the river confluence, namely positions 19, 22 and 25 in Figure 4a, respectively in front of the Chiaravagna River, in the middle of the confluence and in front of the Ruscarolo Creek.

Results
In this section, we present the main results of the experimental campaign, focusing our attention on the flow depth and velocity measurement performed in several stations along the river reaches. The time averaged values of the water depth have been used to produce longitudinal free surface profiles along paths indicated in Figure 4b. The detailed analysis of the signals in terms of Power Spectral Density functions, that help in the identification of the characteristic frequencies, have been left for the Discussion section (Section 4).

Free Surface Local Measurements and Longitudinal Profiles for the Design Discharge
We first present the longitudinal free surface profiles that help to provide an overview of the mean water level along the reach under investigation. In particular, Figure 5a,b report the profiles along the longitudinal paths 3B and 4B, in Figure 4b, for their plan-view locations. Profile 3B starts from upstream the confluence in the centre of the Chiaravagna reach and, following the right bank, it ends downstream the main building. The profile 4B, instead, starts from the confluence and goes through the central aperture below the ELSAG building, ending downstream of the latter, as well as the profile 3B. The profiles, reported in Figure 5, refer to three different configurations, namely 1C, 1R, 2R. Together with the mean value of the water level, indicated with solid symbols, the recorded amplitude of the free surface has been indicated with the error bars. The error bars have been estimated as the mean of the measured water levels that exceeded twice the rms of the signal itself. All configurations share some common features in terms of the mean free surface profiles. Indeed, the two currents upstream are supercritical before entering the confluence, with a Froude number considerably larger than 1; in particular, the Froude number of the Ruscarolo Creek has a typical value of 2.5. The mean profiles are characterised by the presence of two strong hydraulic jumps. The first one is an oblique jump placed just downstream the confluence where the two rivers meet. The fact that the hydraulic jump is oblique is related to the geometry of the confluence. Indeed, the Ruscarolo Creek enters the confluence at an angle of about 30 degrees with respect to the axis of the Chiaravagna River, which is aligned with the stream-wise direction downstream the confluence. Although the noticeable difference in the flow discharges of the two reaches (i.e., the flow discharge of the Ruscarolo is about the 26% or 42% of the one of the Chiaravagna for the configuration C and R, respectively) the position and the intensity of this hydraulic jump are slightly affected by the discharge ratio. In fact, by inspecting the profiles of Figure 5 it is clear that the profiles labelled 1R-3B and 1R-4B are characterised by higher water levels and larger oscillations upstream the Manara Bridge. A second intense hydraulic jump is located downstream the river weir placed downstream of the ELSAG building. This hydraulic jump is less sensitive to the different configurations, since its position is far enough from the confluence where two currents are well mixed. Moreover, also this hydraulic jump is oblique with respect to the longitudinal direction and the reason is the asymmetrical geometry of the weir downstream the ELSAG building.
It is interesting to note that the amplitude of the free surface oscillations, represented by the error bars, tends to decrease moving downstream of the confluence. Moreover, in the configurations 1R, namely when the Ruscarolo has a greater discharge, the oscillations seem to be larger. The black profiles reported in Figure 5 correspond to the configuration 2R, where a longitudinal dyke is introduced to damp the transverse free surface oscillations. In particular, the longitudinal dyke was placed in position B, see Figure 3 bottom panel. This orientation has been selected after a comparison with the free surface measurements obtained testing other two dyke orientation, namely A and C. The effect of the dyke is twofold: on one side it effectively damps the free surface oscillations and, on the other side, it modifies the profiles of the mean water depth, inducing a different redistribution of the flow discharges in the three channels downstream the confluence. In Figure 6 the time signals recorded in positions 3B and 4B along the cross section 12, see Figure 4c, for the three positions of the dyke, namely A, B and C, are compared with the configuration 1R, without the dyke. The water levels recorded in the design configuration show large amplitude oscillations with typical period of few seconds with superimposed high frequency oscillations of small amplitude, the latter are strictly related to the turbulent flow. It is interesting to show the effect of the dyke on the amplitude of the oscillations and on the water levels inside the culverts. Considering the effect on the free surface oscillations damping, the most effective position is the one labelled A. In fact, the oscillation downstream the Chiaravagna River, on the right bank, is almost completely damped, see Figure 6a,b. At the same time, the mean water level is increased under both culverts. On the contrary, the dyke orientation C, Figure 6e,f, is not effective in damping the oscillations and leaves almost unchanged the mean water level. The position B has been selected based on the relative effectiveness in decreasing the amplitude of the oscillations and, at the same time, without modifying the mean water level. For the latter case, in Table 2 we report an estimate of the free surface oscillation amplitudes. In particular, for each profile reported in Figure 5, we have computed the average and maximum oscillation amplitude as a percentage of the measured flow depth. Note that the mean amplitude of the oscillations for the original design configuration ranges between 10.8% and 16.4% of the mean water depth, reaching maximum values around 30% in the configuration 1R. The effect of the dyke (2R) is to strongly reduce the oscillations for the profile 3B, whereas on the left bank (profile 4B) the damping is less pronounced, i.e., from 14.5% of 1R to 11.1% for 2R.

Time Velocity Distributions
The velocity measurements performed in terms of time signals of the longitudinal (v x (x, y, t)) and transversal component (v y (x, y, t)) are presented in this section. The velocity measurements have been performed in several points along the river reach, see Figure 4a

Discussion
The velocity signals show an oscillating character with a clear wide spectrum of frequency. The velocity for the flow under investigation is fully turbulent and, for this reason, the time records of the velocity have been decomposed in a time averaged components plus a fluctuating residual. The velocity residuals have been used to perform a spectral analysis with the aim to describe the characteristic frequency of oscillation of either the free surface instability and the jump toe instability. This approach is very common when the stability of an hydraulic jump is investigated [15][16][17]25,26].
The spectral analysis is presented on stations 19, 22 and 25, in Figure 4a, which lie along a cross section located in the confluence area. This helps to identify the variations of the dominant frequencies along the cross section itself. In the following discussion, we intend to show how the free surface fluctuations observed during the experiments, performed with the design discharge, can be rooted to the occurrence of a specific instability of the hydraulic jump, namely the hydraulic jump toe instability. Moreover, we discuss the role of the longitudinal dyke in decreasing the effect of the jump toe instability.

The Onset of the Free Surface Instability
A dedicated set of experiments on the design configuration, experimental runs from 1 to 23 reported in Table 1, have been performed with the aim of investigating the process that leads to the flow instability already introduced in the previous section.
The critical condition for the onset of the oscillations can be defined based on the spectral analysis of the velocity signals. In fact, the spectral signature of the oscillations should appear quite clearly as a strong peak in the low frequency range, whereas conditions where the flow is undisturbed show a classical scaling of energy spectra for fully three dimensional turbulence. Figure 9 summarises the stability condition of the whole set of experiments performed, using the total discharge as the main parameter. The different colours of the symbols indicate the experiment showing stable and unstable conditions, in blue and red, respectively. Starting from run 1, the onset of the instability has been observed between flow rates equal to 10.00 and 12.11 L/s. The falling branch has been investigated with differential steps much finer with respect to the rising one, looking for the conditions under which the oscillations naturally were damped. Note that the flow disturbances require a much lower discharge to offset with respect to the critical value for their onset. In fact, the flow returns stable for a discharge equal to 5.83 L/s. This behaviour indicates a possible hysteresis in the response of the system to the flow disturbances. Moreover, this aspect has motivated the rest of the experiments, performed with the aim to detect the critical conditions in more details on a new and finer rising branch. Following the runs from number 10 to 23, the instability has been detected more precisely for the total flow rate between 9.69 and 10.24 L/s, confirming the observations of the first experiments. Note that in previous works only unstable conditions have been investigated, as in [12,[15][16][17][23][24][25], without looking specifically at the onset of the instability.
In Figure 10, the spatial spectral analysis of the velocity signals are shown for the runs 11, 23 and 18. The velocities have been recorded along a cross section in the confluence area. The y-axis is oriented from Chiaravagna River towards the Ruscarolo Creek.
The three experiments selected correspond to flow conditions at which no oscillations occur (run 11), strong velocity fluctuations are clearly detectable (run 23) and, finally, at the onset of the main hydraulic jump toe instability (run 18).
Starting from the stable condition (panel a), the spectrum of the velocity signals displays no particular peaks. In panel b, a fully unstable solution is presented, whereby a clear peak is found for a frequency of about 0.24 HZ, which corresponds to a period of about 4.17 s, and spans the entire cross section. It is worth computing the natural frequency of the system using the dispersion relationship in the limit of shallow water [28], which reads: where ω = 2π/T is the wave frequency and T the wave period, k = 2π/L is the wave-number, L is the wavelength, in this case taken as half of the channel width and h is the local water depth and g the acceleration of gravity. In the specific case of run 23, Equation (3) provides a period of 4.13 s that is fairly close to the frequency detected in the velocity signals. This implies that the jump toe instability, which could be considered as a local instability, is able to trigger the oscillation of the entire domain, since its frequency is equal to the natural frequency of the system.  Finally, considering an experiment close to the critical conditions, e.g., run 18 shown in Figure 10c, a weaker peak is shown limited to the central region of the confluence.
Hydraulic jumps at high Froude numbers are well documented in the literature and free-surface instabilities and jump toe instabilities are often observed, see [12,[15][16][17][23][24][25]. In the cited works, the jump toe frequencies F toe or its non dimensional form in terms of the Strouhal number S t = F toe d 1 /V 1 , have been evaluated as a function of the Reynolds number Re = V 1 d 1 /ν or as a function of the supercritical upstream Froude number F r = d 1 / gV 1 .
In the left panel of Figure 11, the present estimation of S t as a function of the upstream Froude number of the Ruscarolo Creek is shown together with the measurements of previous studies. The present data are in fairly good agreement with the measurements of previous studies in terms of Strouhal and Froude numbers. Indeed, the Strouhal number seems to be less affected by varying the Froude number as discussed in [15,25]. Note that only Zhang et al. (2013) examined a range of Froude numbers comparable with the present values, whereas all other authors investigated much higher Froude numbers. Moreover, the Strouhal number, estimated in [16], depends on the measurements position, while the present case corresponds to the position closest to the jump toe. The critical conditions for the onset of the hydraulic jump toe instability have also be sought, in terms of critical Froude number. In the right panel of Figure 11, a stability plot is reported as a function of the Froude numbers of the Chiaravagna River and the Ruscarolo Creek, with the same colour coding of Figure 9 except for the orange dots, which represent experimental runs 1C, 1R and 2R. The data are nicely grouped in clusters: the unstable runs are characterised by Froude numbers of the Ruscarolo Creek (the branch with the observed oscillating hydraulic jumps) higher than 2.6, whereas the stable runs occur for Froude numbers of the Ruscarolo below 2.6. Only two runs do not follow this trend (runs 8 and 9) and, not surprisingly, they are both found on the falling branch discussed in Figure 9 with the occurrence of hysteresis process.
Finally, it is worth noting that the design discharge used to test the hydraulic efficiency of the river works, falls in the unstable branch and, in fact, the free surface presents high amplitude frequency that only with the insertion of a longitudinal dyke has been reduced.

The Effect of the Central Dyke in Damping the Velocity and Free Surface Oscillations
In this section we investigate the influence of the position of the dyke on the flow velocity in terms of oscillation damping. As a general comment, the dyke can influence the dynamics for two main reasons: (I) it represents a rigid longitudinal barrier able to block the transversal flow and damp the oscillations; (II) it segregates the hydraulic jump instability in a region downstream the Ruscarolo reach close to the left bank.
With the aim to understand the effect of changing the position of the dyke (see Figure 3 bottom panel), we performed specific measurements of the velocity along the same cross section used before in the analysis of the frequency of the hydraulic jump instability. Figure 12 reports the comparison between the PSDs calculated in absence of the dyke (panel d) and in the three configurations A, B and C (panels a,b and c, respectively). The effect of the position of the dyke appears quite clearly, comparing panel d (no dyke) with the other panels: with the dyke in position A, peaks in PSDs are detectable both in the Chiaravagna and the Ruscarolo reaches. Moving to position B (panel b), the peaks in the velocity spectrum are strongly reduced in both branches. Position C (panel c) produces a stronger confinement of the disturbances: a less intense peak is still detectable only in front of the Ruscarolo reach.

Conclusions
In this study we have experimentally investigated the hydraulic response of a real river confluence using a scaled physical model of the part of the Chiaravagna River and Ruscarolo Creek including a short undisturbed reach and a long reach after their confluence. To increase the flow capacity of the reach, upstream and downstream the river confluence, river works have been recently designed. The main purpose of the experiments was to test the effectiveness of the design solution. The river confluence has been tested for a flood event with return period equal to 200 years, and the experiments showed the occurrence of intense transversal free surface oscillations, which might dramatically increase the flooding risk. In fact, the amplitude of the oscillations was estimated to be up to 30% of the mean water level. The measurements have demonstrated that the transversal oscillations are triggered by the instability of the toe of the oblique hydraulic jumps of the Ruscarolo Creek in correspondence of the confluence. The typical jump toe instability frequency has been found to be close to the free mode of oscillation of the channel at the confluence. Thus, the jump toe instability is able to excite the transversal mode of instability of the entire channel. The measured non dimensional toe frequency, i.e., the Strouhal number, as a function of the super-critical Froude number have been found in good agreement with published studies. Moreover, the critical conditions have been analysed with the aim to find the critical Froude number for the onset of the hydraulic jump toe instability. The critical Froude number has been found to be around a value of 2.6. Finally, the possibility to add to the original project a longitudinal dyke placed within the river confluence has been investigated. The dyke has been found to be effective in the oscillation damping, confining the instability in a smaller portion of the cross section, namely on the left bank downstream the Ruscarolo Creek.
The results suggest that hydraulic jump instability might play an important role in realistic conditions of river confluence, especially for river reaches in supercritical conditions. This situation in not uncommon in several small basins with relatively high bed slopes. Thus, considering this unsteady phenomena might be important when flood defence river works are planned.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4441/12/7/2027/s1, Figure S1: View of the model confluence from upstream on the Chiaravagna River side, Figure S2: Detail of the model confluence, from right bank, Figure S3: View of the complete model from downstream, Figure S4: Detail of the ELSAG pillars and Via Manara Bridge, view from the left bank, Figure S5: Detail of the ELSAG pillars and Via Manara Bridge, view from downstream, Figure S6: View of the Ruscarolo Creek and Old Arch Bridge from downstream, Figure S7: Detail of the ELSAG pillars and pedestrian bridge, view from downstream: Chiaravagna River is on the left side of the figure, whereas Ruscarolo Creek is on the right, Figure S8: View of the Chiaravagna River from downstream, Figure S9: Detail of the dyke in position B and Ruscarolo weir, Figure S10: View from downstream of the ABB car park, Figure S11: Detail of the model confluence with dyke in position B, view from right bank.
Author Contributions: All authors contributed to the building of this paper in its different aspects. Acknowledgments: The study has been financially supported by ITEC Engineering, responsible for the river works design. We thank Ing. Pinasco of the Comune di Genova, public Authority in charge of the flood risk management of the City of Genova (Italy) and Maurizio Brocchini for the fruitful discussion on the experimental measurements.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: PSD Power Spectral Density F r Foude number S t Strouhal number