User Guide 012: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 133: Line 133:


[[File:fig_2-89.png]]
[[File:fig_2-89.png]]
Figure 2-89. Bed layer block section within the ''Sediment'' tab of the ''CMS-Flow Model Control'' window in SMS 11.1.


A bed layer thickness and composition is defined with a bed layer block. The cards used to define the beginning and end of a bed layer block are described in the table below.  
A bed layer thickness and composition is defined with a bed layer block. The cards used to define the beginning and end of a bed layer block are described in the table below.  
Line 160: Line 162:
The first two options assume a log-normal grain size distribution
The first two options assume a log-normal grain size distribution


{{Equation|<math>\rho_k = \frac{\Delta d_k}{d_k \ ln \ \sigma_g \sqrt{2 \pi}} exp \left [- \frac{1}{2} \left(\frac{ln \ d_k - \ ln \ d_g}{ln \ \sigma_g}  \right)^2 \right]</math>|2-34}}
{{Equation|<math>\rho_k = \frac{\Delta d_k}{d_k \ ln \ \sigma_g \sqrt{2 \pi}} exp \left [- \frac{1}{2} \left(\frac{ln \ d_k - \ ln \ d_g}{ln \ \sigma_g}  \right)^2 \right]</math>|2-33}}




Line 170: Line 172:
It is noted that the geometric standard deviation is dimensionless.  
It is noted that the geometric standard deviation is dimensionless.  


The geometric mean <d<sub>g</sub>  is in given by  
The geometric mean d<sub>g</sub>  is in given by  


{{Equation|<math>d_g = exp \left(\Sigma \rho_k ln \ d_k  \right)</math>|2-35}}
{{Equation|<math>d_g = exp \left(\Sigma \rho_k ln \ d_k  \right)</math>|2-35}}
Line 267: Line 269:
{{Equation|<math>\xi_k = \left(\frac{P_{ek}}{P_{hk}}  \right)^{-m}</math>|2-40}}
{{Equation|<math>\xi_k = \left(\frac{P_{ek}}{P_{hk}}  \right)^{-m}</math>|2-40}}


where ''m''  is an empirical coefficient that varies for each transport formula, approximately equal to 0.6-1.0. P<sub>ek</sub>  and p<sub>hk</sub>  are the total hiding and exposure probabilities and are calculated as  
where ''m''  is an empirical coefficient that varies for each transport formula, approximately equal to 0.6-1.0. P<sub>ek</sub>  and P<sub>hk</sub>  are the total hiding and exposure probabilities and are calculated as  


{{Equation|<math>P_{hk} = \sum_{j=1}^{N}p_{1j} \frac{d_j}{d_k + d_j} \quad \quad P_{ek}= \sum_{j=1}^{N} p_{1j} \frac{d_k}{d_k + d_j} </math>|2-41}}
{{Equation|<math>P_{hk} = \sum_{j=1}^{N}p_{1j} \frac{d_j}{d_k + d_j} \quad \quad P_{ek}= \sum_{j=1}^{N} p_{1j} \frac{d_k}{d_k + d_j} </math>|2-41}}

Revision as of 19:10, 4 May 2015

Fig 2-87.png

Figure 2-87. Example of the total load adaptation length specification within the Adaptation section within the Sediment tab of the CMS-Flow Model Control window in SMS 11.1.

example 2-80

Hard Bottom Specification

Although the hard bottom dataset can be edited (when it’s the active da-taset) by selecting a cell (or group of cells) and changing the scalar (S) value in the Edit Window, an user-friendly window exists which provides specification options. With the Select Grid Cell tool active, make a selection, right click to bring up the tool menu and choose the Specify Hard Bottom... option. This will open the CMS-Flow Hard Bottom Specification window.

The following options are provided in the Hard Bottom Specification window:

Use bathymetric cell depth - Sets the cell hard bottom depth to be the cell geometry value thereby creating an exposed non-erodible condition. If multiple cells were selected, then each cell will use its respective bathymetric depth.
Specified distance below bathymetric cell depth - Sets the cell hard bottom depth to be the cell geometry value plus the specified distance thereby creating a sediment-covered non-erodible condition. The distance is limited to positive values to ensure the hard bottom depth is greater than the geometry value. The cell can provide sediment for transportation, however, the amount of erosion is limited. If multiple cells were selected, then each cell will use its respective bathymetric depth.
Specified depth - Sets the cell hard bottom depth to the specified depth thereby creating a sediment-covered non-erodible condition similar to specified distance. The depth is limited to greater than the geometry value. If multiple cells were selected, then the depth is limited to greater than the largest geometry value and all cells will have the same value.
Unspecified - Resets to an inactive hard bottom condition. The cell hard bottom depth is set to the CMS-Flow null value. If multiple cells were selected, then all cells will be reset.

