Structure of a dense release produced by varying initial conditions

Buoyancy driven flows such as gravity currents, present in nature and human made applications, are conveyors of particles or dissolved substances for long distances with clear implications for the environment. This transport depends on the triggering conditions of the current. Gravity currents are experimentally simulated under varying initial conditions by combining three different initial buoyancies and five volumes of dense fluid released. The horizontal and vertical structures of the gravity currents are analysed and it is shown that the variation on the initial configuration is conditioning for these. Vertical transport through the gravity current is influenced at the bottom by the solid wall over which the current flows, and at the upper interface by the contact with the ambient water. The relative contribution of shear stress at the bottom and at the upper interfaces are estimated and analysed in terms of the initial triggering conditions of the current; these two compete with the buoyancy, the driver of the current, determining mixing and entrainment. By using a proper parametrization, which accounts for both initial volume of release and location of the observation position relative to the lock, a relation between the resistance at the bottom and at the upper interfaces with the initial conditions of release (i.e. the lock-length) has been found; this is found to be independent of the initial density in the lock. The present study shows that the variation of the initial conditions have consequences on (1) the configuration of the currents and on (2) the hydrodynamics of the currents, including mass and momentum exchanges, which are in addition mutually dependent.     

Shoaling internal waves may reduce gravity current transport

Gravity currents descending along slopes have typically been studied in quiescent environments, despite the fact that in many geophysical settings there is significant externally driven motion. Here we investigate how the head of a gravity current is influenced by interfacial internal waves at the pycnocline of a two-layer ambient water column. Our experimental measurements show that larger amplitude internal waves, interacting with the gravity current, reduce both the mass transport by the gravity current and its thickness. These results suggest that the ambient internal wave field should be considered when estimating transport by gravity currents in geophysical settings with strong internal waves, such as lakes and the coastal ocean.     

A numerical model of wave- and current-driven nutrient uptake by coral reef communities

We developed a numerical model capable of simulating the spatial zonation of nutrient uptake in coral reef systems driven by hydrodynamic forcing (both from waves and currents). Relationships between nutrient uptake and bed stress derived from flume and field studies were added to a four-component biogeochemical model embedded within a three-dimensional (3-D) hydrodynamic ocean model coupled to a numerical wave model. The performance of the resulting coupled physical-biogeochemical model was first evaluated in an idealized one-dimensional (1-D) channel for both a pure current and a combined wave-current flow. Waves in the channel were represented by an oscillatory flow with constant amplitude and frequency. The simulated nutrient concentrations were in good agreement with the analytical solution for nutrient depletion along a uniform channel, as well as with existing observations of phosphate uptake across a real reef flat. We then applied this integrated model to investigate more complex two-dimensional (2-D) nutrient dynamics, firstly to an idealized coral reef-lagoon morphology, and secondly to a realistic section of Ningaloo Reef in Western Australia, where nutrients were advected into the domain via alongshore coastal currents. Both the idealized reef and Ningaloo Reef simulations showed similar patterns of maximum uptake rates on the shallow forereef and reef crest, and with nutrient concentration decreasing as water flowed over the reef flat. As a result of the cumulative outflow of nutrient-depleted water exiting the reef channels and then being advected down the coast by alongshore currents, both reef simulations exhibited substantial alongshore variation in nutrient concentrations. The coupled models successfully reproduced the observed spatial-variability in nitrate concentration across the Ningaloo Reef system.     

Gravity currents in two-layer stratified media

An analytical, experimental and numerical study of boundary gravity currents propagating through a two-layer stratified ambient of finite vertical extent is presented. Gravity currents are supposed to originate from a lock-release apparatus; the (heavy) gravity current fluid is assumed to span the entire channel depth, H, at the initial instant. Our theoretical discussion considers slumping, supercritical gravity currents, i.e. those that generate an interfacial disturbance whose speed of propagation matches the front speed, and follows from the classical analysis of Benjamin (J Fluid Mech 31:209?C248, 1968). In contrast to previous investigations, we argue that the interfacial disturbance must be parameterized so that its amplitude can be straightforwardly determined from the ambient layer depths. Our parameterization is based on sensible physical arguments; its accuracy is confirmed by comparison against experimental and numerical data. More generally, measured front speeds show positive agreement with analogue model predictions, which remain strictly single-valued. From experimental and numerical observations of supercritical gravity currents, it is noted that this front speed is essentially independent of the interfacial thickness, ??, even in the limiting case where ?? = H so that the environment is comprised of a uniformly stratified ambient with no readily discernible upper or lower ambient layer. Conversely, when the gravity current is subcritical, there is a mild increase of front speed with ??. Our experiments also consider the horizontal distance, X, at which the front begins to decelerate. The variation of X with the interface thickness and the depths and densities of the ambient layers is discussed. For subcritical gravity currents, X may be as small as three lock lengths whereas with supercritical gravity currents, the gravity current may travel long distances at constant speed, particularly as the lower layer depth diminishes.     

