Temp: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(9 intermediate revisions by the same user not shown)
Line 1: Line 1:
This is a page for setting up pages before loading them.
This is a page for setting up pages before loading them.


The discretized momentum equations are
<math> \int\limits_{A}{\nabla \cdot \left( {{\Gamma }^{\phi }}h\nabla \phi  \right)}\text{d}A=\oint\limits_{S}{{{\Gamma }^{\phi }}h\left( \nabla \phi \cdot \mathbf{n} \right)}\text{d}S=\sum\limits_{f}^{{}}{\bar{\Gamma }_{f}^{\phi }{{{\bar{h}}}_{f}}\Delta {{l}_{f}}{{\left( {{{\hat{n}}}_{i}}{{\nabla }_{i}}\phi  \right)}_{f}}} </math>
<math> \frac{\partial }{\partial t} \int_{A} h U_i dA
 
+ \oint_{F} \frac{\partial (h U_i U_j )}{\partial x_j} dF
First the momentum equations are rewritten as
\int_{A} \epsilon_{ij3} f_c U_j h = - g h \frac{\partial \eta }{\partial x_i} dA
{{Equation| <math> \frac{\partial ( h U_i ) }{\partial t}  
=  
+ \frac{\partial }{\partial x_j} \biggl( (h U_i U_j )- \nu_t  h \frac{\partial U_i }{\partial x_j} \biggr)
- g h \oint_{F} \frac{\partial \eta }{\partial x_i} dF
= - g h \frac{\partial \eta }{\partial x_i} + S_i
  - \oint_{F} \frac{h}{\rho_0} \frac{\partial p_a }{\partial x_i} dF
  </math>|2=1}}
+ \oint_{F} \frac{\partial }{\partial x_j} \biggl ( \nu_t h \frac{\partial U_i }{\partial x_j} \biggr ) dF
 
  + \frac{\tau_i }{\rho}
where <math>S_i</math> includes all other terms. The equation is then integrated over the a control volume as
</math>
{{Equation| <math> \frac{\partial }{\partial t} \int_A h U_i dA
+ \oint_F \frac{\partial }{\partial x_j} \biggl[ (h U_i U_j) -  \nu_t h \frac{\partial U_i }{\partial x_j} \biggr] dF
  = - g h \oint_F \frac{\partial \eta }{\partial x_i} dF + \int_A S_i dA
</math>|2=2}}
 
The resulting di




Line 36: Line 41:


for  <math>  j=1,2  </math>
for  <math>  j=1,2  </math>
{{Equation| <math> \frac{\partial ( h U_i ) }{\partial t} + \frac{\partial (h U_i U_j )}{\partial x_j}
- \epsilon_{ij3} f_c U_j h = - g h \frac{\partial \eta }{\partial x_i}
- \frac{h}{\rho_0} \frac{\partial p_a }{\partial x_i}
+ \frac{\partial }{\partial x_j} \biggl ( \nu_t  h \frac{\partial U_i }{\partial x_j} \biggr )
+ \frac{\tau_i }{\rho}
</math>|2=2}}

Latest revision as of 19:00, 17 March 2011

This is a page for setting up pages before loading them.

First the momentum equations are rewritten as

  (1)

where includes all other terms. The equation is then integrated over the a control volume as

  (2)

The resulting di


The continuity equation is discretized as


where the subscript indicates the cell face, with being the water surface elevation, is equal to the dot product of the velocity unit vector and the cell face unit vector.

The coefficient is equal to

The continuity equation is discretized as

where is the dot product of the cell face unit vector and


The depth-averaged 2-D continuity and momentum equations are given by

  (1)

for