Research Article  Open Access
Xin Chen, Zichao Zhang, Yong Li, Xiaoyan Shi, "Numerical Study for NearBed Variables in VelocitySkewed Oscillatory Sheet Flow Transport", Mathematical Problems in Engineering, vol. 2017, Article ID 2614943, 10 pages, 2017. https://doi.org/10.1155/2017/2614943
Numerical Study for NearBed Variables in VelocitySkewed Oscillatory Sheet Flow Transport
Abstract
The nearbed sediment concentration and vertical sediment flux are important in sediment transport mechanics, but they are known much less than the horizontal sediment flux, especially for fine sediment in unsteady flows. A developed twophase model is applied to study the nearbed sediment concentration, vertical sediment flux, and the relevant total sediment amount for the velocityskewed oscillatory sheet flow transport. With the sediment concentration hindered fall velocity, the classical reference concentration formulas conducted by Engelund and Fredsoe (1976) and Zyserman and Fredsoe (1994) are utilized for the comparison with the twophase model and illustration of the phaselag and sediment size effects in nearbed sediment concentration and vertical sediment flux. The concentration and vertical flux predicted by the twophase model agree well with experimental data and are better than empirical formulas. Furthermore, the sediment size effect for pickup flux function over starved bed is shown to be quite different from that containing sufficient sediment in oscillatory flows.
1. Introduction
The phenomenon of sediment entrainment and settlement is universal in river and coastal environment and extraordinarily important in the subaquatic geomorphology. The vertical sediment fluxes consist of pickup flux due to turbulence, settling flux due to gravity, and their summation denoting the variation rate of total sediment amount. Associated with the sediment concentration, they represent the vertical sediment exchange strength and are important for the study of sediment settling down to deposition and picked up into suspension. Most sediment models need a type of nearbed condition associated with vertical sediment flux for bottom boundary conditions. So the predictions of nearbed concentration, vertical sediment flux, pickup flux, and settling flux are very fundamental for such sediment models.
In fully developed steady flows, the pickup flux and settling flux are in balance and the total vertical flux is zero. Numerous steady flow formulas have been carried out to predict the vertical sediment flux over loose sediment bed in developed equilibrium steady flows. Generally, some are straight driven by the pickup function based on clear water flow over starved bed without considering the influence of suspended sediment [1, 2], so they are not appropriate in the flow already containing large amount of suspended sediment. The others are evaluated by the balance of pickup and settling down according to the concentration and sediment fall velocity [3], and the relevant concentration denotes the vertical sediment flux strength.
The nearbed sediment concentration and vertical sediment fluxes in unsteady flows are much more complicated than that in steady flows due to the phaselag effect, and formulas in unsteady flows are lacking. So the nearbed concentrations denoting vertical flux strength in steady flows [4–6] are applied and widely used in unsteady numerical models for sediment transport [7–9], and even their concentrations drop to zero near the flow reversal when the bottom shear stress is under a threshold. The knowledge of vertical sediment flux in unsteady flow is still necessary for numerical models’ boundary conditions in engineering.
In unsteady flows, the vertical sediment flux can be derived from the variation rate of sediment amount in movement [10], which can be obtained by single phase models or twophase models. The single phase numerical models are dependent on the reference concentration or vertical flux as the bottom boundary conditions and cannot describe the vertical flux below reference height. The twophase models [11–13] can obtain the details of reference concentration and vertical flux as a result of the turbulent upward effect and gravitational settlement. Yu et al. [12] study the vertical flux for medium and coarse sediment, but the model is not valid for fine sediment (diameter mm) when the phaselag is obvious due to the sediment stress. In addition, the vertical sediment fluxes have not been extracted accurately enough due to the noise or error of modelling high concentration. Very little is known about the nearbed concentration and vertical flux for fine sediment and the phaselag effects, which are important in the generation of net sediment transport, especially in velocityskewed oscillatory flows with a sharp wave crest and gentle wave trough.
Present study utilizes a twophase model to study the nearbed sediment concentration and vertical fluxes in velocityskewed oscillatory sheet flows containing sufficient amount of sediment. The results mainly exhibit and explain the phaselag characteristics in concentration and vertical flux at the reference height, nearimmobile bed, and very low concentration area for difference sediment sizes.
2. Applied Models
Chen and Yu [13] twophase model is utilized to obtain the sediment concentration and vertical fluxes. For comparison, several steady flow empirical formulas are also applied.
2.1. TwoPhase Flow Model
The twophase model includes Reynolds averaged continuity and momentum equations:where is the sediment volumetric concentration; is the velocity; is the Cartesian coordinate; is the time; the subscripts and stand for the fluid and the sediment phases; the indices represent horizontal and vertical components and obey summation convention; is the sediment turbulent viscosity confirmed by fluid turbulent viscosity ; is the Schmidt number; is the pressure; is the stress; is the density; is the body force; is the interaction force. , , , and can be seen in Chen and Yu [13]. Equations (1) are closed by  turbulent equations:where is the turbulent kinetic energy (TKE); is the turbulent kinetic energy dissipation rate (DTKE); , /1.33, = 1.44, and = 1.92; is the turbulence generation term [13].
Equations (1) and (2) are applied for sheet flows happening at maximum Shields parameter > 0.8–1.0 when large amounts of sediment are in movement. The computation domain covers the immobile bed and the water tunnel top. The initial bed is located at = 0. The still water condition is assumed to be the initial condition. At the time , all flow variables are zero, while at > 0 and at ≤ 0.
At the tunnel top, , where is the free stream (Figure 1); ; and are the amplitudes of and ; is the period; subscripts and denote acceleration and deceleration stages; scripts and denote the crest and trough durations. Positive ( = 0–0.42) and negative ( = 0.42–1) denote onshore and offshore durations, respectively. The vertical gradient of all variables vanished, and the sediment flux is
Nonslip condition is assumed at the bottom that and the other variables are set to be zero.
In the lateral boundaries, the horizontal gradients of , , , and are all zero. The flow is driven by a horizontal pressure gradient:
2.2. Flux Definition and Empirical Formulas
The sediment amount in movement is defined by Nielsen et al. [10]:
The total vertical flux isand the downward settling flux iswhere is the fall velocity equal to . The pickup flux isThe dimensionless equations (10)–(12) are , where the specific gravity . In the twophase model, and are straight given by (1), and .
For comparison, the classical Engelund and Fredsoe [4] formula (EF76) and Zyserman and Fredsoe [5] formula (ZF94) for reference defined at and are applied. is the single sediment fall velocity in still water [14] and is the same as Yu et al. [12]. The EF76 formula iswhere is the sediment diameter; is the linear concentration; is the Karman constant; is associated with the bed roughness; is the critical Shields parameter; is the dynamic friction coefficient. The ZF94 formula is where . Figure 2 shows the reference against for the two empirical formulas. The EF76 considers the dispersive stress from suspended sediment and drag of moving bed, so the reference is mildly increased when and is hindered mildly to be a constant. ZF94 is based on Garcia and Parker’s [3] prohibiting a particular trend toward a limitation for high sediment concentration. So the ZF94 also has a limitation of and keeps a mild increment after (Figure 2).
For comparison, van Rijn [1] pickup function over a starved bed is applied:where is the fluid molecular viscosity. in EF76 and ZF94 formulas and in van Rijn [1] drop to zero when without phaselag effect.
3. Results
Computations are based on the O’Donoghue and Wright [15] experiments in Table 1, where the phaselag effect increases with the decrement of . Results are shown for the classical reference height first and second the nearimmobile bed and low concentration areas.

