Numerical and wind tunnel modeling on the windbreak effectiveness to control the aeolian erosion of conical stockpiles

There are many industrial sites where open aggregate conical piles exist to store granular materials, like coal, industrial residuals, or other minerals. Usually these storage piles are placed in open areas, making them susceptible to wind erosion, which can create health, environmental, and/or economical concerns. It is common to minimize the dust emission through the placement of windbreaks in the vicinity of the storage piles, which reduces the wind speed in the vicinity of the pile’s surface. In this work, some experimental results from a wind tunnel study on the erosion of a conical sand pile, exposed or protected by a fence with porosities of 0, 70, and 83%, are shown. For the sheltered cases, the windbreak was placed at several distances from the pile’s leading edge, ranging from H to 4H, where H is the initial height of the non-eroded pile. The evolution of the shape of the sand pile is shown, at different instances in time, and the pile deformation quantified, using a novel experimental setup developed for wind erosion studies. This information might be regarded as a useful dataset for the benchmark of computational models aiming to produce the transient simulation of the aeolian erosion of stockpiles. The CFD results are comprised of the modeling of several experimental scenarios. The computational results for the surface wind velocity show a good correlation with the initial deformation of the pile. Based on the results, the isocontours of (u _s /u _r ) presented might be regarded as a good basis for the estimation of the pile shear velocity.     

Experimental assessment of characteristic turbulent scales in two-phase flow of hydraulic jump: from bottom to free surface

A hydraulic jump is a turbulent shear flow with a free-surface roller. The turbulent flow pattern is characterised by the development of instantaneous three-dimensional turbulent structures throughout the air–water column up to the free surface. The length and time scales of the turbulent structures are key information to describe the turbulent processes, which is of significant importance for the improvement of numerical models and physical measurement techniques. However, few physical data are available so far due to the complexity of the measurement. This paper presents an investigation of a series of characteristic turbulent scales for hydraulic jumps, covering the length and time scales of turbulent flow structures in bubbly flow, on free surface and at the impingement point. The bubbly-flow turbulent scales are obtained for Fr = 7.5 with 3.4 × 10⁴ < Re < 1.4 × 10⁵ in both longitudinal and transverse directions, and are compared with the free-surface scales. The results highlight three-dimensional flow patterns with anisotropic turbulence field. The turbulent structures are observed with different length and time scales respectively in the shear flow region and free-surface recirculation region. The bubbly structures next to the roller surface and the free-surface fluctuation structures show comparable length and time scales, both larger than the scales of vortical structures in the shear flow and smaller than the scales of impingement perimeter at the jump toe. A decomposition of physical signals indicates that the large turbulent scales are related to the unsteady motion of the flow in the upper part of the roller, while the high-frequency velocity turbulence dominates in the lower part of the roller. Scale effects cannot be ignored for Reynolds number smaller than 4 × 10⁴, mainly linked to the formation of large eddies in the shear layer. The present study provides a comprehensive assessment of turbulent scales in hydraulic jump, including the analyses of first data set of longitudinal bubbly-flow integral scales and transverse jump toe perimeter integral scales.     

A CFD study of wind patterns over a desert dune and the effect on seed dispersion

In the Namib Desert seed distribution is greatly influenced by wind patterns. Existing literature regarding wind patterns over dunes focuses on two-dimensional simulations of flow over simplified dune structures. The three-dimensional geometries of the sand dunes suggests far more complex flow features exist, which are not captured by two-dimensional simulations. Computational fluid dynamics (CFD) was used to reproduce the three-dimensional near surface wind patterns around a dune with the aim to learn more about seed distribution. Field work included terrain mapping, wind speed, direction and temperature metering. The CFD results show the expected two-dimensional flow features of high pressure at the dune toe, low pressure at the crest and flow acceleration up windward slope. Also observed are some three-dimensional flow features such as a spiral vortex near the crest and transverse flow due to crest-line curvature of the dune. It was also observed how the wall shear stress differs due to the three-dimensional shape of the dune. The wall shear stress suggests that seed accumulation is more likely to occur behind trailing (down-wind) crest edges. Particle tracking showed how seeds tend to move over the dune crest and recirculate towards the crest on the lee-side. The study showed that adding the third dimension makes the simulations more complex, adds to computational requirements and increases simulation time but also provides vital flow information which is not possible with two-dimensional simulations.     

Hydrodynamics and turbulence in emergent and sparsely vegetated open channel flow

