Next Article in Journal
Comparing Topological Partitioning Methods for District Metered Areas in the Water Distribution Network
Next Article in Special Issue
Experimental and Numerical Investigation of Pre-Breaking and Breaking Vorticity within a Plunging Breaker
Previous Article in Journal
Spatiotemporal Variation of Nutrient Concentrations in the Upper Shule River Basin, the Qinghai-Tibetan Plateau, China
Previous Article in Special Issue
Dip Phenomenon in High-Curved Turbulent Flows and Application of Entropy Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Advanced 3D Mapping of Hydrodynamic Parameters for the Analysis of Complex Flow Motions in a Submerged Bedrock Canyon of the Tocantins River, Brazil

1
Federal University of Paraná (UFPR), Post-Graduate Program on Water Resources and Environmental Engineering (PPGERHA), Caixa Postal 19011, Curitiba 81531-990, Brazil
2
Department of Civil Engineering, University of Ottawa, 161 Louis Pasteur St., Ottawa, ON K1N 6N5, Canada
*
Author to whom correspondence should be addressed.
Water 2018, 10(4), 367; https://doi.org/10.3390/w10040367
Submission received: 1 February 2018 / Revised: 13 March 2018 / Accepted: 20 March 2018 / Published: 23 March 2018
(This article belongs to the Special Issue Turbulence in River and Maritime Hydraulics)

Abstract

:
Validation and scaling of sophisticated physical and numerical fluvial hydraulic models to real field conditions are limited by temporal and spatial constraints of field measurement technologies. These limitations increase when analyzing hydraulic properties of complex river forms such as submerged bedrock canyons. The analysis of flow under these conditions has demonstrated non-linear behavior, strong secondary circulation and a high level of turbulence. The objective of this paper is to obtain and analyze turbulence and secondary current information for flow analysis, by evaluating the spatial distribution of bed shear stress and eddy viscosity with acoustic Doppler current profiler (ADCP) measurements. The real field-scale case study was the Lourenço Rock Canyon in Tocantins River, Brazil. A total of 10,531 velocity profiles were measured with an ADCP adapted to a moving boat. The data were post-processed by decomposing the velocity data measurements and calculating the standard deviations. Three different methods are compared for bed shear stress: total kinetic energy (TKE), law of the wall and the depth–slope product. The eddy viscosity was estimated with the Boussinesq approach. The initial data processing confirmed turbulence and flow velocity characteristics described in similar literature findings. The results for the bed shear stress for the canyon area were 8 times greater than the outside and for eddy viscosity, 10 times higher. For the surroundings area, the eddy viscosity showed reasonable results with values around the regular 1.0 m2/s used for rivers. It is concluded that turbulent values can be calculated to better represent physical processes with the intention to improve hydrodynamic numerical models calibrations.

1. Introduction

