CMS-Flow Model Parameters: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
Line 181: Line 181:


=== Adaptation ===
=== Adaptation ===
[[File:CMSFlow Model Parameters GeneralTab Flow SedimentTransport AdaptationV13.2.12.png|500px|thumb|center| Figure 20. Adaptation options]]
[[File:CMSFlow Model Parameters GeneralTab Flow SedimentTransport AdaptationV13.2.12.png|500px|thumb|center| Figure 21. Adaptation options]]


The adaptation coefficient (Figure 20) is an important parameter to consider in setting up the CMS sediment transport model. The sensitivity of results to the adaptation coefficient depends on the spatial and temporal scales of the problem. For example, if a high-resolution grid is used to model short-term dynamics of a nearshore disposal site or a small inlet, then the adaptation coefficient is likely to be important. However, if a relatively coarse grid is used to study sediment pathways in a large estuary entrance, then the adaptation coefficient will not be important, [[https://cirpwiki.info/wiki/User_Guide_011#Adaptation_Coefficient 7]]
The adaptation coefficient (Figure 21) is an important parameter to consider in setting up the CMS sediment transport model. The sensitivity of results to the adaptation coefficient depends on the spatial and temporal scales of the problem. For example, if a high-resolution grid is used to model short-term dynamics of a nearshore disposal site or a small inlet, then the adaptation coefficient is likely to be important. However, if a relatively coarse grid is used to study sediment pathways in a large estuary entrance, then the adaptation coefficient will not be important, [[https://cirpwiki.info/wiki/User_Guide_011#Adaptation_Coefficient 7]]


=== Transport grain sizes classes ===
=== Transport grain sizes classes ===

Revision as of 16:45, 19 December 2022

Figure 1. CMS-Flow Model Control : Model Parameters tab in SMS 13.2.12.

All of the CMS-Flow model parameters, settings, and output options are controlled from the CMS-Flow Model Control window (Figure 1). The window has several tabs including the Flow tab in which most of the general settings are set for CMS-Flow. To open the CMS-Flow Model Control window, click on CMS-Flow | Model Control. The first tab will be the Flow tab. There are several sections within the tab including Time Control and Hot Start options.

General

Figure 2: Time Control parameters for the simulation

Time Control

Below are some of the CMS cards related to the Time Control of the Flow tab. Table 1 provides a brief description of the CMS cards used for time control. Figure 2 shows the time parameters in CMS-Flow.

Card Arguments Default Range Description
STARTING_JDATE REAL none none Julian data in YYDDD with YY being last two digits of the year, and DDD the Julian day of the year.
STARTING_JDATE_HOUR REAL none none Julian hour .
HYDRO_TIME_STEP REAL Calculated based on solution scheme and courant number none Sets to the time step for hydrodynamics in seconds.
DURATION_RUN REAL 48.0 >RAMP>0 Sets the duration of the model simulation in hours.
DURATION_RAMP REAL 1.0 Sets the length of the ramp period in which the model forcing is slowly ramped from zero.

Hot Start

Figure 3: Hot start parameters on CMS-Flow Model Control

The term Hot start (Figure 3) refers to starting a simulation with an initial condition other zero (cold start). Hot starts are used for specifying initial conditions or restarting simulations at intermediate times. The hot start controls are set in the Flow tab of the CMS-Flow Model Control window.

Hot Start File

Figure 4. HDFView showing the structure of the CMS Hot Start File.

The CMS hot start feature CMS lets the user restart simulations that have been stopped due to electric outages, hardware malfunctions, or model crashes. In the case of a model crash the user, may restart the model using larger solver iterations and/or time steps to stabilize the simulation. The user has the option to specify a hot start output time or an interval for outputting a recurring hot start file. Every time the hot start file is written, it overwrites the previous information. The CMS Hot Start file saves information on the water elevation (pressure), and current velocities. If the sediment transport is active, then the water depth and sediment concentrations are also saved for each size class. Only the very last record of information is preserved (no starting from earlier intervals).

The CMS hot start files are written as binary XMDF files by default. Depending on the type of hot start (single file or recurring), the names are as follows are saved in the directory of the CMS-Flow files:

  • SingleHotStart.h5
  • AutoHotStart.h5

After saving a CMS Hot Start file, it is a good idea to rename the file with a different name before using it as an initial conditions file. This way, the file will not be overwritten in future simulations.

Table 1. Hot Start CMS-Flow Cards

Card Arguments Default Range Description
HOT_START_TIME REAL none none Single time after start at which to output a single hot start file.
AUTO_HOT_START_INTERVAL REAL none none Sets the recurring hot start output interval .
  • Initial conditions file

There are several situations where it is convenient to specify a user defined initial condition (hot start) file. For example, if the user forgets to setup the model output a hot start file or when running idealized cases with known initial conditions. A hot start file can easily be created and exported by the user from the SMS interface. The model requires at water levels, current velocities, concentrations, and water depths. Any datasets that are missing from the initial file are assigned a default value which depends on the dataset. If the water level or current velocities are not specified, they are set to zero. If the depth is not specified, then it is set to the input grid depth. If the sediment concentrations are not specified, then they are set to the equilibrium concentrations. It is important to note that the names and paths of the initial condition datasets are important.

Variable Path and name
Water surface elevation Datasets\Water_Elevation
current velocity Datasets\Current_Velocity
Sediment concentrations Datasets\Concentration
Salinity concentrations Datasets\Salinity


Steps for creating a user defined hot start or initial condition file

1. Import CMS-Flow grid and solution file.
2. Sample a time step of the solution datasets for use in the initial condition
  • Click on Data | Data Calculator (Figure 5)
Figure 5 Dataset Toolbox showing a time step sample of the water elevation and current velocity datasets for use in a hot start (initial condition) file.
  • Under the Tools section, select Sample time steps.
  • Under the Datasets section, click on the variable desired to be save (water surface elevation, current velocity, sediment concentrations or salinity concentrations)
3. Export the initial condition datasets to an XMDF file (Figure 6)
Figure 6.a Dataset Toolbox showing a time step sample of the water elevation and current velocity datasets for use in a hot start (initial condition) file.

+

Figure 6.b Saving conditions files
  • Write Hot Start output file [On/Off]
  • Automatic recurring Hot Start file [On/Off]

Solution Scheme

Figure 7 Solution Scheme options

This refers to the temporal discretization of the hydrodynamic, sediment and salinity transport equations. There are two options in CMS: implicit and explicit.

  • Implicit - First order backward Euler scheme. Uses a time step on the order of 5-15 minutes. Appropriate for cases which can be simulated with large computational time steps such as long term morphology change at inlets.
  • Explicit - First order forward Euler scheme. Uses a time step on the order of 0.5-1.0 second. Appropriate for cases that vary quickly in time such as flooding or barrier island breaching.
Matrix Solver
Figure 8. Hot start parameters on CMS-Flow Model Control

Matrix solver (Figure 8):

Threads

Figure 9. CMS Flow General Tab Threads Version 13.2.12
  • Number of threads (Figure 7)

The CMS-Flow is parallelized for PC’s with multi-core processors using OpenMP. The number of threads used for parallel processing. The parallelization works by splitting the computational work into “threads” among several cores. Some cores are hyper-threaded, meaning a single core may support two threads. The number of threads is specified in the CMS-Flow Model Control Window. The number of threads must be equal or greater to 1 and cannot be larger than the number of threads available on the machine. If a number is specified which is larger than the maximum number available on the machine, then the code will default to the maximum number available.

For more detail, see Multi Processor Capability Reference: CMS User Guide 013

Flow

In the Flow tab inside the CMS Flow, users can specify parameters and variables that will be used in the simulation. For example, in this section the Coriolis effect is assigned, friction coefficients and turbulence models are chosen.

Parameters

Figure 10. CMS Flow Parameters Version 13.2.12
  • Hydrodynamic Time Step
  • Wetting and drying*

In the numerical simulation (Figure 10) of the surface water flows with sloped beaches, sand bars and islands, the water edges change with time, with part of the nodes being possibly wet or dry. In the present model, a threshold flow depth (a small value such as 0.02 m in field cases) is used to judge drying and wetting. If the flow depth on a node is larger than the threshold value, this node is considered to be wet, and if the flow depth is lower than the threshold value, this node is dry. Because a fully implicit solver is used in the present model, all the wet and dry nodes participate in the solution. Dry nodes are assigned a zero velocity. On the water edges between the dry and wet nodes, the wall-function approach is applied. [1]

  • Wave Fluxes and Roller Fluxes

User can check the options for wave and roller fluxes if wave outputs are desired.

  • Average latitude for Coriolis

Specifies the average latitude for the grid which is used for the Coriolis parameter. If specified and average latitude is calculated for the grid and used to estimate the Coriolis parameter. By default, the cell latitudes are stored in the XMDF Model Parameters File. For additional details on this file see XMDF CMS-Flow Model Parameters File (*_mp.h5). For more details, see [2] Reference: Geospatial Information More information and tests: [3][4]

Bottom and wall friction

In this section, users may specify wave current bottom friction and roughness dataset (Figure 11).

Figure 11. Bottom and Wall friction options

Bottom roughness dataset

The bottom roughness type (Mannings N, Bottom Friction Coefficient, or Roughness Height) may be specified using a grid dataset referred to as the Bottom Friction Dataset. The dataset contains spatially variable or constant values and is the most general way of specifying the bed roughness. A description of the CMS-Flow cards related to the Bottom Friction Dataset are shown in the table below.

The bottom roughness is specified in the Bottom Friction Dataset section of within the Input tab of the CMS-Flow Model Control window. The Bottom roughness is specified at every computational (ocean) cell and is required for each model simulation. In SMS 11.1, the bottom friction dataset is specified in the Flow tab of the CMS-Flow Model Control window (see Figure 2 23). As in SMS 11.0, the bottom roughness can be specified as a roughness height, bottom friction coefficient, or in Manning's n coefficient. [5]

For the wave-current bottom frictions coefficients, five options are available (Figure 12):

Figure 12. List of wave friction coefficients available
  • Quadratic
  • Soulsby (1995) Data 2
  • Soulsby (1995)Data 13
  • Fredsoe (1984)
  • Huyn-Thanh and Temperville (1991)

For more details, see Bottom Friction See How to create a bottom friction dataset

Turbulence parameters

Figure 13. CMS Flow Turbulance models

Users can choose different Turbulence Model: Subgrid, Falconer, Parabolic or Mixing length (Figure 13). IT is not recommended to edit the default turbulence parameters. For more details, see Subgrid Model [6]

Sediment Transport

Timing

Users can determine from which CMS-Flow will start calculating the morphology change (Figure 14). The default value is 0 hours

Figure 14. Morphology change timing

Formulation

Figure 15. Non-Equilibrium total load

The formulation option available in CMS Flow is the Nonequilibrium total load (Figure 15. For more information, see:

Transport formula

There are five transport formulas to calculate the morphology change (Figure 16).

Figure 16. Transport formulas available
  • Lund-CIRP
  • Van Rijn
  • Soulsby-van Rijn
  • Watanabe
  • C2SHORE

For a detail explanation, see:

Properties

Figure 17. Sediment properties (density and porosity)

Default values for the sediment density and porosity (Figure 17) are given. For more information, see:

Bed Composition

Figure 18. Bed composition options

In these sections (Figure 18), users can choose to run the model with the simplified grain size definitions.

Scaling factors and coefficients

Figure 19. Default values scaling factors and coefficients

Adaptation

Figure 21. Adaptation options

The adaptation coefficient (Figure 21) is an important parameter to consider in setting up the CMS sediment transport model. The sensitivity of results to the adaptation coefficient depends on the spatial and temporal scales of the problem. For example, if a high-resolution grid is used to model short-term dynamics of a nearshore disposal site or a small inlet, then the adaptation coefficient is likely to be important. However, if a relatively coarse grid is used to study sediment pathways in a large estuary entrance, then the adaptation coefficient will not be important, [7]

Transport grain sizes classes

Avalanching

Figure 19. Avalanching parameters

Avalanching is the process of sediment sliding when the critical angle of repose is reached (Figure 19). In CMS, avalanching is simulated using a mass conservative 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.[1] Avalanching is invoked if the bottom slope between two cellsexceeds a specified critical slope. The avalanching algorithm applies an iterative approach to move material down slope until the critical slope is no longerexceeded, while conserving sediment volume. [2]

Hard Bottom

Figure 20. Hardbottom dataset option

The sediment transport and bed change equations assume a loose bottom in which the bed material is available for entrainment. However, hard bottoms may be encountered in practical engineering applications where bed materials are non-erodible, such as bare rocks, carbonate reefs, and concrete coastal structures. Hard-bottom cells in CMS are handled by modifying the equilibrium concentration as in both the sedment transport and bed change equations. The bed-slope term in the bed change equation is also modified so that only deposition (no erosion) may occur at hard-bottom cells. [7]

Although the hard bottom dataset can be edited (when its the active dataset) 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. [8]

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 window, 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.

Salinity/Temperature

The general parameters are the water density and temperature. The CMS-Flow cards for the general parameters are described in the table below.

Water properties

Figure 21. Water properties options
  • Wave information
  • Extrapolation distances

Salinity

Figure 22.a Salinity available inputs: Global or Spatially Varied
Figure 22.b Global salinity option if Calculate Salinity is activated

The characteristics of salinity are important in the coastal environment because salinity can impact marine plants and animals and influence the dynamic behavior of cohesive sediments. Because modifications of coastal inlets, such as channel deepening and widening and rehabilitation or extension of coastal structures, may alter the salinity distribution within estuaries or bays, it is often useful and convenient to simulate the salinity within the scope of an engineering project to determine if a more detailed water quality modeling study is necessary. It is important to emphasize that the CMS is not intended to be used as a water quality model. The CMS solves the depth-averaged (2DH) salinity transport equation and should be used only for cases where the water column is well mixed. If there is flow stratification, a 3D model should be utilized. It is also noted that the salinity is not used to update the water density which is assumed to be constant. Thus any horizontal water density gradients due to varying salinity on the hydrodynamics are assumed to be negligible.2

Temperature

Figure 23. Temperature options

The properties of sea water are a function of temperature, salinity (i.e. dissolved minerals) and pressure (Figure 23). The water density and dynamic viscosity (kinematic viscosity times the water density) are shown in the following two figures as a function of water temperature and salinity. 10

Table 1. CMS-Flow cards related to the general parameters

Card Arguments Default Description
WATER_DENSITY REAL 1025 Water density in kg/m^3.
WATER_TEMPERATURE REAL 15 Water temperature in degrees Celcius.

Wave

Wave information

Figure 24.a Single wave condition
Figure 24.b Inline steering

In this module, users may assign wave values or insert the simulation files from CMS Wave (.sim).

Extrapolation distances

Wind

Users can insert Wind forcing into the model with three different methods:

  • Spatially constant
Figure 25. Wind spatially constant input
  • Meteorological stations
Figure 25.b Wind from meteorological stations
  • Temporally and spatially varying from file
Figure 25.c Wind from different file types

For more details, see CMS Flow Wind Forcing

Output

Figure 26. Output tabs options

Users may specify different time intervals for the multiple outputs of CMS Flow. Different lists with different timeframes can be created and be assigned to specific outputs. For more information, see: CMS Flow Output




Documentation Portal