This present study reports the results of an experimental study characterizing thorough variation of turbulent hydrodynamics and flow distribution in emergent and sparsely vegetated open channel flow. An emergent and rigid sparse vegetation patch with regular spacing between stems along the flow and transverse directions was fixed in the central region of the cross-section of open channel. Experiments were conducted in subcritical flow conditions and velocity measurements were obtained with an acoustic Doppler Velocimetry system. Large variations of the turbulence intensities, Reynolds shear stress, turbulent kinetic energy and vortical motions are found in and around the vegetation patch. At any cross-section through the interior of the vegetation patch, streamwise velocity decreases with increase in streamwise length and the velocity profiles converge from the log-law to a linear profile with increasing slope. Time-averaged lateral and vertical velocities inside the vegetation patch increase with increasing streamwise distance and converge from negative values to positive values. Turbulence intensities interior of the sparse vegetation patch are more than those of without the vegetation patch. Similar to the trend of streamwise velocity profiles inside the vegetation, turbulence intensities and longitudinal-normal Reynolds shear stress profile decreases with streamwise direction. In the interior of the vegetation patch and downstream of the trailing edge, turbulent kinetic energy profiles are exhibiting irregular fluctuations and the maximum values are occurring in the outer layer. Analysis of flow distribution confirms sparse vegetation patch is inducing a serpentine flow pattern in its vicinity. At the leading edge, flow is rushing towards the right hand sidewall, and at the trailing edge, flow is turning to the left hand sidewall. In between the leading and trailing edges, the streamlines are following a zig-zag fashion at varied degree along the streamwise and lateral directions. Immediate upstream of the leading edge and in the interior of the vegetation patch, vortex motion is clearly visible and the vortices are stretched along the width of the channel with streamwise direction.     

Aeolian erosion of storage piles yards: contribution of the surrounding areas

Dust emissions from stockpiles surfaces are often estimated applying mathematical models such as the widely used model proposed by the USEPA. It employs specific emission factors, which are based on the fluid flow patterns over the near surface. But, some of the emitted dust particles settle downstream the pile and can usually be re-emitted which creates a secondary source. The emission from the ground surface around a pile is actually not accounted for by the USEPA model but the method, based on the wind exposure and a reconstruction from different sources defined by the same wind exposure, is relevant. This work aims to quantify the contribution of dust re-emission from the areas surrounding the piles in the total emission of an open storage yard. Three angles of incidence of the incoming wind flow are investigated ( $30^{\circ }, 60^{\circ }$ and $90^{\circ }$ ). Results of friction velocity from numerical modelling of fluid dynamics were used in the USEPA model to determine dust emission. It was found that as the wind velocity increases, the contribution of particles re-emission from the ground area around the pile in the total emission also increases. The dust emission from the pile surface is higher for piles oriented $30^{\circ }$ to the wind direction. On the other hand, considering the ground area around the pile, the $60^{\circ }$ configuration is responsible for higher emission rates (up to 67 %). The global emissions assumed a minimum value for the piles oriented perpendicular to the wind direction for all wind velocity investigated.     

Evaluation of wind environment around a residential complex using a PIV velocity field measurement technique

Wind-tunnel simulations were employed to evaluate the wind environment around a tested residential area located near industrial complexes. The scaled-down geomorphological model of the test area was placed in the test section of a boundary layer wind tunnel. Particle image velocimetry (PIV) measurements were made in five vertical planes and one horizontal plane around the test area for two prevailing wind directions. The results showed that the wind speed decreased in the near surface layer and the velocity fluctuations increased in the upper region due to the presence of hills and high-rise buildings around the test area. Regions of flow separation and low-speed flow were found inside the test area for both the wind directions. The result suggests that the high-rise buildings should be well arranged with respect to the main wind directions to increase the natural ventilation inside the residential complex at the initial design stage.     

Modelling hydraulic jump using the bubbly two-phase flow method