The understanding of complex flow features, as a consequence of complex geometries or hydrodynamic boundary conditions, is very important for the analysis of transport of substances or sediment in rivers and coastal waters, as well as for navigability analysis in waterway projects. The study of those features generally involves sophisticated physical and numerical modeling approaches [1]. Unfortunately, validation with and scaling to real field conditions are still challenging, due to temporal or spatial constraints of field measurement technologies. Whereas lab or numerical data are available in high frequency and high spatial resolution (e.g., particle image velocimetry (PIV) and laser Doppler anemometry (LDA) [2], field data are usually measured at fixed points or along transects, and in the past were limited to mean velocity measurements (no turbulence quantities). Also, regarding extraordinary geometries, such as deep and narrow bedrock canyons, the conventional modeling laws are untested and are most probably not applicable.
Nowadays, acoustic technologies are more accessible and improved, presenting considerable advances in the quality of the field measurements. Compared to intrusive and punctual devices (current meter), acoustic instruments, such as the acoustic Doppler current profiler (ADCP), have been used in conventional, hydrological river measurements (discharge, depth, velocity), yielding increased accuracy and greater spatial and temporal resolution due to hardware improvements providing better accuracy for positioning, but also having several frequencies combined in the same devices, as well as being robust equipment for field deployment, especially compared to optical (e.g., LSPIV) techniques [3]. Recently, higher spatial resolutions were obtained by mapping individual ADCP velocity profiles over large flow areas, as presented by Flener et al. [4] or Jamieson et al. [5], where the latter included measured sections with 5 and 20 m spacing upstream and downstream of two different wing dikes for three different discharges.
Using the Doppler effect of acoustic waves scattered back from particles within the water column, ADCP is proven to measure velocity, depth, and water discharge accurately with high resolution (cm) [6]. Besides the depth and the mean velocity components, additional data inherent to the acoustic measuring method can be analyzed, such as the backscatter intensity and velocity standard deviation. For example, ADCP data have been used to derive other flow characteristics, such as mixing coefficients [7,8,9], sediment transport rates [5], and moving bed velocity [10]. The flow velocity information has also been used for profile analysis and mapping of integrated or characteristic quantities, such as bed shear stress and depth-averaged velocity [6,11]. Moreover, acoustic backscatter signal and its statistics can be used for a suspended sediment transport analysis [12]. The velocity standard deviation information (velocity’s second-moment statistics), can be explored to obtain turbulence quantities [13,14], which will be the focus of this article. This parameter is especially of interest when analyzing and modeling complex flow situations, where strong and large secondary currents exist (3D effects), creating non-conventional turbulence distributions (non-logarithm velocity profiles), thus limiting resistance law-based approaches or hydrostatic modeling concepts.
Bedrock Rivers can present morphology variations of width and depth over small distances caused by climatic variations, abrasion, cavitation or tectonic movements. Relatively narrow and deep canyons emerge in locations of active bedrock incision, sometimes in association with tectonic fault lines [15]. Narrow canyons have flow characteristics (hydraulic radius, velocities, and turbulent effects) vastly different from less confined reaches, with the presence of vertical accelerations and secondary currents that influence the shape of the velocity profiles. Analyzing the flow in canyons of the Fraser River (Canada), Venditti et al. [16] described the flow as nonlinear and with strong secondary circulation and high levels of turbulence. At the entrance to canyon pools, they observed a plunging core of high velocity flow, which suppressed the maximum velocity below the free surface.
The objective of this paper is to obtain turbulence and secondary current information for flow analysis and model calibration in complex flow environments in a peculiar and morphologically unique submerged bedrock canyon of the Tocantins River in Brazil, a tributary of the Amazon estuary. As the canyon reach is considered an obstruction point for navigation within the important Tocantins waterway, recent studies were undertaken to evaluate and plan river engineering projects for Brazil’s National Transport Infrastructure Department (DNIT).
The spatial distributions of bed shear stress (τ0) and eddy viscosity (υt) were assessed in this river reach, based on extensive ADCP measurements collected from a moving boat to better understand the river dynamics in this submerged bedrock canyon. Using this real field-scale case study, three different approaches were used to calculate bed shear stress: the total kinetic energy (TKE) method, law of the wall method and the depth–slope product method. With the results of the TKE method, the eddy viscosity was estimated with the Boussinesq approach.

2. Materials and Methods

2.1. Bedrock Canyon of Tocantins River

The study area is named Lourenço’s Rock Canyon (Figure 1), a reach of Tocantins River, one of the largest and most important rivers in Brazil. Located in northern Brazil, the Tocantins River is the main river of the Tocantins-Araguaia Basin, the second largest hydrographic basin in the country with 967,059 km2, smaller only than the Amazon River Basin.
The Tocantins-Araguaia waterway is one of the main transportation routes of Brazilian products for export, such as minerals and agricultural commodities. Nowadays, the Lourenço Rock Canyon is considered one of the most significant obstructions to navigation in the Tocantins River. Despite the high depth, the canyon path presents a sharp curve and the available width is not adequate for safe navigation during the drought season and for the commonly used vessels. Moreover, through the flood season, the canyon induces strong shear flows, making navigation more difficult and less safe. According to DNIT [17,18], the Brazilian Government is investing more than U$150 million in engineering projects to improve the navigation channel in the study area.
The submerged bedrock canyon is 6 km long, and its depth can reach 60 m (upstream and downstream of the canyon reach the river has a depth of 7 m). The annual mean flow in the reach is 11,800 m3/s, with high seasonal variations ranging from 1898 to 45,717 m3/s. During the flood season, the river width is around 1.5 km, with the presence of a few rocky outcrops, high velocities and strong turbulent flow (large shear layers and longitudinal vortices). During the drought season, the flow is confined along the canyon area (150 m width), exposing the river bedrock at the reach. Although having smaller diameters, the presence of vortices is constant, even under low flow conditions.
According to the Brazilian Mineral Resources Research Company [19], the river reach is located between two tectonic domains, the Couto Magalhães Formation (part of the Tucuruí Group) and the Bacajá Domain (part of the Transamazonas Province). The river bed is composed of sandstone Rock from the Itapecuru (K12it), Pastos Bons (J2pb) and Grájau (K1g) formations on the undifferentiated fractured basement (Fr). The canyon’s origins can be associated with the Tucuruí Fault, caused by tectonic movements that generated a series of structural features in the lithologies and later migmatites, which also induced a dynamic metamorphic event in the rocks under green-schist facies conditions.

2.2. Equipment and Field Data

In April (flood season) and October (drought season) 2015, several flows and velocity measurements were performed in Lourenço’s Rock Canyon using an acoustic Doppler current profiler (ADCP), model M9, series SN-3754 of Sontek Inc. (San Diego, CA, USA) [20]. The equipment was installed with fixed support on a moving boat, and the installation depth of the transducer was 0.47 m below the water surface. The equipment was coupled to a real-time kinematic (RTK) global positioning system providing sub-decimeter accuracy for the measurement coordinates.
As the boat traverses the sampling area, each vertical depth profile (ensemble) is subdivided into several cells (bins), and the velocity profile for each compass axis (easting, northing, and vertical) is calculated for each of the four inclined beams. An internal compass and tilt sensor provide compensation for boat movements. As only three beams are necessary to compute the three velocity components, additional information can be obtained when using ADCPs with four transducers (similar to the one used here), where the fourth beam provides a redundant measurement of vertical velocity, from which velocity heterogeneities and turbulence variations may be quantified. The ADCP can be operated in 3 different modes: NarrowBand, Pulse to Pulse and Broadband [21]. In NarrowBand mode, the transducers emit, in a short period (e.g., 1 s), several pulses (pings) and the water velocity is calculated based on statistical analysis of the received signals. As a result, for every ensemble profile measurement, a short time series is produced along each axis in each depth cell, with a mean velocity (<u>, <v>, <w>) and a standard deviation (σx, σy, σz). Some ADCP brands supply all the velocity samples from the time series, allowing the user to calculate the statistics; however, the M9 only provides the results of mean velocity and standard deviation, which will be used here.
The ADCP M9 is a double frequency instrument, with four transducers (25° from vertical) at both 1 MHz and 3 MHz. The frequency and thus cell sizes are changed automatically on-the-fly depending on the local water depth, if not configured otherwise. Along the sections in the Lourenço’s Rock Canyon, the shallow margins were usually measured with a frequency of 3 MHz, and the deeper canyon with 1 MHz, both in Narrowband mode. The vertical dimension of the measurement cells (ΔZcel) varied between 0.1–0.5 m for the 3 MHz and 0.5–2.0 m for the 1 MHz. The horizontal sampling dimension is determined as a function of the boat’s movement in time, and the data acquisition frequency (f = 1 s). Also, there is a vertical beam of 0.5 MHz in the middle of the instrument head to improve bed level measurements. The Sontek M9 maximum depth for velocity measurements is 40 m, and some parts of the canyon were more than 50 m deep, thus resulting in larger blanking areas near the bottom. Nonetheless, the bed level was measured throughout the canyon using the vertical beam.
Although measurements are available for both flood and drought season, only the data of the flood situation have been processed so far regarding the methods presented in this manuscript. Thus, all following results are based on the flood situation only, being the more interesting one in terms of the observable features and processes.
The mean discharge during the April 2015 flood flow situation was 18,275 m3/s. Several transversal sections and a longitudinal zig-zag type of measurement were undertaken along Lourenço’s Rock Canyon during a short period (short enough to have discharge changes of less than 3%). The measurement resulted in 10,531 velocity profiles (ADCP ensembles) at different points of the canyon, each having on average about 14 cells, thus resulting in approximately 140,000 3D velocity vectors. Figure 2 shows the location of the spatial measurements. Three Elevation Reference (ER) stations were installed and geo-referenced along the study area using the Global Navigation Satellite System (GNSS). The ADCP’s compass was calibrated three times during the measurement. The water surface elevation measured at each point is presented in Figure 2, and the resulting mean water surface slope (SWS) was 2.36 × 10−4 m/km. A detailed bathymetry measurement of the area was available from previous measurements in 2012, with both single and multi-beam echosounder equipment. A DTM (Digital Terrain Model) was created (Figure 1), combining the bathymetry data with SRTM (Shuttle Radar Topographic Mission) data to complete the margin information.

2.3. Initial ADCP Data Processing

The equipment’s software Riversurveyor [20] was used to export all the velocity profiles into Matlab [22] file formats. The data coordinate system was ENU (East-North-Up), the depth reference used was the vertical beam (VB) and the boat velocity used as reference the bottom track (BT) feature. Subsequent post-processing was done using a recently programmed toolbox, VelMap [23]. The toolbox can decompose the velocity data measurements by rotation as a function of an angle chosen by the user. For the study area, an angle equal to 337° clockwise from North (angle = 0°), was considered the longitudinal flow direction (x-axis), and 247° as the transversal (y-axis. This fixed coordinate system has been chosen here due to the large heterogeneities in the flow and thus difficulties in defining a consistent main flow direction. The chosen angle, however, corresponds to the overall mean velocity direction. Moreover, the toolbox uses additional parameters (SNR—signal to noise ratio, STD—standard deviation) for further processing, following instrument specifications.
Many factors affect the uncertainty of ADCP measurements, such as instrument noise and configuration, velocity fluctuations due to flow turbulence, bottom track errors, compass errors and GPS errors [6]. According to Sontek [24], the standard deviation of the velocity data can be attributed to installation and calibration problems, acoustic noise from the equipment, and velocity fluctuations from the local turbulence. Following [25,26], installation and calibration errors were considered negligible, and the standard deviation for each velocity component decomposed as:
σ2 = σt2 + σn2
where σ is the standard deviation from the velocity data (m/s); σt is the velocity fluctuation from the local turbulence (m/s); σn is the acoustic noise effect from the equipment on the standard deviation (m/s). The acoustic noise is a result of the physical process by which sound waves are dispersed by the solid particles in the water. To estimate the turbulent part of the standard deviation, previous approaches calculated the acoustic noise using different assumptions. Vachtman and Laronne [27] considered the acoustic noise equal to the minimum variance of calculated water velocities for all depth cells (σn2 = σ2min). Another approach was presented by Rennie and Church [6], who applied Equation (1), but used the average ADCP error velocity (σerror), instead of the acoustic noise. The average ADCP error velocity considered other uncertainties and can be calculated by the velocity difference (Vd) from 2 vertical velocity measurements and the horizontal velocity error (Ve). According to Sontek [24], although being a random effect, it can be assumed that acoustic noise follows a Gaussian distribution. The magnitude of the noise can vary according to the equipment frequency (F), the vertical dimension of the cell (ΔZcel) and the number of pings (N) emitted in the period measured per cell. For a standard 3D ADCP measurement, the acoustic noise can be estimated with Equation (2).
σn = 235 (F ΔZcel N0.5)−1
Developed by the same manufacturer, Equation (2) was defined for older ADCP models than the M9. Conferring with Sontek [28], there is not yet a specific noise equation for the M9. It is herein assumed that the results from Equation (2) can be used as the best approximation of the noise from the ADCP data.
All results can be visualized as velocity magnitudes along a cross-section for each mean velocity component (<u>, <v>, <w>), related velocity fluctuations from the local turbulence (σtx, σty, σtz) and computed flow parameters (turbulent kinetic energy k and shear stress τ), as shown in Figures 5–7, and explained in the following section. The velocity distributions show the part of the canyon with highest velocities and highest longitudinal velocities occurring in the deeper part of the channel. These effects follow the observations by Venditti et al. [16], where similar flow features were observed for a bedrock canyon of the Fraser River.
The strongest transversal velocity fluctuations occur in the canyon section. The variations of transversal and vertical velocity components indicate large and strong secondary currents induced by the canyon walls. The secondary currents in the horizontal plane were even visible as strong (boat was shaking), and large (several meters) surface vortices (Figure 3) also manifested as strong transversal velocity fluctuations. Normally, these large vortices are difficult to observe in natural rivers. However, these vortices can frequently be viewed along the canyon section during the flood season, reinforcing the river reach peculiarity.

3. Theory and Analytical Methods

This section briefly reviews the different methods applied to calculate the local bed shear stress and eddy viscosity distributions, including discussion of details and assumptions.
ADCPs operated from moving vessels have been used in previous studies to estimate spatial flow characteristics. Flener et al. [4], Rennie and Church [6], Rennie and Millar [10], Sime et al. [11], Szupianyet et al. [29], Guerrero and Lamberti [30], are reference studies that mapped some flow characteristics, such as shear velocity, shear stress, suspended and bedload sediment transport. However, in some river structures, such as meanders, bifurcations and submerged canyons, these characteristics are more difficult to estimate because of the high turbulence levels and the nonlinear velocity profiles along the river area [31]. Three different methods to calculate the bed shear stress were compared here to evaluate the flow characteristic of the measured reach.

3.1. TKE Method

All the velocities can be split into mean components (<u>, <v>, <w>) and fluctuating components (u’, v’, w’) using Reynolds decomposition, where the parentheses denote the time-averaged value, and the prime symbol the turbulent fluctuations. Results of the TKE method relate the intensity of velocity fluctuations (u’, v’, w’) to turbulent kinetic energy (k), as described in Equation (4). Based on the hypothesis that the Reynolds stress and the variances in the region near the bottom are practically constant and proportional to the shear velocity (u*), this method has been applied to ADCP data in studies by Lu and Lueck [21], Nystrom et al. [32], Hurther and Lemmin [33], Vachtman and Laronne [27], Liu and Wu [34]. Therefore, the shear stress was calculated as:
τ0 = C0 TKE = C0 ρ k
k = (<u’2> + <v’2> + <w’2>)/2
where C0 is a dimensionless constant (i.e., C0 = 0.18–0.21) [13]; ρ is the water density (i.e., ρ = 1 g/cm3 at 15 °C). Using the ADCP in Narrowband mode, it is possible to combine the Standard Deviation (STD, Equation (1)) from the velocity data (for each component), and the Reynolds Decomposition to determine the temporal root mean square of the velocity fluctuations (<u’2>).
σtx = (N−1 ∑(u − <u>)2)0.5 = (<(u − <u>)2>)0.5 = (<u’2>)0.5
By rewriting Equation (4) with the other equations and directions, the shear stress for each cell (τ0-TKEcel) can be calculated using the velocity STD components of the ADCP data.
τ0-TKEcel = C0 ρ (σtx2 + σty2 + σtz2)/2
The bed shear stress (τ0-TKE), calculated by the TKE method, can be determined by the integration of the shear stress calculated in each cell (τ0-TKEcel) of the vertical measurement (ensemble). For the vertical average the variance sum law was applied, where the covariance of the velocity time series between each cell should be calculated; however without the samples, this study considered that the sum of the covariance among the ensemble is equal to zero. Togneri et al. [35] used the same consideration to compare the ADCP observations with a 3D model simulation. Besides, Lu and Lueck [21] calculated the distribution (histogram) of the covariance in velocity profiles from an ADCP fixed to the bottom of a tidal channel at a depth of 30 m, during 20-min intervals in 2 different flow events. The results showed a balance between positive and negative covariance, with a mean near to zero.
The numbers of samples in each cell measurement is another relevant point. Depending on the size and duration of the turbulent features to be measured, longer averaging times are required, requiring either anchored boat measurements [11] or slow boat movements. In the present study, the slow boat moving method with large spatial distribution was chosen to capture large-scale turbulent variations.

3.2. Law of the Wall Method

The log law is often used to estimate shear stress in open channels and rivers, as presented in Keulegan [36]. Based on the assumption that the velocity profile in open channel flows has a logarithmic shape, the “Law of the Wall” can be described by the von Karman–Prandtl Equation.
U/u* = 1/κ ln (z/z0)
where κ is the von Karman constant (i.e., κ = 0.40); U is the velocity at height z above the bed; u* is the shear velocity, and z0 is the roughness length. For natural rivers, considered as a hydraulically rough flow, Equation (7) can be described using the characteristic roughness length-scale (ks).
U = (u*/κ) ln(z) + (u*/κ) ln(30/ks)
The shear velocity and the roughness length-scale can be estimated using the velocity profile measurement by the ADCP and fitting the data to the following equation.
U = m ln(z) + b
where m = (u*/κ); b = (u*/κ) ln(30/ks) = m ln(30/ks).
The relationship between shear velocity and bed shear stress (τ0-log), by the log-law, is given by:
τ0-log = ρ (u*)2
According to Pope et al. [37], this method is only valid for steady, uniform, one-dimensional flow and its application to natural streams does have some limitations. Afzalimehr and Rennie [38] comment that the method does not hold in general and can only be considered an approximate representation of experimental data. For measurements at a single point, long time-series data are required for proper averaging. In the present study, the velocity distributions do not seem to follow a logarithmic profile, thus it was expected that this method is not applicable in detail.
Figure 4 shows a typical vertical profile of longitudinal velocities as measured in the canyon. The profile shows the features observed in the cross-section visualizations, with high velocities in deep channel regions. Unfortunately, the depth limitations of the equipment resulted in a large blanking region close to the bed. Regardless, the profile does not visually follow a log law. To minimize the velocity non-linearity problem, each velocity profile was separated in 2 parts: a non-uniform (near the water surface) and a log part (starting at the maximum velocity point in the profile). Using this assumption, only the data below the maximum velocity point were considered for the fitting processes. Several fits also resulted in negative velocities close to the bed, limiting this method further, thus not following it for the subsequent calculations.

3.3. Depth–Slope Product Method

The mean boundary shear stress equation considers a steady and uniform flow and is generally used in sediment transport studies to estimate the channel stability. In large rivers, with small slopes, it is possible to assume that the bed slope (S0) is parallel to the energy slope (Se) and the water level slope (Sws), and the hydraulic radius (Rh) nearly equals the transect-averaged water depth (<h>). The application of these assumptions within the mean boundary shear stress equation is often referred to as the “depth-slope product”, as presented in:
τ0-ds = ρ g Rh Se = ρ g <h> S0
where τ0-ds is bed shear stress, by the depth-slope product method; g is the gravitational acceleration (i.e., g = 9.81 m2/s); <h> is the transect-averaged water depth (m); Se is the local energy slope. The transect-averaged water depth is a standard result of the ADCP measurement, and the bed slope (S0) was considered equal to the water surface slope (Sws) measured in the field. Because of the depth difference, the ADPC cross-section data were subdivided, selecting a mean hydraulic radius of the canyon area and another value for the rest of the cross-section. Unfortunately, the simplifications presented in the depth-slope product method (steady and uniform flow) are an unrealistic representation of the spatial data condition of the measurements in Lourenço’s Rock Canyon. Nonetheless, Sime et al. [11] considered it worthwhile to compare the method with other results because of its high usability. However, it has not been followed further in this manuscript.

3.4. Boussinesq’s Hypothesis

The concept adopted in Reynolds averaging of the Navier–Stokes (RANS) created a non-linear term of the Navier-Stokes equation, called the “Reynolds stress”. The term represents the influence of the turbulent fluctuations in the momentum distribution in the mean flow [39]. The coordinate system of the Reynolds stress is presented in Table 1 [40].
With the extra term, the equations of motion for turbulent flow can never form a closed system in which the number of equations is equal to the number of unknowns. This is called the turbulence closure problem and can be solved by eddy-viscosity models (EVM), which use Boussinesq’s hypothesis.
Boussinesq’s hypothesis was formulated in 1877 and is based on the analogous relationship between the viscous stresses of the laminar regime and the turbulent one, and it is assumed that the turbulent stresses are proportional to the mean velocity gradient of the flow in a similar manner as the viscous stresses (Newton’s Viscosity Law) [41]. This proportionality factor was named the eddy viscosity (υt) and is not a physical property but is dependent on the current flow condition and turbulence levels. The correlation between the turbulent shear stress of the flow and the eddy viscosity (for the x and z coordinate axes) is described by:
τ’xz/ρ = − <u’w’> = υt (∂<u>/∂z + ∂<w>/∂x) − 2/3 δxz k
where τ’xz is the turbulent shear stress (x and z-axes), in kg/ms2; <u’w’> is the Reynolds Stress component (x and z-axes); υt is the eddy viscosity, in m2/s; k is the turbulent kinetic energy, in m2/s2; δ is the dimensionless Kronecker delta.
In natural rivers, the shear stress is the force of moving water against the bed, and the turbulent effects are the principal component (τ ≈ τ’). According to Vachtman and Laronne [32], from the computed turbulent kinetic energy (k), based on the average of the velocity variances for a given depth above the bed, the tridimensional turbulent shear stress (τ’) can be estimated using [13,42]:
τ’ = C1 TKE
where C1 is a dimensionless constant (i.e., C1 ≈ 0.19–0.21). The specific value of C1 may not apply at all levels within the water column because it requires confirmation before using the method universally [13]. Nonetheless, this paper considered the turbulent shear stress for the water column (τ’ensemble) equal to bed shear stress, by the TKE method (C1 = C0 = 0.20).
The velocity gradients and a 3D vorticity field were calculated using the ADCP velocity data. However, the transversal distance between the points was on average 1.5 m, and the longitudinal reached 200 m (distance between transects). Because of the longitudinal distance between the ADCP data, it is reasonable to set the terms with longitudinal distance to zero. Considering the vertical average and applying Boussinesq’s hypothesis, the turbulence eddy viscosity (for each ensemble) can be calculated as:
υt ensemble = τ’ensemble/(ρ ω ensemble)
where υt ensemble is the turbulence eddy viscosity, in m2 s−1, τ’ensemble is shear stress, by the TKE method (Nm−2); ω ensemble is the vertical average of the 3D vorticity along the cross-section (s−1).
Numerical models are generally used to complement field data for further analysis. The most accessible and common models are Reynolds averaged Navier–Stokes (RANS) models, which require either the modeling or the definition of values for eddy viscosity. Usually, constant values are used, and simulations are calibrated using varying roughness values. However, under complex flow situations, the calibration process can result in non-realistic roughness coefficients, thus non-realistic results for simulations beyond the calibrated regime. Under such conditions, varying eddy viscosity values may provide better flow representations, as adopted in Cheng et al. [39].

4. Results

The previously described methods have been applied to all measured velocity profiles to obtain shear stress and eddy viscosity distributions for the whole river section. First, a cross-sectional result will be shown, to illustrate the differences between shallow, regular river sections, compared to the deep canyon section. Those are followed by mapped results for the whole region.

4.1. Single Cross-Section

The section shown in Figure 5 is approximately 1400 m wide and has depths varying between 5 m in the shallow part down to 55 m in the canyon. The canyon itself is only 150 m wide in most of the sections. Figure 6A–F show only a part of the cross-section, between of the distances 700 to 1100 m, focusing on the canyon area.
Applying the TKE method, the velocity data and standard deviations were used to compute the turbulent kinetic energy (k) in all measured cells. For situations where the instrument noise was higher than the turbulence fluctuations, the k value was considered equal to zero. Afterward, all values were interpolated (natural neighbor) on an equidistant vertical grid, to eliminate changing vertical cell sizes due to frequency shifts and to extrapolate the blanking area. The value used was the ADCP maximum vertical cell size in the canyon (2 m). The horizontal grid sizes were maintained (i.e., varying according to boat speed). Afterward, shear stress and eddy viscosity were computed for each cell using Equations (5) and (13). The shear stress results for the TKE method in the canyon area are shown in Figure 7.
For the horizontal analysis, a comparison between the bed shear stress of the three methods is presented in Figure 8. For the TKE method, the vertical-average was used to represent each ensemble.
The velocity information was processed further to identify secondary currents by computing vorticity maps over the interpolated cross-sectional velocity distributions (Figure 9). Those showed large-scale structures at the canyon walls where surface vortices were also mostly present. The large vortices most probably are induced by the shear effects of the strongly inclined canyon walls. Those secondary currents seem to be responsible for pushing the water at high velocities from the surface down into deeper regions of the canyon, creating non-logarithm velocity profiles.
Applying Boussinesq’s hypothesis, the spatial eddy viscosity distribution indicates a strong turbulent effect within the canyon (Figure 10). These high values are different from the homogeneous values obtained for the rest of the cross-section.

4.2. Spatial Distribution

The previously described processing was done for all ADCP profiles over the whole river section. The mapping was done on a regular rectangular grid with cells of 10 m × 30 m (transversal × longitudinal), oriented in the mean flow direction (337° from North). Combining the water level measurements with the bathymetry, the whole wetted region was determined by cutting a water level surface with the digital terrain model of the bathymetry. The data were then extrapolated to regions not measured. The details of the spatial interpolation are presented in Appendix A. The principal results are shown in Figure 11.

5. Discussion and Conclusions

The measured velocity profiles in the submerged canyon showed irregular flow features, characterized by high flow velocities in the deeper part of the canyon, and large secondary currents developing at the canyon side walls. Analyzing the sections velocity field, the higher vorticity values happened around the longitudinal axis. The initial processing of the data showed high velocity fluctuations in the vertical and transversal axis, representing the strong presence of turbulence along the canyon.
The turbulent quantities associated with the ADCP data indicated a large internal shear occurring within the canyon. Those effects have been summarized by computing the bed shear stress and eddy viscosity distributions and showing distinct regions with significant differences. The shear stress values from the depth–slope product method were considered lower than the real flow conditions, reaffirming the method limitations for high turbulence flow. The assumption used in the law of wall method assisted the curve-fitting process for non-logarithm velocity profiles. However, because of the blank area near the bottom, some curves were fitted considering negative velocities near the bottom, resulting in higher values for the bed shear stress along the canyon and some outside parts. The results from the TKE method were considered most realistic for the canyon area, with shear stress values similar to those presented in other studies with similar conditions [13]. The most significant shear stress difference happened in the canyon upstream, with values 8 times greater than outside.
The eddy viscosity distributions show reasonable results with values around the regular 1.0 m2/s used for rivers. In most of the areas of the stretch, the values were between 0.8 and 1.8 m2/s. In the canyon, analyzing each section, some points reached values higher than 20 m2/s. For spatial analysis, the higher value considered was 10 times higher than the other parts of the stretch.
The spatial interpolation was able to show the horizontal variation of shear-stress and eddy viscosity distributions. However, due to the longitudinal distance between the cross-sections, it was not possible to represent the velocity gradients at the x-direction (∂<v>/∂x and ∂<w>/∂x).
The turbulent quantified values can be used with the intention to improve hydrodynamic numerical model calibrations (turbulence models) and better representation of physical processes. However, it is important to highlight the assumptions and the methodology limitations described in this paper, such as the short velocity time series and the interpolation.
To improve future field measurements, we recommend adding extra points with fixed ADCP measurements along each section to evaluate the shear stress results with an extensive time series. To improve the vertical analysis, equipment with a different frequency to check the velocity near the bottom should be used. For longitudinal analysis, the maximum distance should be only 2 times the transversal distance data.
For future analysis, a sensitivity grid test is recommended, comparing different cell sizes or applying a non-regular grid. In addition, different interpolation methods can be tested. The next steps will be to compare the results of the flood and drought season (already measured, but not processed) and to evaluate their applicability to hydrodynamic models.

Acknowledgments

The authors acknowledge the financial support for the field studies from Brazil’s National Transport Infrastructure Department (DNIT), and the continuous support from the Transport Infrastructure Technical Institute (ITTI) of the Federal University of Paraná (UFPR), and Coordination for the Improvement of Higher Education Personnel (CAPES) regarding scholarships. We also acknowledge the participating professor provided by Brazil’s National Council for Scientific and Technological Development (CNPq), as well as additional funding from the Araucaria Foundation for the travel expenses.

Author Contributions

G.T. and T.B. proposed and elaborated the presented idea and performed the field measurements. C.D.R. encouraged G.T. to investigate the law of the wall method and supervised the findings of this work, and provided scientific background information; G.T. performed the data post-processing and the methods; H.G. performed the spatial interpolation. G.T. wrote the manuscript with support from T.B., C.D.R. and H.G.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

For the spatial interpolation of the flow’s parameters, a regular rectangular grid was developed with cells of 10 m × 30 m (transversal × longitudinal). Considering the intrinsic anisotropic patterns of river systems, the grid follows the general longitudinal flow direction of the stretch (337° from North). However, the canyon’s curvature varies along the stretch, resulting in different longitudinal angles of the main flow compared to the canyon’s thalweg. This characteristic results in difficult interpolation inside the canyon with coarse point density.
A similar situation happens with spatial analysis in meandering rivers, where the spatial grid does not follow the longitudinal flow direction. For a proper interpolation of a river’s main channel longitudinal characteristics, Merwade et al. [43] suggested a coordinate transformation from the typical (E, N) coordinate system to one that follows the river’s main channel. The transformation consists of attributing to every point a (s, d) coordinate, where s is the distance of a cross-section to an arbitrary point that starts at the river’s center path or thalweg, and d the distance from a point to the intersection of the central path with the cross-section to which it belongs. The coordinate transformation can be applied to ADCP sample points using the AnyForm QGIS plugin [44]. The algorithm uses the basic concept of the method’s implementation described by Guarneri [45] and described by the Equation (A1).
[s, d] = [Sb + |PB| cos (α), ψ (|PB| sin (α)]
where [] represents a matrix and bold font the vector system; Sb is the longitudinal distance of B from the arbitrary beginning; B is the planimetric coordinate [EB, NB] of the upstream point; P is the planimetric coordinate [EP, NP] of the surveyed point; α is the angle between the vectors |CP| and |CB|; C is the planimetric coordinate [EC, NC] of the downstream point; ψ is a parameter that equals −1 for points on the right side and 1 for points to the left.
For this paper, the spatial interpolation was divided into 2 steps. First, the interpolation occurs over the deepest parts of the canyon, using the samples with bathymetry ADCP data less than 63 m above sea level. The grid was linearized by the canyon’s thalweg (s-line), and the samples were interpolated with kriging interpolation. The results were re-transformed to a regular grid and combined with the other ADCP points as an input for the second step. Figure A1 presents the activity order from the first step. The second step is a simple kriging interpolation with an angle anisotropy of 337° from North and an aperture angle tolerance of 20°. Figure A2 shows the input and output from the second step for the ADCP bathymetry data. The same interpolation procedure was applied to all variables presented in the paper.
Figure A1. (a) Bathymetry data (below 63 m) in a regular coordinate system grid (E-N); (b) Bathymetry data transformed in the s-d grid (by s-line direction); (c) Kriging interpolation using the s-d grid; (d) Interpolated bathymetry grid along the canyon central path in a regular grid.
Figure A1. (a) Bathymetry data (below 63 m) in a regular coordinate system grid (E-N); (b) Bathymetry data transformed in the s-d grid (by s-line direction); (c) Kriging interpolation using the s-d grid; (d) Interpolated bathymetry grid along the canyon central path in a regular grid.
Water 10 00367 g0a1
Figure A2. (a) All the Bathymetry data and the results from the first step at the coordinate system grid (E-N) were used as an input for the second step; (b) Bathymetry results from kriging interpolation after the second step.
Figure A2. (a) All the Bathymetry data and the results from the first step at the coordinate system grid (E-N) were used as an input for the second step; (b) Bathymetry results from kriging interpolation after the second step.
Water 10 00367 g0a2

References

  1. Constantinescu, G.; Kashyap, S.; Tokyay, T.; Rennie, C.D.; Townsend, R.D. Hydrodynamic processes and sediment erosion mechanisms in an open channel bend of strong curvature with deformed bathymetry. J. Geophys. Res. Earth Surf. 2013, 118, 480–496. [Google Scholar] [CrossRef]
  2. Aberle, J.; Rennie, C.D.; Admiraal, D.; Muste, M. Experimental Hydraulics: Methods, Instrumentation, Data Processing and Management: Volume II: Instrumentation and Measurement Techniques; CRC Press: Boca Raton, FL, USA, 2017; ISBN 9781351658034. [Google Scholar]
  3. Pilechi, A.; Mohammadian, A.; Rennie, C.D.; Zhu, D.Z. Efficient method for coupling field data and numerical modeling for the estimation of transverse mixing coefficients in meandering rivers. J. Hydraul. Eng. 2016, 142, 0001116. [Google Scholar] [CrossRef]
  4. Flener, C.; Wang, Y.; Laamanen, L.; Kasvi, E.; Vesakoski, J.M.; Alho, P. Empirical modeling of spatial 3D flow characteristics using a remote-controlled ADCP system: Monitoring a spring flood. Water 2015, 7, 217–247. [Google Scholar] [CrossRef]
  5. Jamieson, E.C.; Rennie, C.D.; Jacobson, R.B.; Townsend, R.D. 3-D flow and scour near a submerged wing dike: ADCP measurements on the Missouri River. Water Resour. Res. 2011, 47, W07544. [Google Scholar] [CrossRef]
  6. Rennie, C.D.; Church, M. Mapping spatial distributions and uncertainty of water and sediment flux in a large gravel bed river reach using an acoustic Doppler current profiler. J. Geophys. Res. 2010, 115, F03035. [Google Scholar] [CrossRef]
  7. Carr, M.L.; Rehmann, C.R. Measuring the dispersion coefficient with acoustic doppler current profilers. J. Hydraul. Eng. 2007, 133. [Google Scholar] [CrossRef]
  8. Shen, C.; Niu, J.; Anderson, E.J.; Phanikumar, M.S. Estimating longitudinal dispersion in rivers using acoustic Doppler current profilers. Adv. Water Resour. 2010, 33, 615–623. [Google Scholar] [CrossRef]
  9. Kim, D. Assessment of longitudinal dispersion coefficients using acoustic Doppler current profilers in large river. J. Hydro-Environ. Res. 2012, 6, 29–39. [Google Scholar] [CrossRef]
  10. Rennie, C.D.; Millar, R.G. Measurement of the spatial distribution of fluvial bedload transport velocity in both sand and gravel. Earth Surf. Process. Landf. 2004, 29, 1173–1193. [Google Scholar] [CrossRef]
  11. Sime, L.C.; Ferguson, R.I.; Church, M. Estimating shear stress from moving boat acoustic Doppler velocity measurements in a large gravel bed river. Water Resour. Res. 2007, 43, W03418. [Google Scholar] [CrossRef]
  12. Wall, G.R.; Nystrom, E.A.; Litten, S. Use of an ADCP to compute suspended-sediment discharge in the tidal Hudson River, New York. U.S. Geological Survey Scientific Investigations Report. 2006. Available online: https://pubs.usgs.gov/sir/2006/5055/ (accessed on 8th July 2017).
  13. Kim, S.C.; Friedrichs, C.T.; Maa, J.P.; Wright, L.D. Estimating bottom stress in tidal boundary layer from acoustic Doppler velocimeter data. J. Hydraul. Eng. 2000, 126, 399–406. [Google Scholar] [CrossRef]
  14. De Serio, F.; Mossa, M. Analysis of mean velocity and turbulence measurements with ADCPs. Adv. Water Resour. 2015, 81, 172–185. [Google Scholar] [CrossRef]
  15. Hancok, G.S.; Anderson, R.S.; Whipple, K.X. Beyond power: Bedrock river incision process and form. In Rivers Over Rock: Fluvial Processes in Bedrock Channels; Tinkler, K., Wohl, E.E., Eds.; American Geophysical Union: Washington, DC, USA, 1998; Volume 3, pp. 35–60. [Google Scholar]
  16. Venditti, J.G.; Rennie, C.D.; Bomhof, J.; Bradley, R.W.; Little, M.; Church, M. Flow in bedrock canyons. Nature 2014, 513, 534–537. [Google Scholar] [CrossRef] [PubMed]
  17. DNIT. Licitação Para Derrocamento do Pedral do Lourenço, na Hidrovia do Tocantins. 2016. Available online: http://www1.dnit.gov.br/editais/consulta/resumo.asp?NUMIDEdital=5682 (accessed on 16th October 2017). (In Portuguese)
  18. Tomas, G.P.; Werner, J.; Ratton, P.; Bleninger, T.; Godoy, P. Modelagem Hidrodinâmica Para o Derrocamento do Rio Tocantins—Pedral do Lourenço; XXI Simpósio Brasileiro de Recursos Hídricos; Associação Brasileira de Recursos Hídricos: Brasília, Brazil, 2015. (In Portuguese) [Google Scholar]
  19. Content Protection for Recordable Media (CPRM). Mapas Geológicos Estaduais—Pará. Companhia de Pesquisa de Recursos Minerais; GEOBANK, 2012. Available online: http://geobank.cprm.gov.br (accessed on 25th January 2016). (In Portuguese)
  20. SonTek. RiverSurveyor S5/M9 System Manual: Firmware Version 1.0; SonTek Inc.: San Diego, CA, USA, 2010. [Google Scholar]
  21. Lu, Y.; Lueck, R.G. Using a broadband ADCP in a tidal channel. Part II: Turbulence. J. Atmos. Ocean. Technol. 1999, 16, 1568–1579. [Google Scholar] [CrossRef]
  22. MATLAB. Matlab and Statistics Toolbox Release 2016a; The MathWorks, Inc.: Natick, MA, USA, 2016. [Google Scholar]
  23. Tomas, G.P.; Hopker, R.; Frigo, A.L.; Bleninger, T. Velocity Mapping Toolbox for Sontek M9 ADCP Data. River Flow; Taylor & Francis Group: Iowa City, IA, USA, 2016; pp. 1571–1576. [Google Scholar]
  24. SonTek. ADP Introductory Documentation. Acoustic Doppler Profiler Technical Documentation; SonTek Inc.: San Diego, CA, USA, 2000. [Google Scholar]
  25. Rennie, C.D. Uncertainty of ADCP spatial velocity distributions. In Proceedings of the 6th International Symposium on Ultrasonic Doppler Methods for Fluid Mechanics and Fluid Engineering, Prague, Czech Republic, 9–11 September; 2008. [Google Scholar]
  26. Greene, A.D.; Hendricks, P.J.; Gregg, M.C. Using an ADCP to Estimate Turbulent Kinetic Energy Dissipation Rate in Sheltered Coastal Waters. J. Atmos. Oceanic Technol. 2015, 32, 318–333. [Google Scholar] [CrossRef]
  27. Vachtman, D.; Laronne, J.B. Flow structure over bed irregularities in a straight cohesive open channel. J. Hydraul. Eng. 2011, 137, 1335–1346. [Google Scholar] [CrossRef]
  28. SonTek; (San Diego, CA, USA). Personal communication, 2017.
  29. Szupiany, R.N.; Amsler, N.L.; Best, J.L.; Parsons, D.R. Comparison of fixed- and moving -vessel flow measurements with an aDp in a large river. J. Hydraul. Eng. 2007, 133, 1299–1309. [Google Scholar] [CrossRef]
  30. Guerrero, M.; Lamberti, A. Flow field and morphology mapping using ADCP and multibeam techniques: Survey in the Po River. J. Hydraul. Eng. 2011, 137. [Google Scholar] [CrossRef]
  31. Muste, M.; Yu, K.; Spasojevic, M. Practical aspects of ADCP data use for quantification of mean river flow characteristics: Part 1. Moving-vessel measurements. Flow Meas. Instrum. 2004, 15, 1–16. [Google Scholar] [CrossRef]
  32. Nystrom, E.A.; Oberg, K.A.; Rehmann, C.R. Measurement of turbulence with acoustic Doppler current profilers—Sources of error and laboratory results. In Hydraulic Measurements and Experimental Methods, Proceedings of the Hydraulic Measurements and Experimental Methods 2002, Estes Park, CO, USA, 28 July–1 August 2002; American Society of Civil Engineers (ASCE): Reston, VA, USA; pp. 346–355.
  33. Hurther, D.; Lemmin, U. A correction method for turbulence measurements with a 3D acoustic Doppler velocity profiler. J. Atmos. Ocean. Technol. 2001, 18, 446–458. [Google Scholar] [CrossRef]
  34. Liu, H.; Wu, J. Estimation of bed shear stresses in the pearl river estuary. China Ocean Eng. 2015, 29, 133–142. [Google Scholar] [CrossRef]
  35. Togneri, M.; Lewis, M.; Neill, S.; Masters, I. Comparison of ADCP observations and 3D model simulations of turbulence at a tidal energy site. Renew. Energy. 2017, 114, 273–282. [Google Scholar] [CrossRef]
  36. Keulegan, G.H. Laws of turbulent flow in open channels. J. Res. Natl. Bur. Stand. 1938, 21, 707–741. [Google Scholar] [CrossRef]
  37. Pope, N.D.; Widdows, J.; Brinsley, M.D. Estimation of bed shear stress using the turbulent kinetic energy approach—A comparison of annular flume and field data. Cont. Shelf Res. 2006, 26, 959–970. [Google Scholar] [CrossRef]
  38. Afzalimehr, H.; Rennie, C.D. Determination of bed shear stress in gravel-bed rivers using boundary-layer parameters. Hydrol. Sci. J. 2009, 54, 147–159. [Google Scholar] [CrossRef]
  39. Cheng, R.T.; Ling, C.H.; Gartner, J.W. Direct Measurements of Turbulence Properties by a BB-ADCP in Bottom Boundary Layer. In Interactions between Estuaries, Coastal Seas and Shelfs Seas; TERRAPUB: Tokyo, Japan, 2000; pp. 37–55. [Google Scholar]
  40. Jamieson, E.C.; Rennie, C.D.; Townsend, R.D. Turbulence and vorticity in a laboratory channel bend at equilibrium clear-water scour with and without stream barbs. J. Hydrol. Eng. 2013, 139, 259–268. [Google Scholar] [CrossRef]
  41. Schmitt, F.G. About Boussinesq’s turbulent viscosity hypothesis: Historical remarks and a direct evaluation of its validity. Comptes Rendus Mec. 2007, 335, 617–627. [Google Scholar] [CrossRef]
  42. Stapleton, K.R.; Huntley, D.A. Seabed stress determination using the inertial dissipation method and the turbulent kinetic energy method. Earth Surf. Process. Landf. 1995, 20, 807–815. [Google Scholar] [CrossRef]
  43. Merwade, V.; Maidment, D.; Hodges, B. Geospatial representation of river channels. J. Hydrol. Eng. 2005, 10, 243–251. [Google Scholar] [CrossRef]
  44. Guarneri, H. Methods for the Analysis and Design of Inland Waterways and Application for the Paraguay River Waterway. Master’s Thesis, Universidade Federal do Paraná, Curitiba, Brazil, March 2017. [Google Scholar]
  45. Guarneri, H. AnyForm—Anisotropic Transformation—QGIS Plugin. GitHub Repos. 2017. Available online: https://github.com/henriqueguarneri/qgis-anyform-plugin (accessed on 21 November 2017).
Figure 1. (A) Geographic location of the Tocantins River and the canyon location at Lourenço’s Rock; (B) Bathymetrical map (above sea level); (C) Image looking downstream during low flow condition; (D) Image looking upstream during low flow conditions and indicating the navigation channel; (E) 3D bathymetrical visualization (vertical exaggeration = 15) overlaid on a satellite image and (C,D) viewing directions.
Figure 1. (A) Geographic location of the Tocantins River and the canyon location at Lourenço’s Rock; (B) Bathymetrical map (above sea level); (C) Image looking downstream during low flow condition; (D) Image looking upstream during low flow conditions and indicating the navigation channel; (E) 3D bathymetrical visualization (vertical exaggeration = 15) overlaid on a satellite image and (C,D) viewing directions.
Water 10 00367 g001
Figure 2. Measurements in Lourenço’s Rock Canyon; (A) ADCP trajectories with colored depth-averaged velocity vectors (April 2015); (B) River banks during high flow are shown in blue, and the location of the three gaging stations with measured values of the water surface elevation (April 2015) is presented.
Figure 2. Measurements in Lourenço’s Rock Canyon; (A) ADCP trajectories with colored depth-averaged velocity vectors (April 2015); (B) River banks during high flow are shown in blue, and the location of the three gaging stations with measured values of the water surface elevation (April 2015) is presented.
Water 10 00367 g002
Figure 3. Large vortices in the canyon shear-layer where steep bathymetrical gradients exist; (A) Flood event—April 2015; (B) Drought event—October 2015.
Figure 3. Large vortices in the canyon shear-layer where steep bathymetrical gradients exist; (A) Flood event—April 2015; (B) Drought event—October 2015.
Water 10 00367 g003
Figure 4. Measured vertical velocity profile of longitudinal velocity (+), and fitted logarithmic curves using a different number of measured points.
Figure 4. Measured vertical velocity profile of longitudinal velocity (+), and fitted logarithmic curves using a different number of measured points.
Water 10 00367 g004
Figure 5. Flow velocity profiles over a full cross-section (3D Velocity magnitude) (April 2015).
Figure 5. Flow velocity profiles over a full cross-section (3D Velocity magnitude) (April 2015).
Water 10 00367 g005
Figure 6. Flow velocity profiles over the canyon part (April 2015); (AC) Magnitude of mean velocity components <u> (longitudinal), <v> (transversal), <w> (vertical) over cross-section; (DF) Magnitude of turbulent velocity fluctuations.
Figure 6. Flow velocity profiles over the canyon part (April 2015); (AC) Magnitude of mean velocity components <u> (longitudinal), <v> (transversal), <w> (vertical) over cross-section; (DF) Magnitude of turbulent velocity fluctuations.
Water 10 00367 g006
Figure 7. Turbulent kinetic energy (k) before (A) and after (B) vertical interpolation; (C) Shear-stress distribution using the TKE method.
Figure 7. Turbulent kinetic energy (k) before (A) and after (B) vertical interpolation; (C) Shear-stress distribution using the TKE method.
Water 10 00367 g007
Figure 8. Bed shear stress distribution along the cross-section, using the TKE method, the log-law method, and the depth–slope product method.
Figure 8. Bed shear stress distribution along the cross-section, using the TKE method, the log-law method, and the depth–slope product method.
Water 10 00367 g008
Figure 9. (A) Distribution plot of the <vw> direction (angle from up) resulting velocity in the vertical plane (<vw> = (<v>2 + <w>2)0.5) and schematized identification of vortex structure; (B) Vorticity (ω) of the canyon section part.
Figure 9. (A) Distribution plot of the <vw> direction (angle from up) resulting velocity in the vertical plane (<vw> = (<v>2 + <w>2)0.5) and schematized identification of vortex structure; (B) Vorticity (ω) of the canyon section part.
Water 10 00367 g009
Figure 10. Eddy viscosity distribution over the cross-section, using the TKE method.
Figure 10. Eddy viscosity distribution over the cross-section, using the TKE method.
Water 10 00367 g010
Figure 11. (A) Velocity magnitude and vector direction within the river section; (B) Shear stress distribution using the log method; (C) Shear stress distribution using the TKE method; (D) Eddy viscosity distribution based on the shear stress computed by TKE method.
Figure 11. (A) Velocity magnitude and vector direction within the river section; (B) Shear stress distribution using the log method; (C) Shear stress distribution using the TKE method; (D) Eddy viscosity distribution based on the shear stress computed by TKE method.
Water 10 00367 g011
Table 1. Reynold Stress notation.
Table 1. Reynold Stress notation.
Measurement Stress ComponentPlaneDirectionVelocity Gradient
−ρ <u’w’>τ’xzyzz∂<w>/∂x
τ’zxyxx∂<u>/∂z
−ρ <u’v’>τ’xyyzy∂<v>/∂x
τ’yxzxx∂<u>/∂y
−ρ <v’w’>τ’yzzxz∂<w>/∂y
τ’zyyxy∂<v>/∂z

Share and Cite

MDPI and ACS Style

Tomas, G.; Bleninger, T.; Rennie, C.D.; Guarneri, H. Advanced 3D Mapping of Hydrodynamic Parameters for the Analysis of Complex Flow Motions in a Submerged Bedrock Canyon of the Tocantins River, Brazil. Water 2018, 10, 367. https://doi.org/10.3390/w10040367

AMA Style

Tomas G, Bleninger T, Rennie CD, Guarneri H. Advanced 3D Mapping of Hydrodynamic Parameters for the Analysis of Complex Flow Motions in a Submerged Bedrock Canyon of the Tocantins River, Brazil. Water. 2018; 10(4):367. https://doi.org/10.3390/w10040367

Chicago/Turabian Style

Tomas, Gustavo, Tobias Bleninger, Colin D. Rennie, and Henrique Guarneri. 2018. "Advanced 3D Mapping of Hydrodynamic Parameters for the Analysis of Complex Flow Motions in a Submerged Bedrock Canyon of the Tocantins River, Brazil" Water 10, no. 4: 367. https://doi.org/10.3390/w10040367

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop