Page 359 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
11:59
0 521 85326 5
May 25, 2005
0521853265appc
CB908/Date
338
ELSE IF(N.EQ.3) THEN APPENDIX C. 2D CARTESIAN CODE
SU(I,J)=TMULT*VWAL+SU(I,J)
SP(I,J)=TMULT+SP(I,J)
ELSE IF(N.EQ.4) THEN
SU(I,J)=TMULT*WWAL+SU(I,J)
SP(I,J)=TMULT+SP(I,J)
ELSE IF(N.EQ.5) THEN
GENR=TMULT*VT/AREA*VT/DELTA
TERM=RHOP**2*CMU*ABS(E(I,J))/VISWAL
IF(YPLUS.GT.11.6)
1 TERM=RHOP*VTTAU*ALOG(EYPLUS)*CMU**0.5/CAPPA/DELTA
SU(I,J)=GENR*VOLP+SU(I,J)
SP(I,J)=TERM*VOLP+SP(I,J)
ELSE IF(N.EQ.6) THEN
TERM=VTTAU**3/CAPPA/DELTA
SU(I,J)=GREAT*TERM
SP(I,J)=GREAT
ELSE IF(N.EQ.7) THEN
IF(TURBUL)THEN
UPLUS=ABS(VT)/VTTAU
STANTON=STAN(UPLUS,YPLUS,PR(7),PRT(7))
TERM=RHOP*ABS(VT)*AREA*STANTON
ELSE
TERM=VISWAL/(PR(N)*DELTA)*AREA
ENDIF
IF(LE.EQ.24)THEN
SU(I,J)=TERM*T(I+1,J)+SU(I,J)
SP(I,J)=TERM+SP(I,J)
QW(I+1,J)=TERM/AREA*(T(I+1,J)-T(I,J))*SPH(I+1,J)
ELSE IF(LE.EQ.26)THEN
SU(I,J)=QW(I+1,J)*AREA/SPH(I+1,J)+SU(I,J)
T(I+1,J)=QW(I+1,J)/TERM*AREA/SPH(I+1,J)+T(I,J)
ENDIF
ENDIF
ENDIF
C *** SOUTH BOUNDARY
200 LS=NTAGS(I,J)
IF(LS.EQ.0)GO TO 300
ASNOW=AS(I,J)
C INLET
IF(LS.EQ.31.OR.LS.EQ.37)THEN
AS(I,J)=0.0