CMS-Flow Hydrodnamics: Variable Definitions: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
(Created page with "Phillips (1977), Mei (1983), and Svendsen (2006) provide a detailed deri-vation of the depth-integrated and wave-averaged hydrodynamic equa-tions. Here, only variable definiti...")
 
No edit summary
 
(16 intermediate revisions by 2 users not shown)
Line 1: Line 1:
Phillips (1977), Mei (1983), and Svendsen (2006) provide a detailed deri-vation of the depth-integrated and wave-averaged hydrodynamic equa-tions. Here, only variable definitions are provided and derivations may be obtained from the preceding references. The instantaneous current velocity u<sub>i</sub>  is split into:


''<math>u_i = \bar{u_i} + \tilde{u_i} + u_i^'</math>''
 
The instantaneous current velocity u<sub>i</sub>  is split into:
{{Equation|
<math>u_i = \bar{u_i} + \tilde{u_i} + u_i^'</math>
|1}}


in which
in which


''<math>\bar{u_i}</math>'' = current (wave-averaged) velocity [m/s]
:<math>\bar{u_i}</math> = current (wave-averaged) velocity [m/s]


''<math>\tilde{u_i}</math>'' =  wave (oscillatory) velocity with wave-average ''<math>\bar{\tilde{u_i}} = 0 [m/s]</math>''
:<math>\tilde{u_i}</math> =  wave (oscillatory) velocity [m/s]with wave-average <math>\bar{{u_i}} = 0</math> below the wave trough


