Temp: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
Deleted (talk | contribs)
No edit summary
Deleted (talk | contribs)
No edit summary
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.
{{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}}


The discretized momentum equations are  
The discretized momentum equations are  
<math> \frac{\partial }{\partial t} \int_{A} h U_i dA
<math> \frac{\partial }{\partial t} \int_{A} h U_i dA
+ \oint_{F} \frac{\partial (h U_i U_j )}{\partial x_j} dF
+ \oint_{F} \biggl{ \frac{ \partial }{\partial x_j} \big[
- \int_{A} \epsilon_{ij3} f_c U_j h = - g h \frac{\partial \eta }{\partial x_i} dA
  (h U_i U_j )
=
- \nu_t h \frac{\partial U_i }{\partial x_j} \big)
- g h \oint_{F} \frac{\partial \eta }{\partial x_i} dF
  - \frac{h}{\rho_0} \frac{\partial p_a }{\partial x_i}
  - \oint_{F} \frac{h}{\rho_0} \frac{\partial p_a }{\partial x_i} dF
 
+ \oint_{F} \frac{\partial }{\partial x_j} \biggl ( \nu_t  h \frac{\partial U_i }{\partial x_j} \biggr ) dF
- \epsilon_{ij3} f_c U_j h = - g h \frac{\partial \eta }{\partial x_i}
  + \frac{\tau_i }{\rho}
  + \frac{\tau_i }{\rho}
  </math>
  </math>
Line 36: Line 45:


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}}

Revision as of 00:05, 30 January 2011

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


  (hUi)t+(hUiUj)xjϵij3fcUjh=ghηxihρ0paxi+xj(νthUixj)+τiρ (2)


The discretized momentum equations are Failed to parse (syntax error): {\displaystyle \frac{\partial }{\partial t} \int_{A} h U_i dA + \oint_{F} \biggl{ \frac{ \partial }{\partial x_j} \big[ (h U_i U_j ) - \nu_t h \frac{\partial U_i }{\partial x_j} \big) - \frac{h}{\rho_0} \frac{\partial p_a }{\partial x_i} - \epsilon_{ij3} f_c U_j h = - g h \frac{\partial \eta }{\partial x_i} + \frac{\tau_i }{\rho} }


Ui,Pn+1=1ai,P(k=1ai,kUi,kn+1+Si)hPai,Pk=1nikΔskpkn+1

The continuity equation is discretized as

hn+1𝐒


where the subscript k indicates the cell face, p=gη with η being the water surface elevation, nik is equal to the dot product of the velocity unit vector and the cell face unit vector.

The coefficient ai,P is equal to ai,P=ai,k+aP0

The continuity equation is discretized as pPn+1=pPngΔtΔAPk=1nkFkn+1

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


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

  ht+(hUj)xj=S (1)

for j=1,2