Two-Dimensional Numerical Modeling of Flow in Physical Models of Rock Vane and Bendway Weir Configurations

: While there are a wide range of design recommendations for using rock vanes and bendway weirs as streambank protection measures, no comprehensive, standard approach is currently available for design engineers to evaluate their hydraulic performance before construction. This study investigates using 2D numerical modeling as an option for predicting the hydraulic performance of rock vane and bendway weir structure designs for streambank protection. We used the Sedimentation and River Hydraulics (SRH)-2D depth-averaged numerical model to simulate flows around rock vane and bendway weir installations that were previously examined as part of a physical model study and that had water surface elevation and velocity observations. Overall, SRH-2D predicted the same general flow patterns as the physical model, but over- and underpredicted the flow velocity in some areas. These over- and underpredictions could be primarily attributed to the assumption of negligible vertical velocities. Nonetheless, the point differences between the predicted and observed velocities generally ranged from 15 to 25%, with some exceptions. The results showed that 2D numerical models could provide adequate insight into the hydraulic performance of rock vanes and bendway weirs. Accordingly, design guidance and implications of the study results are presented for design engineers.


Introduction
Meandering river channels are complex and dynamic systems with planimetric flow path variability confined by the valley floor or meander belt width. Infrastructure, agricultural lands, and buildings are commonly placed adjacent to and within the historical meander migration zone. Bank erosion due to lateral migration of meandering channels can encroach upon infrastructure. Engineers and scientists have developed methodologies to inhibit outer bank erosion to combat deleterious channel bend migration. Flow redirecting approaches, such as transverse instream structures, are used for stream bank protection rather than riprap revetments because these structures cost less and have more habitat benefits. These structures do not directly increase bank erosion resistance but rather alter flow patterns along the bank, thereby indirectly reducing applied flow shear stress. Examples of transverse features include rock vanes (also known as barbs) and bendway weirs. Rock vanes are structures oriented upstream that extend from the top of the bank with a sloping crest and intersect the bed after crossing the thalweg. Bendway weirs also are oriented upstream and have flat tops positioned above the low water surface and below the bank-full water surface elevation (WSE). These features produce a weir effect at lower flows. At high flows, bendway weirs redirect near-bed flow away from a bank and change secondary currents near banks, all of which help reduce bank-toe erosion.
Physical-model and numerical-model studies, therefore, are an important resource for design engineers and offer invaluable insight into how structure designs are likely to perform in an actual field setting.
Physical model studies, while providing valuable measurements and results, are time consuming, can be expensive when major configuration alterations are needed, and sometimes do not allow rapid assessment of alternative designs based on multiple transverse feature configurations or channel configurations. Numerical modeling offers an attractive alternative, and 2D studies have been performed to study design structures [28]. 3D models, however, are also limited by computational expenses and have mainly been applied to examine only a single structure or a couple of structures [15,20]. Thus, using 2D depth-averaged models is attractive from an engineering design standpoint. Depthaveraging of the flow field does not account for the vertical velocity component in 3D flow fields generated by these structures. Therefore, identifying and quantifying deviations of the predicted flow field from the actual flow field can examine these limitations and help ascertain whether a 2D model can adequately be used to examine structure performance.
To examine the predictive capability of a 2D model of flow in a channel fitted with configurations of rock vanes or bendway weirs, we compared data from the numerical model with data from a physical model. This paper describes the physical model, which provided measurements from baseline and three rock vane and four bendway weir configurations. The numerical model results section presents results of the 2D numerical model for one rock vane and one bendway weir configuration. Lastly, information is provided on how SRH-2D can be effectively applied during the design process.