''<math>u_i^'</math>'' = turbulent fluctuation with ensemble average ''<math>\langle u_i^' \rangle</math>'' = 0 and wave average ''<math>\bar{u_i^'}</math>'' = 0 [m/s]
:<math>u_i^'</math> = turbulent fluctuation [m/s] with ensemble average <math>\langle u_i^' \rangle</math> = 0 and wave average <math>\bar{u_i^'}</math> = 0  


The wave-averaged total volume flux is defined as
The wave-averaged total volume flux is defined as


''<math>hV_i</math>'' = ''<math>\bar{{\int_z^\eta} {u_i dz }}</math>''
{{Equation|
<math> h V_i = \overline{{\int_{z_b}^\eta} {u_i dz }}</math>
|2}}


where
where


''<math>h</math>'' = wave-averaged water depth ''<math>h=\bar{\eta} - z_b </math>'' [m]
:<math>h</math> = wave-averaged water depth <math>h=\bar{\eta} - z_b </math> (Figure 1) [m]
 
''<math>V_i</math>'' = total mean mass flux velocity or simply total flux velocity for short [m/s]
 
''<math>u_i</math>'' = instantaneous current velocity [m/s]
 
''<math>\eta</math>''  = instantaneous water level with respect to the Still Water Level (SWL) [m]
 
''<math>z_b</math>''  = bed elevation with respect to the SWL [m]
 
For simplicity in the notation, the over bar in subsequent wave-averaged variables is omitted. The total flux velocity is also referred to as the mean transport velocity (Phillips 1977) and mass transport velocity (Mei 1983). The current volume flux is defined as
 
''<math>hU_i = \int^\eta_{z} \bar{u_i}dz</math>''    (2-3)
 
where ''<math>U_i</math>'' is the depth-averaged current velocity. Similarly, the wave volume flux is defined as by
 
''<math>Q_{wi} = hU_{wi} = \bar{\int \tilde{u_i} dz}</math>''
 
where ''<math>U_{wi}</math>'' is the depth-averaged wave flux velocity [m/s], and ''<math>\eta_t</math>'' = wave trough elevation [m]. Therefore the total flux velocity may be written as
 
''<math>V_i = U_i + U_{wi}</math>''
 
On the basis of the above definitions, and assuming depth-uniform cur-rents, the general depth-integrated and wave-averaged continuity and momentum equations may be written as (Phillips 1977; Mei 1983; Svendsen 2006)
 
''<math>\frac{\partial h}{\partial t} + \frac {\partial(hV_j)} {\partial x_j} = S_M</math>''        (2-6)
 
''<math>\frac{\partial(hV_i)}{\partial t} + \frac {\partial(hV_jV_i)}{\partial x_j} - \varepsilon_{ij}f_chV_j = -gh\frac{\partial \bar{\eta}}{\partial x_i} - \frac{h}{\rho} \frac{\partial p_{atm}}{\partial x_i}</math>''
 
''<math>+ \frac {\partial}{\partial x_j} {\left(v_{t}h \frac {\partial V_i} {\partial x_j} \right)} - \frac{1}{\rho} \frac{\partial} {\partial x_j} \left(S_{ij} + R_{ij} - \rho h U_{wi}U_{wj} \right) + \frac{\tau_{si}}{\rho} - m_{b}\frac{\tau_{bi}}{\rho}</math>''    (2-7)
 
where
 
t = time[s]
 
''<math>x_j</math>'' = Cartesian coordinate in the ''<math>j^{tj}</math>'' direction [m]
 
''<math>f_c</math>'' = Coriolis parameter [rad/s] ''<math> = 2\Omega sin\phi</math>'' where ''<math>\Omega = 7.29 x 10^{-5} </math>'' rad/s is the earth’s angular velocity of rotation and ''<math>\phi</math>'' is the latitude in degrees
 
h = wave-averaged total water depth ''<math> h = \zeta + \bar{\eta}</math>'' [m]
 
''<math>\bar{\eta} = </math>''wave-averaged water surface elevation with respect to refer-ence datum [m]
 
''<math>S_M = </math>'' water source/sink term due to precipitation, evaporation and structures (e.g. culverts) [m/s]
 
''<math>V_i = </math>''total flux velocity defined as ''<math> V_i = U_i + U_{wi}</math>'' [m/s]
 
''<math>U_i = </math>''wave- and depth-averaged current velocity [m/s]


''<math>U_{wi} = </math>''mean wave mass flux velocity or wave flux velocity for short [m/s]
:<math>V_i</math> = total mean mass flux velocity or simply total flux velocity [m/s]


''<math>g = </math>'' gravitational constant (~9.81 m/s<sup>2</sup>)
:<math>\eta</math> = instantaneous water level with respect to the Still Water Level (SWL) [m]


''<math>p_{atm} = </math>''atmospheric pressure [Pa]
:<math>\bar{\eta}</math> = wave-averaged water surface elevation with respect to the SWL (Figure 2-1) [m]


''<math>\rho = </math>'' water density (~1025 kg/m<sup>3</sup>)
:<math>z_b</math> = bed elevation with respect to the SWL (Figure 1) [m]


''<math>v_t = </math>'' turbulent eddy viscosity [m<sup>2</sup>/s]
<center>
[[File:fig_2_1.png|400px]]


''<math>\tau_{si} = </math>'' wind surface stress [Pa]
'''Figure 1. Vertical conventions used for the bed and mean water surface elevation.'''
</center>


''<math>S_{ij} = </math>'' wave radiation stress [Pa]
The total flux velocity is also referred to as the mean transport velocity (Phillips 1977) and mass transport velocity (Mei 1983). The current volume flux is defined as
{{Equation|
<math>hU_i = \int^\bar{\eta}_{z_b} \bar{u_i}dz</math>|3}}


''<math>R_{ij} = </math>'' surface roller stress [Pa]
where <math>U_i</math> is the depth-averaged current velocity. Similarly, the wave volume flux is defined as by
{{Equation|
<math>Q_{wi} = hU_{wi} = \overline{\int_{\eta_t}^\eta \tilde{u_i} dz}</math>
|4}}


''<math>m_b = </math>'' bed slope coefficient [-]
where


''<math>\tau_{bi} = </math>'' combined wave and current mean bed shear stress [Pa]
:<math>U_{wi}</math> = depth-averaged wave flux velocity [m/s]
:<math>\eta_t</math> = wave trough elevation [m]


The equations above are similar to those derived by Svendsen (2006), ex-cept for the inclusion of the water source/sink term in the continuity equation and the atmospheric pressure and surface roller terms in the momentum equation. It is also noted that the horizontal mixing term is formulated differently as a function of the total flux velocity, similar to the Generalized Lagrangian Mean (GLM) approach (Andrews and McIntyre 1978; Walstra et al. 2000). This approach is arguably more physically meaningful and also simplifies the discretization.


Figure 2-1 shows the vertical conventions used for mean water surface elevation and bed elevations. The arrows indicate the reference elevation for each variable. Note that the total water depth is always positive.
Therefore the total flux velocity may be written as
{{Equation|
<math>V_i = U_i + U_{wi}</math>
|5}}


=References=


Figure 2-1 here.
* Mei, C. 1989. The applied dynamics of ocean surface waves. New York: John Wiley.
[[File:fig_2_1.png]]


Figure 2.1 Vertical conventions used for mean water surface elevation and bed elevations.
* Phillips, O. M. 1977. The dynamics of the upper ocean. (2nd Edition). Cambridge, England: Cambridge University Press.
----
[[CMS#Documentation_Portal | Documentation Portal]]

Latest revision as of 15:26, 18 February 2015


The instantaneous current velocity ui is split into:

 

(1)

in which

= current (wave-averaged) velocity [m/s]
= wave (oscillatory) velocity [m/s]with wave-average below the wave trough
= turbulent fluctuation [m/s] with ensemble average = 0 and wave average = 0

The wave-averaged total volume flux is defined as

 

(2)

where

= wave-averaged water depth (Figure 1) [m]
= total mean mass flux velocity or simply total flux velocity [m/s]
= instantaneous water level with respect to the Still Water Level (SWL) [m]
= wave-averaged water surface elevation with respect to the SWL (Figure 2-1) [m]
= bed elevation with respect to the SWL (Figure 1) [m]

Fig 2 1.png

Figure 1. Vertical conventions used for the bed and mean water surface elevation.

The total flux velocity is also referred to as the mean transport velocity (Phillips 1977) and mass transport velocity (Mei 1983). The current volume flux is defined as

 

(3)

where is the depth-averaged current velocity. Similarly, the wave volume flux is defined as by

 

(4)

where

= depth-averaged wave flux velocity [m/s]
= wave trough elevation [m]


Therefore the total flux velocity may be written as

 

(5)

References

  • Mei, C. 1989. The applied dynamics of ocean surface waves. New York: John Wiley.
  • Phillips, O. M. 1977. The dynamics of the upper ocean. (2nd Edition). Cambridge, England: Cambridge University Press.

Documentation Portal