Abstract
A thorough understanding of the mixing and diffusion of turbulent jets released in porous obstructions is still lacking in literature. This issue is undoubtedly of interest because it is not strictly limited to vegetated flows, but also includes outflows which come from different sources and which spread among oyster or wind farms, as well as aerial pesticide treatments sprayed onto orchards. The aim of the present research is to analyze this process from a theoretical point of view. Specifically, by examining the entrainment coefficient, it is deduced that the presence of a canopy prevents a momentum jet from having an entrainment process, but rather promotes its detrainment. In nature, detrainment is usually associated with buoyancydriven flows, such as plumes or density currents flowing in a stratified environment. The present study proves that detrainment occurs also when a momentumdriven jet is issued in a notstratified obstructed current, such as a vegetated flow.
Introduction
The interaction of aquatic vegetation and, more generally, of porous obstructions in natural environments with a stream flow remains insufficiently explored and represents a challenge for scientists. In fact, multiple ecological, morphological and physical aspects characterizing plant canopies, such as shoot density, leaf length, plant stiffness, standing biomass, as well as plant submergence ratio, size and location, influence the flow hydrodynamics^{1,2}. On the other hand, hydrodynamics affect the structure of the underwater vegetation fields by dispersing spores, by mediating the availability of nutrients^{2,3}, and by exerting drag forces and associated turbulence^{4,5}, which could be responsible for vegetation development and survival^{6}.
There is a continuous and mutual interaction between the canopy and its surrounding flow. It is well known that seagrasses are able to significantly influence the hydrodynamic environment by reducing current velocity, dissipating wave energy and increasing deposition or retention of finer sediments, to which benthic invertebrates are sensitive^{7,8,9,10}. They have important ecological consequences, reducing turbidity and increasing light penetration, thus enhancing primary production and photosynthesis, which in turn guarantee their growth and reproduction.
Therefore, as also pointed out by Nepf^{11}, aquatic vegetation is of tremendous significance to many ecosystem functions: (i) promoting biodiversity and creating different habitats with spatial heterogeneity in the stream velocity; (ii) improving water quality due to the uptake of nutrients and production of oxygen^{12}; (iii) reducing coastal erosion and enhancing bank stability, as in the case of marshes mangroves and riparian vegetation^{13}.
For all the above mentioned reasons, vegetation hydrodynamics should not be strictly studied from a hydraulic perspective. Nepf^{11} observed that the presence of vegetation modifies the velocity field across several scales, relevant to different processes. For example, the uptake of nutrients by an individual blade depends on the boundary layer on that blade, i.e. on the bladescale flow. Similarly, the capture of pollen is mediated by the flow structure generated around individual stigma^{14}. On the contrary, the retention or release of organic matter, sediments, seeds and pollen from a meadow or patch depends on the flow structure at the meadow or patch scale^{15}. Furthermore, spatial heterogeneity in the canopyscale parameters and architecturally varying components may originate complex flow patterns, which is difficult to interpret, even if using data collected from real channels with live vegetation.
In this extremely complex context, the current research^{16,17}, carried out at the DICATECh of the Technical University of Bari (Italy), adds a turbulent jet as a new element interacting with a streamflow through a canopy. The aim is to provide a more thorough understanding of the hydrodynamic behavior of jets released in obstructed channels. In fact, it is worth to noting that vegetated channel flows have been studied in depth by many researchers^{11,18,19,20}. Because of their practical applications, ranging from effluents discharged in atmosphere and water bodies, to combustion and thrust control, also turbulent jets have been widely analytically, computationally and experimentally studied^{16,21,22,23,24}.
Nevertheless, studies on turbulent jets interacting with porous obstructions are still rare in literature. Moreover, the issue can also be of interest for its applicability not strictly limited to vegetated flows (Fig. 1a), but extended to a wider range of real cases. In fact, we can consider that also the planting of random or regular arrays of trees, used as protection and management systems of floodplains and banks, as well as solar power plants, offshore wind farms^{25} and oyster or mussel farms (Fig. 1b) could be generally considered as canopies. Hence, the release of water of boats between oyster farms, the airflow produced by spraying during pesticide treatments on orchards, submarine pollination can be considered as typical examples of this process.
In this paper, firstly a thorough study of the theory governing the phenomenon is carried out, secondly, some experimental results of momentum jets issued in a vegetated flow are examined. The theoretical study highlights an innovative aspect in the jet diffusion process, due to the interaction with the canopy. In fact, jet dynamics are generally associated with the entrainment process, i.e. the oneway transport process from the ambient fluid to the flowing turbulent fluid. As pointed out by Mc Climans^{26}, the entrainment process was defined in pioneering studies as the erosion by turbulence of surrounding nonturbulent fluid. It can be intended as the net transport of fluid from a less to a more turbulent fluid, considering that turbulence intensity can be conceived as a measurement of turbulence.
Nevertheless, in some specific configurations, jets experience detrainment, which is the process by which fluid is expelled from a turbulent flow. Detrainment is an important mechanism in some fluid flows in nature, often observed in cloud dynamics, when cloudy air is transferred outside of the cloud volumes, with consequent influence on the vertical heat and moisture fluxes in the atmosphere^{27}. Also in the ocean, detrainment is an important process, which explains the sediment transport from large density currents on the continental slope as well as the fate of hydrothermal plumes from submarine volcanic rift zones in the middle of the ocean^{27}.
In any case, to the authors knowledge, previous studies point out that detrainment can occur when a jet or a plume impinges on a stratified interface^{23,27,28,29}, without investigating into the necessary conditions for detrainment. Previous research studies^{27,30} analyzed plumes and thermals in a uniform environment, observing that a type of detrainment could be produced by buoyancy reversal, affecting the trajectory of the flow. It should be remarked that transport across a stratified interface is an essential aspect of many geophysical processes, such as pollution transport and dispersion by winds in urban areas or dispersion of effluents from sewage in coastal waters (Fig. 2). For this reason, to quantify the physical processes occurring at the interface, a profound understanding of turbulent mixing, entrainment and vortex dynamics at the interface has always been pursued, also considering that the more accurate these deductions are, the more transport models can be developed accurately^{27}.
In any case, all these previous research studies principally focused on the entrainment rate or the behavior of turbulence at the interface, rather than on the conditions for detrainment. This brief review shows that the physics of jet detrainment is still not fully understood, and even its occurrence in different flows is an unresolved question. While detrainment in nature is usually associated with buoyancydriven flows, like plumes or density currents flowing in a stratified environment, the present study demonstrates that detrainment occurs also when a momentumdriven jet is issued in a notstratified obstructed current, such as a vegetated flow.
Theoretical Framework and Results
Considering that generally jets are issued in an ambient fluid with the presence of currents^{31}, without loss of generality in terms of effects of the vegetation on the jets, below we will consider the case of a turbulent jet with a circular nozzle whose fluid is discharged into a large body of water with equal density, i.e. no buoyancy effects are present. Considering Fig. 3, the diameter of the jet at the exit is equal to D = 2b_{0} and the exit uniform velocity is equal to U_{0}. The ambient flow velocity is U_{e}. As wellknown, the fully developed flow starts downstream of the core region, whose length has an order of magnitude of 12b_{0}. Below we will consider only the fully developed region, where the longitudinal velocity distribution in each cross section, i.e. the distribution of the uvelocity in the radial direction, has the same wellknown Gaussian shape shown in Fig. 3. At each cross section the maximum longitudinal velocity component is u_{m} and b is the typical length scale, which is generally assumed as the distance from the jet centerline where the longitudinal velocity is u = u_{m}/2.
In the present paper we consider the case of jets issued in an ambient flow with a periodic square array of cylinders of uniform diameter d and distance s. Other key parameters of the cylinder array used in the present paper are the solid volume fraction ϕ, i.e. the volume within the canopy occupied by solid elements, which is the complement of the canopy porosity 1 − ϕ, the frontal area per unit volume of the canopy a = nd, which is equal to d/s^{2} in the case of a periodic square array, where n is the number of elements per unit planar area. For the sake of brevity, the theoretical development will report only some equations. Appendix 1 shows the complete theoretical development of a plane turbulent jet issued in a still obstructed fluid. The case of a compound circular jet has been preferred in the present analysis because, even if more complex from a theoretical point of view, generally it is much more common in nature. In any case, the mathematical development of a circular compound jet starts from the analysis of the simpler case of a plane turbulent jet in still water. Therefore, Appendix 1 enables the reader to better understand that reported below.
Within the array the flow is spatially heterogeneous at the scale of the individual elements and often unsteady in time. Generally, in the case of obstructed flows the doubleaveraging method is used to remove the temporal and elementscale spatial heterogeneity of the current^{32,33}. In other words, the instantaneous equations of the vegetated current are first averaged over a time longer than the time scale of turbulence or instabilities in the flow and then averaged over an infinitesimally thin area that spans many cylinders, including only the area occupied by the fluid. In the case of the present study, only the timeaverage operator will be used, since a spatial average of the jet flow would lack information on the typical spatial variations of the jet longitudinal and transversal sections. For the sake of simplicity, the ambient current U_{e} is considered uniform.
Using cylindrical coordinates, the motion equation in the longitudinal direction x is
where (referring to Fig. 3)
τ being the turbulent shear stress and D_{x} the drag force in the x direction, i.e. the resistance due to the solid medium, sum of form and viscous drag over the stem. As shown by Nepf^{34}, various resistance laws for flows in porous media can be derived. Particularly, in open channel or atmospheric vegetated flow, the quadratic form is used, which, in the case of a compound turbulent circular jet, can be written using local velocity scales as follows
Koch and Ladd^{35} verified that a quadratic drag law describes the resistance at moderatetohigh Reynolds numbers in cylinder arrays.
Integrating along the radial direction the motion eq. (1) and using eq. (3) for the drag force, we get
assuming that C_{D}a is constant with the radial direction. This hypothesis is rigorous only in the case of obstructions formed by circular stems concentric with the jet axis, which generally is not common in nature. Nevertheless, the assumed hypothesis is reasonable considering the stem configuration, the ratios between the jet scale (b), the stem diameter (d) and distance (s) which will be considered and shown below in the present paper. The integral of the left side of eq. (4) is the integral momentum of a compound pure circular jet. Therefore, eq. (4) shows an interesting result, since in an obstructed flow the integral momentum of the jet is not preserved along the longitudinal direction x, as in the analogous case of an unobstructed flow, when C_{D} is nil.
For the similarity analysis, we can write
being η = r/b.
Rewriting eq. (4), with consideration of eqs (2) and (5), we get
The similarity analysis allows to write also
Therefore, assuming
and inserting eqs (5), (7) and (8) in eq. (1), the motion equation becomes
Furthermore, we can assume simple forms for U_{m} and b, such as:
The case of strong jets
Referring to the case of strong jets, i.e. U_{e}/U_{m} ≪ 1, eq. (6) can be simplified and becomes
Therefore, taking into account eq. (10), we derive
In the same hypothesis of strong jets, eq. (9) can be simplified too, providing
Since the left hand side of eq. (13) is independent of x, the same must be for the right hand side, i.e.
In order to consider the last term of the right side of eq. (13) independent of x, some considerations are reported below. For eqs (10) and (14) we can write that
In the present study, some simplifications can be done, considering that we refer to jets with
where n = O(10–100).
In other words we are taking into account the case most frequently observed in nature, where d/s = O(10^{−1}–1), b/d = O(10–10^{2}), b/s = O(10–10^{2}), as shown in Fig. 4.
Therefore, with the above mentioned orders of magnitude, it is possible to write that, if
the order of magnitude of b changes when it becomes
i.e. when x^{q}increases of an order of magnitude or more. Therefore, b has the same order of magnitude between two values of x, i.e. from x_{1} to x_{2} > x_{1}, when
Since in the analyzed case q = 1, as derived from eq. (14), it is possible to write that
Therefore, along a longitudinal distance between x_{1} and x_{2} satisfying eq. (19), the last term of eq. (13) can be considered approximately constant in the limits above described.
Therefore, using eq. (12) we finally get
The case of weak jets
When the jets are weak, i.e. U_{e}/U_{m} ≫ 1, eq. (9) can be simplified and rewritten as follows
The first two terms of the right side of eq. (22) are independent of x. In order to consider also the last term of the right side of eq. (22) independent of x we can assume b with the same order of magnitude between two values of x, i.e. from x_{1} to x_{2} > x_{1}, for which eq. (19) is verified.
Furthermore, assuming
it is possible to assert that, considering two distance x_{1} and x_{2} > x_{1}, the order of magnitude of U_{m} does not change if
Considering the most restrictive hypothesis between (19) and (24), it is possible to conclude that, assuming that the right side of the eq. (22) is independent of x, also the left side must be independent of x.
In order to consider the left side of eq. (22) independent of x, we have
Furthermore, eq. (6) becomes
Therefore, eq. (26) can be written as follows
Considering eq. (8), from eq. (27) we can obtain
Therefore, using eqs (25) and (28) we get
For the intermediate case of jets with U_{m}/U_{e} ≈ 1 it is not possible to disregard some terms of the derived equations, therefore a simple exponential relation, as in the cases analyzed above, cannot be deduced.
Entrainment coefficient analysis
Assume
where x_{0} is the distance from the nozzle where the jet flow starts to be fully developed and C_{b} and C_{um} are the values of b and U_{m}, respectively, for x equal to x_{0}.
In the case of strong compound circular jets, the flow rate becomes equal to
where C_{Q} is a dimensionless coefficient which takes into account both the geometry of each cross section of the jet and the ratio between the average longitudinal velocity U_{b} and the maximum longitudinal velocity of each analyzed cross section. In other words, it is possible to write
where is the nominal outer boundary of the jet where u is close to U_{e}. Therefore, we assume the hypothesis that the jet cross sections preserve their original circular shape, which of course is rigorous only in the case of obstructions formed by circular stems concentric with the jet axis. The assumed hypothesis is generally reasonable and, in any case, it is always possible to consider a rectified circular cross section equivalent to the real cross section.
The variation of the jet flow rate along the longitudinal axes is equal to
where α_{e} is the entrainment coefficient, i.e. the ratio between the velocity of the entrained fluid v_{e} and U_{m}, with v_{e} the velocity transversal to the nominal outer boundary of the jet, which is oriented towards the jet centerline in the case of a positive entrainment coefficient and vice versa in the case of a negative entrainment coefficient, i.e. when a detrainment flow is present.
Considering eq. (34) and the derivative of eq. (32)
we get the following equation for the entrainment coefficient
A different result is obtained in the case of weak compound jet, where
Therefore
Considering eqs (34) and (38), it is possible to obtain the following equation of the entrainment coefficient
It is possible to highlight that in the case of a weak compound jet the variation of the jet flow rate is null as well as the entrainment coefficient, as shown by eq. (39), when obstructions are absent, i.e. with a nil bulk drag coefficient C_{D}. On the contrary, in the presence of obstructions, the jet experiences a detrainment process as shown by eq. (39).
Method
Experimental procedure
The experimental runs were carried out in a smooth horizontal rectangular channel in the hydraulic laboratory of the DICATECh of the Technical University of Bari (Italy). The channel was 25.0 m long, 0.40 m wide and 0.50 m deep. The lateral walls and the bottom surface of the channel were constructed of Plexiglas. The outlet and the inlet structures of the channel were connected to a hydraulic circuit, allowing a continuous recirculation of steady discharges. The channel was equipped with a series of stilling grids and a sidereservoir spillway with adjustable height in order to maintain a constant and uniform water head. In addition, it was equipped with two movable gates, placed at the inlet and the outlet of the channel, in order to regulate the channel flow rate and tailwater. At the downstream end of the channel, the water flow was intercepted by a rectangular reservoir, which was 3 m long, 1 m wide and 1 m deep, equipped with a Vnotch sharp crested weir, in order to measure the channel flow rate. For further details on the experimental apparatus, please refer to Ben Meftah et al.^{17}.
To simulate vegetation stems, a square array of rigid circular steel cylinders was used. The stem diameter, d, was equal to 0.003 m. The stems were inserted into a plywood board 3.0 m long, 0.398 m wide and 0.02 m thick, which in turn was fixed along the channel bottom, forming the experimental area. In order to reduce the effect of plywood board thickness on the experimental area, two other 3.0 m × 0.398 m × 0.02 m plywood boards, without vegetation stems, were attached to the channel bottom at both the upstream and the downstream side of this area. Stems were spaced longitudinally and transversally, with the same distance s = 0.05 m, so that the stem density was 400 stems/m^{2} and the vegetation density per meter was a = nd = dH/s^{2}H = d/s^{2} = 1.2 m^{−1}, with ϕ = nπd^{2}/4 = 0.00283 and n = 1− ϕ = 0.997.
The jet source was placed at the center of the experimental area, 15.0 m and 0.2 m from the inlet and the sidewalls of the channel, respectively. It consisted of a circular metallic pipe with a diameter, D, of 0.003 m. The jet was connected to a rectangular fiberglass tank by means of a plastic pipe. The tank was 1.0 m long, 0.5 m wide and 0.5 m deep and was positioned at a height of 3.6 m over the channel bottom surface. In order to maintain the jet discharge constant, the water flow was pumped continuously into the fiberglass tank by an electropump with a discharge larger than that of the jet. The water excess, distributed by the sidetank spillway, was driven via a pipeline to rereach the reservoir from where the electropump absorbed the water.
The instantaneous three velocity components were measured accurately using a 3D Nortek Acoustic Doppler Velocimeter (ADV) system. The ADV was used with a velocity range equal to ±0.30 m/s, a velocity accuracy of ±1%, a sampling rate of 25 Hz and a sampling volume of 27 mm^{3}.
Discussion
The main characteristics of the analyzed runs are shown in Table 1, where U_{0} is the initial jet velocity, Re is the channel Reynolds number and Re_{0} is the initial jet Reynolds number. Therefore, both channel and jet flows were turbulent. The set of experiments were limited to the case of weak jets, because of peculiar interest for many natural cases, such as the pollen jets released into the wind. Each run was performed twice, firstly without obstructions and successively with obstructions. The geometry of the obstructions remained unchanged for all runs.
As shown by Nepf^{34}, the bulk drag coefficient C_{D} is a function of the array density and, in the case of the present experiments, it was assumed equal to 1.2.
Figure 5 shows the experimental values of U_{m} along the channel, compared with the theoretical law of eq. (31), for the case of weak jets, with the error bars. Some points of the diagrams are outside of the error bars, because of local velocity variations from the theory due to the relative distance of the measurement points from the stems. Considering all these aspects, we can conclude that the model presented in this paper can satisfactorily reproduce the flow field and, furthermore, we can infer other conclusions on the jet entrainment or detrainment. The analysis of Fig. 5 allows us to deduce that the coefficient C_{um} of eq. (31) is not a universal constant, changing with the velocity of the ambient fluid, and therefore with its Reynolds number.
The entrainment process can be defined as a oneway transport process from the ambient fluid to the flowing turbulent fluid of a jet. The implication of the entrainment process is an increased volume flux, which has to be taken into account in the continuity equation for volume. Of course, the opposite process occurs in the case of a detrainment process, which the jets experience much more rarely and in presence of peculiar situations, not always known.
Figure 6a shows the trend of the flow rates of compound strong jets with obstructions, using eq. (32). Figure 6b shows the trend of the entrainment coefficients of compound strong jets using eq. (36). The figures show the values of the fully developed flow region, i.e. x/x_{0} > 1, where the analysis of the present paper is valid. The analogous values for the case of jets in unobstructed flows can be derived from eq. (36) assuming a nil bulk drag coefficient. As already written, the figures show a trend, since the coefficients of eqs (32) and (36) have been assumed of O(1), apart from C_{D} and a = d/s^{2}. It is important to highlight that all these coefficients can be correctly evaluated using experimental data. The entrainment coefficient of a compound strong jet without obstructions has a positive value, which does not change with x, as shown by eq. (36) with a nil bulk drag coefficient. On the contrary, in some configurations of obstructed flows, which depend on the canopy characteristics, the entrainment coefficients have positive values in a limited zone, close to the jet origin, with, consequently, an increasing jet flow rate. Beyond this zone, the entrainment coefficients become negative, revealing a reduction of the jet flow rates (Fig. 6b). In other configurations the entrainment coefficients are always negative. When the entrainment coefficient is negative, the jet flow rate decreases, revealing a release of the transported solute towards the ambient fluid.
Figure 7a,b display the same results for the case of a compound weak jet. For the coefficients of eqs (37) and (39) the same hypothesis as the strong jets has been assumed. In any case, it is important to highlight that the curve trend, as well as the conclusions reported below, does not depend on the coefficient values. As shown by eq. (39), the entrainment coefficient of a compound jet in an unobstructed flow is nil, i.e. the jet flow rate is preserved. On the contrary, the same jet issued in an obstructed flow is characterized by a negative entrainment coefficient, with the consequence of a flow rate release towards the ambient fluid. Particularly, the entrainment coefficient reaches a minimum negative value at a certain value of x/x_{0}, which depends on the canopy characteristics. The entrainment coefficient then starts to increase, remaining always negative. This conclusion is particularly interesting considering its effects on the ambient fluid, which experiences a release of the jet flow rate and, consequently, a release of its transported solute.
Conclusions
This paper investigates the behavior of strong and weak compound jets, issued in obstructed flows, from a theoretical point of view. For both strong and weak jets issued in obstructed flows, the paper presents new equations of: (i) the flow rate; (ii) the entrainment coefficient. It is worth noting that these equations deeply depend on the bulk drag coefficient and canopy density. Some results have been also validated by laboratory experiments.
The obtained equations prove that a detrainment flow occurs when the jets are issued in obstructed flows, analogously to the cases, generally shown in literature, of jets released in stratified flows. The principal findings of our study can be summarized as follows.

