Page 357 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
May 25, 2005
0 521 85326 5
CB908/Date
0521853265appc
336
APPENDIX C. 2D CARTESIAN CODE
1 TERM=RHOP*VTTAU*ALOG(EYPLUS)*CMU**0.5/CAPPA/DELTA 11:59
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(LW.EQ.14)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(LW.EQ.16)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 *** EAST BOUNDARY
100 LE=NTAGE(I,J)
IF(LE.EQ.0)GO TO 200
AENOW=AE(I,J)
C INLET
IF(LE.EQ.21.OR.LE.EQ.27)THEN
AE(I,J)=0.0
IF(N.EQ.2)SU(I,J)=AENOW*U(I+1,J)+SU(I,J)
IF(N.EQ.3)SU(I,J)=AENOW*V(I+1,J)+SU(I,J)
IF(N.EQ.4)SU(I,J)=AENOW*W(I+1,J)+SU(I,J)
IF(N.EQ.5)SU(I,J)=AENOW*E(I+1,J)+SU(I,J)
IF(N.EQ.6)SU(I,J)=AENOW*D(I+1,J)+SU(I,J)
IF(N.EQ.7)SU(I,J)=AENOW*T(I+1,J)+SU(I,J)
SP(I,J)=AENOW+SP(I,J)
ENDIF
C SYMMETRY
IF(LE.EQ.22)THEN