190 likes | 276 Views
o. c. v. +. u. z. dsu ( i,j ). o. +. i ,j. dsz ( i,j ). o. +. i ,j. dsv ( i,j ). o. +. i ,j. dsc ( i,j ). o. +. i ,j. dnu ( i,j ). o. +. i ,j. dnz ( i,j ). o. +. i ,j. d nv ( i,j ). o. +. i ,j. dnc ( i,j ). o. +. i ,j. dsdnu ( i,j ). o. +.
E N D
o c v + u z
dsu(i,j) o + i,j
dsz(i,j) o + i,j
dsv(i,j) o + i,j
dsc(i,j) o + i,j
dnu(i,j) o + i,j
dnz(i,j) o + i,j
dnv(i,j) o + i,j
dnc(i,j) o + i,j
dsdnu(i,j) o +
Momentum balance in u direction • Only consider advection of velocity component in the local grid direction
Continuity equation • Multiply continuity equation with u
Subtract continuity eq. times u from momentum equation and divide by V
Procedure • Go around the cell centered at u-point • Compute q across cell boundary by averaging qx resp. qy • If q is inward • Compute Qin by multiplying with dnzrespdsc • Compute uin as • Add Qin (uin–u) to advection