Physical Model Facility
A two-bend physical model that yielded extensive water-surface elevation and velocity measurements was compared with numerical-model results from the numerical model SRH-2D. The Engineering and Research Center at Colorado State University constructed the physical model at 1:12 length scale. The overall model gradient was about 0.000863 V:H. The model was constructed using cross-sectional plywood templates with steel flashing and filled with sand material. A brushed concrete surface was placed between the plywood templates. Flows were delivered from the model sump to the headbox in a recirculating system. Tailwater conditions were established to provide approximately normal flow conditions throughout the model.
The physical model was based on a reach of the Middle Rio Grande between Cochiti Dam downstream to Angostura Diversion Dam near Algodones, New Mexico ( Figure 2). In this reach, the historic low-flow, braided, sand-dominated channel had been transitioned to an incising single-thread, gravel-dominated, lateral migrating channel as a result of installation of Cochiti Dam [29]. The river's top-width (TW) and radius of curvature (RC) were identified for channel bends in the prototype reach. In the physical model, two representative bend geometries (Figure 3) spanned the ranges of RC/TW observed from the prototype ( Table 1).
The channel configurations described here represent fixed-bed native topography. Even though Cochiti bend is upstream of the San Felipe bend in the middle Rio Grande (Figure 2), the physical model switched these two to represent the natural topography in the San Felipe site as the upstream bend and the topography in the Cochiti site as the downstream bend in the physical model ( Figure 4). This was done to accommodate the physical model within the laboratory configuration. The RC/TW of both bends are in Table  1.

Physical Model Transverse Feature Configurations
This section describes the physical model configurations with fixed bed native topography [30]. To measure and evaluate flows, [22,24,27]. We designed three rock vane and four bendway weir configurations, based upon the most common design values from the available design guidance [1][2][3][4][5][6][7][8][9][10][11]. Parameters are based on design recommendations. Table 2 provides the vane parameters and Table 3 provides the bendway weir parameters. Planimetric views of the vane design configurations (V01, V02, and V07) are presented in Figure 5 and bendway weir designs (BW01-BW04) are presented in Figure 6. The vanes and bendway weirs were constructed using plywood to create impermeability and with tight-fitting angular stones ranging from 91 to 152 mm along the intermediate axis, in accordance with the guidelines (Tables 2 and 3). Stones were sized to remain stable while testing physical models in the laboratory and for structures to be constructed according to the design recommendations in Tables 2 and 3. d100 is the largest rock diameter; L is the projected structure length; * crest length is reported. d100 is the largest rock diameter; L is the projected structure length.    [21]. Plywood in the center of each structure is used for impermeability.

Figure 8.
Bendway weir configuration looking downstream [23]. Plywood in the center of each structure is used for impermeability.

Physical Model Measurements
Two Acoustic Doppler Velocimeters (ADV) were used to measure velocity: a Son-Tek ® ADV with ±1% measured accuracy, and a Nortek Vectrino ® ADV with ±0.5% accuracy. Velocity was measured to the nearest 0.001 m per second (m/s), with measurement values ranging from approximately 0.003 to 1.3 m/s. The associated measurement error for these values (based on reported precisions) is negligible. For each configuration, 2500 velocity data points were measured [22,24,27]. To directly compare between laboratory configurations and between laboratory and numerical model results, measurements were made at consistent locations for baseline, rock vanes, and bendway weirs. To measure depth-averaged velocity (approximated as the velocity at 0.6 of the flow depth below the water surface), these data points were along the thalweg, along seven cross-sections, placed in a more dense grid around one structure in each bend for each test, and placed throughout the physical model (Table 4 and Figure 9).
The water-surface elevation for each configuration was measured at the bankline at approximately 0.15 m intervals using a Leica Total Station laser-surveying unit with a 3D point accuracy of 0.2 mm. Bed topography was measured with a Leica ScanStation 2 ® (Leica Geosystems AG, Heerbrugg, Switzerland), a terrestrial light detection and ranging unit (LiDAR) with 6-mm position accuracy, 4-mm distance accuracy, and 2-mm modeled surface precision.

SRH-2D Numerical Model
For this study, we used the SRH-2D (v 3.0) model [31], a two-dimensional (2D) depthaveraged numerical model specifically designed for river-flow hydraulics. SRH-2D adopts a zonal approach for coupled modeling of channels and floodplains. In this approach, a river system is broken down into modeling zones (delineated based on natural features such as topography, vegetation and bed roughness), each with unique parameters including flow resistance, which is specified in a Manning's roughness coefficient. SRH-2D uses an unstructured hybrid mixed element mesh that is based on the arbitrarily shaped element method [32] for geometric representation. This meshing strategy is flexible enough to model zones, enables greater modeling detail in areas of interest-and ultimately leads to increased modeling efficiency by balancing solution accuracy and computing demands.