Turbulent entrainment into fluid mud gravity currents

The entrainment of ambient water into non-Newtonian fluid mud gravity currents was investigated in this study. Constant volume release gravity currents were generated in a lock-exchange tank for a wide range of experimental conditions. A technique similar to the so-called light attenuation technique was used to find the boundary of the current, allowing for the calculation of both temporal and bulk entrainment parameters (in terms of the temporal and bulk entrainment velocities, respectively). It was found that the temporal entrainment velocity is dependent on different parameters in the different propagation phases. The slumping phase begins with an adjustment zone (henceforth, non-established zone) in which the temporal entrainment velocity is not a function of the current front velocity, followed by the established zone in which the temporal entrainment velocity is a function of the current front velocity. This dependence of the temporal entrainment velocity on the current front velocity carries through to the inertia-buoyancy phase. As expected, temporal entrainment velocity in the viscous-buoyancy phase was negligible in comparison to average entrainment velocity in the other phases. It is observed that the temporal entrainment characteristics in the non-established zone is governed by the competition between the entrainment-inhibiting density stratification effects and the entrainment-favouring effects of the Kelvin–Helmholtz billows that are quantified by the Richardson number and the Reynolds number of the gravity current, respectively. In the established zone, Reynolds number effects were observed to dominate over Richardson number effects in dictating temporal entrainment characteristics. A parameterization for the temporal entrainment velocity for non-Newtonian fluid mud gravity currents is developed based upon the experimental observations. This study also found that the bulk entrainment characteristics for the non-Newtonian fluid mud gravity currents can be parameterized by the Newtonian bulk entrainment parameterizations that rely solely on a bulk Richardson number. Interestingly, it was found that the non-Newtonian characteristics of the gravity current have little to no effect on the entrainment of the Newtonian ambient fluid.     

Mass-based depth and velocity scales for gravity currents and related flows

Gravity driven flows on inclines can be caused by cold, saline or turbid inflows into water bodies. Another example are cold downslope winds, which are caused by cooling of the atmosphere at the lower boundary. In a well-known contribution, Ellison and Turner (ET) investigated such flows by making use of earlier work on free shear flows by Morton, Taylor and Turner (MTT). Their entrainment relation is compared here with a spread relation based on a diffusion model for jets by Prandtl. This diffusion approach is suitable for forced plumes on an incline, but only when the channel topography is uniform, and the flow remains supercritical. A second aspect considered here is that the structure of ET’s entrainment relation, and their shallow water equations, agrees with the one for open channel flows, but their depth and velocity scales are those for free shear flows, and derived from the velocity field. Conversely, the depth of an open channel flow is the vertical extent of the excess mass of the liquid phase, and the average velocity is the (known) discharge divided by the depth. As an alternative to ET’s parameterization, two sets of flow scales similar to those of open channel flows are outlined for gravity currents in unstratified environments. The common feature of the two sets is that the velocity scale is derived by dividing the buoyancy flux by the excess pressure at the bottom. The difference between them is the way the volume flux is accounted for, which—unlike in open channel flows—generally increases in the streamwise direction. The relations between the three sets of scales are established here for gravity currents by allowing for a constant co-flow in the upper layer. The actual ratios of the three width, velocity, and buoyancy scales are evaluated from available experimental data on gravity currents, and from field data on katabatic winds. A corresponding study for free shear flows is referred to. Finally, a comparison of mass-based scales with a number of other flow scales is carried out for available data on a two-layer flow over an obstacle. Mass-based flow scales can also be used for other types of flows, such as self-aerated flows on spillways, water jets in air, or bubble plumes.     

Thin-layer models for gravity currents in channels of general cross-section area,a review