Hydraulic jumps have complex flow structures, characterised by strong turbulence and large air contents. It is difficult to numerically predict the flows. It is necessary to bolster the existing computer models to emphasise the gas phase in hydraulic jumps, and avoid the pitfall of treating the phenomenon as a single-phase water flow. This paper aims to improve predictions of hydraulic jumps as bubbly two-phase flow. We allow for airflow above the free surface and air mass entrained across it. We use the Reynolds-averaged Navier–Stokes equations to describe fluid motion, the volume of fluid method to track the interface, and the k–ε model for turbulence closure. A shear layer is shown to form between the bottom jet flow and the upper recirculation flow. The key to success in predicting the jet flow lies in formulating appropriate bottom boundary conditions. The majority of entrained air bubbles are advected downstream through the shear layer. Predictions of the recirculation region’s length and air volume fraction within the layer are validated by available measurements. The predictions show a linear growth of the shear layer. There is strong turbulence at the impingement, and the bulk of the turbulence kinetic energy is advected to the recirculation region via the shear layer. The predicted bottom-shear-stress distribution, with a peak value upstream of the toe of the jump and a decaying trend downstream, is realistic. This paper reveals a significant transient bottom shear stress associated with temporal fluctuations of mainly flow velocity in the jump. The prediction method discussed is useful for modelling hydraulic jumps and advancing the understanding of the complex flow phenomenon.     

Large-eddy simulation of turbulent flow and dispersion over a complex urban street canyon

Turbulent flow and dispersion characteristics over a complex urban street canyon are investigated by large-eddy simulation using a modified version of the Fire Dynamics Simulator. Two kinds of subgrid scale (SGS) models, the constant coefficient Smagorinsky model and the Vreman model, are assessed. Turbulent statistics, particularly turbulent stresses and wake patterns, are compared between the two SGS models for three different wind directions. We found that while the role of the SGS model is small on average, the local or instantaneous contribution to total stress near the surface or edge of the buildings is not negligible. By yielding a smaller eddy viscosity near solid surfaces, the Vreman model appears to be more appropriate for the simulation of a flow in a complex urban street canyon. Depending on wind direction, wind fields, turbulence statistics, and dispersion patterns show very different characteristics. Particularly, tall buildings near the street canyon predominantly generate turbulence, leading to homogenization of the mean flow inside the street canyon. Furthermore, the release position of pollutants sensitively determines subsequent dispersion characteristics.     

Evaluation of a fast response pressure solver for flow around an isolated cube

This work describes and evaluates a pressure solver that has been incorporated into a fast response three-dimensional building-resolving diagnostic wind modeling system. The solver computes the three-dimensional pressure field around buildings and on exterior walls in terms of a coefficient of pressure by solving a simplified pressure Poisson equation (that neglects turbulence stresses in the Navier-Stokes) for incompressible flow. The input to the solver is the three-dimensional mean wind field obtained from a fast response empirical-diagnostic urban wind model. The present study is an evaluation of the pressure solver using wind-tunnel data for flow normal to and at a 45° angle to an isolated cubical building. Results for the normal incident wind angle case indicate that the model satisfactorily reproduces the general spatial patterns and the magnitude of the pressure difference around much of the cube. Details of the flow field that are not satisfactorily predicted include the spatial distribution of pressure on the roof and the lower half of the front side of the building and the magnitude along the sidewalls where pressures are over predicted. The results for the 45° case show reasonable agreement between the model and experiments on the front and the back walls, but over predict pressures on the leading edge of the rooftop. Regions with poor pressure predictions appear to be a result of unsatisfactory mean wind modeling.     

Effect of the topography of an industrial site on dust emissions from open storage yards

The erosion, transport and subsequent deposition of particles by wind has been identified as an important environmental problem worldwide. Investigations assessing fugitive dust emissions represent a relevant subject for various applications such as prevention of desert expansion, farmland erosion, quantification of dust emissions due to wind erosion from mining zones or open storage yards of granular materials. The present study, within this latter framework, aims to improve the methods of quantification of fugitive emissions from granular materials stockpiles on industrial sites. The methodology mostly used by the industrial operators to quantify dust emissions from their plants is based on the wind flow analysis and application of emission factors into such pattern. In this paper, three-dimensional numerical simulations were done to simulate wind flow over a real configuration of a power plant to highlight the influence of the surroundings buildings on the real exposure of the granular materials stockpiles. This study emphasizes the necessity to take into consideration the presence of surrounding buildings for the calculation of the particle emissions rate on stockpiles in order to obtain a more accurate and relevant evaluation of fugitive dust emissions from open storage systems on industrial sites.     

Turbulence investigation in the roughness sub-layer of a near wall flow

