To investigate the influences of river bend geometry on flow structures, three-dimensional simulations were carried out. The influence of aspect ratio, channel sinuousness and curvature ratio were investigated. Simulation results show that with an increasing of channel aspect ratios, the streamwise variation of flow velocities decreases significantly, which is not the case for the spanwise direction. Simulation results also reveal the significant influence of the curvature ratio on the redistribution of primary flows. The streamwise variation of free-surface and bottom-wall velocities was compared and the phase differences between them was identified, which was described as anti-phase and out-of-phase patterns. Deep-water channels with an aspect ratio of 0.15 show significant phase-lag of flow velocities compared with channel geometry, this newly observed phenomenon implies the shift of the peak scouring point in fluvial river meanders. For the first time to our knowledge, the streamwise and spanwise variation of primary and secondary flow strength, as well as their spatial derivatives, are quantitatively provided, which are valuable for theoretical studies. It is also found that the widely adopted assumption that the streamwise velocity component is one order magnitude larger than the spanwise, is only valid for river bends with small curvatures and aspect ratios.
INTRODUCTION
As one of the most common river patterns in nature, meandering rivers have flow structures much more complicated than those observed in straight rivers. They are characterized by a traverse-inclined free-surface, secondary currents and spiral flows (Johannesson & Parker 1989a, 1989b). These complicated flow structures contribute to sediment transport, bed deformation, bank erosion and river meander development (Kikkawa et al. 1976; Olsen 2003; Hooke & Yorke 2010, 2011; Xu et al. 2011; Kasvi et al. 2013).
A deep understanding on the flow structures in meandering rivers is of great importance in theoretical studies by using the perturbation method (Ikeda et al. 1981; Liu et al. 2009). For the perturbation method, it is crucial to correctly estimate the magnitude of each variable in the governing equations of flow motion, and choose small parameters rationally. Different estimations may lead to completely different theoretical results (Smith & McLean 1984). To the authors’ knowledge, however, a comprehensive understanding on the flow structures in river bends with diverse geometric shapes is still lacking.
Two-dimensional (2D) numerical modeling of river bend flows also have some limitations. To investigate the flow and the resultant landform development in river bends, many 2D numerical models have been developed based on the shallow-water equations (Liu et al. 2009; Hardy et al. 2011; Parker et al. 2011). However, it should be noted that most of these models are based on certain assumptions, such as hydrostatic pressure, small curvature, large aspect ratio (wide-to-depth ratio) etc., which limit their applications.
Three-dimensional (3D) numerical simulations of flows in open channels using computational fluid dynamics (CFD) is nowadays becoming more and more powerful and efficacious owing to the significant development of high-performance computers and numerical methods. Using modern CFD techniques like direct numerical simulation and large eddy simulation (LES), even the most complex turbulent structures in straight open channels can be resolved (Ji et al. 2012, 2013; San 2014; Li et al. 2015). Recently, LES of turbulent flows in curved open channels were carried out and the results agreed well with experimental measurement (Booij 2003; Van Balen et al. 2009; Xu et al. 2013). However, as indicated by Xu et al. (2013), LES of turbulent flows requires extremely high-resolution grids especially near the channel walls and even the most powerful parallel computing facility cannot afford the simulation of turbulent channel flows with high Reynolds numbers. This greatly limits its application in systematic investigations of flow structures in river bends with a variety of geometric shapes. With regards to this, traditional CFD models based on the Reynolds averaged Navier–Stokes (RANS) equations and certain turbulent models are effective options for simulation of turbulent flows in channel bends due to their high computing efficiency and reasonably good results, especially where the time-averaged flow velocities are concerned (Nguyen et al. 2007; De Marchis & Napoli 2008; Kashyap et al. 2012; Ramamurthy et al. 2013). Kashyap et al. (2012) investigated the influences of the curvature ratio (bend radius/channel width, R/B) and aspect ratio (channel width/flow depth, B/H) using the RANS method. Ramamurthy et al. (2013) also achieved good simulation results on flow in a 90 degree bend using both k-ε renormalization group (RNG) and Reynolds stress model in a RANS framework. However, the available studies only focus on a specified channel shape, for example circular shape in Kashyap et al. (2012) and 90 degree bend in Ramamurthy et al. (2013). Considering the diversity of river bend shapes in nature, more effort is needed for a comprehensive understanding on the complicated flow structures in river bends.
In this paper, 3D RANS-based numerical simulations of turbulent flow in a curved open channel are carried out. The RNG k-ɛ turbulent model is adopted for the closure of the governing equations and the finite volume method (FVM) is applied for spatial discretization. The rest of the paper is organized as follows. In the next, the adopted methodology is introduced. After that, the numerical model is verified by comparing the numerical results with three groups of experimental data. In the next section, the influences of geometric shapes of channel bends on flow structures are discussed for channels with different curvatures, channel widths and water depths. Streamwise and spanwise velocities, together with their spatial gradients, are presented. The conclusions are presented in the final section.
3D CFD SOLVER AND VERIFICATION CASES
Governing equations





