#### Abstract

The near-bed 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 two-phase model is applied to study the near-bed sediment concentration, vertical sediment flux, and the relevant total sediment amount for the velocity-skewed 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 two-phase model and illustration of the phase-lag and sediment size effects in near-bed sediment concentration and vertical sediment flux. The concentration and vertical flux predicted by the two-phase model agree well with experimental data and are better than empirical formulas. Furthermore, the sediment size effect for pick-up 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 pick-up 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 near-bed condition associated with vertical sediment flux for bottom boundary conditions. So the predictions of near-bed concentration, vertical sediment flux, pick-up flux, and settling flux are very fundamental for such sediment models.

In fully developed steady flows, the pick-up 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 pick-up 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 pick-up and settling down according to the concentration and sediment fall velocity [3], and the relevant concentration denotes the vertical sediment flux strength.

The near-bed sediment concentration and vertical sediment fluxes in unsteady flows are much more complicated than that in steady flows due to the phase-lag effect, and formulas in unsteady flows are lacking. So the near-bed 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 two-phase 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 two-phase 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 phase-lag 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 near-bed concentration and vertical flux for fine sediment and the phase-lag effects, which are important in the generation of net sediment transport, especially in velocity-skewed oscillatory flows with a sharp wave crest and gentle wave trough.

Present study utilizes a two-phase model to study the near-bed sediment concentration and vertical fluxes in velocity-skewed oscillatory sheet flows containing sufficient amount of sediment. The results mainly exhibit and explain the phase-lag characteristics in concentration and vertical flux at the reference height, near-immobile bed, and very low concentration area for difference sediment sizes.

#### 2. Applied Models

Chen and Yu [13] two-phase model is utilized to obtain the sediment concentration and vertical fluxes. For comparison, several steady flow empirical formulas are also applied.

##### 2.1. Two-Phase Flow Model

The two-phase 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 pick-up flux isThe dimensionless equations (10)–(12) are , where the specific gravity . In the two-phase 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] pick-up 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 phase-lag effect.

#### 3. Results

Computations are based on the O’Donoghue and Wright [15] experiments in Table 1, where the phase-lag effect increases with the decrement of . Results are shown for the classical reference height first and second the near-immobile 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 phase-lag effect. The results of two-phase 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 phase-lag. Some differences between two-phase 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 CA5010The EF76 and ZF94 are strictly in phase with without any phase-shift. 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 phase-lags. 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 two-phase model results at the flow reversal and flow crest. RMS errors in Table 2 for two-phase model are much less than EF76 and ZF94 formulas.

Corresponding to the concentration due to phase-lag, 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 two-phase 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 pick-up 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 pick-up 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 two-phase model for (Table 2).

Against the settling down induced by gravity, the pick-up effect induced by turbulence maintains the suspension of sediment, and the total vertical sediment flux () is the balance of pick-up flux () and settling flux (). Figure 6 shows the two-phase 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 phase-lag effect. For comparison, is shown considering theoretical in steady flow, where is the maximum . The difference between and denotes the phase-lag. There is a time phase-shift between and . For fine 0.13 mm, the phase-lag effect is the largest with the most obvious difference between and , while the phase-lag 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 two-phase 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] pick-up flux function is shown in Figure 7 under the same = 1.5 m/s. van Rijn [1] pick-up 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 pick-up function (Figure 7) with given by two-phase 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] pick-up 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 quasi-steady saturated and hinders the pick-up process. The vertical fluxes display the sediment exchange in vertical direction, but they can not be used as the index for the sediment-laden capacity of oscillatory sheet flow. The sediment amount is in a quasi-steady status in the developed stage and can be used to estimate the sediment-laden 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 phase-lag, 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 near-bed flux at the locations above or below the reference height is essential for better understanding of the near-bed sediment exchange. Figure 9 shows the two-phase 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 mmAlso, 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 near-bed conditions in velocity-skewed oscillatory sheet flows with a two-phase model and classical empirical formulas. The phase-lag 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 phase-lag, 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 two-phase model computed reference concentration and settling flux by two-phase 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 pick-up flux and settling flux are in periodic averaged balance in oscillatory sheet flow, and their magnitudes are around each other with opposite direction. The pick-up 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 pick-up flux also increases first and then decreases with the increment of sediment size.

The sediment amount magnitude increases with the decrement of sediment size. Phase-lag 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 pick-up flux and varies almost following the pick-up 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 Science-Technology Support Plan of China (Grant no. 2015BAD20B01).