The turbulence behaviour along a wall roughened by pyramidal elements was analysed in the region extending from the apex of the roughness elements up to the external limit of the roughness sub-layer. The data used for the analysis were obtained by particle image velocimetry technique. The rough wall turbulent boundary layer flow is characterized by a relatively low Reynolds number. All the results on the rough wall were compared with data referring to the canonical flow on a smooth wall turbulent boundary layer. Mean values and turbulence quantities for the two flows collapse when approaching the external limit of the roughness sublayer. The quadrant analysis of the Reynolds shear stress, in the region near the surface, shows that the contribution of the sweep motions is about equivalent for the two flows (except for wall distances lower than 40 viscous units). The contribution of the ejection motions appears to be more important over the smooth wall than over the rough wall with increasing differences approaching the wall. The probability density functions of the streamwise fluctuating velocity field for the rough wall case appear to be positively skewed in the zone very close to the pyramid apex, in contrast with the behavior observed for the smooth wall case at corresponding distances from the wall. The integral and Taylor scales for the rough wall case appear to be strongly reduced by the presence of the roughness, while the Kolmogorov microscale shows higher values.     

Compound channel flow with a longitudinal transition in hydraulic roughness over the floodplains

Flows in a compound open-channel (two-stage geometry with a main channel and adjacent floodplains) with a longitudinal transition in roughness over the floodplains are experimentally investigated in an 18 m long and 3 m wide flume. Transitions from submerged dense vegetation (meadow) to emergent rigid vegetation (wood) and vice versa are modelled using plastic grass and vertical wooden cylinders. For a given roughness transition, the upstream discharge distribution between main channel and floodplain (called subsections) is also varied, keeping the total flow rate constant. The flows with a roughness transition are compared to flows with a uniformly distributed roughness over the whole length of the flume. Besides the influence of the downstream boundary condition, the longitudinal profiles of water depth are controlled by the upstream discharge distribution. The latter also strongly influences the magnitude of the lateral net mass exchanges between subsections, especially upstream from the roughness transition. Irrespective of flow conditions, the inflection point in the mean velocity profile across the mixing layer is always observed at the interface between subsections. The longitudinal velocity at the main channel/floodplain interface, denoted \(U_{int}\), appeared to be a key parameter for characterising the flows. First, the mean velocity profiles across the mixing layer, normalised using \(U_{int}\), are superimposed irrespective of downstream position, flow depth, floodplain roughness type and lateral mass transfers. However, the profiles of turbulence quantities do not coincide, indicating that the flows are not fully self-similar and that the eddy viscosity assumption is not valid in this case. Second, the depth-averaged turbulent intensities and Reynolds stresses, when scaled by the depth-averaged velocity \(U_{d,int}\) exhibit two plateau values, each related to a roughness type, meadow or wood. Lastly, the same results hold when scaling by \(U_{d,int}\) the depth-averaged lateral flux of momentum due to secondary currents. Turbulence production and magnitude of secondary currents are increased by the presence of emergent rigid elements over the floodplains. The autocorrelation functions show that the length of the coherent structures scales with the mixing layer width for all flow cases. It is suggested that coherent structures tend to a state where the magnitude of velocity fluctuations (of both horizontal vortices and secondary currents) and the spatial extension of the structures are in equilibrium.     

Spatial evolution of coherent motions in finite-length vegetation patch flow

A number of experimental studies on submerged canopy flows have focused on fully-developed flow and turbulent characteristics. In many natural rivers, however, aquatic vegetation occurs in patches of finite length. In such vegetated flows, the shear layer is not formed at the upstream edge of the vegetation patch and coherent motions develop downstream. Therefore, more work is neededz to reveal the development process for large-scale coherent structures within vegetation patches. For this work, we considered the effect of a limited length vegetation patch. Turbulence measurements were intensively conducted in open-channel flows with submerged vegetation using Particle Image Velocimetry (PIV). To examine the transition from boundary-layer flow upstream of the vegetation patch to a mixing-layer-type flow within the patch, velocity profiles were measured at 33 positions in a longitudinal direction. A phenomenological model for the development process in the vegetation flow was developed. The model decomposed the entire flow region into four zones. The four zones are the following: (i) the smooth bed zone, (ii) the diverging flow zone, (iii) the developing zone and (iv) the fully-developed zone. The PIV data also confirmed the efficiency of the mixing-layer analogy and provided insight into the spatial evolution of coherent motions.     

Instantaneous transport of a passive scalar in a turbulent separated flow