If no cells are selected when opening the Hard Bottom Specification win-dow, then all computational (ocean) cells will be used. If a selection of only non-computational cells, then specification cannot occur. If a selection contains computational and non-computational cells, then the specification will only apply to the computational cells.

If multiple computational cells with differing specifications are selected, the window will not display a selected specification type and the OK button will be disabled. This is to protect the previous specifications from being overwritten by mistake. The OK button will be enabled when an option is selected. The minimum hard bottom depth of the multiple computational cells selected will be displayed in the Depth edit field and the minimum hard bottom depth minus the maximum geometry depth of the multiple computational cells selected will be displayed in the Distance edit field.

Display Options

The hard bottom dataset (when it’s the active dataset) will only display the cells with hard bottom specified if the Ocean cell display option is turned on. Inactive hard bottom cells are not displayed.

CMS-Flow includes hard bottom symbols to differentiate specifications. On the Cartesian Grid page of the Display Options window (when CMS-Flow is the active model), the Hard bottom symbols check box controls the display of symbols that will appear in hard bottom cells (even if the hard bottom dataset is not active). If this is turned on, then the user must be aware of the individual symbol settings accessed by clicking on the Options... button. The Options... button displays the CMS-Flow Hard Bottom Symbols window.

Hard bottom symbols can be selected for three hard bottom specification types:

Non-erodible - Displayed in exposed hard bottom cells (cell hard bottom depth is equal to cell bathymetric depth).
Erodible to specified depth - Displayed in sediment-covered hard bottom cells (cell hard bottom depth is greater than cell bath-ymetric depth).
Invalid specification - Displayed in hard bottom cells where the hard bottom depth is less than cell bathymetric depth (the geometry is below the erosion limit).

If the Hard bottom symbols check box is turned off, no symbols will be displayed and the individual settings cannot be accessed, however, the individual settings will not be changed.

Bed Slope Term

The bed slope term accounts for the effect of gravity on sloped beds. The larger the bed slope coefficient, the more sediment tends to move downslope, thus smoothing the solution. The CMS-Flow used to specify the slope coefficient is described in the table below. The bed slope coeffi-cient is set in the Sediment tab of the CMS-Flow Model Control window in SMS 11.0.

table 2-84

example 2-81

=Note=:

• A constant bed slope coefficient is used spatially, temporally and for all sediment size classes.
• Increasing the bed slope coefficient will increase the smoothing of the bed.
• The measured morphology change can be used to calibrate this pa-rameter. However, its effect is less significant than the capacity transport formula and scaling factors.

Total-Load Correction Factor

The total-load correction factor accounts for the nonuniform vertical pro-file of sediment concentration and current velocity and produces temporal lag between the flow and sediment transport. The factor is used in the nonequilibrium total load sediment transport formula.

table 2-85

example 2-82

=Note=:

• The total load correction factor is only used for the total-load non-equilibrium transport model (NET).
• Computing the total-load correction factor based on the sediment concentration and current velocity profiles increases slightly the computational time but is almost negligible.

=Recommendations=:

• It is NOT recommended to use the total-load correction factor as a calibration factor.  

Bed Composition

Single-sized Sediment Transport

In the case of single-sized sediment transport the bed is usually assumed to be uniform and variations in the bed composition are ignored. In CMS, the option is provided to use a unique approach which the bed composi-tion is allowed to be different from the single transported grain size and the bed composition is considered in hiding and exposure corrections to the sediment transport. The approach works well for cases where the bed is uniform and well sorted with local patches of coarse material such as shell hash. In this approach the bed composition is specified as a median grain size and optionally the 90th percentile diameter. Since only one sed-iment size class is transported, the bed composition cannot change and the input bed composition is assumed to be constant over time. Therefore, the approach does not consider armoring or downstream fining. In addition, since only the surface bed composition is used, the there is no need to specify the bed layering. The advantage of the single-sized sediment transport mode is that it is much less computationally and data intensive compared to the multiple-sized sediment transport model.

table 2-86

example 2-83

example 2-84

example 2-85

example 2-86

=Note=:

• In the case of the single-sized sediment transport model, the bed composition (specified using the d50 and d90 datasets) is optional. The input bed composition is only used to correct the sediment transport for hiding and exposure and does NOT change during simulation.

Multiple-sized Sediments