It should be noticed that the standard k-ɛ model is based on the hypothesis of isotropic eddy viscosity and turbulent diffusion. Therefore, it cannot be applied in the simulation of anisotropic turbulent flows in its original form. However, the k-ɛ model is still one of the most popular turbulence models, owing to its simplicity and good performance. Fortunately, its modified version, RNG k-ɛ model is capable of handling anisotropic turbulent flow structures, which have been widely verified in the field of river hydrodynamics (Olsen 2003; Abad 2005; Stoesser et al. 2010; Bai et al. 2014). Therefore, the RNG k-ɛ model was adopted in the present study.
Discretization, grids and boundary conditions
Computational grid (a) structured grid indexing scheme (b) schematic of 3D mesh.
Computational grid (a) structured grid indexing scheme (b) schematic of 3D mesh.
The no-slip boundary condition was applied on the channel bottom and side-walls, while the free-slip boundary condition was adopted on the free-surface. Note that the rigid wall approximation of the free-surface adopted in this study suppresses the development of surface waves. However, the introduced errors on the water level variations is negligible (less than 5% of the water depth (Stoesser et al. 2010)) due to the small Froude number (smaller than 0.71 for all cases) in this study. A uniform inflow was imposed on the inlet boundary, and the non-reflecting outflow condition was imposed on the outlet boundary. A standard wall function was applied on the no-slip walls. The Reynolds number based on mean velocity and water depth ranges from 888.9 to 8,571.4. A uniform grid was adopted in all three directions with a total element number around 106. There were 64 grids in the vertical direction, which yields ranging from 26 to 97 for different cases. A constant time step of 0.01 s was adopted, which yields a CFL (Courant–Friedrichs–Lewy) number ranging from 0.01 to 0.18. Each case ran around 10,000 time steps to reach a quasi-steady state. The running time for each case is around 48 hours on a 4-core 2.4 GHz CPU.
Verification cases
To validate the numerical model, three groups of laboratory experiments with different curvatures were adopted, namely Group I (sharply curved channel), Group II (mildly curved channel) and Group III (sine-generated channel).
Group I: sharply curved channel
Verification of vertical profiles of streamwise velocity – u and spanwise velocity – v at the cross-section at θ = 90° (x is the coordinate along the channel width with zero locates at the inner bank).
Verification of vertical profiles of streamwise velocity – u and spanwise velocity – v at the cross-section at θ = 90° (x is the coordinate along the channel width with zero locates at the inner bank).
Group II: mildly curved channel
Setup of channel bend: (a) experimental flume (Booij 2003); (b) computational domain.
Setup of channel bend: (a) experimental flume (Booij 2003); (b) computational domain.

