RPM2: Difference between revisions
No edit summary (change visibility) |
No edit summary (change visibility) |
||
| (4 intermediate revisions by one other user not shown) | |||
| Line 1: | Line 1: | ||
<div class="title">Ripple Predictor Mark 2</div> | <div class="title">Ripple Predictor Mark 2</div> | ||
Time-dependent ripple geometry is implemented in subroutine RPM2 using the approach described in [[Bibliography#SoulsbyR2005a | Soulsby and Whitehouse (2005)]]. | Time-dependent ripple geometry is implemented in subroutine RPM2 using the approach described in [[Bibliography#SoulsbyR2005a | Soulsby and Whitehouse (2005)]]. Presently, this routine is in the version of SSW_BBL.h located in [https://www.myroms.org/projects/cstm/browser/branches/crs/cstm_trunk Sherwood's branch]. | ||
Ripple geometry is described by ripple height <math>\eta</math> (m), wavelength <math>\lambda</math> (m), and orientation <math>\theta</math> (positive clockwise from the <math>y</math>-axis). The orientation is defined as follows. Under unidirectional flow (no waves), ripple orientation is the same as the current direction. Under pure waves, ripple orientation is in the direction of wave propagation. For linear ripples (assumed here), ripple orientation is perpendicular to the ripple crest. In a slight modification of the original algorithm, the ripple orientation is actually stored as <math>x</math>- and <math>y</math>- components <math>\theta_x</math> and <math>\theta_y</math>. | |||
\eta | |||
The [[Bibliography#SoulsbyR2005a | Soulsby and Whitehouse (2005)]] approach is to calculate equilibrium ripple geometry under existing conditions, and then nudge the existing ripple geometry toward equilibrium at a rate that depends on sediment-transport conditions. Ripples remain unchanged if the transport threshold is not exceeded. The formula for updating ripple height is as follows. | |||
</ | <math display="block">\eta _{t} = \eta _{t-1} + (\eta _{eq} - \eta _{t-1} )(1 - \exp ( -\frac{\beta}{T_e} \Delta t) )</math> | ||
where <math>\eta_t</math> is the ripple height to be calculated for the present time, <math>\eta_{t-1}</math> is ripple height at the previous step, | |||
<math>\eta_{eq}</math> is the equilibrium ripple height under present conditions, <math>\beta</math> is the evolution rate of ripple height under present conditions, <math>T_{e}</math> is the time scale for ripple evolution, and <math>\Delta t</math> is the time step. Similar equations are used to nudge <math>\lambda</math> and components of ripple direction <math>\theta_x</math> and <math>\theta_y</math> toward the equilibrium associated with present sediment-transport conditions. | |||
The Shields parameters for waves alone and for currents alone are calculated, and the equilibrium geometry and evolution rate is determined by the larger of the two. | |||
One potential modification to this routine might be to calculate wash-out criteria and evolution rates based on the combined wave-current stress. | |||
Latest revision as of 11:52, 18 May 2016
Time-dependent ripple geometry is implemented in subroutine RPM2 using the approach described in Soulsby and Whitehouse (2005). Presently, this routine is in the version of SSW_BBL.h located in Sherwood's branch.
Ripple geometry is described by ripple height (m), wavelength (m), and orientation (positive clockwise from the -axis). The orientation is defined as follows. Under unidirectional flow (no waves), ripple orientation is the same as the current direction. Under pure waves, ripple orientation is in the direction of wave propagation. For linear ripples (assumed here), ripple orientation is perpendicular to the ripple crest. In a slight modification of the original algorithm, the ripple orientation is actually stored as - and - components and .
The Soulsby and Whitehouse (2005) approach is to calculate equilibrium ripple geometry under existing conditions, and then nudge the existing ripple geometry toward equilibrium at a rate that depends on sediment-transport conditions. Ripples remain unchanged if the transport threshold is not exceeded. The formula for updating ripple height is as follows.
where is the ripple height to be calculated for the present time, is ripple height at the previous step, is the equilibrium ripple height under present conditions, is the evolution rate of ripple height under present conditions, is the time scale for ripple evolution, and is the time step. Similar equations are used to nudge and components of ripple direction and toward the equilibrium associated with present sediment-transport conditions.
The Shields parameters for waves alone and for currents alone are calculated, and the equilibrium geometry and evolution rate is determined by the larger of the two.
One potential modification to this routine might be to calculate wash-out criteria and evolution rates based on the combined wave-current stress.