In the case of where multiple sediment size classes are specified, the bed material is discretized vertically into a number of layers and horizontally using the 2DH computational grid. The bed composition refers to the grain size distribution in each layer and cell.

The basic bed layer settings are specified in the Bed Composition section of the Sediment tab in the CMS-Flow Model Control window (see figure below). The basic parameters for the bed layers are the maximum number of bed layers, and the minimum and maximum bed layer thicknesses.

Fig 2-88.png

Figure 2-88. Bed Composition section within the Sediment tab of the CMS-Flow Model Control window.

The CMS-Flow cards used to specify the maximum number of bed layers, and the minimum and maximum bed layer thickness are described in the table below.

table 2-87

Mixing Layer Thickness

The mixing layer is the first layer from the surface which is allowed to ex-change (or mix) sediments with the sediment transport. The mixing layer thickness is calculated based on the median grain size and bed form size. The option is also available to set the mixing layer to a constant. The CMS-Flow cards used for setting the mixing layer thickness to a constant is described in the table below.


table 2-88

=Notes=:

• Using a larger mixing layer will slow down the temporal change in bed material composition and enhance model stability, but will re-duce the accuracy of the bed composition estimates.
• Erosional cases tend to have smaller mixing layer thicknesses and depositional cases tend to have larger mixing layer thicknesses.
• The mixing layer thickness should NOT be used as a calibration parameter.

Bed Layer Thickness and Composition

The initial bed layer thickness and composition are specified at least one layer for the whole grid. The Layer ID corresponds to the layer number from the surface downwards. If the Layer ID is not specified in the card file, then it is assumed to be equal to the bed layer block sequential num-ber. In order to avoid having to specify all of the bed layers, the bed layer ID’s are allowed to have gaps in which case the missing bed layers or cop-ied from the first specified bed layer above. For example if bed layers 1 and 3 are specified and there is a maximum of 10 layers, then bed layer 2 is copied from bed layer 1, and bed layers 4 through 10 are copied from bed layer 3.

The bed composition is specified using percentile diameters which indicate the percentage of diameters smaller than a specific diameter. The percentile diameters are denoted by ,dper where the subscript per indicates the percentage of diameters smaller than d.

=Notes=:

• There are several advantages to using percentile diameter datasets instead of the fraction of each sediment size class. The percentile diameters are independent of the sediment size classes and there-fore do not need to be changed if the size classes are changed. In addition the percentile diameters are easier to manipulate spatially than sediment fractions because the sum of the later must equal to one.

Fig 2-89.png

Figure 2-89. Bed layer block section within the Sediment tab of the CMS-Flow Model Control window in SMS 11.1.

A bed layer thickness and composition is defined with a bed layer block. The cards used to define the beginning and end of a bed layer block are described in the table below.

table 2-89

Figure 2-89. Bed layer block section within the Sediment tab of the CMS-Flow Model Control window in SMS 11.1.

There are two methods for specifying each bed layer thickness:

1. Constant value
2. Spatially variable (dataset)

The CMS-Flow cards related to the bed layer thickness are described in the table below.

table 2-90

There are four methods for specifying the bed layer composition:

1. D16, D50, and D84 Datasets
2. D35, D50, D90 Datasets
3. Cumulative Grain Size Distribution

The first two options assume a log-normal grain size distribution

  (2-33)


where is the fraction, is the size class bin width, is the size class diameter, and is the median grain size. The geometric standard deviation can be defined using the method of moments

  (2-34)


It is noted that the geometric standard deviation is dimensionless.

The geometric mean dg is in given by

  (2-35)


An alternate approach to the above equation for obtaining geometric standard deviation is by using the graphical measures (cumulative percentile values) (Folk and Ward 1957)

  (2-36)

where is a graphical measure of the geometric standard deviation.

Sediments with a geometric standard deviation close to one are classified as well sorted and poorly graded (see table below).

table 2-91

=Recommendations=:

• Areas with sediments with greater then approximate 1.41 should be simulated with multiple grain sizes.


D16, D50, and D84 Datasets

In this option, 2DH spatial datasets for the diameters corresponding to the 16th, 50th, and 84th percentiles are specified for the whole computational domain. The bed composition is calculated using the log-normal distribution (Equation 2.28) with and

  (2-37)

where the subscript indicates the percentile of the grain size diameter.

D30, D50, and D90 Datasets

In this option, 2DH spatial datasets for the diameters corresponding to the 30th, 50th, and 90th percentiles are specified for the whole computational domain. The bed composition is calculated using the log-normal distribution (Equation 2.28)with and

  (2-38)