Model Setup and Calibration
We used SRH-2D to compute channel hydraulics in the fixed-bed mode. Inputs included the physical model's topographies, discharge at the upstream boundary, and the measured downstream water surface elevation. We used the Surface-water Modeling System (SMS) software (v 13.1) [33] to generate a 2D mesh of quadrilateral elements in the main channel, to delineate model roughness areas and to assign the model boundary conditions. The pre-processed mesh as well as the roughness and boundary variables served as inputs to the SRH-2D solver.
The SRH-2D model was calibrated for baseline configuration (corresponding with the physical model's baseline configuration), one dataset for a rock vane configuration and one dataset for the bendway weir configuration. During the calibration, we adjusted the bed roughness values iteratively, and we compared the predicted flow fields with the measured flow fields using water-surface elevations and depth-average resultant velocities.
To calibrate the baseline condition, we applied a discharge rate of 0.34 cubic meters per second (m³/s) at the upstream boundary in the physical model. However, as the physical model operations had a small but unquantified amount of seepage, we needed to adjust the upstream boundary inflow rate and bed roughness in SRH-2D as part of the calibration process. We compared SRH-2D computed velocity results with lab-measured depth-averaged velocities and calculated the corresponding Root Mean Square Error (RMSE), Mean Absolute Error (MAE) and the Mean Absolute Percentage Error (MAPE). For this calibration, "error" is defined as the difference between computed and measured fields. Based on the RMS and MAE results, a discharge value of 0.325 m³/s (seepage of 0.015 m³/s) and a Manning's roughness coefficient of 0.018 were determined to be the values that best matched the measured data (Table 5). We used the velocity field for the RMSE, MAE, and MAPE analyses rather than the water-surface elevation (WSE) as velocity field had a larger variance in the velocity distribution than the WSEs did. This larger variance, which can be readily seen in the "filled" contour plots presented in the results section, made velocity a better discriminator of model performance between the configurations examined. Moreover, velocity is a more important parameter than water-surface elevation when evaluating the efficacy of a transverse feature design. Qi is ith discharge simulated; ni is ith Manning's n coefficient simulated.
The mesh resolution used in a 2D model can influence the model's results if not selected carefully. Before the roughness calibration, we tested three mesh sizes to establish grid independence (i.e., the most computationally efficient mesh where smaller cell sizes did not change the results significantly) ( Table 6 and Figure 10). The results are shown in Figure 11. While the three mesh sizes did not appreciably affect velocity results, the coarse mesh size did not adequately resolve the flow separation zones associated with the upstream bend contraction ( Figure 11). The medium-sized mesh ( Figure 10) was therefore selected for all the model runs.
After establishing grid independence and calibrating the roughness under the baseline configuration, the next step was calibrating SRH-2D against model measurements with structures in place. Because the LIDAR had a high accuracy rating, the surface topography of the rock structures had enough accuracy for the numerical model. Because the physical model replications of rock vanes and bendway weirs used angular rock approximately 91 to 152 millimeters (mm), at least 15 to 25 LiDAR data points were measured over the diameter for each angular rock-providing good representation of the rock structure roughness with the topography. The medium mesh cell sizes adequately captured rock forms as the cells were generally smaller than the rock diameter. Manning's roughness values representing skin friction were successively changed (n = 0.018-0.05) over the footprint of transverse features. Setting the structure roughness at the same value of the baseline calibrated n = 0.018 produced an ideal match to measured depth-averaged velocity. Scurlock et al. [34] also found that setting the physical model structure's roughness value equal to the numerical model's baseline roughness value resulted in the most ideal measured velocity match for models using FLOW 3D, Flow Science ® 's 3D numerical model.

Velocity Zone Comparison
Flow field comparisons between the physical model and the numeric model are categorized by the distinct characteristic flow zones that define rock vane and bendway weir hydraulics. These separate regions of identifiable velocity trends from rock vane and bendway weir structures (described in Figure 12 and  [12,13,15,23]. These regions include: the outer-bank reduction zone, the shifted-conveyance zone at the center of the channel and inner-bank, the outer-bank increased zone, and zones of convective acceleration at the structure crest and tip. Although vanes slope from the top of bank (crest zone) down to the bed at the free end (tip) they have essentially the same flow zones as bendway weir (Figure 12)..

Flow Zone Definition
Shifted Conveyance Zone that extends from near the structure tip to the point bar where velocity increases above baseline as longitudinal flows contract and a new thalweg location is created. Crest Zone where there is convective acceleration of flow passing over the crest of the structure.
Tip Turbulent zone at the tip where flow contracts horizontally. This zone is frequently associated with vortex shedding.
Outer Bank Increase Zone Zone along the outer bank shortly downstream of the crest of each structure as flow vertically expands after passing over the structure crest.

Outer Bank Reduction
Zone along the outer bank covering most of the bank between structures and downstream of the most downstream structure, where bank-line velocity is decreased from baseline conditions. This zone generally extends laterally from the bankline to near structure tips.

Depth-Averaged Velocity Verification
The physical modeling took measurements at the same flow locations for the baseline, rock vanes and bendway weirs as described in Section 2. To approximate depth-averaged velocities, we took these measurements at 60% of the flow depth from the water surface (Table 4 and Figure 9). To determine if the 60% depth measurements were appropriate to directly compare with SRH-2D output and to measure the flow field, we measured a denser grid at 5% depth intervals around one structure in each rock vane and bendway weir configuration, and measured 20% depth intervals along the outer bank. To compare the depth-averaged velocity at each of these measurement locations with the 60% depth measurements, we integrated the 5% and 20% depth interval data over the depth using the trapezoidal rule. Overall, the integrated velocities from the 5% and 20% depth compared well with the velocities measured at 60% of the flow depth (Figures 13 and 14). This result indicates that the 60% depth velocity measurements are a good approximation of the depth-averaged velocity. Figure 13. Point velocity (60% of depth from water surface) compared to the integrated depth average velocity using either 5% or 20% depth intervals for bendway weir configuration BW02. Figure 14. Point velocity (60% of depth from water surface) compared to the integrated depth average velocity using either 5% or 20% depth intervals for V01 configuration. Figure 15 shows the physical model's measured and SRH-2D model's predicted water-surface elevation (WSE) ,and Figure 16 shows the velocity fields for the upstream bend (Bend 1) and the downstream bend (Bend 2). In each of these figures, the colored shading indicates the contours, which were derived from the data using inverse distance weighting to create a Raster surface. The "difference" plot derived by subtracting the physical model's measured WSEs from the SRH-2D-predicted WSEs is also shown in Figure 15. Because there are larger variations in the velocity field, the "difference" plot for velocity is shown in a separate figure (Figure 17). To avoid distortions introduced by interpolation, only results at the measurement locations are shown in Figure 17.

Baseline Configuration: Physical and SRH-2D Numerical Model Comparisons
Overall, the WSEs from the physical and numerical model results show good agreement, with a slight underprediction in Bend 1 (predicted values are generally 0 to 0.006 meters (m) lower than the measured values) and a slight overprediction in Bend 2 (predicted values are generally 0 to 0.006 m higher than the measured values). The differences between the measured and predicted velocity contours are more pronounced than the differences in WSE. The remainder of this section, therefore, focuses on comparing the measured and predicted velocity contours.
At the upstream section of the physical model, velocity is highest in the area on the inside of Bend 1 (Figure 16). The high-velocity zone extends through and immediately downstream of the contraction, with the highest velocity occurring on the outside of the bend. Upstream of the contraction, velocity is lower as the channel experiences a backwater that develops the head for the flow to pass through the contraction with associated energy losses. Downstream of the contraction, the next-highest velocity area is along the outside of the Bend 1, starting from just upstream of the bend apex and stretching through the downstream side of the bend. Through the crossing (i.e., transition between Bend 1 and Bend 2), the high-velocity zone is along the left bank until the zone reaches the inside of Bend 2. Through the bendway, the zone of high velocity crosses over the channel, such that the high velocity reaches the outside bank of Bend 2 near the bend apex. The highvelocity zone stays along the outside bank through to the physical model exit. SRH-2D model results for the baseline configuration (Figure 16) also show a higher velocity area on the inside of Bend 1 in the upstream section. There is a high-velocity area at the contraction, with the highest velocity near the inside of the bend, which appears to be due to the contraction effect rather than entrance conditions. Throughout Bend 1, the high-velocity area is nearly symmetrical with a smaller higher-velocity area near the outer bank on the downstream of the bend. A small skew in the high-velocity area towards the outer bank at the downstream boundary of Bend 1 remains on the left bank as the flow transitions into Bend 2. As the flow enters Bend 2, high-velocity flows are on the inside of Bend 2. Downstream of the transition, the high-velocity area becomes more symmetrical within the channel width, with a slight skew towards the outer bank and high-velocity flows remain in this alignment until the downstream model boundary. The highest-velocity area is just downstream of the Bend 2 apex and extends for a short distance.
The velocity differences are compared in Figure 16. The overprediction is mostly about 0 to 0.09 m/s with a few areas of 0.09 to 0.18 m/s. Overpredictions are also within 20% of the measured velocities, while most differences are much less. These differences are relatively small, and it may be concluded that there is good correspondence between the measured and predicted through both bends.
There is general agreement between the SRH-2D WSE and velocity predictions and the measurements.

Rock Vane Configuration: Physical and SRH-2D Numerical Model Comparisons
This section presents results for the rock vane configuration V01 installation WSE ( Figure 18) and depth-average velocity ( Figure 19). Like the baseline configuration, the measured WSEs compared well with the predicted elevations, with most of the differences from −0.015 to 0 m, and some small areas with differences between 0 and 0.015 m, and the rest of this section, therefore, focuses on velocity.
Physical model measurements show that the highest velocity area (0.63 to 0.79 m/s) is in the shifted conveyance zone between the bend apex and near the downstream model boundary. The remainder of the shifted conveyance zone is a high-velocity area, with magnitudes ranging from 0.50 to 0.62 m/s. There is also high velocity in the crest and tip zones for the first five vanes. The vanes induce low-velocity (0.15 to 0.23 m/s) banklines out to near the tips of the vanes, in the outer bank reduction zone. The highest velocity region is in the shifted conveyance zone from about the bend apex downstream to near the second-to-last vane.
SRH-2D-predicted velocity patterns are close to the physical model's measured flow patterns. The outer bank reduction zone is well defined from the bank to near the vane tips, with low flow velocities ranging from 0 to 0.23 m/s. The high velocity in the shifted conveyance zone has similar flow regions and velocity ranges in the physical model's measurements and the SRH-2D predictions. The SRH-2D-predicted highest velocity is in the same location as the physical model's measured highest velocity in the streamwise direction, but the predicted velocity area is farther away from the vane tips in the transverse direction than the physical model's measured area.
The velocity differences are compared in (Figure 20

Bendway Weir Scenario Configuration: Physical and SRH-2D Numerical Model Comparisons
This section presents results for bendway weir configuration BW02 WSE ( Figure 22) and velocity (Figures 23 and 24).
Physical model measurements show a high-velocity zone that moves away from the outside of Bend 1 in the outer bank reduction zone, in contrast to the baseline configuration without the structures, where the high velocity was along the outer bank. The high velocity shifts away from the weirs toward the inner bank in the shifted conveyance zone. There is, however, a local high-velocity area between Bendway Weirs 7 and 8 within the outer bank reduction zone, but it does not extend to the bankline [23]. This local area could be from the interactions between flow eddies and the structures. Within the shifted conveyance zone, local pockets of highest velocities are seen close to the tips of Bendway Weirs 6 and 8. Upstream, at the contraction, there is high flow acceleration over the first bendway weir that is not observed over the other weirs. This high flow acceleration of the first bendway weir is also predicted by SRH-2D-suggesting that it is driven largely by the contraction effects. SRH-2D was able to predict the shift in high velocity away from the outer bank reduction zone. The highest-velocity zone was, however, a little farther away from the weir tips than the measured zone in the physical model. SRH-2D did not predict the local highvelocity zone between Bendway Weirs 7 and 8-suggesting that the presence of the zone in the physical model likely caused flow features with significant vertical velocity components ( Figure 25). In design, underprediction along the bankline could lead to a failure to recognize the potential for accelerated bank erosion that could flank Bendway Weir 8. Approaches to provide a stable design configuration are included in Section 4.3. Overall, SRH-2D predicted regions of both higher and lower velocities in the shifted conveyance zones as well as between the weirs. Velocity vectors illustrate the high velocity along the crest zone of the bendway weirs (Sections A-A and C-C Figure 25). The velocity downstream of the bendway weir near the bankline (Section C-C in Figure 25) is nearly as great as the velocity in the shifted conveyance zone (Section D-D in Figure 25).
The WSEs at Bendway Weir 6 (counted from bend entrance) and an area between Bendway Weirs 7 and 8 were overpredicted by SRH-2D. However, the measured WSE is about 0.15 m lower than nearby measurements-suggesting these locations were inaccurately measured.

SRH-2D Model Performance
A summary of velocity results for the baseline, rock vane and bendway weir configurations is presented in Table 8. Velocity zones, defined by distinct hydraulic conditions ( Figure 12 and Table 7), are identified where differences occur between computed and measured velocities.
Measured and predicted results from SRH-2D have the same general velocity patterns, indicating that SRH-2D can be used for design applications with a factor of safety as described in Section 4.3. SRH-2D overpredicts velocity in the crest zone, tip zone and sometimes in the shifted conveyance zone-generally where there is a vertical velocity component and turbulent kinetic energy. These increases are often on the order of 0 to 0.09 m/s and occasionally as high as 0.18 m/s. SRH-2D also shows both over-and underpredictions compared to measurements in zones where it cannot adequately resolve flow conveyance shifts due to the neglect of vertical-velocity components. In these zones, the increased and decreased velocities are necessary for mass conservation. Nonetheless, the predicted velocities are within ±25% of the measured velocities, with some localized exceptions. The difference is much less in most locations, which suggests good agreement between the measured and predicted results. In the shifted conveyance zone (Section D-D in Figures 21 and 25), the velocity vectors are nearly uniform.
The structures induce secondary currents in the opposite direction of typical bend flow in the shifted conveyance zone (Section A-A in Figures 21 and 25). In general, however, the longitudinal velocity magnitude is much greater than the lateral or vertical velocity components. This higher longitudinal velocity is the underlying reason for the velocity measured at 60% of the depth below the water surface matching the integrated/computed depth average velocity (Figures 13 and 14).
Papanicolaou et al. [28] compared the FESWMS 2D model with a physical model. They reported differences of up to 60% between measured and predicted flow velocity using a single set of model input (i.e., the same Manning's n and eddy viscosity) for the whole domain. When Papanicolaou et al. [28] changed the model input variables for each flow zone, the model performance improved, predicting differences comparable to those noted in this study (i.e., velocity differences up to 25%). SRH-2D was able to reach an accuracy of within ±25% using a single Manning's n input value for the whole domain, which is easier for design application.

Baseline
The highest velocity area starts on the inside of upstream part of Bend 1, moves to outside of Bend 1 downstream of constriction and stays on that bank to the upstream part of Bend 2, then switches to outside of Bend 2 about the bend apex and continues to the downstream boundary. SRH-2D predicted less variability in the lateral velocity area, which continues until the wide section near model outlet.

V01
Both measured and predicted results show that velocity in outer bank reduction zone is lower in the rock vane configuration than in the baseline configuration. Within the rock vane configuration, there is a higher velocity in tip zone, and the highest velocity is in the shifted conveyance zone near bend apex downstream of vane downstream to Vane 6 (predicted), and Vane 7 (measured). SRH-2D overpredicts in crest and tip zone and between Vanes 1 and 2, and 2 and 3, and downstream of Vanes 3 and 5, and in shifted conveyance zone downstream of Vane 5, but by less than ±25%.

BW02
Both measured and SRH-2D results show that velocity in outer bank reduction zone is lower and is higher in the shifted conveyance zone in the bendway weir configuration than in the baseline configuration.SRH-2D overpredicts in the crest, tip, and outer bank increase zone downstream of contraction. SRH-2D both shows increase and decrease in shifted conveyance zone (overpredicted and underpredicted regions in zone). Over-and under-predictions are less than ±25%.

Modeling Summary
Transverse features or indirect bank-stabilization methods such as rock vanes and bendway weirs, based on experience or anecdotal design assumptions, have been used for many decades. We examined the ability of a two-dimensional depth-averaged model to predict the flow field around these structures in a two-bend 1:12-scale physical model developed from two Middle Rio Grande river bends. These bends represent conditions downstream of Cochiti Dam where the river channel has transitioned from a wide, shifting sand substrate and low-flow braided condition to a gravel-dominated bed, single thread, slightly sinuous channel.
We measured depth and velocity in the native topography simulated in the physical model without the rock structures (baseline) and then in three rock vane configurations and four bendway weir configurations. The measured data comprised about 2500 velocity measurements within the physical model at the design flow of 0.34 m³/s. LiDAR topography was measured for baseline and each configuration.
The physical model was not impervious and leaked. Leakage was not measured but was determined by evaluating combinations of reduced flow rates and Manning's roughness coefficient. This calibration resulted in selection of water flowrate = 0.325 m³/s and Manning's n = 0.018. The depth-averaged SRH-2D developed by [31] used the adjusted model inflow and calibrated roughness together with each set of LiDAR (converted to a surface) and a computational mesh to evaluate each physical model configuration.

Design Implications of Model Results
The good correspondence between the SRH-2D predictions and the physical model measurements suggests that SRH-2D can be a useful tool for modeling hydrodynamics of rock vanes and bendway weirs for design. Key SRH-2D model observations derived from evaluating all the model runs for the configurations in Tables 2 and 3-and how they may  influence design-are:  SRH-2D model results matched the overall flow patterns found in physical model measurements for rock vanes and bendway weirs. Thus, SRH-2D models can be used to compare the effects of spacing, length, horizontal angle (rock vanes) and orientation angle of different structure configurations.  SRH-2D did not compute the high-velocity zone between Bendway Weirs 7 and 8 ( Figure 25). Higher velocities exist on the downstream portion of river bends ( Figure 16). Additional bank erosion countermeasures such as bioengineering [35] could be employed between bendway weirs located in the downstream section of river bends, as well as between other weirs.  SRH-2D computes increased resultant velocity when compared to physical model measurements over transverse feature crests and right at the tips where there is a vertical-velocity component and turbulent kinetic energy. Overprediction varies between approximately 15 and 25%, with a few locations over this range. This overprediction provides conservative design. And we recommend that model results be used over weir crests when designing erosion countermeasures such as bioengineering [35].  SRH-2D may not fully resolve conveyance shifts. This failure to fully resolve can lead to both overprediction and underprediction of velocity in certain zones of the flow field by up to 25%, depending on the configuration. The highest velocity found in the shifted conveyance zone or at the tip is conservative and should be used for design tip scour and riprap sizing. Transverse feature rock sizing procedures from large-scale flume experiments [36] should be used, along with an appropriate factor of safety [37]. These higher velocities might affect sediment transport analysis, but this study did not determine the magnitude of this effect.
 SRH-2D underpredicts bankline velocity in some locations by up to about 25%. Increase bankline velocity by at least this amount during design.

Conclusions
The velocity differences between SRH-2D predictions and physical-model measurements are more pronounced (Figures 17, 20 and 24) than WSE differences (Figures 15, 18  and 22). Results show a useful agreement between SRH-2D estimates and the measurements. Therefore, SRH-2D is a suitable tool to compare combinations of transverse feature length, spacing, horizontal and vertical angles (rock vanes), and crest elevation (bendway weirs). SRH-2D generally predicted increased velocity (between 15 and 25%) in flow regions with a vertical velocity component and turbulence kinetic energy (tip and crest zones). These numerical model predictions are suitable for conservative designs. There also were local regions where SRH-2D underpredicted the velocity (up to 25%) at the outer bank zone and at the tip zone of select features. Increasing bankline velocity for design up to 25% is recommended. The crest extending from the top of bank and sloping down to past the thalweg in rock vanes allows rock vanes to provide much better bank-line velocity reduction than bendway weirs. As such, bendway weirs are not recommended for bank protection-and should only be used for thalweg management unless additional erosion countermeasures are also employed.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to not having a web site available for public access.