Page 361 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
0 521 85326 5
0521853265appc
340
APPENDIX C. 2D CARTESIAN CODE
VTTAU=CMU**0.25*SQRT(ABS(E(I,J))) May 25, 2005 11:59
YPLUS=VTTAU*DELTA*RHOP/VISWAL
EYPLUS=ELOG*YPLUS
TMULT=VISWAL*AREA/DELTA
TAUW(I,J-1)=VISWAL*VT/DELTA
IF(TURBUL)TAUW(I,J-1)=RHO(I,J-1)*VTTAU**2
IF(YPLUS.GT.11.6)TMULT=RHOP*CAPPA*VTTAU*AREA/ALOG(EYPLUS)
IF(N.EQ.2) THEN
SU(I,J)=TMULT*UWAL+SU(I,J)
SP(I,J)=TMULT+SP(I,J)
ELSE IF(N.EQ.3) THEN
SU(I,J)=ASNOW*VWAL+SU(I,J)
SP(I,J)=ASNOW+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(LS.EQ.34)THEN
SU(I,J)=TERM*T(I,J-1)+SU(I,J)
SP(I,J)=TERM+SP(I,J)
QW(I,J-1)=TERM/AREA*(T(I,J-1)-T(I,J))*SPH(I,J-1)
ELSE IF(LS.EQ.36)THEN
SU(I,J)=QW(I,J-1)*AREA/SPH(I,J-1)+SU(I,J)
T(I,J-1)=QW(I,J-1)/TERM*AREA/SPH(I,J-1)+T(I,J)