where the subscript indicates the percentile of the grain size diameter.

Cumulative Grain Size Distribution

This option should be used in cases where the initial bed grain size distri-bution does not follow a log-normal distribution such as the case of a bi-modal distribution. In this option more than 3 percentile diameter da-tasets need to be specified with at least one below and above the 5oth percentile. Currently the percentile diameters which are allowed in SMS 11.1 are: 5, 10, 16, 20, 30, 50, 64, 84, 90, and 95. The cumulative grain size distribution is linearly interpolated from the percentile diameters at the size class bounds and then converted to sediment fractions.

The CMS-Flow cards used to specify the bed composition are described in the table below

table 2-92

=Notes=:

• When creating a percentile diameter dataset in SMS 11.1, the block number is appended to the end of the dataset name as for example “D50 2”. The number is used to distinguish datasets from different bed layer blocks and is not necessarily the same as the bed layer ID number.
• If the percentile diameters are specified outside of a bed layer block, then they are assumed to belong to the first and only layer.

Advanced Bed Layer Composition

Two additional options are available for specifying the bed composition:


1. D50 dataset and constant geometric standard deviation
2. Size class fractions

A description of the CMS-flow card used to select each bed material com-position option is described in the table below.

table 2-93

=Notes=:

• The above two options are not available in SMS and are therefore considered advanced.
• The geometric standard deviation must always be equal or larger than 1.0.
• When the user specifies the sediment fractions, the sum of the fractions must equal 1.0. If they do not equal 1.0, the fractions are divided by the sum so that they equal 1.0.

Examples

In the case of single-sized sediment transport, the option is available to specify a median grain size dataset which is used to calculate a hiding and exposure correction to the single transport grain size.

example 2-87

example 2-88

Hiding and Exposure

The hiding and exposure is an important parameter for sediment transport. Currently, there are no well established coastal sediment transport formulas for nonuniform sediments. The most common ap-proach is to adapt single-size transport formulas by applying correction factors. The most common of these factors is the hiding and exposure correction. The implementation of this factor varies depending on the transport formula and usually involves a coefficient (exponent) which needs to be calibrated. Because all of the transport formulas used in CMS were originally developed for uniform sediments, there is still a large un-certainty in the value of the hiding and exposure coefficient. It is expected that the hiding and exposure coefficient will be different for different transport formula especially since the implementation of the hiding and exposure correction may be different for different formula.

Single-sized Sediment Transport

For single-sized sediment transport, the correction function for hiding and exposure is calculated following Parker et al. (1982) as

  (2-39)

where m is the hiding and exposure coefficient. The approach has been successfully applied to Shinnecock Inlet, NY to simulate morphology change at a coastal inlet (Sánchez and Wu 2011a). If the sediment is not well sorted, then a multiple sediment fractions need to be simulated (see section below)

Multiple-sized Sediment Transport

The hiding and exposure of the each sediment size class is considered by modifying the critical shields parameter for each sediment size class based on Wu et al. (2000)

  (2-40)

where m is an empirical coefficient that varies for each transport formula, approximately equal to 0.6-1.0. Pek and Phk are the total hiding and exposure probabilities and are calculated as

  (2-41)

where N is the number of grain size classes.

The hiding and exposure coefficient used for both the single- and multiple-sized sediment transport is specified with the same CMS-Flow card (see table below).

table 2-94

=Recommendations=:

• Due to the uncertainty in the hiding and exposure correction coeffi-cient. It is recommended to always calibrate the hiding and exposure coefficient. If no or limited data is available for calibration then a sensitivity analysis should be done


=Note=:

• The hiding an exposure coefficient may have a different “best” value for different transport formula.

Avalanching

Avalanching is the process of sediment sliding when the critical angle of repose is reached. In CMS, avalanching is simulated using a mass con-servative relaxation method which limits the bed slope to the critical angle of repose. For most coastal applications, the critical angle of repose is never reached, so it is not needed. The CMS-Flow cards used for specifying avalanching, and its options, are described in table below.

Fig 2-90.png

Figure 2-90. Avalanching section within the Sediment tab of the CMS-Flow Boundary Conditions window.

The table below describes the CMS-Flow cards used to specify the ava-lanching options.

table 2-95

example 2-89

=Notes=:

• Turning on avalanching increases slightly the computational time.
• For most coastal morphodynamic applications, the avalanching can be turned off since the angle of repose is almost never reached. If a simulation is run without avalanching and steep angles are found in the calculated morphology change, then the simulation should be rerun with avalanching turned on.