We present a brief review of the recent investigations on gravity currents in horizontal channels with non-rectangular cross-section area (such as triangle, \(\bigvee \)-valley, circle/semi-circle, trapezoid) which occur in nature (e.g., rivers) and constructed environment (tunnels, reservoirs, canals). To be specific, we discuss the propagation of a gravity current (GC) in a horizontal channel along the horizontal coordinate x, with gravity g acting in the \(-z\) direction, and y the horizontal–lateral coordinate. The bottom and top of the channel are at \(z=0,H\). The “standard” problem is concerned with 2D flow in a channel with rectangular (or laterally unbounded) cross-section area (CSA). Recent investigations have successfully extended the standard knowledge to the channels of CSA given by the quite general \(-f_1(z)\le y \le f_2(z)\) for \(0 \le z \le H\). This includes the practical \(\bigvee \)-valley, triangle, circle/semi-circle and trapezoid; these geometries may be in “up” or “down” setting with respect to gravity, e.g., \(\bigtriangleup \) and \(\bigtriangledown \). The major objective of the extended theory is to predict the height of the interface \(z=h(x,t)\) and the velocity (averaged over the CSA) u(x, t), where t is time; the prediction includes the speed and position of the nose \(u_N(t), x_N(t)\). We show that the motion is governed by a set of simplified equations, called “model,” that provides versatile and insightful solutions and trends. The emphasis in on a high-Reynolds-number current whose motion is dominated by buoyancy–inertia balance; in particular a GC released from a lock, which also contains general effects such as front and internal jumps (shocks), and reflected bore. We discuss two-layer, one-layer, and box models; Boussinesq and non-Boussinesq systems; compositional and particle-driven cases; and the effect of stratification of the ambient fluid. The models are self-contained, and admit realistic initial and boundary conditions. The governing equations are amenable to analytical solutions in some special circumstances. Some salient features of the buoyancy-viscous regime, and the estimate for the length at which transition to this regime takes place, are also presented. Some experimental support to the theory, and open questions for further investigations, are also mentioned. The major conclusions are (1) The CSA geometry has significant influence on the motion of the GC; and (2) The new theory is a useful, very significant, extension of the standard two-dimensional GC problem. The standard current is just a particular case, \(f_{1,2} =\) constants, among many other covered by the new theory.     

Curvature-induced secondary flow in submarine channels

The curvature-driven secondary flow in sinuous submarine channels has been a subject of considerable interest and controversy. Here, results from numerical model studies involving saline flow in laboratory-scale channels are presented. A 3D finite volume model of density and turbidity currents is used and simulations are run with different inflow discharges and channel-axis slopes. The simulation results show strong influence of bend wave length, channel gradient, confinement and cross sectional shape on the structure of secondary flow in submarine channels. Major findings are: (i) reversal of secondary flow in submarine channels is strongly associated with a tight bend characterized by a smaller wave length to width ratio or larger wave number, (ii) for the same inflow condition and planform characteristics, a trapezoidal channel cross section is more favorable to secondary flow reversal than a rectangular cross section, (iii) lateral convection resulting from the interaction between in-channel and overbank flows leads to the reversal of secondary flow in an unconfined channel at a lower channel slope than in a confined channel with the same dimensions, (iv) flow discharge has only nominal effect on the secondary flow in submarine channels.     

Experimental and numerical analysis of flow instabilities in rectangular shallow basins

Free surface flows in several shallow rectangular basins have been analyzed experimentally, numerically and theoretically. Different geometries, characterized by different widths and lengths, are considered as well as different hydraulic conditions. First, the results of a series of experimental tests are briefly depicted. They reveal that, under clearly identified hydraulic and geometrical conditions, the flow pattern is found to become non-symmetric, in spite of the symmetrical inflow conditions, outflow conditions and geometry of the basin. This non-symmetric motion results from the growth of small disturbances actually present in the experimental initial and boundary conditions. Second, numerical simulations are conducted based on a depth-averaged approach and a finite volume scheme. The simulation results reproduce the global pattern of the flow observed experimentally and succeed in predicting the stability or instability of a symmetric flow pattern for all tested configurations. Finally, an analytical study provides mathematical insights into the conditions under which the symmetric flow pattern becomes unstable and clarifies the governing physical processes.     

Pigeon homing in relation to geomagnetic,gravitational, topographical,and meteorological conditions