1
In the case of a compound strong jet, the entrainment coefficient does not change with x in unobstructed flows. On the contrary, in some configurations of obstructed flows, which depend on the canopy characteristics, the entrainment coefficients have positive values in a limited zone, close to the jet origin, with, consequently, increasing flow rates of the jet. Beyond this zone, the entrainment coefficients become negative, revealing a reduction of the flow rates. In other configurations the entrainment coefficients are always negative. When the entrainment coefficient is negative, the jet flow rate decreases, revealing a release of the transported solute towards the ambient fluid.

2
In the case of a compound weak jet issued in an unobstructed flow, the entrainment coefficient is nil, as the bulk drag coefficient is nil, and consequently the jet flow rate is preserved. On the contrary, in obstructed flows, the jet flow rate decreases since the jet experiences a detrainment process. Specifically, the jet reaches the maximum detrainment coefficient at a specific value of x/x_{0}. For greater values of x/x_{0}, the detrainment coefficient increases, without becoming positive. This conclusion is particularly interesting, meaning that the ambient fluid experiences a release of the jet flow rate and, consequently, of its transported solute.
Additional Information
How to cite this article: Mossa, M. and De Serio, F. Rethinking the process of detrainment: jets in obstructed natural flows. Sci. Rep. 6, 39103; doi: 10.1038/srep39103 (2016).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
Bouma, T. J. et al. Spatial flow and sedimentation patterns within patches of epibenthic structures. Continental Shelf Research 27(8), 1020–1045 (2007).
Mendez, F. J. & Losada, I. J. An empirical model to estimate the propagation of random breaking and nonbreaking waves over vegetation fields. Coastal Engineering 51, 103–118 (2004).
Koehl, M. A. R. How do benthic organisms withstand moving water? Am. Zool. 24, 57–70 (1984).
De Serio, F. & Mossa, M. Streamwise velocity profiles in coastal currents. Environmental Fluid Mechanics 14(4), 895–918 (2014).
De Serio, F. & Mossa, M. Analysis of mean velocity and turbulence measurements with ADCPs. Advances in Water Resources 81, 172–185 (2015).
Pujol, M. I. & Larnicol, G. Mediterranean Sea eddy kinetic energy variability from 11 years of altimetric data. J. Mar. Syst. 58, 3–4 (2005).
Gambi, M. C., Nowell, A. R. M. & Jumars, P. A. Flume observations on flow dynamics in Zostera marina (eelgrass) beds. Marine Ecology Progress Series, 61, 159–169 (1990).
Fonseca, M. S., Fisher, J. S., Zieman, J. C. & Thayer, G. W. Influence of the seagrass Zostera marina (L.) on current flow. Estuarine, Coastal and Shelf Science 15, 351–364 (1982).
Fonseca, M. S. & Koehl, M. A. R. Flow in seagrass canopies: The influence of patch width. Estuarine, Coastal and Shelf Science 67, 1–9 (2006).
Ondiviela, B. et al. The role of seagrasses in coastal protection in a changing climate. Coastal Engineering 87, 158–168 (2014).
Nepf, H. M. Hydrodynamics of vegetated channels. Journal of Hydraulic Research 508(3), 262–279 (2012).
Wilcock, R., Champion, P., Nagels, J. & Crocker, G. The influence of aquatic macrophytes on the hydraulic and physicochemical properties of a New Zealand lowland stream. Hydrobiologia 416(1), 203–214 (1999).
Pollen, N. & Simon, A. Estimating the mechanical effects of riparian vegetation on stream bank stability using a fiber bundle model. Water Resour. Res. 41, W07025 (2005).
Ackerman, J. Submarine pollination in the marine angiosperm Zostera marina. Am. J. Bot. 84(8), 1110–1119 (1997).
Zong, L. & Nepf, H. Spatial distribution of deposition within a patch of vegetation. Water Resour. Res. 47, W03516 (2011).
Ben Meftah, M., De Serio, F. & Mossa, M. Hydrodynamic behavior in the outer shear layer of partly obstructed open channels. Physics of Fluids 26(6), 065102 (2014).
Ben Meftah, M., De Serio, F., Malcangio, D., Mossa, M. & Petrillo A. F. Experimental study of a vertical jet in a vegetated crossflow. Journal of Environmental Management 164, 19–31 (2015).
Carollo, F., Ferro, V. & Termini, D. Flow Velocity Measurements in Vegetated Channels. J. Hydraul. Eng. 128(7), 664–673 (2002).
Poggi, D., Porporato, A., Ridolfi, L., Albertson, J. D. & Katul, G. G. The effect of vegetation density on canopy sublayer turbulence. BoundaryLayer Meteorology 111, 565–587 (2004).
Armanini, A., Righetti, M. & Grisenti, P. Direct measurement of vegetation resistance in prototype scale. Journal of Hydraulic Research 43(5), 481–487 (2005).
Jirka, G. H. & Harleman, D. R. F. Cooling Impoundments: Classification and Analysis. Journal of Hydraulic Division 105, 291–309 (1979).
Eroglu, A. & Breidenthal, R. E. Exponentially Accelerating Jet in Cross ow. AIAA Journal 36(6) (1998).
Morton, B. R. & Ibbetson, A. Jets deflected in a crossflow. Exp. Therm. Fluid Sci. 12, 112–133 (1996).
Muppidi, S. & Mahesh, K. Direct numerical simulation of round turbulent jets in crossflow. J. Fluid Mech. 574, 59–84 (2007).
Kinzel, M., Araya, D. B. & Dabiri, J. O. Turbulence in vertical axis wind turbine canopies. Physics of Fluids 27, 115102 (2015).
McClimans, T. A. Transport of pollutants in the Arctic by drift ice. IAHR Ice Symp. 3, 966–971 (1994).
Cotel, A. J. A review of recent developments on turbulent entrainment in stratified flows. Phys. Scr. T142, 014044 (2010).
Kumagai, M. Turbulent buoyant convection from a source in a confined twolayered region. Journal of Fluid Mechanics 147, 105–131 (1984).
Cotel, A. J. & Breidenthal, R. E. A model of stratified entrainment using vortex persistence. Appl. Sci. Res. 57 349–66 (1997).
Johari, H. & Paduano, R. Dilution and Mixing in an Unsteady Turbulent Jet. Experiments in Fluids 23(4), 272–280 (1997).
Rajaratnam, N. Turbulent Jets, Elsevier, 304 pp (1976).
Raupach, M. & Shaw, R. Averaging procedures for flow within vegetation canopies. Boundary Layer Meteorol. 22, 79–90 (1982).
Finnigan, J. J. Turbulent transport in flexible plant canopies. In ForestAtmosphere Interaction (ed. B. A. Hutchison & B. B. Hicks ), 443–480 (1985).
Nepf, H. M. Drag, Turbulence and diffusion in flow through emergent vegetation. Water Resource Research 35(2), 479–489 (1999).
Koch, D. l. & Ladd, A. J. C. Moderate Reynolds number flows through periodic and random arrays of aligned cylinders. Journal of Fluid Mechanics 349, 31–66 (1997).
Author information
Affiliations
Contributions
The theoretical analysis investigated in the paper was carried out by M.M. Both M.M. and F.D. discussed the theoretical and experimental results, wrote and revised the paper.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Electronic supplementary material
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Mossa, M., De Serio, F. Rethinking the process of detrainment: jets in obstructed natural flows. Sci Rep 6, 39103 (2016). https://doi.org/10.1038/srep39103
Received:
Accepted:
Published:
Further reading

Predictive model of bulk drag coefficient for a naturebased structure exposed to currents
Scientific Reports (2021)

The recent 500th anniversary of Leonardo da Vinci’s death: a reminder of his contribution in the field of fluid mechanics
Environmental Fluid Mechanics (2021)

Evaluating the hydrodynamics of a round jet in a vegetated crossflow through large eddy simulation
Environmental Fluid Mechanics (2019)

Vertical dense jet in flowing current
Environmental Fluid Mechanics (2018)

How vegetation in flows modifies the turbulent mixing and spreading of jets
Scientific Reports (2017)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.