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


<math> U_{i,P}^{n+1} = \frac{1}{a_{i,P}} \biggl( \sum_{k=1} a_{i,k} U_{i,k}^{n+1} + S_i \biggr)  
<math> U_{i,P}^{n+1} = \frac{1}{a_{i,P}} \biggl( \sum_{k=1} a_{i,k} U_{i,k}^{n+1} + S_i \biggr)  
- \frac{h_P}{a_{i,P}} \sum_{k=1} n_i n_k  \Delta s_k p_k^{n+1}
- \frac{h_P}{a_{i,P}} \sum_{k=1} n_{ik} \Delta s_k p_k^{n+1}
</math>
</math>


where the subscript <math>k</math> indicates the cell face, <math>p = g \eta</math> with <math>\eta</math> being the water surface elevation, <math>n_i</math> is the unit vector in the <math>i</math> direction, and <math>n_k</math> is the unit vector normal to the cell face.  
where the subscript <math>k</math> indicates the cell face, <math>p = g \eta</math> with <math>\eta</math> being the water surface elevation, <math>n_{ik}</math> is equal to the dot product of the velocity unit vector and the cell face unit vector.  


The coefficient a_{i,P} is equal to <math> a_{i,P} = \sum a_{i,k} + a_P^0 </math>
The coefficient <math>a_{i,P}</math> is equal to <math> a_{i,P} = \sum a_{i,k} + a_P^0 </math>


The continuity equation is discretized as
The continuity equation is discretized as

Revision as of 23:29, 29 January 2011

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

The discretized momentum equations are

Ui,Pn+1=1ai,P(k=1ai,kUi,kn+1+Si)hPai,Pk=1nikΔskpkn+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

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

  ht+(hUj)xj=S (1)

for j=1,2

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