TR0813:Chapter2
Model Description
Waveaction balance equation with diffraction
Taking into account the effect of an ambient horizontal current or wave behavior, CMSWave is based on the waveaction balance equation as (Mase 2001)

(1) 
where

(2) 
is the waveaction density to be solved and is a function of frequency and direction . E() is spectral wave density representing the wave energy per unit watersurface area per frequency interval. In the presence of an ambient current, the waveaction density is conserved, whereas the spectral wave density is not (Bretherton and Garrett 1968; Whitham 1974). Both wave diffraction and energy dissipation are included in the governing equation. Implementation of the numerical scheme is described elsewhere in the literature (Mase 2001; Mase et al. 2005a). C and C_{g} are wave celerity and group velocity, respectively; x and y are the horizontal coordinates; C_{x}, C_{y}, and C_{$\theta$} are the characteristic velocity with respect to x, y, and, respectively; N_{y} and N_{yy} denote the first and second derivatives of N with respect to y, respectively; is an empirical parameter representing the intensity of diffraction effect; _{b} is the parameterization of wave breaking energy dissipation; S denotes additional source S_{in} and sink S_{ds} (e.g., wind forcing, bottom friction loss, etc.) and nonlinear wavewave interaction term.
Wave diffraction
The first term on the right side of Equation 1 is the wave diffraction term formulated from a parabolic approximation wave theory (Mase 2001). In applications, the diffraction intensity parameter (values 0) needs to be calibrated and optimized for structures. The model omits the diffraction effect for = 0 and calculates diffraction for > 0. Large (values > 15) should be avoided as it can cause artificial wave energy losses (Mase 2001). In practice, values of between 0 (no diffraction) and 4 (strong diffraction) have been determined in comparison to measurements. A default value of = 2.5 was used by Mase et al. (2001, 2005a, 2005b) to simulate wave diffraction for both narrow and wide gaps between breakwaters. In CMSWave, the default value of assigned by SMS is 4, corresponding to strong diffraction. For wave diffraction at a semiinfinite long breakwater or at a narrow gap, with the opening equal or less than one wavelength, = 4 (maximum diffraction allowed in the model) is recommended. For a relatively wider gap, with an opening greater than one wavelength, = 3 is recommended. The exact value of in an application is dependent on the structure geometry and adjacent bathymetry, and should to be verified with measurements.
Wavecurrent interaction
The characteristic velocties C_{x}, C_{y}, and C_{$\theta$} in Equation 1 can be expressed as:
(3) 
(4) 
(5) 
where U and V are the depthaveraged horizontal current velocity components along the x and y axes, k is the wave number, and is the water depth. The dispersion relationships between the relative angular frequency , the absolute angular frequency , the wave number vector , and the current velocity vector are (Jonsson 1990)
(6) 
and
(7) 
where
is the Dopplershifting term, and g is the acceleration due to gravity. The main difference between the wave transformation models with and without ambient currents lies in the solution of the intrinsic frequency. In treatment of the dispersion relation with the Doppler shift, there is no solution corresponding to wave blocking, if intrinsic group velocity C_{g} is weaker than an opposing current (Smith et al. 1998; Larson and Kraus 2002):
(8) 
Under the wave blocking condition, waves cannot propagate into a strong opposing current. The wave energy is most likely to dissipate through breaking with a small portion of energy either reflected or transformed to lower frequency components in the wave blocking condition. In CMSWave, the wave action corresponding to the wave blocking is set to zero for the corresponding frequency and direction bin.
Wave reflection
The wave energy reflected at a beach (or upon the surface of a structure) is calculated under assumptions that the incident and reflected wave angles, relative to the shorenormal direction, are equal in magnitude and that the reflected energy is a given fraction of the incident wave energy. The reflected wave action N_{r} is assumed to be linearly proportional to the incident wave action N_{i }:

(9) 
where K_{r} is a reflection coefficient (0 for no reflection and 1 for full reflection) defined as the ratio of reflected to incident wave height (Dean and Dalrymple 1984).
CMSWave calculates the wave energy reflection toward the shore, e.g., reflection from a sidewall or jetty, within the wave transformation routine. It can also calculate reflection toward the sea boundary, e.g., wave reflection off the beach or detached breakwater, with a backward marching calculation routine (Mase et al. 2005a). Users should be aware that although the computer execution time calculating forward reflection is relatively small, the time is almost double for the backward reflection routine.
Wave breaking formulas
The simulation of depthlimited wave breaking is essential in nearshore wave models. A simple wave breaking criterion that is commonly used as a first approximation in shallow water, especially in the surf zone, is a linear function of the ratio of wave height to depth. For random waves, the criterion is (Smith et al. 1999)

