View source for SSW BBL

<div class="title">SSW bottom boundary layer formulation</div>
__TOC__
==Wave-orbital calculations==
Near-bed wave-orbital characteristics, including representative orbital velocity <math>u_{br}</math>, representative period <math>T_r</math>, and average direction of wave propagation <math>\theta_w</math> (degrees, nautical convention, which is positive clockwise from north) are defined according to Madsen (1994). When SWAN results are used, these correspond to UBOT, PWAVE, and DWAVE. If surface-wave statistics (e.g., <math>H_s</math>, <math>T_d</math>, and <math>\theta_w</math>) are provided, they can be converted to bottom orbital velocity externally (using, for example, the routines suggested in [[Bibliography#WibergP_2008a | Wiberg and Sherwood (2008)]] and provided as UBOT in a SWAN input file. Alternatively, if SSW_CALC_UB is defined, orbital velocity <math>u_{br}</math> is calculated according to linear wave theory as follows:

<math display="block">u_{br} = \frac{H_s}{2\sinh (kh)}</math>

where <math>kh</math> is wavenumber x depth, and <math>k</math> can be approximated using one of the methods described [[SoulsbyKH | here]].
==Ripple Geometry==
Ripple height <math>\eta_r</math> and wavelength <math>\lambda_r</math> are calculated using information from the previous time step and the Malarkey and Davies (2003) implementation of the Wiberg and Harris (1994) formulation, which is valid for wave-dominated conditions. They approximate ripple wavelength as <math>535 D_{50}</math> and ripple steepness as:

<math display="block">\frac{{\eta _r }}
{{\lambda _r }} = \exp \left[ { - 0.095\left( {\ln \left( {\frac{{d_0 }}
{{\eta _r }}} \right)} \right)^2  + 0.442\left( {\ln \left( {\frac{{d_0 }}
{{\eta _r }}} \right)} \right) - 2.28} \right]</math>

where <math>d_0 = u_{br}T/\pi</math> is the wave-orbital diameter. When transport stage is below the threshold for sediment transport (<math>T_*  = frac{\tau_{wc}}(\tau_{ce} < 1</math>), ripple dimensions from the previous time step are retained.

An alternative formulation for time-dependent ripple roughness is in development branches and is described [[RPM2 | here]].

==Bottom Roughness==
Roughness lengths associated with grain roughness <math>z_{0N}</math>, sediment transport roughness <math>z_{0ST}</math>, and bedform roughness (ripples) <math>z_{0BF}</math> are estimated as:

<math display="block">z_{0N}  = 2.5D_{50} /30</math>


<math display="block">z_{0ST}  = \alpha D_{50}  \cdot a_1 \frac{{T_* }}
{{1 + a_2 T_* }}</math>


<math display="block">z_{0BF}  = a_r \eta _r^2 /\lambda _r</math>

where the sediment transport coefficients are <math>\alpha = 0.056</math>, <math>a1 = 0.068</math>, and <math>a_2  = 0.0204\ln (100D_{50}^2 ) + 0.0709\ln (100D_{50} )</math> (Wiberg and Rubin, 1989) with the bedform roughness <math>D_50</math> expressed in meters, and where <math>a_r</math> is a coefficient that may range from 0.3 to 3 (Soulsby, 1997). Grant and Madsen (1982) proposed <math>a_r = 27.7/30</math>, but we use a default value of <math>a_r = 0.267</math>, as suggested by Nielsen (1992). The roughness lengths are additive, so subsequent BBL calculations use <math>z_0  = \max \left[ {z_{0N}  + z_{0ST}  + z_{0BF} ,\quad z_{0MIN} } \right]</math>, where <math>z_{0MIN}</math> allows setting a lower limit on bottom drag (default <math>z_{0MIN} = 5e^{-5}</math>m).

==Wave-current combined stress and roughness==
Initial estimates of (kinematic) bottom stresses based on pure currents <math>\tau_c</math> ( = <math>\tau_b</math> ) and pure waves <math>\tau_w</math> (<math>\tau_b = 0</math>) are made as follows.

<math display="block">\tau _c  = \frac{{\left( {u^2  + v^2 } \right)\kappa ^2 }}
{{\ln ^2 (z/z_0 )}}</math>

and <math>\tau _w  = 0.5f_w u_b^2</math>, where <math>f_w</math> is the Madsen (1994) wave-friction factor, which depends on the ratio of the wave-orbital excursion amplitude to the bottom roughness length <math>A_b/k_b</math>, where <math>A_b = u_{br}T /(2\pi)</math> and <math>k_b = 30z_0</math>. <math>f_w =</math>

<math display="block">0.3,  A_b /k_b  \leq 0.2</math>
 
<math display="block">\exp ( - 8.82 + 7.02(A_b /k)^{ - 0.078} ), 0.2 < A_b /k_b  \leq 100</math>

<math display="block">\exp ( - 7.30 + 5.61(A_b /k)^{ - 0.109} ),  A_b /k_b  > 100</math>

The pure-current and pure-wave limits are used as initial estimates for calculations towards consistent profiles for eddy viscosity and velocity between <math>z_0</math> and <math>z_r</math>, using either the model of Madsen (1994) or Styles and Glenn (2000). Both of these models assume eddy viscosity profiles scaled by <math>u_{*wc}  = \sqrt {\tau _{wc} }</math> in the wave boundary layer (WBL) and <math>u_{*c}  = \sqrt {\tau _b }</math> in the current boundary layer, calculated as <math>K_M  =</math>

<math display="block">\kappa u_{*wc} z, z < \delta _{wbl}</math>

<math display="block">\kappa u_{*c} z, z > \delta _{wbl}</math>

where <math>\delta _{wbl}</math> is the thickness of the WBL, which scales as <math>u_{*wc}T/(2\pi)</math>. <math>\tau_{wc}</math> represents the maximum vector sum of wave- and current-induced stress, but the <math>\tau_b</math> is influenced by the elevated eddy viscosity in the WBL, and must be determined through an iterative process. The shape and elevation of the transition between these profiles and other details differ among the two models, but both the models of Madsen (1994) or Styles and Glenn (2000) return values for the horizontal vectors <math>\tau_b</math>, <math>\tau_w</math>, and <math>\tau_{wc}</math>. The parameter <math>\tau_b</math> is the mean (over many wave periods) stress used as the bottom boundary condition in the momentum equations, and <math>\tau_{wc}</math> is the maximum instantaneous stress exerted over the bottom by representative waves and currents.

==Skin friction - form drag partitioning==
When ripples are present, <math>\tau_{wc}</math> is a combination of form drag, which does not directly contribute to sediment transport, and skin friction, which does. The next step in the BBL calculations is to estimate the skin-friction component of <math>\tau_{wc}</math> using the ripple dimensions and a bedform drag-coefficient approach (Smith and McLean, 1977; Wiberg and Nelson, 1992), as follows.

<math display="block">\tau _{sfm}  = \tau _{wc} \left[ {1 + 0.5C_{dBF} \frac{{\eta _r }}
{{\lambda _r \kappa ^2 }}\left( {\ln \frac{{\eta _r }}
{{\left( {z_{0N}  + z_{0ST} } \right)}} - 1} \right)^2 } \right]^{ - 1}</math>

where <math>C_{dBF} \approx 0.5</math> is a bedform drag coefficient for unseparated flow (Smith and McLean, 1977).

==Maximum shear stress==
Because shear stress varies between ripple crests and troughs, an estimate of the maximum shear stress at the crests <math>\tau_{sfm}</math> is calculated for use in sediment-transport algorithms as: 

<math display="block">\tau _{sf}  = \tau _{sfm} (1 + 8\frac{{\eta _r }}
{{\lambda _r }})</math>