3.1. Reference Vertical Flux
The references α and are shown in Figure 3. The absolute locations of are slightly different and are very close to each other, so they can be treated as an approximation to the initial bed. The distribution pivot, around which is almost constant in a wave period, exists at the location a little below the initial bed [15]. In Figures 3(a), 3(c), and 3(e), the location is very close to the pivots, and both the maximum free stream flow and flow reversal would probably lead to peaks. For the experiment (dot), the difference between flow crest () and flow trough () is not obvious and even the crest flow strength is much larger than that at flow trough (Figure 3), because many sediments picked up at the flow crest are still in movement at the flow trough. Large corresponds to small concentration magnitude due to the large fall velocity. never decrease to zero in a period due to the phaselag effect. The results of twophase model exhibit clearly the decrement of concentration magnitude when increases, and the computed is also always above a proper nonzero value (Figures 3(a), 3(c), and 3(e)) due to the phaselag. Some differences between twophase model and experimental data exist from a qualitative perspective, probably due to the high concentration measurement uncertainties [15] without mass conservation.
(a) of FA5010
(b) of FA5010
(c) of MA5010
(d) of MA5010
(e) of CA5010
(f) of CA5010
The EF76 and ZF94 are strictly in phase with without any phaseshift. Mild increment at large shear stress (Figure 2) leads to a small concentration variation near onshore flow peak around (Figures 3(a) and 3(c)). There are also small concentration variations for EF76 and ZF94 near the offshore flow peak ( = 0.6~0.8) caused by the nearly constant U (Figure 1). When is less than a critical level, EF76 and ZF94 predicted are unrealistic zero (Figure 2), because the formulas are derived without the phaselags. Predicted concentrations of EF76 and ZF94 increase rapidly in the acceleration stage and decrease fast in the deceleration stage (Figures 3(a), 3(c), and 3(e)). Furthermore, their prediction between onshore and offshore are quite different because of the large different flow strength and bottom shear stress. Although both EF76 and ZF94 can predict the experiments near the flow trough ( = 0.6~0.8), they much overestimate the experiments near the flow crest (). Generally, the EF76 and ZF94 are weaker than the twophase model results at the flow reversal and flow crest. RMS errors in Table 2 for twophase model are much less than EF76 and ZF94 formulas.