Summary Experienced homing pigeons were released at sites unfamiliar to them and with magnetic and gravity anomalies as well as in areas with rather normal fields throughout the FRG (41 releases when sunny, 14 when overcast; Figs. 1–3). The second-order release data were subjected to both univariate and multivariate statistical analysis (stepwise regression, factor analysis). The magnetic field strength and its gradients within the 1-km circle around the release site were determined from aeromagnetic maps of the anomalies of magnetic total intensity. Analogous variables were derived from gravity anomaly maps. It was tested whether the pigeons fly along that gradient to minimize the difference between the magnetic or gravity field at the release site and the loft at maximum rate. Further independent variables described magnetic K index, day-to-day variations of the magnetic components, topography, meteorological conditions, the number of the releases the pigeons had done, and the distance.Over magnetic anomalies widely varying in strength (departing – 250 nT to 300 nT from normal 600 m above ground), extent, and distance from loft, the pigeons vanished with less deviation from the homeward direction and faster than they did in areas with less irregular fields under sunny conditions; this is in contrast to other studies on magnetic anomalies, except one. At sites of gravity anomalies (15–49 mgal), the pigeons were significantly less homeward oriented and homed slower than at less anomalous sites (–9 to 14 mgal).Variables related to gravity were best predictors in 8 and and second predictors in 3 out of 15 regression analyses of the navigational parameters for the releases under sun. Six times the (absolute) amount of the gravity difference between the release site and the loft was selected first (Figs. 7B, C, 8B, C). The results suggest gravity to be involved in navigation as the pigeons' distance measure. Homeward directedness declined with increasing amount of the gravity gradient in the first 12 releases under sunny skies as well as when overcast (Fig. 8A). A preferred compass direction towards north-northeast was determined, being closest to the grand mean vector of the ascending gravity gradient (Fig. 6). The analyses failed to show directional preferences as assumed by the hypotheses tested. Temperature and degree of cloud cover provided some information for predicting mean vector lengths and mean vanishing times, respectively. The surface wind component in the homeward direction was correlated with median homing performances.     

Processesof Tamarix invasion and floodplain development along the lower Green River, Utah.

Significant ecological, hydrologic, and geomorphic changes have occurred during the 20th century along many large floodplain rivers in the American Southwest. Native Populus forests have declined, while the exotic Eurasian shrub, Tamarix, has proliferated and now dominates most floodplain ecosystems. Photographs from late 19th and early 20th centuries illustrate wide river channels with largely bare in-channel landforms and shrubby higher channel margin floodplains. However, by the mid-20th century, floodplains supporting dense Tamarix stands had expanded, and river channels had narrowed. Along the lower Green River in eastern Utah, the causal mechanism of channel and floodplain changes remains ambiguous due to the confounding effects of climatically driven reductions in flood magnitude, river regulation by Flaming Gorge Dam, and Tamarix invasion. This study addressed whether Tamarix establishment and spread followed climate- or dam-induced reductions in annual peak flows or whether Tamarix was potentially a driver of floodplain changes. We aged 235 Tamarix and 57 Populus individuals, determined the hydrologic and geomorphic processes that controlled recruitment, identified the spatial relationships of germination sites within floodplain stratigraphic transects, and mapped woody riparian vegetation cohorts along three segments of the lower Green River. The oldest Tamarix established along several sampling reaches in 1938, and 1.50-2.25 m of alluvium has accreted above their germination surfaces. Nearly 90% of the Tamarix and Populus samples established during flood years that exceeded the 2.5-year recurrence interval. Recruitment was most common when large floods were followed by years with smaller peak flows. The majority of Tamarix establishment and Green River channel narrowing occurred long before river regulation by Flaming Gorge Dam. Tamarix initially colonized bare instream sand deposits (e.g., islands and bars), and most channel and floodplain changes followed the establishment of Tamarix. Our results suggest that Tamarix recruitment was triggered by large annual floods that were followed by years with lower peak flows, not by periods of low flow alone. Tamarix appears to have actively invaded floodplains, while Populus colonization has been limited. Thus, Tamarix invasion may have greatly influenced floodplain development and riparian vegetation composition along the lower Green River since the early 20th century.     

Numerical simulation of particle-driven gravity currents

