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




{{Equation| <math> \frac{\partial ( h U_i ) }{\partial t} + \frac{\partial (h U_i U_j )}{\partial x_j}
{{Equation| <math> \frac{\partial ( h U_i ) }{\partial t}  
- \epsilon_{ij3} f_c U_j h = - g h \frac{\partial \eta }{\partial x_i}
+ \frac{\partial }{\partial x_j} \biggl( (h U_i U_j )-  \nu_t  h \frac{\partial U_i }{\partial x_j} \biggr) dS
- \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{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}
  + \frac{\tau_i }{\rho}
  </math>|2=2}}
  </math>|2=2}}

Revision as of 00:12, 30 January 2011

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


  (hUi)t+xj((hUiUj)νthUixj)dSϵij3fcUjh=ghηxihρ0paxi+τ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