Total-Load Adaptation Length

There are several options to calculate the total-load adaptation length. The simplest option is to utilize a constant adaptation length calibrate the parameter using field measurements. Experience has shown that for most field applications, this approach is sufficient. However, several other methods are available with varying degrees of complexity. The total-load adaptation length may also be determined.
The advection coefficient considers several processes which affect the velocity and direction of the sediment transport including the vertically non-uniform distribution of the horizontal current velocity and sediment concentration profiles and bed-slope effects.
There are four methods for calculating the total load adaptation coefficient in HEC-RAS:

  1. Constant adaptation length
  2. Weighted averaged of bed- and suspended load adaptation lengths

In the case that a constant adaptation length is specified, the total-load adaptation coefficient is computed as

1) \alpha_t = \frac{hU}{L_t\omega_s}

where

h: water depth [L]

U: Current velocity [L/T]

ω_{s}: sediment settling velocity [L/T]

L_t: total-load adaptation length [L]


The fall velocity is equal to the transport grain size, for single size sediment transport, or the median grain size, in the case of multiple-sized sediment transport.

A more realistic approach to use a transport weighted average of the bed- and suspended-load adaptation lengths (Wu 2007)

2) \alpha_t = r_s \alpha_s + (1-r_s) \frac{hU}{L_b\omega_s}

where  

\alpha_s: suspended-load adaptation coefficient [-]

L_b : bed-load adaptation length [L]

 r_s : fraction of suspended load of the total load.


The methods for determining  and  are described in subsequent sections.

Bed-load Adaptation Length

There are two options for the bed-load adaptation length. The first is for the user to specify a constant value. There are two methods for calculating the total load adaptation coefficient in HEC-RAS:

  1. Constant bed-load length
  2. Depth-dependent method

In the case of the depth-dependent method, the adaptation length is computed as

3) L_b = f_{bL} h

where  

f_{bL}: depth-dependent factor approximately between [-]

h: water depth [L]

Suspended-load Adaptation Coefficient

There are three options for computing the suspended-load adaptation coefficient:

  1. Constant coefficient
  2. Zhou and Lin (1998)
  3. Armanini and Di'Silvio

The figure below shows two examples of the two formulations for the suspended-load adaptation coefficient.

Figure 1. Suspended-load adaptation coefficient. 

In the case of the Zhou and Lin (1998) formulation, the adaptation coefficient is computed by solving a 2DV advection-diffusion equation of suspended load with constant diffusivity and uniform flow. A concentration bottom boundary condition is specified for the erosion case, and a gradient boundary condition for the deposition case. The analytical solution is given by

4) \alpha_s = \frac{R}{4} + \frac{\sigma_1^2}{R}

where  

R=6 \omega_s / (\kappa u_*): suspendesion parameter [-]

ω_{s}: sediment settling velocity [L/T]

u_*: total bed shear velocity [L/T]

\sigma_1: first positive root of the following equations:

\tan \sigma = - \sigma / R   for erosion

2 \cot \sigma = 2 \sigma - R/2  for deposition

A plot of the suspended-load adaptation coefficient following Zhou and Lin (1998) is shown in the figure below. The adaptation coefficients for erosion and deposition can be very different, especially for low suspension parameter values. In practice, a weighted average of the erosion and deposition values is used in order to avoid stability issues. It is also noted that the adaptation coefficient of Zhou and Lin (1998) is always greater than 1.

Armanini and Di’Silvio (1986) proposed the following formula for the adaptation coefficient

5) \frac{1}{\alpha_s} = \frac{a}h} + \left ( 1 - \frac{a}{h} \right ) \exp \left [ -1.5 \left ( \frac{a}{h}\right )^{1/6} \frac{\omega}{u_*} \right ]

where

 ω_{s}: sediment particle settling velocity [L/T]

a=33 z_0 : thickness of the bottom layer [L]

 z_0: roughness length [L]

 u_*: total bed shear velocity [L/T]

 h : water depth [L]