Personal tools

Module Waves

From MohidWiki

Revision as of 16:20, 29 March 2011 by Francisco (talk | contribs) (Equations)
Jump to: navigation, search

Overview

Module Waves is responsible to simulate surface waves in MOHID. It can act as a wave database, reading wave data from files, or it can predict wave parameters through a surface wave generation model which is responsible for computing surface wave properties and its related parameters, therefore communicating directly with the water-air interface module, responsible for all the fluxes between the water column and the atmosphere. The water-air interface module provides wind velocity, direction and shear stresses, receiving from the surface wave module the necessary information to compute water surface rugosity, which can ultimately be used by the turbulence and hydrodynamic module. As surface waves depend and influence water elevations and currents, the wave module also communicates with the hydrodynamic module, providing it with radiation stresses and receiving water depths. Also wave parameters as orbital velocity and excursion are used by the sediment-water interface module to compute bottom shear stresses, which will be used both in hydrodynamics and in sediment transport, controlling erosion and deposition processes.

Fetch-based wave generation model

The surface wave generation model predicts significant wave height and wave period at the end of generation zone based on fetch, water depth and wind speed and direction. The model is restricted to surface wave generation by local wind and propagation is not considered explicitly. As a consequence, this model can only be applied to areas where wave propagation from the exterior is limited (e.g. ocean swell). Examples of such areas are coastal lagoons, lakes or estuaries with limited wave interaction with the ocean.

Algorithm

In areas where the fetch is not restricted (e.g. open sea), there is a dominant fetch and small changes in wind direction don’t significantly change fetch length or the direction of the generated wave. On the other hand, for wind directions parallel to the shore and mainly for lakes, lagoons or estuaries, a small variation in the wind direction can generate significant changes in the wind’s length of action. In these areas there is not a predominant direction for fetch. In order to calculate the effective fetch it is then necessary to take into account the local morphology of the study area considering several directions, each one with a different weight. Thus, for the same wind direction, a point with more obstacles in the near directions has assigned a smaller fetch than a point with more water length in the same directions. These distances are then used to calculate the effective fetch on each cell, in only 8 cardinal directions. Thus, this calculus is only made once, having each cell 8 values of effective fetch. The result is stored in memory during the entire simulation. During the simulation period, in each time step and in each and every cell, depending on wind direction, the effective fetch is selected from the list of 8 values, and it’s used together with wind speed and water depth to calculate the wave parameters. There are some limitations computing effective fetch when simulating domains with large open boundaries, because there is no information to evaluate fetch distances in the boundaries when wind is blowing from the exterior. Effective fetch is computed considering the highest water level in the simulation domain. This means that when if there are areas that during the simulation become uncovered (e.g. forming a temporary tidal flat island), effective fetch is not adjusted. In those cells, all calculus is switched off until they are flooded again.

Equations

As propagation is not considered, the model does not require any initial conditions and there’s no time dependency in the wave parameters evolution. Thus, significant wave height and mean wave period can be computed directly using eq. 1 and eq. 2, being only required to provide the necessary input variables: wind speed, water depth and fetch. The equations were adapted from the CE-QUAL-W2 model (Cole, 2003) and were previously developed by Seymour (1977) and Kang (1982).


<mathtex>Hs=\frac{W^2}g0.283tanh\left[0.53\left(\frac{gH}{W^2}}\right)^{0.75}\right]tanh\left[\frac{0.0125\left(\frac{gF}{W^2}\right)^{0.42}}{tanh\left[0.53\left(\frac{gH}{W^2}\right)^{0.75}\right]}\right]</mathtex>

       Wave height equation 

<mathtex>Ts=\frac{2piW}g1.2tanh\left[0.833\left(\frac{gH}{W^2}}\right)^{0.375}\right]tanh\left[\frac{0.077\left(\frac{gF}{W^2}\right)^{0.25}}{tanh\left[0.833\left(\frac{gH}{W^2}\right)^{0.375}\right]}\right]</mathtex> Wave period equation

Where, Hs is Significant Wave Height, Ts is Significant Wave Period, W is the wind velocity modulus (m/s), H the water depth (m) and F, the fetch distance (m).

