CMS-Wave Model Parameters: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
(Created page with ==CMS-Wave Setup== ===Grid Generation=== Now that all modifications to the CMS-Flow grid alternatives have been finalized, the wave grid can be generated from the Flow grid. CM...)
 
(72 intermediate revisions by 3 users not shown)
Line 1: Line 1:
[[File:CMSWave Model Parameters GeneralTab V13.2.12.png|thumb|right|400px|CMS-Wave Model Control window Version 13.2.12.]]
= Parameters =
[[File:CMS Wave Model Control V. 13.2.12.png|500px|thumb|center]]
==CMS Wave plane mode ==


==CMS-Wave Setup==
===Full half plane with reverse spectra===
===Grid Generation===
===Full plane===
Now that all modifications to the CMS-Flow grid alternatives have been finalized, the wave grid can be generated from the Flow gridCMS-Wave grids are typically generated from scatter datasets, but can be a mirror image of the flow grid if the boundary condition (BC) locations are close enough to force both the tide and waves from the same BC area.  In the near future, CMS will have the option to be run fully inline, where only one grid will contain all of the information needed to run both CMS-Wave and CMS-Flow.
In this mode, CMS-Wave performs two half-plane runs in the same grid.  The first run is in the half-plane with the principle wave direction toward the shoreThe second run is in the seaward half-plane as opposed to the first run.  Upon the completion of the second run, two half-plane results are combined to one full-plane solution (Lin et al. 2012).  Because the run time for the full-plane is approximately twice of the regular half-plane, users shall consider the full-plane mode only if the full-plane features like wave generation and propagation in a bay or around an island.
===Half Plane===


The CMS-Wave IJ triad should be adjusted where the I-direction is pointing in the dominant wave direction, or more specifically from the open ocean. This means that it will be 180 deg off as compared to the CMS-Flow grid. Therefore, after the CMS-Flow grid is generated, the model domain should be set to include all major reaches of the bay, extending laterally a couple kilometers from the inlet and extending seaward well beyond the anticipated reach of the ebb jet. A good rule of thumb for setting the ocean boundary for inlet modeling is to place the boundary three or four times the length of the ebb jet (leaving provision for extreme discharges, such as during storm surge). The given bathymetry ends at approximately four times the length of the ebb jet at Shark River.
== Source terms ==
CMS-Wave is a phase-averaged model for propagation of directional irre-gular waves over complicated bathymetry and nearshore where wave re-fraction, diffraction, reflection, shoaling, and breaking simultaneously act at inlets. Wave diffraction terms are included in the governing equations following the method of Mase et al. (2005). Four different depth-limiting wave breaking formulas can be selected as options including the interaction with a current. The wave-current interaction is calculated based on the dispersion relationship including wave blocking by an opposing current (Larson and Kraus 2002). Wave generation and whitecapping dissipation are based on the parameterization source term and calibration using field data (Lin and Lin 2004a and b, 2006a). Bottom friction loss is estimated from the classical drag law formula (Collins 1972).