The results of large-eddy simulations of flow and transient solute transport over a backward facing step and through a 180° bend are presented. The simulations are validated successfully in terms of hydrodynamics and tracer transport with experimental velocity data and measured residence time distribution curves confirming the accuracy of the method. The hydrodynamics are characterised by flow separation and subsequent recirculation in vertical and horizontal directions and the solute dispersion process is a direct response to the significant unsteadiness and turbulence in the flow. The turbulence in the system is analysed and quantified in terms of power density spectra and covariance of velocity fluctuations. The injection of an instantaneous passive tracer and its dispersion through the system is simulated. Large-eddy simulations enable the resolution of the instantaneous flow field and it is demonstrated that the instabilities of intermittent large-scale structures play a distinguished role in the solute transport. The advection and diffusion of the scalar is governed by the severe unsteadiness of the flow and this is visualised and quantified. The analysis of the scalar mass transport budget quantifies the mechanisms controlling the turbulent mixing and reveals that the mass flux is dominated by advection.     

Initial contact,exploration and attachment of barnacle (Balanus amphitrite) cyprids settling in flow

Settlement responses of barnacle (Balanus amphitrite) cyprids to boundary-layer flows were examined in laboratory flume-experiments. The leading-edge configuration of flat plates was altered in order to manipulate flows without changing surface topography or freestream velocity. Settlement along the plates correlated strongly with downstream gradients in shear stress. Analyses of video images taken during the experiments indicate that cyprids first contact plates in regions where plate-ward advection is high, and subsequent exploratory movement along the plate is oriented with flow direction at the plate surface. After exploration, cyprids reject a surface more frequently in a fast flow (10 cm s^–1 freestream velocity) than in a slow flow (5 cm s^–1), but rejection occurs in shear stresses well below the threshold that would prevent attachment and exploration. A higher rejection rate does not result in lower settlement, however, since contact rate is higher in fast than slow flows. The movement of cyprids in flow thus appears to be a passive transport process during the initial contact stage of settlement, but an active behavioral response to flow direction and shear stress during later stages of exploration and attachment.     

Experimental investigation of bubbly flow and turbulence in hydraulic jumps

Many environmental problems are linked to multiphase flows encompassing ecological issues, chemical processes and mixing or diffusion, with applications in different engineering fields. The transition from a supercritical flow to a subcritical motion constitutes a hydraulic jump. This flow regime is characterised by strong interactions between turbulence, free surface and air–water mixing. Although a hydraulic jump contributes to some dissipation of the flow kinetic energy, it is also associated with increases of turbulent shear stresses and the development of turbulent eddies with implications in terms of scour, erosion and sediment transport. Despite a number of experimental, theoretical and numerical studies, there is a lack of knowledge concerning the physical mechanisms involved in the diffusion and air–water mixing processes within hydraulic jumps, as well as on the interaction between the free-surface and turbulence. New experimental investigations were undertaken in hydraulic jumps with Froude numbers up to Fr = 8.3. Two-phase flow measurements were performed with phase-detection conductivity probes. Basic results related to the distributions of void fraction, bubble frequency and mean bubble chord length are presented. New developments are discussed for the interfacial bubble velocities and their fluctuations, characterizing the turbulence level and integral time scales of turbulence representing a “lifetime” of the longitudinal bubbly flow structures. The analyses show good agreement with previous studies in terms of the vertical profiles of void fraction, bubble frequency and mean bubble chord length. The dimensionless distributions of interfacial velocities compared favourably with wall-jet equations. Measurements showed high turbulence levels. Turbulence time scales were found to be dependent on the distance downstream of the toe as well as on the distance to the bottom showing the importance of the lower (channel bed) and upper (free surface) boundary conditions on the turbulence structure.     

Influence of planform geometry and momentum ratio on thermal mixing at a stream confluence with a concordant bed

