Computational Options

The General tab of the Sediment Computation Options and Tolerances editor contains the 1D computational options but also includes some options which the 1D and 2D sediment models share (see figure below). The 2D sediment transport model uses the Bed Roughness Predictor, the Computation Multiplier and Warmup Periods.


Figure 1. General tab of the Sediment Computation Options and Tolerances editor.

Bed Roughness Predictor

The Bed Roughness Predictor computes bed roughness dynamically based on the hydrodynamics and bed materials. The sediment transport model and hydraulic models use the computed roughness, then update it based on the new hydraulics and bed materials. The 1D model allows users to define sub-reaches to apply the bed roughness predictor, but the 2D model will apply the selected algorithm to the all wet cells. When using bed-roughness predictor it is recommended to also output the cell bed roughness's and evaluate the computed bed roughness to maker sure results are realistic. The Bed Roughness Predictor is specified in the General Tab of the Sediment Computation Options and Tolerances editor. 

The Bed Roughness Predictor options include three methods:

  1. Limerinos
  2. Brownlie
  3. Van Rijn

Limerinos

The Limerinos (1970) formula is a relatively simple equation based on the d84 and hydraulic radius, with no consideration for bed form mechanics. Therefore, it should be applied mainly to coarse, gravel and cobble systems where the grain roughness is the primary source of channel roughness.

Brownlie

Brownlie (1983) computed bed roughness based on bed form mechanics in large rivers. In particular, he tried to capture the non-linear drop in roughness when bed form dominated transport shifts from lower to higher regimes (EM 4000). Brownlie, evaluates the bed form regime based on hydraulic parameters, grain size, and the gradational distribution, then applies separate equations for low and high regime transport.

Van Rijn

The Van Rijn (1984) bed roughness formula is based on flume and field data. Van Rijn computes bed form dimensions from flow and sediment parameters and converts these into equivalent bed roughness. This method was designed to compute bed roughness in both dune and plane bed regimes. 

Sediment Computation Multiplier

The Sediment Computation Multiplier is used to compute sediment computations at a larger time step than the hydraulic computations. The Sediment Computation Multiplier is an integer which defines the number of hydraulic time steps before computing a sediment time step. A 2D sediment time step calculation can be significantly more computationally expensive than the hydraulic time, especially if there are many grain classes. It is important to note that the actual number of hydraulic time steps between sediment time steps may differ from the Sediment Computation Multiplier at certain times during the simulation in order for the model to synchronize time steps with output intervals. When the multiplier is utilized, the face flows and areas are averaged during the sediment time step. The parameter must be used with caution since it can reduce the model accuracy or produce instabilities. Typical values that work well are in the range of 5 to 20. Testing should be be done to confirm that the the parameter does not significantly affect the computed results. The Sediment Computation Multiplier has a default value of 1 meaning the sediment and hydraulics are computed at the same computational time step. The Sediment Computation Multiplier is specified in the General Tab of the Sediment Computation Options and Tolerances editor. 

Warmup Method

HEC-RAS applies sediment warmup period slightly differently in 1D and 2D sediment transport. The 1D model adds the warmup period to the simulation, as an additional (negative) time period computed before the simulation starts. In 2D the warmup period is considered as part of the simulation time period.
The 2D model also computes three sediment warmup periods, in sequence (after the hydraulic warmup). First, HEC-RAS computes equilibrium concentrations in each cell based on the first specified flow or stage for the first time step at each hydraulic boundary, the boundary sediment fluxes, and the bed gradations. The Concentration warmup is only available in 2D sediment transport, and avoids the initial erosion associated with starting a model with clear water throughout the domain.


Figure 2. Specifying the warmup periods for 2D simulations.

After the concentration warmup, the 2D model warms up the bed gradations. Starting with the cell concentrations at the end of the concentration warmup, HEC-RAS will run the unsteady, 2D, sediment model with constant flow, stage, and sediment boundary conditions (those associated with the first time step in the simulation window). During the gradation warmup, HEC-RAS will not adjust the bed, the elevation-volume curves for the cells, or the elevation-area curves of the cell faces. It will only adjust the bed gradations, allowing them to fine or (more often) coarsen in response to the flow field. The 2D sediment model will also begin to develop vertical gradational stratigraphy in this warmup period, even bed layers were not initially specified.

After the Gradation warmup, the 2D model will warmup the bathymetry. In this warmup period, the model will still hold they hydraulic and sediment boundary conditions constant, but will allow the bed elevations, and the elevation-volume/elevation-area curves to adjust in response to transport. This phase will tend to fill local depressions, scour local raised areas, plane out bed forms, highlight problems with the terrain, and, generally smooth the terrain.