Temp

From CIRPwiki
Revision as of 00:16, 30 January 2011 by Deleted (talk | contribs)
Jump to navigation Jump to search

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

First the momentum equations are rewritten as

  (hUi)t+xj((hUiUj)νthUixj)=ghηxi+Si (2)

where Siincludesallotherterms. The equation is then integrated over the a control volume tAhUidA+Fxj[(hUiUj)νthUixj]dF=ghFηxidF+ASidA


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