CMS-Flow:Hydro Eqs: Difference between revisions
No edit summary |
No edit summary |
||
Line 1: | Line 1: | ||
<big> | <big> | ||
= Governing Equations = | |||
The depth-averaged 2-D continuity equation may be written as | The depth-averaged 2-D continuity equation may be written as | ||
{{Equation|<math> \frac{\partial h}{\partial t}+\nabla \cdot (h\mathbf{U})=S </math>|2=1}} | {{Equation|<math> \frac{\partial h}{\partial t}+\nabla \cdot (h\mathbf{U})=S </math>|2=1}} | ||
Line 10: | Line 10: | ||
where <math>g</math> is the gravitational constant, <math> \mathbf{B}=\left( \begin{matrix} 0 & {{f}_{c}} \\ -{{f}_{c}} & 0 \\ \end{matrix} \right) </math> where <math>f_{c}</math> is the Coriolis parameter, <math>\nu_t</math> is the eddy viscosity, <math> \tau_{wi} </math> is the wind stress, <math> \tau_{Si} </math> is the wave stresses, and <math> \tau_{bi} </math> is the combined wave-current mean bed shear stress. | where <math>g</math> is the gravitational constant, <math> \mathbf{B}=\left( \begin{matrix} 0 & {{f}_{c}} \\ -{{f}_{c}} & 0 \\ \end{matrix} \right) </math> where <math>f_{c}</math> is the Coriolis parameter, <math>\nu_t</math> is the eddy viscosity, <math> \tau_{wi} </math> is the wind stress, <math> \tau_{Si} </math> is the wave stresses, and <math> \tau_{bi} </math> is the combined wave-current mean bed shear stress. | ||
= Numerical Methods = | |||
== Temporal Term == | |||
The temporal term of the momentum equations is discretized using a first order implicit Euler scheme | |||
{{Equation| <math> \int\limits_{A}{\frac{\partial (h\phi )}{\partial t}}\text{d}A=\frac{\partial }{\partial t}\int\limits_{A}{(h\phi )\text{d}A}=\frac{{{h}^{n+1}}\phi _{{}}^{n+1}-{{h}^{n}}\phi _{{}}^{n}}{\Delta t}\Delta A </math>|2=3}} | |||
where <math> \Delta A </math> is the cell area, and <math> \Delta t </math> is the hydrodynamic time step. | |||
== Advection Term == | |||
The advection scheme obtained using the divergence theorem as | |||
where is the outward unit normal on cell face f, is the cell face length and is the total water depth linearly interpolated to the cell face. Here the overbar indicates a cell face interpolation operator described in the following section. For Cartesian grids the cell face unit vector is always aligned with one of the Cartesian coordinates which simplifies the calculation. Defining the cell face normal velocity as the above equation simplifies to | |||
{{Equation| <math> \int\limits_{A}{\nabla \cdot (h\mathbf{U}\phi )}\text{d}A=\oint\limits_{L}{h\phi \left( \mathbf{U}\cdot \mathbf{n} \right)}\text{d}L=\sum\limits_{f}^{{}}{{{{\bar{h}}}_{f}}\Delta {{l}_{f}}{{\left( {{{\hat{n}}}_{i}}{{U}_{i}} \right)}_{f}}{{{\tilde{\phi }}}_{f}}} </math>|2=4}} | |||
where <math> \mathbf{n}={{\hat{n}}_{i}}=({{\hat{n}}_{1}},{{\hat{n}}_{2}}) </math> is the outward unit normal on cell face f, <math> \Delta {{l}_{f}} </math> is the cell face length and <math> {{\bar{h}}_{f}} </math> is the total water depth linearly interpolated to the cell face. Here the overbar indicates a cell face interpolation operator described in the following section. For Cartesian grids the cell face unit vector is always aligned with one of the Cartesian coordinates which simplifies the calculation. Defining the cell face normal velocity as the above equation simplifies to | |||
---- | |||
{| border="1" | {| border="1" |
Revision as of 16:42, 6 April 2011
Governing Equations
The depth-averaged 2-D continuity equation may be written as
(1) |
where is the total water depth , is the water surface elevation, is the still water depth, is the depth-averaged current velocity, is a source term due to precipitation and evaporation, and is the divergence operator.
The momentum equation can be written as
(2) |
where is the gravitational constant, where is the Coriolis parameter, is the eddy viscosity, is the wind stress, is the wave stresses, and is the combined wave-current mean bed shear stress.
Numerical Methods
Temporal Term
The temporal term of the momentum equations is discretized using a first order implicit Euler scheme
(3) |
where is the cell area, and is the hydrodynamic time step.
Advection Term
The advection scheme obtained using the divergence theorem as where is the outward unit normal on cell face f, is the cell face length and is the total water depth linearly interpolated to the cell face. Here the overbar indicates a cell face interpolation operator described in the following section. For Cartesian grids the cell face unit vector is always aligned with one of the Cartesian coordinates which simplifies the calculation. Defining the cell face normal velocity as the above equation simplifies to
(4) |
where is the outward unit normal on cell face f, is the cell face length and is the total water depth linearly interpolated to the cell face. Here the overbar indicates a cell face interpolation operator described in the following section. For Cartesian grids the cell face unit vector is always aligned with one of the Cartesian coordinates which simplifies the calculation. Defining the cell face normal velocity as the above equation simplifies to
Symbol | Description | Units |
---|---|---|
Time | sec | |
Total water depth | m | |
Still water depth | m | |
Water surface elevation with respect to the still water elevation | m | |
Current velocity in the jth direction | m/sec | |
Sum of Precipitation and evaporation per unit area | m/sec | |
Gravitational constant | m/sec2 | |
Water density | kg/m3 | |
Atmospheric pressure | Pa | |
Turbulent eddy viscosity | m2/sec |