Verification of vertical profiles of streamwise velocity – U, spanwise velocity – V, vertical velocity – W at the cross-section at θ = 135° (y is the coordinate along the channel width with zero locates at the inner bank).
Verification of vertical profiles of streamwise velocity – U, spanwise velocity – V, vertical velocity – W at the cross-section at θ = 135° (y is the coordinate along the channel width with zero locates at the inner bank).
Group III: sine-generated channel
Experimental facility for validation on surface flow in channel bends (unit: mm).
Experimental facility for validation on surface flow in channel bends (unit: mm).
RESULTS AND DISCUSSION
In order to investigate flow structures in channel bends with a wide range of geometries, in total 10 simulation cases were performed, including three cases (RUN F30, F60, F110) with different deflection angles (channel sinuousness) ranging from 3 to 110°, four cases (RUN H15, H30, H60, H100) with different aspect ratios (B/H) ranging from 0.15 to 1.0 and three cases (RUN B2, B5, B15) with different curvature ratios ranging from 1.0 to 7.5. The Reynolds numbers ranged from 888.9 to 8,571.4 and Froude numbers ranged from 0.03 to 0.20.
Influences of channel bend sinuousness on flow structures
In order to study the influences of channel sinuousness on flow structures, three maximal deflection angles of 30, 60 and 110° were adopted in the numerical simulations. Table 1 lists other geometry and simulation parameters.
Geometric parameters and simulation parameters of bends with various sinuousness (Group I)
RUN . | Geometric parameters . | Reynolds number, Re . | Water depth, H/m . | Mean velocity, U0/m/s . | ||
---|---|---|---|---|---|---|
Maximal deflection angle,° . | Wave length along centerline, M/m . | River width, B/m . | ||||
F30 | 30 | 2 | 0.15 | 8,571.4 | 0.10 | 0.2 |
F60 | 60 | 2 | 0.15 | 8,571.4 | 0.10 | 0.2 |
F110 | 110 | 2 | 0.15 | 8,571.4 | 0.10 | 0.2 |
RUN . | Geometric parameters . | Reynolds number, Re . | Water depth, H/m . | Mean velocity, U0/m/s . | ||
---|---|---|---|---|---|---|
Maximal deflection angle,° . | Wave length along centerline, M/m . | River width, B/m . | ||||
F30 | 30 | 2 | 0.15 | 8,571.4 | 0.10 | 0.2 |
F60 | 60 | 2 | 0.15 | 8,571.4 | 0.10 | 0.2 |
F110 | 110 | 2 | 0.15 | 8,571.4 | 0.10 | 0.2 |
Note: Reynolds number Re is based a characteristic length scale of hydraulic radius.
Secondary current at bend apex (a) and near bank velocity vector (b).
Streamwise profiles of flow velocities in bends with various sinuousness.
It can be seen that for all cases, all three velocity components show quasi-periodic variations along the curved channel, while the streamwise velocity is consistently larger than the others. The flow transition length, i.e. the distance from the inlet of the channel to the position where the quasi-periodic flow developed, decreases with an increasing sinuousness. For example, obvious quasi-periodic variations can be found in the cross-flow velocities beyond s >0.8 m for RUN F30, while large variations of the cross-flow velocities develop from the channel inlet for RUN F110. For the case with small curvature (RUN F30), the streamwise velocity near the free-surface is significantly larger than that near the bottom-wall due to the shear drag effects. However, this is not necessarily the truth for the cases with a large sinuousness. The strong secondary flow developed in the channel bends with a large curvature deflects the ‘flow core’ (the area with large streamwise velocity) from the center of the channel to the side-walls, as shown in Figure 8. The vertical velocity is close to zero, which verifies the assumption of adopted in many theoretical studies (Ikeda et al. 1981). The streamwise velocity near the side-walls has a larger variation amplitude and a two-time longer wave length than that on the center lines. The spanwise velocities near the bottom-wall and free-surface show similar variation amplitude and wave length but in an anti-phase synchronized pattern due to the presence of the secondary currents.
To compare the magnitude of characteristic variables of flow in channel bends, the maximum values of the following variables are computed in Table 2. In the following variables, subscribe ‘s’ refers to the variation along a streamwise curve and ‘n’ along a spanwise direction:
Range of streamwise variation of streamwise velocity and spanwise velocity;
Range of spanwise variation of streamwise velocity and spanwise velocity;
Mean gradient of streamwise variation of streamwise velocity and spanwise velocity;
Mean gradient of spanwise variation of streamwise velocity and spanwise velocity.
Magnitude of flow velocities in bends with various sinuousness
RUN . | ||||||||
---|---|---|---|---|---|---|---|---|
F30 | 1.00 | 0.35 | 0.75 | 0.25 | 0.150 | 0.052 | 0.75 | 0.25 |
F60 | 0.88 | 0.40 | 0.60 | 0.25 | 0.131 | 0.060 | 0.60 | 0.25 |
F110 | 1.10 | 0.50 | 1.00 | 0.40 | 0.165 | 0.075 | 1.00 | 0.40 |
RUN . | ||||||||
---|---|---|---|---|---|---|---|---|
F30 | 1.00 | 0.35 | 0.75 | 0.25 | 0.150 | 0.052 | 0.75 | 0.25 |
F60 | 0.88 | 0.40 | 0.60 | 0.25 | 0.131 | 0.060 | 0.60 | 0.25 |
F110 | 1.10 | 0.50 | 1.00 | 0.40 | 0.165 | 0.075 | 1.00 | 0.40 |
It can be seen from Table 2 that when the maximal reflection angle ranges from 30 to 110°, the sinuousness has trivial effects on the magnitude of flow velocities in channel bends. The streamwise variation ranges of velocities are equivalent to their spanwise counterpart. However, the spanwise variation gradients of velocities are over five times larger than their streamwise counterpart, because the streamwise length of the channel is much larger than the channel width. Considering the feature of vertical profile of lateral velocity distributed over the depth, see Figure 6(b), can be roughly estimated as twice of the maximum lateral velocity
, namely
. Here,
represents the characteristic velocity of secondary current and can be approximated as
. ‘U’ refers to the mean flow velocity in the streamwise direction, which can be viewed as the primary flow velocity. Therefore,
/2U can be taken as a ratio of the strength of secondary current to the strength of primary flow. RUN F110 has a ratio of
/U = 0.4, see Table 2, which means that the ratio of secondary flow to primary flow strength is about 20%. The velocity magnitude of secondary current in highly sinuous channel bends can be as large as 20% of the primary flow, which may invalidate many theoretical assumptions (Bai & Yang 2011; Bai et al. 2012).
Influences of aspect ratio on flow structures
As one of the key features which significantly affect flow structures, the aspect ratio of the cross-section (channel width/flow depth, B/H) diverges drastically from one river to another. For instance, in China, the averaged aspect ratio of the Yellow River is over 40 times larger than that of the Changjiang River (Bai et al. 2008). In fact, even for the same river, the aspect ratio of cross-section varies within a quite wide range in different sections and in different seasons. To study the influences of aspect ratio on flow structures, numerical simulations were carried out for sine-generated river bends with a fixed channel width of B = 15 cm and a varying water depth H from 15 to 60 cm. The corresponding aspect ratio of the channel ranges from 0.15 to 1. The maximum deflection angle of 110° was adopted on the purpose of strengthening the secondary flows. See Table 3 for other geometric and simulation parameters.
Geometric parameters and simulation parameters of bends with various aspect ratios (Group II)
RUN . | Geometric parameters . | Reynolds number, Re . | Water depth, H/m . | Mean velocity, U0/m/s . | |||
---|---|---|---|---|---|---|---|
Maximal deflection angle /° . | Wave length along centerline M/m . | River width B/m . | Aspect ratio (B/H) . | ||||
H15 | 110 | 2 | 0.15 | 1.0 | 2,500.0 | 0.015 | 0.2 |
H30 | 110 | 2 | 0.15 | 0.5 | 4,285.7 | 0.030 | 0.2 |
H60 | 110 | 2 | 0.15 | 0.25 | 6,666.7 | 0.060 | 0.2 |
H100(F110) | 110 | 2 | 0.15 | 0.15 | 8,571.4 | 0.100 | 0.2 |
RUN . | Geometric parameters . | Reynolds number, Re . | Water depth, H/m . | Mean velocity, U0/m/s . | |||
---|---|---|---|---|---|---|---|
Maximal deflection angle /° . | Wave length along centerline M/m . | River width B/m . | Aspect ratio (B/H) . | ||||
H15 | 110 | 2 | 0.15 | 1.0 | 2,500.0 | 0.015 | 0.2 |
H30 | 110 | 2 | 0.15 | 0.5 | 4,285.7 | 0.030 | 0.2 |
H60 | 110 | 2 | 0.15 | 0.25 | 6,666.7 | 0.060 | 0.2 |
H100(F110) | 110 | 2 | 0.15 | 0.15 | 8,571.4 | 0.100 | 0.2 |
Note: RUN H100 is RUN F110 in simulation of Group I.
Streamwise profiles of flow velocities in bends with various aspect ratios (case F110 is equivalent to H100).
Streamwise profiles of flow velocities in bends with various aspect ratios (case F110 is equivalent to H100).
Table 4 further quantitatively shows the variations of flow velocities in channel bends with different aspect ratios. It can be seen that, with the decreasing of the aspect ratio, the variation amplitude of U in the streamwise direction decreases monotonously, while the variation amplitude of V in the streamwise direction increases monotonously. This can be attributed to the stronger secondary current in the channel bend with a larger water depth. Similar situations can be found in the spanwise variation of velocity components.
Magnitude of flow velocities in bends with various aspect ratios
RUN . | ||||||||
---|---|---|---|---|---|---|---|---|
H15 | 1.60 | 0.63 | 1.25 | 0.35 | 0.24 | 0.090 | 1.25 | 0.35 |
H30 | 1.50 | 0.75 | 1.25 | 0.50 | 0.225 | 0.113 | 1.25 | 0.50 |
H60 | 1.25 | 0.80 | 0.75 | 0.60 | 0.188 | 0.120 | 0.75 | 0.60 |
H100(F110) | 1.10 | 0.50 | 1.00 | 0.40 | 0.165 | 0.075 | 1.00 | 0.40 |
RUN . | ||||||||
---|---|---|---|---|---|---|---|---|
H15 | 1.60 | 0.63 | 1.25 | 0.35 | 0.24 | 0.090 | 1.25 | 0.35 |
H30 | 1.50 | 0.75 | 1.25 | 0.50 | 0.225 | 0.113 | 1.25 | 0.50 |
H60 | 1.25 | 0.80 | 0.75 | 0.60 | 0.188 | 0.120 | 0.75 | 0.60 |
H100(F110) | 1.10 | 0.50 | 1.00 | 0.40 | 0.165 | 0.075 | 1.00 | 0.40 |
Influences of channel curvature on flow structures
Channel curvature, usually denoted by the ratio of curvature radius to channel width (R/B), is a geometric parameter which describes the channel's non-straight morphology. The difference between channel curvature and channel sinuousness is that channel sinuousness only represents the curve degree of the channel centerline, while channel curvature combines the effects of channel shape and channel width. Emmett & Leopold (1964) measured 50 meandering rivers and found that the value of (R/B)2/3 ranged from 1.5 to 4.2 with an average of 2.7. They concluded that for alluvial meandering rivers, although their flow rates and geometric dimensions differed drastically, the value of (R/B)2/3 was close to a constant. Therefore, meandering rivers in alluvial plain areas show similar appearances on satellite images. However, for mountain rivers, the value of (R/B)2/3 may be far beyond that range due to the geological conditions. To study the influences of the channel curvature on flow structures, three sine-generated channel bends with the same maximum deflection angle of 110° but different channel widths were adopted in the simulations of Group III. The corresponding curvature ratio ranges from 1.0 to 7.5 (with (R/B)2/3 ranges from 1.0 to 20.5). The aspect ratios of the channels were set to 0.25. Table 5 shows the geometric and simulation parameters.
Geometric parameters and simulation parameters of bends with various curvature ratios (Group III)
RUN . | Geometric parameters . | Reynolds number, Re . | Water depth, H/m . | Mean velocity, U0/m/s . | |||
---|---|---|---|---|---|---|---|
Maximal deflection angle /° . | Wave length along centerline M/m . | River width B/m . | Curvature ratio (R/B) . | ||||
B2 | 110 | 2 | 0.02 | 7.5 | 888.9 | 0.008 | 0.2 |
B5 | 110 | 2 | 0.05 | 3.0 | 2,222.2 | 0.02 | 0.2 |
B15(H60) | 110 | 2 | 0.15 | 1.0 | 6,666.7 | 0.06 | 0.2 |
RUN . | Geometric parameters . | Reynolds number, Re . | Water depth, H/m . | Mean velocity, U0/m/s . | |||
---|---|---|---|---|---|---|---|
Maximal deflection angle /° . | Wave length along centerline M/m . | River width B/m . | Curvature ratio (R/B) . | ||||
B2 | 110 | 2 | 0.02 | 7.5 | 888.9 | 0.008 | 0.2 |
B5 | 110 | 2 | 0.05 | 3.0 | 2,222.2 | 0.02 | 0.2 |
B15(H60) | 110 | 2 | 0.15 | 1.0 | 6,666.7 | 0.06 | 0.2 |
Note: RUN B15 is RUN H60 in simulation of the Group II.
Contour of velocity magnitude in bends with various curvature ratios.
Contour of velocity magnitude in bends with various curvature ratios (0.1H below free surface).
Contour of velocity magnitude in bends with various curvature ratios (0.1H below free surface).
Streamwise profiles of flow velocities in bends with various curvature ratios.
The strengthening of the secondary current with increasing channel curvature can also be demonstrated by the variations of spanwise velocity in the streamwise and spanwise directions, i.e. and
, as shown in Table 6.
Magnitude of flow velocities in bends with various curvature ratios
RUN . | ||||||||
---|---|---|---|---|---|---|---|---|
B2 | 0.85 | 0.25 | 0.75 | 0.10 | 0.09 | 0.005 | 5.625 | 0.75 |
B5 | 0.70 | 0.25 | 0.60 | 0.10 | 0.07 | 0.005 | 1.80 | 0.3 |
B15(H60) | 1.25 | 0.80 | 0.75 | 0.60 | 0.188 | 0.120 | 0.75 | 0.60 |
RUN . | ||||||||
---|---|---|---|---|---|---|---|---|
B2 | 0.85 | 0.25 | 0.75 | 0.10 | 0.09 | 0.005 | 5.625 | 0.75 |
B5 | 0.70 | 0.25 | 0.60 | 0.10 | 0.07 | 0.005 | 1.80 | 0.3 |
B15(H60) | 1.25 | 0.80 | 0.75 | 0.60 | 0.188 | 0.120 | 0.75 | 0.60 |
Notice that the Reynolds number in the present study is relatively low, which ranges from 2,500 to 8,571.4 for most run cases, except RUN B2, which is 888.9. The flow falls in a regime of hydraulic smooth turbulent or laminar-turbulent transition (RUNB2). These conditions are close to laboratory experiments in mini-flumes. Although the Reynolds stress may not play an important role compared to the inertia force and centrifugal force in channel bends, the quantitative influence of the Reynolds number need to be further investigated. Considering the fact that the same uniform inlet flow velocity was adopted for all run cases in this study, the Froude number may change with the variation of aspect ratio of channel. The influence of Froude number on flow structures was also neglected in the study.
Further, in a natural meandering river, the water elevation near the concave sidewall is usually higher than that near the convex sidewall – the well-documented water super-elevation transverse slope. The super-elevation usually favors the development of the secondary current and produces a higher velocity and lower pressure region near the concave bank (Nguyen et al. 2007; De Marchis & Napoli 2008; Abad & Garcia 2009). However, this effect is usually thought to be relatively weak in the flows with a small Froude number (Stoesser et al. 2010), the rigid-lid approximation of the free-surface adopted in this study was also used by many other similar numerical studies (Van Balen et al. 2009; Stoesser et al. 2010; Bai et al. 2014) which yield acceptable agreement with experimental data.
CONCLUSIONS
By using a RANS-based 3D CFD solver for incompressible viscous flows, numerical simulations for open channel flow in curved channel bends were carried out. A series of run cases were performed and the influences of river channel sinuousness, aspect ratio and curvature ratio were studied. Simulation results show the following:
For all of the channel bends in this study (maximum deflection angle ranges from 30 to 110°, aspect ratio ranges from 0.25 to 1.0, channel curvature ratio ranges from 0.42 to 3.14), the streamwise velocity is 4–10 times larger than the spanwise velocity and the spanwise velocity is 2–10 times larger than the vertical velocity. The streamwise variation of velocity magnitude along the banks are around 50% of the mean flow velocity, while it is less than 25% along the centerline of the channel.
Highly sinuous channel bends with a maximum deflection angle of 110° show much larger cross-sectional variation of velocity compared with 60 and 30° bends. The streamwise variation ranges of velocities are of the same order of magnitude of their spanwise counterpart. However, the spanwise variation gradients can be five times higher than those in the streamwise direction. The velocity magnitude of secondary current in highly sinuous channel bends can be as large as 20% of the primary flow, which may invalidate many theoretical assumptions that secondary current is much weaker than the primary flow.
The spanwise gradient of flow velocities decreases when the aspect ratio of channel bends increases from 0.15 to 1.0. Channels with higher aspect ratios tend to present more uniform velocity distributions on cross-sections. Simulation results show that the lower the water depth, the closer the dynamic flow axis to the channel banks. The variation of streamwise velocity along channel centerlines decreases 70% when the aspect ratio increases from 0.15 to 1.0. Deep-water channels with an aspect ratio of 0.15 show significant phase-lag between flow velocities and channel geometry. The spanwise velocities near the free-surface and the bottom-wall on the center plane are found to be out-of-phase when the water depth is large, and can be attributed to the large Dean number.
When the channel curvature ratio decreases from 7.5 to 1.0, the spanwise and streamwise variation amplitude of velocities increases. For channel bends with a curvature ratio 7.5 and 3.0, the spanwise velocities near the free-surface and bottom-wall vary in the anti-phase pattern due to the Dean numbers being smaller than the critical value. However, when the curvature ratio is 1.0, the out-of-phase pattern is observed.
Note that the channel bend shapes presented in this study are controlled by three factors, i.e. sinuousness, aspect ratio and curvature ratio, and may not cover the practical meandering river shapes in nature due to their diversity. A major contribution of the present study is the flow features in channel bends with various geometries. Such knowledge is necessary for understanding the fundamental mechanism of river meandering dynamics in fluvial environment. Moreover, a systematic analysis of flow velocity magnitudes in various channel bends are presented, which is especially informative and meaningful to theoretical studies.
Nowadays, turbulent models remain a huge challenge for RANS-based CFD simulations, including the RNG k-ɛ model adopted in the present study. In this regard, further investigations on the behavior of various turbulent models and their influences on secondary current are necessary, especially for sharply curved channel bends with a strong secondary current. A more advanced numerical technique, such as LES, could be a better option. However, considering LES is computationally expensive, investigation can firstly focus on one geometric factor, such as the deflection angle, before extending to more complex channel shapes.
ACKNOWLEDGEMENTS
The study is financially supported by the Science Fund for Creative Research Groups of the National Natural Science Foundation of China (51621092), the National Natural Science Foundation of China (51279124, 51009105, 51579175).