The effects of planform geometry and momentum flux ratio on thermal mixing at a stream confluence with concordant bed morphology are investigated based on numerical simulations that can capture the dynamics of large-scale turbulence. In two simulations, the bathymetry and asymmetrical planform geometry are obtained from field experiments and the momentum flux ratio is set at values of one and four. These two conditions provide the basis for studying differences in thermal mixing processes at this confluence when the wake mode and the Kelvin–Helmholtz mode dominate the development of coherent structures within the mixing interface (MI). The effects of channel curvature and angle between the two incoming streams on thermal mixing processes are investigated based on simulations conducted with modified planform geometries. Two additional simulations are conducted for the case where the upstream channels are parallel but not aligned with the downstream channel and for the zero-curvature case where the upstream channels are parallel and aligned with the downstream channel. The simulations highlight the influence of large-scale coherent structures within the MI and of streamwise-oriented vortical (SOV) cells on thermal mixing processes within the confluence hydrodynamics zone. Simulation results demonstrate the critical role played by the SOV cells in promoting large-scale thermal mixing for cases when such cells form in the immediate vicinity of the MI and in modifying the shape of the thermal MI within cross sections of the downstream channel—predictions consistent with empirical measurements of thermal mixing at the confluence. The set of numerical simulations reveal that the degree of thermal mixing occurring within the confluence hydrodynamic zone varies dramatically with planform geometry and incoming flow conditions. In some cases thermal mixing at the downstream end of the confluence hydrodynamic zone is limited to the MI and its immediate vicinity, whereas in others substantial thermal mixing has occurred over most of the cross-sectional area of the flow. Overall, the simulations highlight the flow conditions and the controls of these conditions that influence mixing within the immediate vicinity of a confluence.     

Self-aeration in the rapidly- and gradually-varying flow regions of steep smooth and stepped spillways

In high-velocity chute flows, free-surface aeration is often observed. The phenomenon is called self-aeration or white waters. When the turbulent shear stresses next to the free-surface are large enough, air bubbles are entrained throughout the entire air–water column. A rapidly-varied flow region is observed immediately downstream of the inception point of free-surface aeration. An analytical solution of the air diffusion equation is proposed and the results compare well with new experimental data. Both experiments and theory indicate that the flow bulking spans over approximately 3–4 step cavities downstream of the inception point of free-surface aeration on a stepped chute. Further downstream the void fraction distributions follow closely earlier solutions of the air diffusion equation. The application of the diffusion equation solution to prototype and laboratory data shows air bubble diffusivities typically larger than the momentum transfer coefficient. The result highlights however a marked decrease in the ratio of air bubble diffusivity to eddy viscosity with increasing Reynolds number. The finding might indicate some limitation of laboratory investigations of air bubble diffusion process in self-aerated flows and of their extrapolation to full-scale prototype applications.     

Canopy-wake dynamics and wind sheltering effects on Earth surface fluxes

The atmospheric boundary layer adjustment at the abrupt transition from a canopy (forest) to a flat surface (land or water) is investigated in a wind tunnel experiment. Detailed measurements examining the effect of canopy turbulence on flow separation, reduced surface shear stress and wake recovery are compared to data for the classical case of a solid backward-facing step. Results provide new insights into the interpretation for flux estimation by eddy-covariance and flux gradient methods and for the assessment of surface boundary conditions in turbulence models of the atmospheric boundary layer in complex landscapes and over water bodies affected by canopy wakes. The wind tunnel results indicate that the wake of a forest canopy strongly affects surface momentum flux within a distance of 35–100 times the step or canopy height, and mean turbulence quantities require distances of at least 100 times the canopy height to adjust to the new surface. The near-surface mixing length in the wake exhibits characteristic length scales of canopy flows at the canopy edge, of the flow separation in the near wake and adjusts to surface layer scaling in the far wake. Components of the momentum budget are examined individually to determine the impact of the canopy wake. The results demonstrate why a constant flux layer does not form until far downwind in the wake. An empirical model for surface shear stress distribution from a forest canopy to a clearing or lake is proposed.     

Near-critical free-surface flows: real fluid flow analysis

An open channel flow with a flow depth close to the critical depth is characterised by a curvilinear streamline flow field that results in steady free surface undulations. Near critical flows of practical relevance encompass the undular hydraulic jump when the flow changes from supercritical (F > 1) to subcritical (F < 1), and the undular weir flow over broad-crested weirs where the flow changes from subcritical (F < 1) to supercritical (F > 1). So far these flows were mainly studied based on ideal fluid flow computations, for which the flow is assumed irrotational and, thus, shear forces are absent. While the approach is accurate for critical flow conditions (F = 1) in weir and flumes, near-critical flows involve long distances reaches, and the effect of friction on the flow properties cannot be neglected. In the present study the characteristics of near-critical free-surface flows are reanalysed based on a model accounting for both the streamline curvature and friction effects. Based on the improved model, some better agreement with experimental results is found, thereby highlighting the main frictional features of the flow profiles.