br style=clear:both /
== [https://cirpwiki.info/wiki/CMS-Flow:Wave-current_Interaction Current interaction] ==
[[Image:Shark_Fig26.PNG|thumb|right|400px|Figure 26.''''''''Data'' ''''''Options, ''''''''Switch Current Model'''''''': Available Models.]]


* Transforming the CMS-Flow grid to a CMS-Wave Grid
== Bottom friction ==
Bottom friction could be assigned as constant or by a dataset inside SMS. To see more details, see the [https://cirpwiki.info/wiki/CMS-Flow_Model_Parameters#Bottom_and_wall_friction Bottom and wall friction of the CMS Flow].
Usually, the same roughness method used on CMS Wave is used on the CMS Flow parameters bottom friction.
For a detail explanation of how Bottom Friction variable is used on the CMS Wave model, see [https://cirpwiki.info/wiki/Bottom_Friction Bottom Friction]


1. Open one of the generated grids rt-click the Grid in the SMS list, and click on ''Duplicate,'' select the copied grid (copy) to activate it, and go to ''Data'', ''Switch Current Model'' (Figure 28),
==Surge Fields==
==Wind Fields==
If the spatial wind field input is required, users shall prepare a wind.dat file or *.wind (in the same format as *.cur) to provide the x- and y-component wind data corresponding to the incident wave conditions in the model grid.


2. Change the Model from ''CMS-Flow'' to ''CMS-Wave'', click ''OK'', and rename the grid by rt-clicking on the copy name and clicking ''Rename'',
==Matrix Solver==
https://cirpwiki.info/wiki/CMSFlow_Matrix_Solver


br style=clear:both /
= Boundary control =
[[Image:Shark_Fig27.PNG|thumb|right|400px|Figure 27. Setting structure cells in CMS-Wave.]]
[[File:CMS Wave Model Control Boundray ControlV. 13.2.12.png|500px|thumb|center]]
Water level and wind information are optional source as specified under
Wave Source in addition to the spectral input data.
==Source==
Spatially varied spectral input – This is simply the case as in a child grid that spatially varied wave spectra are permitted to assign at user specified locations along or near the seaward boundary of the child grid. To apply spatially varied spectra for wave input without a parent grid, users will need to prepare the wave input file with the format as described in the child grid run.


3. Set the jetty and groin cells to rubble-mound type structures by selecting the cells (Figure 27), and rt-click, ''Cell Attributes..'', and select a ''Rubble-mound'' from the ''Structure Type'', click ''OK''.
https://cirpwiki.info/wiki/CMS-Wave_Input_Spectra


4. ''File'', ''Save As'', ''Wave.sim'' (selecting the ''Save As Type'' as a .sim for simulation) in the folder with the CMS-Flow grid.
==Interpolation==
*Inverse distance weighting


br style=clear:both /
The inverse-distance interpolation also referred to as Shepard interpola-tion is given by (Shepard 1968)


{{Equation|<math>\phi(\overrightarrow{x}) = \sum_{i=1}^N w_i \phi_i</math>|17-1}}


===Providing Ocean Buoy Data to CMS-Wave ===
where the interpolation weights are given by


After all modifications to the CMS-Wave grid are complete, the spectral waves or wave
{{Equation|<math>w_i = \frac{d_i ^{-\rho_s}}{ \sum_{j} d_j^{-\rho_s}}</math>|17-2}}
parameters can be generated for the wave grid. Full (directional) spectra can be imported into the SMS for the CMS-Wave, as well as simplified wave parameters (angle, wave height, and period, etc). Directional spectral data collected by the National Data Buoy Center (NDBC) or Coastal Data Information Program (CDIP) buoys, available from the
National Oceanographic Data Center, can be processed (transformed to the model domain)
and used as a source for wave input to CMS-Wave.


br style=clear:both /
where 
[[Image:Shark_Fig28.PNG|thumb|right|400px|Figure 28. NODC buoy data access website.]]


==== Importing NDBC bouy data ====
<math>\rho_s</math>  = real and positive power parameter [-]
*Run ndbc-spectra.exe to read the NDBC standard directional wave file and prepare the CMS-Wave input spectral *.eng file.


1. Download the NDBC standard monthly directional wave spectral file from
d = distance between the known points <math>\overrightarrow{x_i}</math>  and the unknown interpolation points <math>\overrightarrow{x}</math>  equal to the Euclidean norm <math>d = ||\overrightarrow{x} - \overrightarrow{x}_i||</math> .  
http://www.nodc.noaa.gov/BUOY/buoy.html (e.g., 44025_200912) (Figures 28 to 31), this is provided under Hands On Materials - Shark River Inlet/Waves.


[[Image:Shark_Fig29.PNG|thumb|right|400px|Figure  29. NODC buoy data access worldmap.]]
In this interpolation, the weight of each point decreases with distance from the interpolated point. One advantage of the inverse-distance interpolation is the interpolation weights are independent of the interpolation function, and therefore only need to be calculated once and can be saved for computational efficiency.[[User_Guide_028|[1]]]
br style=clear:both /
[[Image:Shark_Fig30.PNG|thumb|right|400px|Figure  30. NODC buoy data access regional map.]]
br style=clear:both /
[[Image:Shark_Fig31.PNG|thumb|right|400px|Figure  31. NDBC buoy spectral data download website.]]
br style=clear:both /


2. In the DOS mode (Command Prompt), run ndbc-spectra.exe (Figure 32),
== Computational spectral grid==
Spectral waves or wave parameters can be generated for the wave grid forcing, or wind direction and speeds can provide the necessary information for wind- wave generation. Full (directional) spectra can be imported into the SMS for the CMS-Wave, as well as simplified wave parameters (angle, wave height, and period, etc).


3. Responding to the on-screen input, type the NDBC spectral filename,
==Sides==


4. Type the starting time stamp (default value is 0) for saving output file(s),
==Case data==


5. Type ending time stamp (default is 99999999) for saving output file(s),
===Wind direction angle conversion===
*Cartesian
*Meteorologic
*Oceanographic
*Shore normal


6. Type the time interval (hr) for saving output data,
===Populate from Spectra===
===Set Reference Time===


7. Type 2 to save the CMS-Wave *.eng and *.txt files,
= Output control =
[[File:CMS Wave Model Control OutputControl. 13.2.12.png|500px|thumb|center]]
== Limit observation output ==
== Radiation stresses ==
== Sea/swell ==
==Breaking type==


8. Type the local shoreline orientation (actually the CMS-Wave grid y axis)
= Options =
in Polar Coordinates (deg, positive from North covering the sea, e.g., 180
[[File:CMS Wave Model Control Options 13.2.12.png|500px|thumb|center]]
deg for a north-south oriented project on the Atlantic coast),
To include (trigger) either of wave run-up, infra-gravity wave, nonlinear wave-wave interaction, binary (xmdf or *.h5) output, multiple processors, muddy bed, and spatial wind field input is just a one-click step in the SMS interface. Additional files are required for the muddy bed and spatial wind field input.
==[https://cirpwiki.info/wiki/CMS-Flow_NUmerical_Methods:_Wetting_and_Drying Allow wetting and drying]==


9. Type the NDBC buoy location water depth (m) and then the CMS-Wave
==Infragravity wave effect==
seaward boundary mean water depth (m), e.g. Buoy 44025 has a nominal
==[https://cirpwiki.info/wiki/CMS-Wave:Diffraction Diffraction intensity]==
depth of 36.3 m,


[[Image:Shark_Fig32.PNG|thumb|right|400px|Figure 32. Run NDBC-spectra in DOS window.]]
==Nonlinear wave effect==
==Run up==
==Fast-mode run==
==Roller effects==
The wave roller parametric formulation is commonly applied in the wave spectral model to modify breaking wave energy dissipation nearshore to mimic better the surf zone dynamicsThe wave roller effect used in CMS-Wave is based on the roller model developed by Zhang et al. (2014).


10. Type 1 to include wind or 0 to skip the wind, and to complete the run
==Forward reflection==
(Figure 32),
==Backward reflection==
==Muddy bed==
If the muddy bed calculation is required, users shall prepare a mud.dat file or *.mud (in the same format as *.dep) to list the spatial varying max-imum kinematic viscosity for the entire grid (recommended maximum kinematic viscosity for mud is 0.04 m2/sec).
==Wave breaking formula==
==Date format==


The fort.15 and fort.18 files are the parameter and spectra, respectively, to be imported in
To setup the model parameters for CMS-Wave:
to CMS-Wave. Two additional files are generated: fort.12 and fort.19, which are the
tabulated wave (and wind) parameters for the full-plane and half-plane wave propagation
coverage, respectively. Matlab fig1set.m and fig2sets.m can read these files to plot the
wave and wind data time series.


br style=clear:both /
1. Go to CMS-Wave, Model Control, Options and turn on Allow wetting and drying
and Bed friction


11. Rename fort.15 to *.txt and fort.18 to *.eng (half-plane wave energy),
2. Users can also specify constant or varied forward and backward reflection
which was renamed to wave_shark_2009.eng,
coefficients in Settings,
 
[[Image:Shark_Fig33.PNG|thumb|right|400px|Figure  33. Importing wave spectra.]]


12. In SMS, go to CMS-Wave, then Spectral Energy, and select Import, and
select the *.eng file just created (Figure 33),


br style=clear:both /
<br style="clear:both" />
[[Image:Shark_Fig34.PNG|thumb|right|400px|Figure    34. Generated wave spectra parameters with a snapshot of spectral output.]]
To read more about CMS Wave defintions, see [https://cirpwiki.info/wiki/CMS-Wave_Model_Control CMS Wave Model Control Definitions]


13. Save as a new CMS-Wave file (.sim) and all necessary files including the
= CMS-Wave Model Control File formats=
spectral energy file will be saved with that name (Figure 34).
The CMS-Wave Options File (*.std) can have one of 3 input formats. Click one of the options below for the format description.
 
If these steps (from pages 27-31) are skipped, there is a grid provided under Hands On
Materials - Shark River Inlet/CMS-Wave that has the spectra included.
 
br style=clear:both /
 
==== Importing CDIP buoy data ====
Importing CDIP buoy data is similar to the method of importing NDBC buoy data.
 
1. Download the CDIP wave file from http://cdip.ucsd.edu/?nav=historicsub=data (e.g., sp154-200912) (Figures 35 to 37), also provided under Hands On Materials - Shark River
Inlet/Waves,
 
Note that the only difference is the processed CDIP wave spectra do not need to be
transformed to the CMS-Wave seaward boundary. The data can be directly input to
CMS-Wave. Examples given are CDIP 154 and NDBC 44025 standard spectral files for
Dec 2009.
 
br style=clear:both /
[[Image:Shark_Fig35.PNG|thumb|right|400px|Figure 35. CDIP buoy data access website.]]
 
br style=clear:both /
[[Image:Shark_Fig36.PNG|thumb|right|400px|Figure 36. CDIP buoy data access website.]]
 
br style=clear:both /
[[Image:Shark_Fig37.PNG|thumb|right|400px|Figure  37. CDIP buoy spectral data download website.]]
 
2. Run cdip-spectra.exe and follow the same steps as those from pages 30-31
(from processing the NDBC buoy data with the ndbc-spectra.exe program),
to read the CDIP standard directional wave file and prepare the data for
input to CMS-Wave as a *.eng file.
 
 
br style=clear:both /
 
===Model Control===
br style=clear:both /
[[Image:Shark_Fig38.PNG|thumb|right|400px|Figure 38. CMS-Wave Model Control window.]]
To setup the model parameters for CMS-Wave:
 
1. Go to CMS-Wave, Model Control, and turn on Allow wetting and drying
and Bed friction (Figure 38),
 
2. Users can also specify constant or varied forward and backward reflection
coefficients in Settings,


3. Water level and wind information are optional source as specified under
{|
Wave Source in addition to the spectral input data,
|-
|1. [[Wave_NoCards|No cards - order of values on one line is very important.]] || '''(Used with SMS 11.1 and previous)'''
|-
|2. [[Wave_CardFormat1|Card-based - [value]  [!card name]]] || '''(Used with SMS 11.2 to 13.1)'''
|-
|3. [[Wave_CardFormat2|Card-based - [card name] [value(s)]]] || '''(Used with SMS 13.2 and after)'''
|}


4. File, Save As, Wave.sim (selecting the Save As Type as


br style=clear:both /
[[CMS-Wave_File_Formats | Back to CMS-Wave File Formats]]

Latest revision as of 19:52, 20 December 2022

CMS-Wave Model Control window Version 13.2.12.

Parameters

CMS Wave Model Control V. 13.2.12.png

CMS Wave plane mode

Full half plane with reverse spectra

Full plane

In this mode, CMS-Wave performs two half-plane runs in the same grid. The first run is in the half-plane with the principle wave direction toward the shore. The second run is in the seaward half-plane as opposed to the first run. Upon the completion of the second run, two half-plane results are combined to one full-plane solution (Lin et al. 2012). Because the run time for the full-plane is approximately twice of the regular half-plane, users shall consider the full-plane mode only if the full-plane features like wave generation and propagation in a bay or around an island.

Half Plane

Source terms

CMS-Wave is a phase-averaged model for propagation of directional irre-gular waves over complicated bathymetry and nearshore where wave re-fraction, diffraction, reflection, shoaling, and breaking simultaneously act at inlets. Wave diffraction terms are included in the governing equations following the method of Mase et al. (2005). Four different depth-limiting wave breaking formulas can be selected as options including the interaction with a current. The wave-current interaction is calculated based on the dispersion relationship including wave blocking by an opposing current (Larson and Kraus 2002). Wave generation and whitecapping dissipation are based on the parameterization source term and calibration using field data (Lin and Lin 2004a and b, 2006a). Bottom friction loss is estimated from the classical drag law formula (Collins 1972).

Current interaction

Bottom friction

Bottom friction could be assigned as constant or by a dataset inside SMS. To see more details, see the Bottom and wall friction of the CMS Flow. Usually, the same roughness method used on CMS Wave is used on the CMS Flow parameters bottom friction. For a detail explanation of how Bottom Friction variable is used on the CMS Wave model, see Bottom Friction

Surge Fields

Wind Fields

If the spatial wind field input is required, users shall prepare a wind.dat file or *.wind (in the same format as *.cur) to provide the x- and y-component wind data corresponding to the incident wave conditions in the model grid.

Matrix Solver

https://cirpwiki.info/wiki/CMSFlow_Matrix_Solver

Boundary control

CMS Wave Model Control Boundray ControlV. 13.2.12.png

Water level and wind information are optional source as specified under Wave Source in addition to the spectral input data.

Source

Spatially varied spectral input – This is simply the case as in a child grid that spatially varied wave spectra are permitted to assign at user specified locations along or near the seaward boundary of the child grid. To apply spatially varied spectra for wave input without a parent grid, users will need to prepare the wave input file with the format as described in the child grid run.

https://cirpwiki.info/wiki/CMS-Wave_Input_Spectra

Interpolation

  • Inverse distance weighting

The inverse-distance interpolation also referred to as Shepard interpola-tion is given by (Shepard 1968)

  (17-1)

where the interpolation weights are given by

  (17-2)

where

= real and positive power parameter [-]

d = distance between the known points and the unknown interpolation points equal to the Euclidean norm .

In this interpolation, the weight of each point decreases with distance from the interpolated point. One advantage of the inverse-distance interpolation is the interpolation weights are independent of the interpolation function, and therefore only need to be calculated once and can be saved for computational efficiency.[1]

Computational spectral grid

Spectral waves or wave parameters can be generated for the wave grid forcing, or wind direction and speeds can provide the necessary information for wind- wave generation. Full (directional) spectra can be imported into the SMS for the CMS-Wave, as well as simplified wave parameters (angle, wave height, and period, etc).

Sides

Case data

Wind direction angle conversion

  • Cartesian
  • Meteorologic
  • Oceanographic
  • Shore normal

Populate from Spectra

Set Reference Time

Output control

CMS Wave Model Control OutputControl. 13.2.12.png

Limit observation output

Radiation stresses

Sea/swell

Breaking type

Options

CMS Wave Model Control Options 13.2.12.png

To include (trigger) either of wave run-up, infra-gravity wave, nonlinear wave-wave interaction, binary (xmdf or *.h5) output, multiple processors, muddy bed, and spatial wind field input is just a one-click step in the SMS interface. Additional files are required for the muddy bed and spatial wind field input.

Allow wetting and drying

Infragravity wave effect

Diffraction intensity

Nonlinear wave effect

Run up

Fast-mode run

Roller effects

The wave roller parametric formulation is commonly applied in the wave spectral model to modify breaking wave energy dissipation nearshore to mimic better the surf zone dynamics. The wave roller effect used in CMS-Wave is based on the roller model developed by Zhang et al. (2014).

Forward reflection

Backward reflection

Muddy bed

If the muddy bed calculation is required, users shall prepare a mud.dat file or *.mud (in the same format as *.dep) to list the spatial varying max-imum kinematic viscosity for the entire grid (recommended maximum kinematic viscosity for mud is 0.04 m2/sec).

Wave breaking formula

Date format

To setup the model parameters for CMS-Wave:

1. Go to CMS-Wave, Model Control, Options and turn on Allow wetting and drying and Bed friction

2. Users can also specify constant or varied forward and backward reflection coefficients in Settings,



To read more about CMS Wave defintions, see CMS Wave Model Control Definitions

CMS-Wave Model Control File formats

The CMS-Wave Options File (*.std) can have one of 3 input formats. Click one of the options below for the format description.

1. No cards - order of values on one line is very important. (Used with SMS 11.1 and previous)
2. Card-based - [value] [!card name] (Used with SMS 11.2 to 13.1)
3. Card-based - [card name] [value(s)] (Used with SMS 13.2 and after)


Back to CMS-Wave File Formats