Particle-driven gravity currents frequently occur in nature, for instance as turbidity currents in reservoirs. They are produced by the buoyant forces between fluids of different density and can introduce sediments and pollutants into water bodies. In this study, the propagation dynamics of gravity currents is investigated using the FLOW-3D computational fluid dynamics code. The performance of the numerical model using two different turbulence closure schemes namely the renormalization group (RNG) ${k-\epsilon}$ scheme in a Reynold-averaged Navier-Stokes framework (RANS) and the large-eddy simulation (LES) technique using the Smagorinsky scheme, were compared with laboratory experiments. The numerical simulations focus on two different types of density flows from laboratory experiments namely: Intrusive Gravity Currents (IGC) and Particle-Driven Gravity Currents (PDGC). The simulated evolution profiles and propagation speeds are compared with laboratory experiments and analytical solutions. The numerical model shows good quantitative agreement for predicting the temporal and spatial evolution of intrusive gravity currents. In particular, the simulated propagation speeds are in excellent agreement with experimental results. The simulation results do not show any considerable discrepancies between RNG ${k-\epsilon}$ and LES closure schemes. The FLOW-3D model coupled with a particle dynamics algorithm successfully captured the decreasing propagation speeds of PDGC due to settling of sediment particles. The simulation results show that the ratio of transported to initial concentration C _o /C _i by the gravity current varies as a function of the particle diameter d _s . We classify the transport pattern by PDGC into three regimes: (1) a suspended regime (d _s is less than about 16 μm) where the effect of particle deposition rate on the propagation dynamics of gravity currents is negligible i.e. such flows behave like homogeneous fluids (IGC); (2) a mixed regime (16 μm < d _s <40 μm) where deposition rates significantly change the flow dynamics; and (3) a deposition regime (d _s ?> 40 μm) where the PDGC rapidly loses its forward momentum due to fast deposition. The present work highlights the potential of the RANS simulation technique using the RNG ${k-\epsilon}$ turbulence closure scheme for field scale investigation of particle-driven gravity currents.     

A review of gravity currents formed by submerged single-port discharges in inland and coastal waters

In this paper, the authors review the current state of the science on the dynamics of gravity currents generated by positively and negatively buoyant jet discharges from submerged round outfalls (i.e., a point source) in inland and coastal waters. Specifically, this article focuses on describing gravity currents occurring at both the bottom boundary and the free surface of the receiving fluid. The manmade discharge operations generating both types of gravity currents and their significance to sustainability of the surrounding hydro-environment are first described. The authors then summarize the flow regimes characteristics of these discharges before becoming gravity currents and how those flow regimes influence the dynamics of the gravity currents. The gravity current dynamics in the calm receiving waters are then analyzed. This analysis is followed by an analysis of the influence of the hydrodynamic forces (e.g., currents, turbulence, waves) on the dynamics of gravity currents. Finally, the authors review quantitative modeling approaches for different forms of gravity current, and identify the current knowledge gaps and research needs.     

Experimental investigation of cross-shore sandbar volumes

Sandbars are critical to the cross-shore movement of sediment. Prediction of cross-shore sandbar volumes requires knowledge about the functional relationship of sediment transport rate conditions with waves, currents, base slope, sediment property and water depth. In this study, experiments on cross- shore sediment transport were carried out in a laboratory wave channel for initial base slopes of 1/8, 1/10 and 1/15. Using regular waves with different deep-water wave steepness generated by a pedal-type wave generator, bar volumes caused by cross-shore sediment transport are investigated for beach materials with the medium diameter of d₅₀?=?0.25, 0.32, 0.45, 0.62 and 0.80 mm. A non-dimensional equation for sandbar volume was obtained by using linear and non-linear regression methods through the experimental data and was compared with previously developed equations in the literature. The results have shown that the experimental data fitted well to the proposed equation with respect to the previously developed equations.     

Experimental study of recirculating flows generated by lateral shock waves in very large channels

Due to the lack of data on hydraulic-jump dynamics in very large channels, the present paper describes the main characteristics of the velocity field and turbulence in a large rectangular channel with a width of 4 m. Although a hydraulic jump is always treated as a wave that is transversal to the channel wall, in the case of this study it has a trapezoidal front shape, first starting from a point at the sidewalls and then developing downstream in an oblique manner, finally giving rise to a trapezoidal shape. The oblique wave front may be regarded as a lateral shockwave that arises from a perturbation at a certain point of the lateral wall and travels obliquely toward the centreline of the channel. The experimental work was carried out at the Coastal Engineering Laboratory of the Water Engineering and Chemistry Department of the Technical University of Bari (Italy). In addition to the hydraulic jump formation, a large recirculating flow zone starts to develop from the separating point of the lateral shock wave and a separate boundary layer occurs. Intensive measurements of the streamwise and spanwise flow velocity components along one-half width of the channel were taken using a bidimensional Acoustic Doppler Velocimeter (ADV). The water surface elevation was obtained by means of an ultrasonic profiler. Velocity vectors, transversal velocity profiles, turbulence intensities and Reynolds shear stresses were all investigated. The experimental results of the separated boundary layer were compared with numerical predictions and related work presented in literature and showed good agreement. The transversal velocity profiles indicated the presence of adverse pressure gradient zones and the law of the wall appears to govern the region around the separated boundary layer.     