Corresponding to the concentration due to phaselag, the experimental never decreases to zero near the flow reversal in Figure 3, and the value near the flow crest is close to that near the flow trough. The twophase model predicted computed agree very well with the experiments for nonzero duration and the magnitude increment with . Although the magnitude decreases with the increment of D, the magnitude of increases with the increment of in Figure 3, because the flux is hindered by high concentration at small sediment size. Based on the same = 1.5 m/s, the pickup flux () against for ZF94 and ZF76 in steady flow is shown in Figure 4, where the dimensionless is represented by /. The flux increases first and then decreases with D (/). The magnitude would decrease with small when is large enough as the balanced pickup flux () in Figure 4, and finally it decreases to zero when reference magnitude reduces to zero.
Near the flow crest, there are hollows in EF76 and ZF94 that predicted which are far from the experiments (Figures 3(b), 3(d), and 3(f)) with overestimation of (Figures 3(a), 3(c), and 3(e)), because the maximum and optimal exist in the two empirical formulas when we applied , as seen in Figure 5. The maximum approximately corresponds to optimal , after which decreases with the increment of by the hindered fall velocity. Near the flow trough, both EF76 and ZF94 can predict well due to the agreement of concentration. At flow reversal, both the EF76 and EF94 predicted drop to zero because of unrealistic zero (Figures 3(a), 3(c), and 3(e)) when is less than a critical value. Predicted of EF76 and ZF94 changes rapidly near flow reversal with rapid change of . Generally, the RMS errors for EF76 and ZF94 are much larger than twophase model for (Table 2).
Against the settling down induced by gravity, the pickup effect induced by turbulence maintains the suspension of sediment, and the total vertical sediment flux () is the balance of pickup flux () and settling flux (). Figure 6 shows the twophase model computed and at and the immobile bed surface level denoted by (). To illustrate their relationships availably and validate indirectly, in Figure 3 is included again. Magnitudes of and are around each other with opposite direction, and they are in periodic averaged balance with averaged . As the difference of and , is almost one magnitude smaller than and . increases and reaches a local maximum value before flow crest () or flow trough () and decreases to a local minimum value before flow reversal both in the onshore and in the offshore durations. Correspondingly, almost follows the same variation tendency as . Similar process can be seen in Ribberink et al. [16], Nielsen et al. [10], Li et al. [17], and Bakhtyar et al. [11]. Onshore magnitudes of and are larger than that at offshore, and the onshore variations are more rapid than offshore.
(a) at , FA5010
(b) , FA5010
(c) at , MA5010
(d) , MA5010
(e) at , CA5010
(f) , CA5010
almost varies together with which decided the bed morphology in the present three cases in Figure 6. The bed morphology repeats periodically in oscillatory sheet flow conditions without horizontal variation. When , decreases because sediment below the initial bed is picked up. decreases to minimum (t/T = 0.25–0.3) almost after decreases to zero. While , increase because sediment above the initial bed settles down. increases to local maximum (t/T = 0.45–0.5) almost after increases to zero. never drops to zero due to the phaselag effect. For comparison, is shown considering theoretical in steady flow, where is the maximum . The difference between and denotes the phaselag. There is a time phaseshift between and . For fine 0.13 mm, the phaselag effect is the largest with the most obvious difference between and , while the phaselag effect for 0.46 mm is the smallest because present almost coincides with .
Corresponding to , the magnitude also increases with the increment of for the three sizes. But it does not increase with the increment of monotonously, which is similar to Figure 4. In the twophase model, is obtained from the turbulence diffusion term ) modelled by the classical gradient transport assumption. There are two reasons for the smallest amplitude with fine D = 0.13 mm in Figure 6: the concentration gradient and the turbulence intensity. With D = 0.13 mm, much sediment is picked up for suspension and at becomes smaller than D = 0.27–0.46 mm. When much sediment is already in suspension in the boundary layer, much turbulence energy is dissipated [15, 18] corresponding to small and restrains other sediments to be picked up.
For comparison without the effect of suspended hindered sediment, van Rijn [1] pickup flux function is shown in Figure 7 under the same = 1.5 m/s. van Rijn [1] pickup function decreases with the increment of , which is opposite to the flows containing sufficient sediment amount in Figures 4 and 5. The onshore magnitude is much larger than offshore magnitude without the effect of suspended hindered sediment. In the experiment [15], did not drop to zero due to the shear (velocity overshoot) and turbulence that exist near flow reversal. But van Rijn [1] function based on steady flow drops to zero around flow reversal. Comparing the van Rijn pickup function (Figure 7) with given by twophase model (Figure 6), their magnitudes are the closest only for = 0.46 mm when the suspension sediment amount is the smallest (Figure 8). Over all, the van Rijn [1] pickup function is not valid for the oscillatory sheet flow when the maximal > 0.8–1.0 and sufficient sediment is in suspension.
In developed oscillatory sheet flow, the averaged is zero, and suspended sediment is quasisteady saturated and hinders the pickup process. The vertical fluxes display the sediment exchange in vertical direction, but they can not be used as the index for the sedimentladen capacity of oscillatory sheet flow. The sediment amount is in a quasisteady status in the developed stage and can be used to estimate the sedimentladen capacity. above is shown in Figure 8, where = 0.5 m is the water depth. increases with the decrement of . The variation rate of is (see (6)), so reaches a local maximum value when positive decreases to zero near the flow peak and reaches a local minimum value when negative increases to zero near the flow reversal. Due to large phaselag, the magnitude’s ratio between onshore ( = 0–0.42) and offshore ( = 0.42–1) for 0.13 mm case is 0.87, much larger than 0.41 for 0.46 mm case.
3.2. Other Vertical Fluxes
The nearbed flux at the locations above or below the reference height is essential for better understanding of the nearbed sediment exchange. Figure 9 shows the twophase model computed vertical fluxes at = 10 mm and = −3.5 mm. −3.5 mm is in the main eroded area with high concentration and almost fully hindered . 10 mm is out of the sheet flow layer where particle interactions can be negligible and is still in the oscillatory boundary layer with small concentration and unhindered . Because the absolute locations of in the main eroded area or low concentration area for three sizes are far from each other, the vertical locations in Figure 9 are in absolute millimetres rather than as a function of the diameter. At mm, , , and decrease with the increment of . Two peak values exist in Figure 9(a) and onshore magnitude is larger than offshore, but only one peak value exists near the onshore flow crest in Figure 9(c), while in Figure 9(e) almost zero fluxes exist because −3.5 mm is nearly located at the immobile bed for D = 0.46 mm. At = 10 mm, obvious phases for all the fluxes are shifted after . Both the maximum and are around 0.01, which are much larger than . For fine and medium sediment (Figures 9(b) and 9(d)), their magnitudes of , , and are similar to and larger than coarse sediment (Figure 9(f)). For coarse sediment, both and only exist near the onshore flow peak, and almost disappear compared to the fine and medium sediment cases. Comparing the results in Figure 9 to that in Figure 6, all , , and decrease from the initial bed toward the immobile bed and the top of tunnel. The decrement is also seen in Li et al. [17] and Bakhtyar et al. [11]. All the fluxes finally decrease to zero when the location is down into the immobile bed or up to the top of tunnel.
(a) FA5010, = −3.5 mm
(b) FA5010, = 10 mm
(c) MA5010, = −3.5 mm
(d) MA5010, = 10 mm
(e) CA5010, = −3.5 mm
(f) CA5010, = 10 mm
Also, above the locations −3.5 mm and 10 mm are shown (Figure 10). At = −3.5 mm and = −10 mm, the variation amplitudes of each sediment sizes are much smaller than that at (Figure 8) due to smaller in Figure 9. Finally, the variation amplitudes decrease to zero in the immobile bed or tunnel top corresponding to (Figure 9).
(a) Above = −3.5 mm
(b) Above = 10 mm
4. Conclusions
The present work studies the nearbed conditions in velocityskewed oscillatory sheet flows with a twophase model and classical empirical formulas. The phaselag and sediment size effects in concentration, vertical sediment flux, and sediment amount are focused on the area near the initial bed, the immobile bed, and the low concentration area.
Due to the phaselag, the reference concentration and the settling flux never drop to zero in a period even at the flow reversal, and the relevant magnitudes between onshore and offshore are similar. The twophase model computed reference concentration and settling flux by twophase model agree with experiment very well and are better than empirical formulas. Under the same flows, reference concentration magnitudes decrease with the increment of sediment size due to the increment of fall velocity. While, for settling flux, as the product of concentration and fall velocity, the magnitude increases first and then decreases with the increment of sediment size.
The pickup flux and settling flux are in periodic averaged balance in oscillatory sheet flow, and their magnitudes are around each other with opposite direction. The pickup flux increases and reaches a local maximum value before flow peak and decreases to a local minimum value before flow reversal both in the onshore and in the offshore durations. Corresponding to the settling flux, the pickup flux also increases first and then decreases with the increment of sediment size.
The sediment amount magnitude increases with the decrement of sediment size. Phaselag effect increases the sediment amount magnitude’s ratio between onshore and offshore. As the sediment amount variation rate, the total vertical flux is much smaller than the pickup flux and varies almost following the pickup flux. The vertical sediment flux magnitude decreases from the initial bed toward immobile bed and top of tunnel.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Grants no. 51609244 and no. 11472156) and National ScienceTechnology Support Plan of China (Grant no. 2015BAD20B01).
References
 L. C. van Rijn, “Sediment pickup functions,” Journal of Hydraulic Engineering, vol. 110, no. 10, pp. 1494–1502, 1984. View at: Publisher Site  Google Scholar
 D. Zhong, G. Wang, and Y. Ding, “Bed sediment entrainment function based on kinetic theory,” Journal of Hydraulic Engineering, vol. 137, no. 2, pp. 222–233, 2011. View at: Publisher Site  Google Scholar
 M. Garcia and G. Parker, “Entrainment of bed sediment into suspension,” Journal of Hydraulic Engineering, vol. 117, no. 4, pp. 414–435, 1991. View at: Publisher Site  Google Scholar
 F. Engelund and J. Fredsoe, “A sediment transport model for straight alluvial channels,” Nordic Hydrology, vol. 7, pp. 295–324, 1976. View at: Google Scholar
 J. A. Zyserman and J. Fredsoe, “Data analysis of bed concentration of suspended sediment,” Journal of Hydraulic Engineering, vol. 120, no. 9, pp. 1021–1042, 1994. View at: Publisher Site  Google Scholar
 S. Wright and G. Parker, “Flow resistance and suspended load in sandbed rivers: Simplified stratification model,” Journal of Hydraulic Engineering, vol. 130, no. 8, pp. 796–805, 2004. View at: Publisher Site  Google Scholar
 M. Li, P. T. Fernando, S. Pan, B. A. O'Connor, and D. Chen, “Development of a quasi3d numerical model for sediment transport prediction in the coastal region,” Journal of HydroEnvironment Research, vol. 1, no. 2, pp. 143–156, 2007. View at: Publisher Site  Google Scholar
 B. G. Ruessink, T. J. J. Van Den Berg, and L. C. Van Rijn, “Modeling sediment transport beneath skewed asymmetric wavesabove a plane bed,” Journal of Geophysical Research: Oceans, vol. 114, no. 11, article C11021, 2009. View at: Publisher Site  Google Scholar
 W. N. M. Hassan and J. S. Ribberink, “Modelling of sand transport under wavegenerated sheet flows with a RANS diffusion model,” Coastal Engineering, vol. 57, no. 1, pp. 19–29, 2010. View at: Publisher Site  Google Scholar
 P. Nielsen, K. Van Der Wal, and L. Gillan, “Vertical fluxes of sediment in oscillatory sheet flow,” Coastal Engineering, vol. 45, no. 1, pp. 61–68, 2002. View at: Publisher Site  Google Scholar
 R. Bakhtyar, A. YeganehBakhtiary, D. A. Barry, and A. Ghaheri, “Twophase hydrodynamic and sediment transport modeling of wavegenerated sheet flow,” Advances in Water Resources, vol. 32, no. 8, pp. 1267–1283, 2009. View at: Publisher Site  Google Scholar
 X. Yu, T.J. Hsu, J. T. Jenkins, and P. L.F. Liu, “Predictions of vertical sediment flux in oscillatory flows using a twophase, sheetflow model,” Advances in Water Resources, vol. 48, pp. 2–17, 2012. View at: Publisher Site  Google Scholar
 X. Chen and X. Yu, “A numerical study on oscillatory flowinduced sediment motion over vortex ripples,” Journal of Physical Oceanography, vol. 45, no. 1, pp. 228–246, 2015. View at: Publisher Site  Google Scholar
 L. C. van Rijn, Principles of Sediment Transport in Rivers, Estuaries and Coastal Seas, Aqua Publications, 1993.
 T. O'Donoghue and S. Wright, “Concentrations in oscillatory sheet flow for well sorted and graded sands,” Coastal Engineering, vol. 50, no. 3, pp. 117–138, 2004. View at: Publisher Site  Google Scholar
 J. S. Ribberink, C. M. DohmenJanssen, D. M. Hanes, S. R. McLean, and C. Vincent, “Nearbed sand transport mechanisms under waves,” in Proceedings of the 27th International Conference on Coastal Engineering, ICCE 2000, pp. 3263–3276, Sydney, Australia, July 2000. View at: Publisher Site  Google Scholar
 M. Li, S. Pan, and B. A. O'Connor, “A twophase numerical model for sediment transport prediction under oscillatory sheet flows,” Coastal Engineering, vol. 55, no. 12, pp. 1159–1173, 2008. View at: Publisher Site  Google Scholar
 K. C. Wilson, “Friction of waveinduced sheet flow,” Coastal Engineering, vol. 12, no. 4, pp. 371–379, 1989. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Xin Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.