Turbulent jet through porous obstructions under Coriolis effect: an experimental investigation

The present study has the main purpose to experimentally investigate a turbulent momentum jet issued in a basin affected by rotation and in presence of porous obstructions. The experiments were carried out at the Coriolis Platform at LEGI Grenoble (FR). A large and unique set of velocity data was obtained by means of a Particle Image Velocimetry measurement technique while varying the rotation rate of the tank and the density of the canopy. The main differences in jet behavior in various flow configurations were assessed in terms of mean flow, turbulent kinetic energy and jet spreading. The jet trajectory was also detected. The results prove that obstructions with increasing density and increased rotation rates induce a more rapid abatement of both jet velocity and turbulent kinetic energy. The jet trajectories can be scaled by a characteristic length, which is found to be a function of the jet initial momentum, the rotation rate, and the drag exerted by the obstacles. An empirical expression for the latter is also proposed and validated.


Introduction
Rivers and wastewater discharge into lakes and coastal waters often form turbulent jets. Long-shore and cross-shore coastal currents can also be intended as jet-like flows in some cases when constrained by local topography or manmade structures. These types of environmental flows all have a profound effect on aquatic ecosystems due to the transport of turbulence, tracers, and sediment particles. At the same time, the surrounding environment affects their mixing and dispersion processes. Our ability to predict their behavior strongly depends on our capacity to describe how jets interact with the real environment they discharge into.
For this reason, the action played by Coriolis force cannot be disregarded. Similarly, the presence of obstructions in the receiving water body must also be given due consideration. Examining these mutual interactions could certainly be of help in mitigating undesired impacts and in developing best environmental management practices.
Mixing processes are rather well understood for unobstructed jets and are mainly controlled by the initial jet characteristics, the boundary conditions, and the hydrodynamic features of the ambient current (Smith and Mungal 1998;Nepf 2012;Liberzon and Fernando 2014). When such jets are discharged in the presence of rotation, the Coriolis force changes the orientation of turbulent eddies, especially affecting the energy cascading process (Lin and Atkinson 2000;Moisy et al. 2011). For density-driven jets, mixing and entrainment depends on bottom slope, Froude and Reynolds number (Fernando and Ching 1993;Praud et al. 2006; Thomas and Linden 2007;Cossu et al. 2010).
Our knowledge of the effect of vegetation on flows, in general, is rather extensive. The vegetation canopy, characterized particularly by its density and geometry (Oldham and Sturman 2001;Nepf et al. 1997;Ghisalberti and Nepf 2005;Poggi et al. 2009;Nepf 2012;Termini 2015Termini , 2018, interacts with the mean and turbulent flow, perturbing both advection and dispersion (Raupach and Thom 1981;Ben Meftah et al. 2015;De Serio et al. 2018;Righetti 2008). For canopy flows, the formation of coherent structures has been studied also numerically, e.g. Bailey and Stoll (2016).
Momentum jets obstructed by vegetation have been investigated only recently (Ben Meftah et al. 2015;Ben Meftah and Mossa 2016). It is well known that, when released in an unobstructed flow, the momentum jet experiences entrainment, while detrainment is observed only for density jets. Results by Mossa and De Serio (2016) and Mossa et al. (2017) demonstrated that the vegetation canopy is capable of inducing detrainment for strong momentum jets, which is unexpected behavior.
Most previous experimental studies of both obstructed jets in absence of rotation and unobstructed jets in rotating background have been conducted mainly in small-scale facilities (Etling and Fernando 1994;Lin and Atkinson 2000;Nepf 2012;Termini 2015). Even when large-scale facilities have been used the experiments have never investigated the joint influence of jets, obstacles, and rotation (Moisy et al. 2011;Ben Meftah et al. 2015;Davarpanah et al. 2020).
Although there is a good understanding of how turbulent jets interact either with rotation alone or with vegetation alone, to the best of our knowledge no one has so far experimentally investigated jets that are affected simultaneously by vegetation and the Coriolis force, always present in the field. Hence, this is the main aim and the novel aspect of the present paper. Although in real sea, due to real scales, the direct effect of rotation on turbulence induced by an obstructed pattern is expected quite negligible, the mean flow is modified by rotation and consequently, also the transport and spreading of turbulent kinetic energy and scalars (tracers) by the mean flow are modified accordingly.
As known, geophysical flows are generally characterized by time scales τ such that τ ≥ 1/ω, where ω is the Earth's rotation rate. Being geophysical flows often shallow, with the horizontal length scale L much larger than the vertical one H, for the horizontal and vertical velocities we have U >> W, where U and W are the horizontal and vertical velocity scales, respectively. Consequently, such flows are usually considered bi-dimensional, developing in the horizontal plane. They are characterized by a Rossby number Ro = U/(2ω L) which is less than unity, indicating that the Coriolis effect, i.e. rotation, prevails on the advective motion, and thus can effectively modify the flow dynamics.
To evaluate some effects induced by the Coriolis force on a turbulent flow, a turbulent Rossby number can be defined analogously as Ro' = u t '/(2ωl), being u t ' a turbulent velocity scale associated with an integral length scale l. Considering the typical order of magnitude for u t ' and ω, we generally get Ro' >> 1. Signifying the turbulent flow is unperturbed by Earth rotation (Godeferd and Moisy 2015), as the rotation timescale of 1/ω > > l/u t ' which is the timescale over which the turbulent structure evolves.
Addressing mesoscale flows, such as coastal flows, with typical horizontal length scales of O(10-100)km and velocity scales of O(0.1-1)m/s, the Rossby number tends to approach unity Ro≈1. In this case, the inertial and the rotation time scales become comparable so that a mutual influence between the mean flow and the rotation may become important. When such flows are in presence of obstructions, the turbulent velocity scale has an order of magnitude u t ' = (0.01 ÷ 0.1)U (De Serio et al. 2018). As for the integral length scale l, it can vary significantly (Nepf 2012;De Serio et al. 2018), depending on the density and location of the obstructions, spanning over O(0.1-100)m. Consequently, also Ro' can vary greatly. In such cases, when Ro' exceeds unity, the turbulent structures are unaffected by rotation directly. Instead, it is the transport and spreading by the mean flow of the turbulent fluctuations that are expected to be modified by rotation. Based on these scales, the examined conditions can be considered representative of rivers or wastewater discharges flowing into coastal waters through reeds and marsh beds, as well as offshore wind plants located in estuarine areas.
The paper outline is as follows. The experimental apparatus, the executed tests, and the methods used for processing experimental data are presented in Sect. 2. Section 3 includes examination and discussion of the obtained results in terms of the mean flow, jet width, turbulent kinetic energy, and transport of turbulent kinetic energy by the mean flow. Finally, some considerations about the jet trajectory are provided in Sect. 4.

Experimental setup and measurement technique
The experiments, carried out in the large-scale Coriolis rotating Platform, were performed using two configurations: in the first, the jet was only affected by background rotation; in the second, the jet flows through a canopy field consisting of rigid rods, mimicking rigid and emergent vegetation. To allow a comparison, the unobstructed and obstructed jets were examined under numerous identical clockwise rotation velocities of the tank Ω = 2 ∕T , T being the rotation period.
Besides the rotation rate, an additional variation parameter was the spacing between the rods, s. The Coriolis Platform has a diameter of D = 13 m, the working water depth was set at H = 0.80 m. The momentum jet was released from a dedicated pipe, rigidly fixed at 0.40 m depth ( Fig. 1), fed from a constant head tank. The outlet ensured a horizontal release of the jet with an initial diameter of d j = 0.08 m (internal pipe diameter) and initial flowrate Q j . A specially designed canopy was placed in the tank for the obstructed flow configuration. The rigid, emergent vegetation stems were modeled by an array of plastic transparent rods (made of Plexiglass®) d = 0.02 m in diameter, arranged on a 2 m × 2 m panel fixed at the bottom of the tank. The rods were manually mounted on the panel inside pre-drilled holes and constituted a regular pattern with the center-to-center distance s. The panel was placed in the tank as shown in Fig. 1. The center of the jet outlet (O) was 1 m away from the upstream edge of the panel and 0.77 m from its external edge.
The instantaneous measurements of the velocity field were made using a Particle Image Velocimetry (PIV) system. The laser source of the PIV (continuous Nd Yag laser, 532 nm, 25 W) was mounted in the center of the tank on a vertical axis. The laser sheet was set in a horizontal plane and was 5 mm thick. It covered a surface of 3 × 3 m 2 by means of an optical lens of 60° angle (Fig. 1a,b). Polyamide particles (Orgasol®) with density of 1016 kg/m 3 and diameter of 60 µm were used to seed the flow. They were added and mixed to the water in the head tank feeding the jet. The images were acquired using three synchronized CMOS cameras (two PCO® filming at 2560 × 2160 pixels resolution and one Falcon® filmings at 2432 × 1728 pixels resolution) mounted above the water surface, fixed on rigid supports located at 4.29 m, 4.74 m, and 4.87 m from the bottom, correspondingly. They filmed at a frequency of 33 Hz for 60 s, resulting in a total of 2000 images per measurement. This PIV acquisition frequency was properly selected to avoid possible disturbances of the free surface deformations on the measurements, that is the time delay between two successive PIV images was set to be much smaller than the typical time scale of the free surface disturbances induced by the jet.
The position of the cameras ensured the three fields of view (FOV1 corresponding to PCO1, FOV2 corresponding to PCO2, and FOV3 corresponding to FALCON) partially overlapped, to ensure full coverage of the jet propagation path (Fig. 2). The reference system has origin at the jet outlet O, with the X and Y axes in the horizontal plane, X directed along the jet longitudinal axis, and Y in the transverse direction, oriented towards the center of the tank; the Z-axis is vertical and positive upward (Fig. 1c). To map the velocity fields in the examined volume of water, the 2D velocity fields were acquired at six equispaced horizontal layers (Fig. 1c). Nevertheless, for the scope of the present study, only data assessed in layer Z = 0 m, i.e. at mid-depth, are examined.
The control number characterizing our experiments is the jet-based Rossby number Ro j = u j /(2ΩD), where u j is the jet velocity at the outlet (along X). We also define a jet Reynolds number Re j = u j d j /ν, with ν being the water kinematic viscosity. In the analyzed experiment, the jet flow rate Q j = 0.006 m 3 /s was used, corresponding to u j = 1.144 m/s setting the jet Reynolds number is in all cases at Re j = 91514, indicating a fully turbulent jet flow. Experiments executed at two different rotation rates are considered in this paper, Ω = 0.070 rad/s and Ω = 0.105 rad/s, respectively corresponding to T = 90 s and T = 60 s. For each value of T, two different densities of the obstructed array were examined. The low rods density was n = 25/m 2 for rods distanced s = 0.20 m from each other, resulting in a solid volume fraction occupied by obstacles of ϕ = 0.79%. The high rods density was n = 100/m 2 for rods distanced s = 0.10 m, resulting in solid volume fraction of ϕ = 3.14% (Nepf 2012). Such values of ϕ in the range [0.001 ÷ 0.04] can be considered representative of reeds and marsh beds (Liu et al. 2021) in a coastal environment. Table 1 lists all the analyzed experiments (EXP), for which similar analysis procedure was implemented. The jet-based Rossby numbers are both of order of magnitude O(10 -1 ), allowing us to evaluate the effect of rotation on the jet mean flow, expected to be more pronounced in the experiments characterized by T = 60 s.
For the sake of brevity, experiments without rotation carried out as control tests, are not described in the present study, as they showed typical and expected trends (   modus operandi for processing the PIV data from all examined experiments was as follows. The PIV images acquired by each camera were successively coupled and correlated, to obtain the PIV measured velocity fields in the horizontal plane. As for the PIV images, velocities along the X direction are named u, while those in the Y direction are named v. The open-source ParaPIV software (Wang 2018) running in Matlab® was used to process the acquired images. Extensive parallelization (16 cores) was implemented to reduce computational times. The interrogation window size was set at 32X32 pixels with a step size of 16 pixels. Two passes were used. The first pass was 64X64 pixels with a step size of 32, and the second pass was 32×32 pixels with an overlap of 16 pixels. The uniform distribution of the canopy base panel was used as the real-world-coordinates translation grid. The equidistant holes (in the bottom panel) were exposed completely in EXP15 and their centers were detected. The centers were used to perform the geometric transformation of all three cameras' images, aligning them on the same grid and correcting minor optical distortions.

Results and discussion
Reaching a steady flow condition of the jet was verified first, before implementing time-averaged analysis of the assessed data. Figure 3 shows, as an example, the cumulative average and cumulative rms of u velocity component in one selected point of the jet core, respectively for EXP15, EXP19, and EXP23. It is quite evident that an asymptotic value is achieved in all cases.

Jet velocity, width, and turbulent kinetic energy
A local reference system is conveniently placed to better follow the jet behavior, that is a curvilinear abscissa Xs along the jet trajectory with its origin at the jet outlet. The jet behavior was studied in the zone of established flow (Sami et al. 1967) for all the experiments, starting at Xs = Xs 0 . The value of Xs 0 is detected as the first distance from the outlet, greater than 10d j (Sami et al. 1967 (New et al. 2006) as the streamline starting from the center of the jet. In the present work, the trajectory is identified as the streamline starting from the central point of the crosssection at Xs = Xs 0 (Mossa et al. 2021). In the local reference system, at each point of the curvilinear abscissa Xs , the radial axis Ys is also placed (positive towards the center of rotation). Us and Vs are the streamwise and spanwise velocity components along Xs and Ys , respectively. The timeaveraged velocity fields and the detected jet trajectories are displayed in Fig. 4. Specifically, for all the experiments investigated (Table 1), some selected transversal profiles are plotted in the left panels of Fig. 4, while the overall velocity fields together with the jet trajectory are plotted in the right panels of Fig. 4. As previously written, the possible misalignment of the cylinders along the edges of the FOVs, essentially due to the different camera perspectives, was adequately taken into account in the data processing. The selected cross-sections are equidistant, at different Xs i , named with increasing values of -i suffix. In all cases, it is evident that the rotation affects the jet deflecting its trajectory. As expected, the jet velocity reduces with increasing distance from the outlet. We observe the maxima values of the velocity generally along the jet trajectory, especially in the fields of view FOV1 and FOV2. In the portion of the jet covered by FOV3, especially for EXP15, EXP16, and EXP19, the maxima velocities are slightly shifted towards the center of rotation.  Compared with the unobstructed cases, the presence of rods array determines a more rapid decrease of the velocity magnitude along the jet trajectory. The rate of velocity decrease is higher when the array density is higher (EXP23 and EXP24), as expected due to the increased blocking effect of the rods and the reduced drag effect. The comparison between experiments with identical rods array density (such as EXP19 and EXP18) shows higher velocity reduction rate at faster rotation rates, thus manifesting jet inertia limitation by Coriolis effect.
In Fig. 5 we plot the time-averaged streamwise velocities normalized by their maxima values, i.e. Us/Us max , varying along Ys (rated by Xs − Xs 0 ). Figure 5 shows steep gradients near the jet edges in all the experiments. For EXP15 and EXP16, the curves Us/Us max tend to collapse on a single curve as expected for free jets (Hussein et al. 1994) with a peak close to the center of the jet, thus following the jet self-similarity behavior. The Gaussian shape of the curve typical for free jets is also observed in the obstructed experiments for sections not intercepting the canopy (magenta marked trends in Fig. 5). On the contrary, when a section intercepts the canopy, Us/Us max has an undulating trend, with lower values corresponding to the location of obstructions and higher values to the free space between two consecutive obstructions (De Serio et al. 2018). Such undulating behavior is clearly observed in EXP19 and EXP18 for sections Xs 1 ÷ Xs 7 , as well as in EXP23 and EXP24 especially for sections Xs 1 ÷ Xs 5 (green marked trends in Fig. 5). Nevertheless, the envelope curve resembles the Gaussian shape. The symmetry with respect to the jet axis is generally preserved in EXP15 and EXP16. In presence of obstructions, the jet symmetry is almost lost inside the canopy. Moreover, in the obstructed cases the jet has the tendency to spread more towards the center of rotation, as revealed by the position of the Gaussian toe; such shift also intensifies with the increasing obstruction density.
Referring to the spanwise velocities Vs, they are one order of magnitude lower than the streamwise ones for all cases. In Fig. 6 they are plotted along each examined profile Ys∕(10d j ) as dimensionless Vs/0.1u j . Figure 6 demonstrates that the spanwise velocity components are mainly antisymmetric with respect to the jet axis, supporting previous observations by Mossa and De Serio (2016). A more detailed examination of the profiles gives information on entrainment or detrainment processes. In the experiments without obstructions (EXP15 and EXP16), for sections Xs 1 ÷ Xs 9 , the spanwise velocities in the core region are always directed towards the center of the jet, with similar entrainment as in free jets. Outside the core region, mainly for external points at Ys∕(10d j ) < − 0.25, Vs values are negative, i.e. directed outwards. For reference, see the trends marked in orange in Fig. 6 for EXP15 and EXP16. Starting from section Xs 9 and up to section Xs 13 , the jet trajectory represents a specific reference line. In fact, especially for EXP15, internal points ( Ys∕(10d j ) > 0) are characterized by positive Vs values, while in both EXP15 and EXP16 external points ( Ys∕(10d j ) < 0) are characterized by negative Vs (azure marked trends in Fig. 6). These observations show a tendency of the jet to a transversal spreading, mainly towards the external edge, suggesting that rotation also induces a detrainment at a distance from the outlet where the Coriolis effect prevails over the inertial one. Such distance is investigated in more detail in Sect. 3.2. In the obstructed cases (EXP19 and EXP18), the Vs values are mainly positive for Ys∕(10d j ) > 0, with exception for the initial sections of the set, and are generally negative for Ys∕(10d j ) < 0, thus showing a substantial detrainment towards the ambient flow. Apart from the last sections of the set, where the jet tends to close on itself, an analogous detrainment is also observed for the velocities in the obstructed cases EXP23 and EXP24 (see orange marked trends in Fig. 6 for EXP19, EXP18, EXP23, and EXP24). Such behavior agrees well with observations of Mossa and De Serio (2016), where jets interacting with rigid obstructions in absence of rotation were characterized by detrainment.
In each investigated cross-section, the jet transversal width B was obtained, defined as the total distance from the two points along Ys where the streamwise velocity tends to become negligible. We assume such condition is reached when Us = Us max /e, e being the Euler number. The obtained values of B are plotted in Fig. 7 along Xs∕d j , showing a linear trend with increasing distance Xs . The linear relation corresponds to typical free jets released in quiescent fluid, where the coefficient of proportionality between the jet radius and the downstream distance from the outlet is 0.1.
For same rotation rate T = 90 s, the jet widths are comparable for the cases EXP15 and EXP19 along the same distances Xs∕d j = 25 ÷ 40, rather for EXP23 the jet starts dispersing closer to the outlet. Moreover, for equal distances, the jet width is greater in EXP23 than in EXP19, consistently with the higher density of the rods, a condition for which the jet exhibits larger transversal deviation. Aanalogous behavior is observed also for the corresponding experiments with T = 60 s. B values are comparable for same distances Xs∕d j in EXP16 and EXP18, while they are greater in EXP24 and are reached at lower distances from the origin. The slopes of these linear trends (i.e. the coefficients of proportionality between B and Xs ) are 0.25 for high-density cases and 0.35 for both low-density obstruction array and unobstructed cases. Consequently, we can deduce that the jet dispersion (i) begins closer to the outlet for higher density array cases, (ii) is less rapid, over an equal distance along Xs , for unobstructed and low density obstruction array cases. Figure 8 displays the trend of Us/u j varying along the jet trajectory Xs∕Xs 0 . The exponential decreasing trend for Us is most evident for EXP15 and EXP16, similarly to free jets, starting from the region where the jet is well defined (Fig. 4), i.e. Xs∕Xs 0 > 2. Exponentially decreasing trends for Us values along Xs∕Xs 0 are also observed for all other cases, with generally high values of the correlation coefficient (Fig. 8). It is worth noting that while in the unobstructed cases a loss close to 50% of the initial value of Us occurs along a range on the trajectory equal to 2Xs∕Xs 0 , in the obstructed cases it is much more rapid, occurring inside the canopy, within 0.5Xs∕Xs 0 . This is especially evident for EXP19, where the initial values of Us are comparable with those of EXP15 and EXP16, but it is also valid for EXP23 and EXP24.
Since the vertical velocity component was not measured, we define a time-averaged turbulent kinetic energy in the horizontal plane as K = 1.33 2 u �2 + v �2 , following Ting and Kirby (1996). The prime symbol indicates the turbulent fluctuations of the considered velocity component, and the overbar indicates time-averaging over the measurement period. As displayed in Fig. 9, the reduction of Us is accompanied by the decay of K, with faster decay rates in presence of obstacles. For the unobstructed experiments, we note a clear exponential decreasing trend for K along the jet trajectory.
In the experiments with the higher obstruction array density, the decay is twice more rapid along the jet axis. For the same obstruction array density, along the jet path we observe higher values of K at higher rotation rate, indicating higher turbulent activity. This implies that, in our experiments, rotation has an impact also on the turbulent velocities, in the region where the jet is well developed.
The transport of K by the mean flow, that is by the streamwise and spanwise velocity (UsK and VsK, respectively), is analyzed along the jet trajectory, for each examined section. Based on the lower order of magnitude of Vs with respect to Us, the spanwise transport is generally lower than the streamwise one, despite the transversal jet dispersion. In Fig. 10, the transport UsK/u j 3 is displayed for some selected sections along the trajectory, showing peak values in the jet core, close to the trajectory. Moreover, it reduces in presence of obstructions, with a faster decay rate for higher density arrays, confirming that in such cases the jet mixes faster with the ambient fluid.
All the obtained results show that neglecting the presence of porous obstructions in the environment where the jet is issued leads to overestimation of spreading of both the jet streamwise velocity and its turbulent kinetic energy.

Jet trajectory and maximum penetration length of the jet
We apply the dimensional analysis to better investigate the trends of jet trajectory deviations, in the simplified hypothesis of a not vertically nor horizontally confined jet (i.e. sufficiently deep water and large tank). Thus, the main variables governing the phenomenon are the kinematic momentum of the jet M j , the rotation rate Ω , the water density ρ and molecular viscosity μ, a characteristic horizontal length scale of the jet L, and the drag force per unit planar area exerted by the obstructions . In this approach, in fact, the obstructions are accounted for by the drag force exerted in the X and Y directions, being the sum of the form and viscous drag components over the rod. As shown by Nepf (1999), various resistance laws for flow in porous media can be derived for such cases. Particularly, in open channel or in atmospheric vegetated flow, the quadratic form can be assumed to express the drag force per unit planar area: with i = 1, 2 and u 1 = u; u 2 = v in the current case. In Eq. (1) the drag coefficient C D is shown, which depends on density and diameter of the obstacles (Nepf 1999). The unit planar area addressed in Eq. (1) is a = nd , also equal to d/s 2 in the case of a periodic square array, with n being the density of the rods, i.e. the number of obstructions per unit planar area.
Consequently, we can express the studied phenomenon as a function f (M, Ω, , , L, ) = 0 . If we define the kinematic momentum of the jet M j as the product u j Q j , we can consider M j , the water density ρ, and the length scale L as three independent variables, thus arriving at equivalent function of the following non-dimensional numbers: The first no-dimensional number in Eq. Therefore, rewriting the first non-dimensional number in terms of L 0 , Eq. (2) can be rewritten as Trend of the turbulent kinetic energy along the jet trajectory, with fitting exponential curve for unobstructed cases or analogously as We take into consideration that when the obstructions are absent C D = 0 and thus Eq. (4) must reduce to L = L 0 , which is formally analogous to what was already found by Etling and Fernando (1994) in the case of unobstructed jets, with being a coefficient depending on the experimental conditions. Consequently, we now propose to assume for the function f 1 the exponential form f 1 = e C D , thus providing In our experiments, the characteristic horizontal scale L in the rotating tank is taken to be the maximum jet penetration distance from its outlet, named L C . Its coordinates in the fixed reference system are (Xc, Yc). Its values, for all the experiments, confirm that rotation effect prevails on the jet inertial one. Figure 11 displays all the experimental trajectories scaled by L C , demonstrating their collapse on a single curve. In this figure, we also plot the experimental results of Sobey (1972) for an unobstructed momentum jet of 0.79 cm 4 /s 2 , issued in a tank of 5 m diameter and rotating with Ω = 0.0751 rad/s. Nevertheless, in Sobey's tests (1972) the value of M 1∕2 j (2ΩL 2 ) is in the range similar to that of experiments, allowing a comparison of jet behavior. Our current results, compared with those of Sobey's (1972), confirm jets trajectories are homothetic.
To evaluate the coefficient C D we follow suggestion of Nepf (1999) for our values of ad = 0.01 and ad = 0.04, that correspond, respectively to our cases with sparse and dense obstacles array. We thus obtain C D = 0.7 for EXP 19 and EXP 23, and C D = 0.4 for EXP 18 and EXP 24. Figure 12 illustrates the comparison between the experimental Lc values and the theoretical ones obtained by the dimensional analysis with Eq. (5), revealing a high correlation (R 2 = 0.9), supporting the validity of the proposed dimensional analysis approach. The calibration of the coefficient provides = 5.55 for the unobstructed cases and = 2.25 for the obstructed cases.
This result is of high significance, allowing estimation of the possible maximum jet penetration distance in both cases of absence or presence of obstacles, based only on known parameters of the jet (i.e. Mj), of the rotating frame (i.e. Ω ), and of the obstruction array (diameter and density of the rods). Such information may, for example, allow evaluation of the effective deltaic area influenced by a river issuing in a vegetated canopy, or evaluation of the actual distance at which cages or pens in a coastal farm could be affected by any polluting coastal discharge, thus rendering our findings important in planning and managing engineering activities in coastal areas.

Conclusion
The present study means to broaden a more thorough knowledge of momentum jets issued in obstructed environments, under the effect of Coriolis force. The experiments were carried out in the large-scale Coriolis platform at LEGI Grenoble (FR), where a large set of velocity data was obtained by means of a PIV system while varying the rotation rate of the tank and the density of the obstructions. The main results of this experimental analysis can be summed up as in the following.
The rotation evidently affects the jet, inducing its propagation trajectory curving; at the same time, the jet velocity reduces with increasing distance from the outlet. Presence of obstructions array determines a more rapid decrease of the velocity magnitude along the jet trajectory. The rate of velocity decrease is higher when the array density is higher; while for equal array density, it is higher at greater tank rotation rates, thus manifesting the effect of Coriolis on the jet inertia.
The profiles of the time-averaged streamwise and spanwise velocity along the examined cross-sections show that the former has higher rates of decrease along jets trajectories. For unobstructed cases, a Gaussian shape, typical for free jets, is observed in the profiles of the streamwise component. As well, the jet resembles the self-similarity behavior. The spanwise velocities are antisymmetric with respect to the jet axis and specifically in the obstructed experiments show a jet detrainment towards the ambient flow inside the array, confirming previous results of Mossa and De Serio (2016) derived in absence of rotation. The variability of the jet width along the trajectory path highlights that the jet dispersion begins closer to the outlet for higher density array cases; nevertheless, over equal distances, it is less rapid for unobstructed and low-density obstruction array cases.
The exponential decreasing trend of the streamwise velocity along the jet trajectory is evident for unobstructed cases, similarly to free jets. Its abatement rate is higher in the obstructed cases than in absence of obstructions. The turbulent kinetic energy was also observed to decay, at exponential rate in the unobstructed cases. Such decay is more rapid with higher obstruction array density. These results are relevant, as they show that neglecting the presence of porous obstructions in the environment where the jet is issued means overestimating the spreading of both jet streamwise velocity and jet turbulent kinetic energy. Finally, in all our experiments, a characteristic horizontal scale is observed, that is the maximum jet penetration distance from its outlet, which scales all the experimental trajectories, being them homothetic. The study proves that this scaling length is a function of jet momentum, rotation rate, and drag coefficient of the canopy. A simple expression is proposed and validated for such maximum length, thus being a useful tool in the planning and managing engineering activities of a coastal area.