These equations predict wave height and wave period at the end of the generation zone, therefore the duration of the wind action it is not accounted for, because it is admitted that the total development of the waves is achieved. The hyperbolic function in eq. 1 and eq. 2 simulates the effects that, for high depths, waves are not influenced by the bottom and is also used to avoid the indefinite increasing of wave height and wave period with increasing fetch. A multiplicative factor for both wave height and period are included in the equations (unitary default value) enabling the calibration of the model results against available waves data.

The broader equation (eq. 3 and Figure 4) for effective fetch (Fe) was obtained from Rogala (1997) and Howes (1997):

eq. 3


Figure 5. Fetch calculation with 3 and 5 directions.

Where Fi is the fetch distance (m) in the i direction and αi = angle (º) between the wind direction and the i direction. For the fetch directions, 3 or 5 distances (eq. 4 and eq. 5 respectively) can be used (Figure 5). This is called the effective modified fetch (Fm): eq. 4


eq. 5


Where, F45L e F22L are the distances to land in the directions 22.5º and 45º left of the wind direction (m); F45R e F22L are the distances to land in the directions 22.5º and 45º right of the wind direction (m); and F0 is distance to land along wind direction (m). The weight given to each distance is based on the cosine function, taking a unitary weight for the wind direction and decreasing weights to angles away from this direction. In order to compute the effective fetch, distances to land need to be computed. As described above the distances can be computed in two different ways: i) in a grid-based method with an algorithm that moves in the grid along the 8 or 16 directions until it finds a land cell or the end of the grid; ii) in a “graphical” method that computes distances between cell centre and land polygons.

Wave orbital velocity and excursion

Linear wave theory is generally applied to determine the near-bed velocities. In case of symmetrical (sinusoidal) small-amplitude waves in relatively deep waters this theory yields good results. When waves approach shallower waters, the waves will be distorted leading to asymmetrical wave profiles and higher order wave theories are necessary to determine the near-bed velocities. In this model only linear wave theory is considered, therefore applying it, the peak value of the orbital excursion (Aδ) and velocity (Uδ) at the edge of the wave boundary layer can be expressed as:

User manual

Main options

In this chapter it is briefly described how to use Module Waves. The different options provided by the module can be defined through an input data file, similarly to other modules in MOHID. In order to couple this module with the rest of the model simulation, it has to be activated in the Module Model input data file, by defining the following keyword:

WAVES                                     : 1

ModuleWaves can be used as a database, reading the solution for the different parameters (radiation stresses, significant wave height, mean wave period and mean wave direction) from a file, or it can be used to compute wave parameters based on external information. Depending on the options chosen, ModuleWaves might be dependent on the activation of other modules, namely the Atmosphere and InterfaceWaterAir modules, which provide wind information. In order to use the wave module so that it influences hydrodynamics (water level and currents), the user must activate the radiation stress options and activate in the Hydrodynamic input data file the following keyword:

WAVESTRESS                           : 1

To use Module Waves so that the waves effect contributes to the calculation of bottom shear stress it is necessary to turn on the following keyword in InterfaceSedimentWater data file:

WAVETENSION                       : 1

This will also enable to consider waves influence on resuspension and deposition of particulate matter in the water-sediment interface.

Input data file

Outputs

Time series

Maps (HDF5 format)

References

  • Braunschweig, F., P. Chambel, L. Fernandes, P. Pina and R. Neves, 2004. The object-oriented design of the integrated modelling system MOHID. Computational Methods in Water Resources International Conference, Chapel Hill, North Carolina, USA.
  • Cole, T.M. and S.A. Wells, 2003. CE-QUAL-W2: A two-dimensional, laterally averaged, Hydrodynamic and Water Quality Model, Version 3.2. Instruction Report EL-03-1, US Army Engineering and Research Development Center, Vicksburg, MS.
  • Howes, 1997. British Columbia Estuary mapping Systems – Appendix A: Wave Exposure Calculation. Ministry of Sustainable Resource Management, British Columbia U.S.A.
  • Rogala, J.T., 1997. Estimating Fetch for Navigation Pools in the Upper Mississippi River Using a Geographic Information System. United States Geological Survey -Project Status Report 97-08 .
  • Seymour, R.J., 1977. Estimating Wave Generation in Restricted Fetches. J. ASME WW2, May 1977 pp251-263.