(10) 
where H_{b} denotes the significant breaking wave height. A more comprehensive criterion is based on the limiting steepness by Miche (1951) for random waves as

(11) 
where k_{p} is the wave number corresponding to the spectral peak. In the shallow water condition (k_{p}h small), Equation 11 reduces asymptotically to Equation 10. Iwagaki et al. (1980) verified that Miche's breaker criterion could replicate laboratory measurements over a sloping beach with a current present, provided that the wavelength was calculated with the current included in the dispersion equation.
In CMSWave, the depthlimited spectral energy dissipation can be selected from four different formulas: (a) Extended Goda formulation (Sakai et al. 1989), (b) Extended Miche (Battjes 1972; Mase et al. 2005b), (c) Battjes and Janssen (1978), and (d) Chawla and Kirby (2002). These formulas, considered more accurate for wave breaking on a current, can be divided into two generic categories (Zheng et al. 2008). The first class of formulations attempts to simulate the energy dissipation due to wave breaking by truncating the tail of the Rayleigh distribution of wave height on the basis of some breaker criterion. The Extended Goda and Extended Miche formulas belong to this class. The second category of wave breaking formulas uses a bore model analogy (Battjes and Janssen 1978) to estimate the total energy dissipation. The Battjes and Janssen formula and Chawla and Kirby formula are in this class. The spectral energy dissipation is calculated based on one of these four wave breaking formulas, and the computed wave height is limited by both Equations 10 and 11.
Extended Goda formula
Goda (1970) developed a breaker criterion, based on laboratory data, taking into account effects of the bottom slope and wave steepness in deep water. This criterion is used widely in Japan. Goda's formula was modified later by Sakai et al. (1989) to include the action of an opposing current in which a coefficient accounted for the combined effects of current, depth, bottom slope angle , and deepwater wavelength L_{o}
(12) 
where A = 0.17 is a proportional constant, and
(13) 
with
(14) 
in which L_{o}g, and T_{p} is the spectral peak period. The change in breaker height with respect to the cell length dx is defined as
(15) 
In CMSWave, the breaking heights at the seaward and landward sides, denoted as H_{bi} and H_{bo}, respectively, of a grid cell are given by
(16) 
(17) 
and the rate of wave breaking energy is (Mase et al. 2005b)
(18) 
where H_{1/3 } is the significant wave height defined as the average of the highest onethird waves in a wave spectrum.
Extended Miche formula
Battjes (1972) extended Miche's criterion, Equation 11, to water of variable depth as
(19) 
where L_{b} is the wavelength at the breaking location including the current, is an adjustable coefficient varying with the beach slope, and a = 0.14. This formula reduces to a steepness limit in deep water and to a depth limit in shallow water, thus incorporating both wave breaking limits in a simple form. The coefficient has been treated as a constant value of 0.8 in application for random waves (Battjes and Janssen 1978). Based on field and laboratory data, Ostendorf and Madsen (1979) suggested that
(20) 
which is applied in CMSWave. The change of breaker height with respect to the cell length dx can be obtained as follows (Mase et al. 2005b):
(21) 
Equations 19 to 21 are incorporated in Equations 16 to 18 to calculate the spectral energy dissipation rate _{b}.
Battjes and Janssen formula
Battjes and Janssen (1978) developed a formula for predicting the mean energy dissipation in a bore of the same height as a depthinduced breaking wave as
(22) 
where is an empirical coefficient of order one, is the sea water density, is the spectral mean frequency, and Q_{b} is the probability that at a given location the wave is breaking. By assuming the wave height has a Rayleigh distribution, the probability of wave breaking can be determined from the following expression
(23) 
where is the rootmeansquare wave height. Battjes and Janssen (1978) calculated the maximum possible height from Equation 19 using a constant breaker value = 0.8. Booij et al. (1999) and Chen et al. (2005) investigated Battjes and Janssen formula and obtained a better wave breaking estimate with = 0.73. In CMSWave, Equation 23 is adapted to parameterize the wave breaking energy dissipation by applying the Battjes and Janssen formula with = 0.73. The calculation of the wave breaking dissipation rate is from:
(24) 
Chawla and Kirby formula
Chawla and Kirby (2002) proposed an alternative expression for the bulk dissipation in random waves assuming the probability of wave breaking is dependent on the wave slope and a bore type of dissipation. Their modified bore dissipation formula worked well for wave breaking under a strong opposing current. The rate of energy dissipation was defined as
(25) 
where is the wave number corresponding to the spectral mean frequency , and scaling parameters b and are equal to 0.4 and 0.6, respectively. The rate of wave breaking energy dissipation is calculated by Equation 24.
Wind forcing and whitecapping dissipation
The evolution of waves in the largescale, open coast is more affected by windoceanwave interactions than on the nearshore wavecurrentbottom processes. The result is a nonlinear wave field that is balanced between wind forcing, whitecapping, and wave growth. The surface wind can feed energy into the existing waves and can also generate new waves. On the other hand, the energy can dissipate through whitecapping from turbulencewave interactions and airwavewater interactions. In CMSWave, these wind forcing and whitecapping processes are modeled as separate sink and source terms (Lin and Lin 2004a and b).
Wind input function
The windinput source S_{in} is formulated as functions of the ratio of wave celerity C to wind speed W, the ratio of wave group velocity to wind speed, the difference of wind speed and wave celerity, and the difference between wind direction _{wind} and wave direction (Lin and Lin 2006b):
(26) 
where
(27) 
(28) 
(29) 
and
(30) 
is the functional form of the PiersonMoskowitz (PM) spectrum, _{o }= g/W is the Phillips constant, and
(31) 
is a normalized directional spreading. The function F_{1} presents the wind stress effect, F_{2} designates Phillips' mechanisms (Phillips 1957) and F_{3} accounts for the wave age effect. For swell or long waves, the wave group velocity C_{g} is generally large and F_{3} < 1. If C_{g} W, then F_{3} = 0. For short waves, the phase velocity is generally small and F_{3} > 1.
Whitecapping dissipation function
The wave energy dissipation (sink) S_{ds} (Lin and Lin 2006b) for whitecapping including current and turbulent viscous effect is
(32) 
with
(33) 
and
(34) 
where c_{ds} is a proportionality coefficient, and is for the turbulent viscous dissipation. The wave amplitude is calculated at each grid cell. To avoid numerical instability and considering the physical constraint of energy loss for the dissipation, the function F_{4} is set to 1 if the computed value is greater than 1.
Wave generation with arbitrary wind direction
In the case of wind forcing only, with zero wave energy input at the sea boundary, CMSWave can assimilate the fullplane wave generation. The model will execute an internal grid rotation, based on the given wind direction, to calculate the wave field and map the result back to the original grid. This feature is convenient for the local wave generation by wind in a lake, bay, or estuary, neglecting swell from the ocean.
Bottom friction loss
The bottom friction loss (sink) S_{ds} is calculated by a drag law model (Collins 1972)
(35) 
with
(36) 
where _{b} presents the ensemble mean of horizontal wave orbital velocity at the sea bed, E_{total} is the total energy density at a grid cell, and c_{f }is the DarcyWeisbach type friction coefficient. The relationship between c_{f} and the DarcyWeisbach friction factor f_{DW} is c_{f} = f_{DW}/8. Typical values of c_{f} for sandy bottoms range from 0.004 to 0.007 based on the JONSWAP experiment and North Sea measurements (Hasselmann et al. 1973; Bouws and Komen 1983). Values of c_{f} applied for coral reefs range from 0.05 to 0.40 (Hardy 1993; Hearn 1999; Lowe et al. 2005). Application of this model capability to a specific site requires validation to field data. If the Manning friction coefficient n is used instead of the DarcyWeisbach type coefficient, the relationship between the two drag coefficients is
(37) 
Estimates of Manning coefficient n are available in most fluid mechanics reference books (e.g., 0.01 to 0.05 for smooth to rocky/weedy channels).
Wave runup
Wave runup is the maximum shoreward wave swash on the beach face for engineering structures such as jetties and breakwaters by wave breaking at the shore. Wave runup is significant for beach erosion as well as wave overtopping of seawalls and jetties. The total wave runup consists of two components: (a) rise of the mean water level by wave breaking at the shore, known as the wave setup, and (b) swash of incident waves. In CMSWave, the wave setup is computed based on the horizontal momentum equations, neglecting current, surface wind drag and bottom stresses
(38) 
(39) 
where is the water density and S_{xx}, S_{xy,} and S_{yy} are radiation components from the excess momentum flux caused by waves. By using the linear wave theory (Dean and Dalrymple 1984), S_{xx}, S_{xy,} and S_{yy} can be expressed as
(40) 
(41) 
(42) 
where . Equations 38 and 39 also calculate the water level depression from the stillwater level resulting from waves known as wave setdown outside the breaker zone. Because CMSWave is a halfplane model, Equation 38 controls mainly wave setup and setdown calculations, whereas Equation 39 acts predominantly to smooth the water level alongshore. The swash oscillation of incident natural waves on the beach face is a random process. The most landward swash excursion corresponds to the maximum wave runup. In the engineering application, a 2% exceedance of all vertical levels, denoted as R2, from the swash is usually estimated for the wave runup (Komar 1998). This quantity is approximately equal to the local wave setup on the beach or at structures such as seawalls and jetties, or the total wave runup is estimated as
(43) 
In CMSWave, R2 is calculated at the landwater interface and averaged with the local depth to determine if the water can flood the proceeding dry cell. If the wave runup level is higher than the adjacent land cell elevation, CMSWave can flood the dry cells and simulate wave overtopping and overwash at them. The feature is useful in coupling CMSWave to CMSFlow (Buttolph et al. 2006) for calculating beach erosion or breaching. Calculated quantities of S_{xx}/x, S_{xy}/x, S_{xy}/y, and S_{yy}/y are saved as input to CMSFlow. CMSWave reports the calculated fields of wave setup and maximum water level defined as
Maximum water level = Max (R2, + H1/3/2)  (44) 
Wave transmission and overtopping at structures
CMSWave applies a simple analytical formula to compute the wave transmission coefficient K_{t} of a rigidly moored rectangular breakwater of width B_{c} and draft D_{c} (Macagno 1953)
(45) 
Wave transmission over a structure or breakwater is caused mainly by the fall of the overtopping water mass. Therefore, the ratio of the structure crest elevation to the incident wave height is the prime parameter governing the wave transmission. CMS calculates the rate of overtopping of a vertical breakwater based on the simple expression (Goda 1985) as
(46) 
where h_{c} is the crest elevation of the breakwater above the stillwater level, and H_{i} is the incident wave height. Equation 46 is modified for a composite breakwater, protected by a mound of armor units at its front, as
(47) 
For rubblemound breakwaters, the calculation of wave transmission is more complicated because the overtopping rate also depends on the specific design of the breakwater (e. g., toe apron protection, front slope, armor unit shape and size, thickness of armor layers). In practice, Equation 47 still can be applied using a finer spatial resolution with the proper bathymetry and adequate bottom friction coefficients to represent the breakwater.
Grid nesting
Grid nesting is applied by saving wave spectra at selected locations from a coarse grid (parent grid) and inputting them along the offshore boundary of the smaller fine grid (child grid). For simple and quick applications, a singlelocation spectrum saved from the parent grid can be used as the wave forcing for the entire sea boundary of the child grid. If multilocation spectra were saved from the parent grid, they are then interpolated as well as extrapolated for more realistic wave forcing along the sea boundary of the child grid. Multiple grid nesting (e.g., several coexisting child grids and grandchild grids) is supported by CMSWave. The parent and child grids can have different orientations, but need to reside in the same horizontal coordinate system. Because CMSWave is a halfplane model, the difference between grid orientations between parent and child grids should be small (no greater than 45 deg) for passing sufficient wave energy from the parent to child grids.
Variablerectangularcell grid
CMSWave can run on a grid with variable rectangular cells. This feature is suited to largedomain applications in which wider spacing cells can be specified in the offshore, where wave property variation is small and away from the area of interest, to save computational time. A limit on the shorenormal to shoreparallel spacing ratio in a cell is not required as long as the calculated shoreward waves are found to be numerically stable.
Nonlinear wavewave interaction
Nonlinear wavewave interactions are a conserved energy transfer from higher to lower frequencies. They can produce transverse waves and energy diffusion in the frequency and direction domains. The effect is more pronounced in shallower water. Directional spreading of the wave spectrum tends to increase as the wavelength decreases. The exact computation of the nonlinear energy transfer involves sixdimensional integrations. This is computationally too taxing to be used in practical engineering nearshore wave transformation models. Mase et al. (2005a) have shown that calculated wave fields differ with and without nonlinear energy transfer. Jenkins and Phillips (2001) proposed a simple formula as an approximation of the nonlinear wavewave interaction. Testing of this formula in CMSWave is underway.
Fastmode calculation
CMSWave can run in a fast mode for simple and quick applications. The fast mode calculates the halfplane spectral transformation on either five directional bins (each 30deg angle for a broadband input spectrum) or seven directional bins (each 5deg angle for a narrowband input spectrum or 25deg angle with wind input) to minimize simulation time. It runs at least five times faster than the normal mode, which operates with 35 directional bins. The fastmode option is suited for a long or timepressing simulation if users are seeking preliminary solutions. The wave direction estimated in the fast mode is expected to be less accurate than the standard mode because the directional calculation is based on fewer bins.
Table of Contents  Chapter 3  CMSWave Interface 