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 11: Line 11:


The discretized momentum equations are  
The discretized momentum equations are  
<math> \frac{\partial }{\partial t} \int_{A} h U_i dA </math>
<math> \frac{\partial }{\partial t} \int_A h U_i dA  
<math> + \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) \biggr} </math>
+ \oint_A \frac{\partial (h U_i U_j )}{\partial x_j} dS
<math> - \frac{h}{\rho_0} \frac{\partial p_a }{\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 )
- \epsilon_{ij3} f_c U_j h = - g h \frac{\partial \eta }{\partial x_i}
- \epsilon_{ij3} f_c U_j h = - g h \frac{\partial \eta }{\partial x_i}
  + \frac{\tau_i }{\rho}
  + \frac{\tau_i }{\rho}

Revision as of 00:09, 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 tAhUidA+A(hUiUj)xjdShρ0paxi+xj(νthUixj)ϵij3fcUjh=ghηxi+τiρ


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