Experimental study of a positive surge. Part 2: comparison with literature theories and unsteady flow field analysis

A positive surge is a unsteady open channel flow resulting from the rapid rise of the free-surface. The phenomenon may be observed in water supply canals and channels as well as in some estuaries during spring tidal conditions. The formation and development of positive surges can be predicted using the method of characteristics and shallow water equations. The paper is the second part of a study presenting the results from new experimental investigations conducted in a large rectangular channel. Detailed unsteady velocity measurements were performed with a high temporal resolution using acoustic Doppler velocimetry and non-intrusive free-surface measurement devices. Several experiments were conducted with the same initial discharge (Q = 0.060 m³/s) and six different gate openings after closure resulting in both non-breaking undular and breaking bores. A comparison between main features of the undular surges with literature theories demonstrated that the experimental data were mostly in agreement with Andersen’s theory. The analysis of unsteady flow field including Reynolds stresses confirmed and extended previous findings about positive surge hydrodynamics.     

Direct numerical simulations of boundary condition effects on the propagation of density current in wall-bounded and open channels

The propagation of density current under different boundary conditions is investigated using high resolution direct numerical simulations (DNS). A revised Kleiser and Schumann influence-matrix method is used to treat the general Robin type velocity boundary conditions and the related “tau” error corrections in the numerical simulations. Comparison of the simulation results reveals that the boundary conditions change the turbulent flow field and therefore the propagation of the front. This paper mainly focuses on the effects of boundary conditions and initial depth of the dense fluid. The differences in energy dissipation and overall front development in wall-bounded and open channels are examined. Through DNS simulations, it is evident that with the decrease of initial release depth ratio ( $D/H$ ), the effect of the top boundary becomes less important. In wall-bounded channels, there are three distinctive layers in the vertical distribution of energy dissipation corresponding to the contributions from bottom wall, interface, and top wall, respectively. In open channels, there are only two layers with the top one missing due to the shear free nature of the boundary. It is found that the energy dissipation distribution in the bottom layer is similar for cases with the same $D/H$ ratio regardless the top boundary condition. The simulation results also reveal that for low Reynolds number cases, the energy change due to concentration diffusion cannot be neglected in the energy budget. To reflect the real dynamics of density current, the dimensionless Froude number and Reynolds number should be defined using the release depth $D$ as the length scale.     

Analysis of the velocity field in a large rectangular channel with lateral shockwave

In this work the authors describe the main characteristics of the velocity field of hydraulic jumps in a very large channel where lateral shockwaves occur. Experiments were carried out at the Coastal Engineering Laboratory of the Water Engineering and Chemistry Department of the Technical University of Bari (Italy). Extensive flow velocity measurements were investigated in order to have a clearer understanding of both hydraulic jump development and lateral shockwave formation in a very large channel. Eight experiments were performed in a 4m wide rectangular channel; the experiments differed in the inlet Froude number F ₀ and the jump type. Seven tests were carried out with undular jumps and one with a roller jump. The flow velocity and the flow free surface measurements were taken using a two-dimensional Acoustic Doppler Velocimeter (ADV) and an ultrasonic profiler, respectively. The experimental results can be summarized as follow: (i) the formation of well developed lateral shockwaves similar to those of oblique jumps were observed; (ii) the comparison of the experimental and theoretical data shows that the classic shockwave theory is sufficiently confirmed in the analyzed range of Reynolds number, taking into account the experimental errors and the difference between the theoretical and experimental assumptions; (iii) the transversal flow velocity profiles in the recirculating zone show a good agreement with the numerical simulations presented in literature in the case of a separated turbulent boundary layer over a flat plate. This conclusion enables us to confirm the hypothesis that the lateral shockwaves in the channel are the result of a boundary layer which, as observed, forms on the channel sidewalls.