CMS-Flow:Features

Boundary Conditions

CMS-Flow has multiple types of boundary conditions which are listed and discussed below. All CMS-Flow boundary conditions are forced at the edges of the domain by use of cellstrings defined with the Surfacewater Modeling System.

Flow Rate-Forcing

User creates a cellstring for the given boundary, chooses a Flow Rate type for the cellstring, then creates a time-series curve of flow rates. NOTES:

• Total flow rate specified is divided between the total number of cells in the cellstring with each carrying a portion of the total.
• The sign of the flow rate curve is dependent on the direction of flow with respect to the origin (always lower-left hand corner of the grid). This guide should assist in proper assignment.
  • Flow rate from the East - Negative value
  • Flow rate from the West - Positive value
  • Flow rate from the North - Negative value
  • Flow rate from the South - Positive value

Water Surface Elevation Forcing

Two types of Water Surface Elevation Forcing exist for CMS-Flow. Once a water surface elevation curve (or series of curves) is applied, the user is able to display the curve information graphically.

1. Single value for all cells on a cellstring
2. Multiple values on a cellstring (one for each cell)

Single Value

User creates a cellstring for the given boundary and defines a time-series curve. The value for each time on this curve is applied to all cells along the designated cellstring.

Multiple Values

User creates a cellstring for the given boundary and extracts multiple time-series curves from a dataset or database. Each cell along the cellstring is given its own time-series curve information. Examples are:

• Extraction of water surface elevation values from a larger domain solution (ie. Larger CMS-Flow or ADCIRC grid)
• Extraction of tidal constituent information from a tidal database, from which a water surface elevation curve can be generated.

Water Surface Elevation and Velocity Forcing

Users are able to extract both water surface elevations and velocity components from a larger domain solution (ie. Larger CMS-Flow or ADCIRC grid).

Salinity Concentration Forcing

If salinity transport is active for the simulation, the user has the ability to use existing hydrodynamic cellstrings in the interface in order to provide a time-series curve of salinity concentrations.

Global Forcing

There are two main types of global forcing available in CMS-Flow, wind and wave. Global forcing means that the forcing is applied on a cell-by-cell nature, rather than forcing along a boundary.

Wind Forcing

Temporally varying, spatially constant (ie. one value of wind applied for a given time to every computational cell in the domain.

• In the near future a variable wind will be allowed and an interface within the SMS will be provided.

Wave Forcing

Temporally and spatially varying.

• Wave forcing is generally provided by the user selecting to use the Steering Module within SMS. The mapping of wave data from wave grid to flow grid is automated during the course of the steering process.
• If a choice is made not to use the steering process, the user must provide several datasets of information which has been mapped to the flow grid geometery.
  • Radiation stress gradient
  • Wave height
  • Wave period
  • Wave direction
  • Wave dissipation

Transport Options

Sediment Transport

CMS-Flow includes three sediment transport models:

1. Equilibrium Total load (Exner Equation)
2. Equilibrium Bed load plus Non-Equilibrium Suspended Load
3. Non-equilibrium Total load

All of the previously mentioned models account for hard bottom and effect of the bed slope.

The nearbed sediment concentation or concentration capacity are calculated with one of the following transport formula:

1. Lund-CIRP (2006)
2. Van Rijn (1998)
3. Watanabe (1987)
4. Soulsby-van Rijn (1997) (V4.0 and higher)

Equilibrium Total Load

In this model, both the bed load and suspended load are assumed to be in equilibrium. The bed change is solved using a simple mass balance equation known as the Exner equation.

More information on the this model can be found here.

Equilibrium Bed load plust Advection-Diffusion Suspended Load

Calculations of suspended load and bed load are conducted separately. The bed load is assumed to be in equilibrium and is included in the bed change equation while the suspended load is solved through the solution of an advection-diffusion equation.

More information on the this model can be found here.

Non-equilibrium Total Load

The non-equilibrium sediment transport algorithm (NET) simulates non-cohesive, single size sediment transport and bed change using a Finite Volume method and includes advection, diffusion, hiding and exposure, and avalanching. NET sediment transport is calculated with a non-equilibrium bed-material (total load) formulation. In this approach, the suspended- and bed-load transport equations are combined into a single equation and thus there is one less empirical parameter to estimate (adaptation length).

Additional information on NET can be found here.

Salinity Transport

In many estuaries, the density gradients caused by spatial variations in salinity can be an important driving force in the circulation. Salinity is also a key water quality variable in estuaries, since it affects the chemical and biological processes. Salinity is simulated in the Coastal Modeling System (CMS) in a depth-averaged sense. This means that the estuary or body of water is assumed to be well mixed vertically and the salinity is constant over the water column.

The salinity transport equation is solved with an explicit, finite volume method. The advection term is discretized with upwind scheme, and the diffusion term is discretized with the standard central difference scheme.

Additional information on Salinity Transport in CMS-Flow can be found here.

Other Processes

Bottom Friction

The bottom friction parameter (related to Manning $n$) is spatially varying (cell-specific) over the grid domain. The default value upon grid creation is 0.025. At times a user may desire to represent locations where added friction is needed due to structures or increased turbulence due to sharp changes in current speed. More information on using this feature of CMS-Flow can be found here.

Hard Bottom

Hard Bottom is a morphologic constraint that provides the capability to simulate mixed bottom types within a single simulation. This cell-specific feature limits the erodibility of the constrained cells down to a specified depth below the water surface. More information on the use of hard bottom within the SMS can be found here.

Variable D50

When sediment transport and morphology change are activated, sometimes it may be desired to designate areas of the grid to have a certain grain size. This cell-specific parameter allows each cell to take on the characteristics of a certain sediment grain size. More information on the use of variable grain size (D50) can be found here.

Eddy Viscosity

Eddy Viscosity in CMS-Flow is not a user definable parameter. It is presently calculated with a subgrid turbulence model. More information can be found here.

Other Features

Parallelization with OpenMP

Both Intel and AMD processors now are shipping chips with multiple cores/processors (henceforth referred to as "processors") available. CMS-Flow is now configured to make use of these extra processes that are available on newer machines.

Additional information on using Multiple Processors with CMS-Flow can be found here.

Units of Measurement

Variable	Units	Symbol
Water Surface Elevation	meters	($m$)
Current Velocity	meters per second	($m/sec$)
Flow Rate	cubic meters per second	($m^{3}/sec$)
Salinity Concentration	parts per thousand	($ppt$)
Sediment Concentration	kilogram per meter cubed	($kg/m^{3}$)
Sediment Transport	meter squared per second	($m^{2}/sec$)

---

CMS-Flow