Page 356 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
0 521 85326 5
0521853265appc
335
APPENDIX C. 2D CARTESIAN CODE
IF(LW.EQ.13.OR.LW.EQ.15) THEN May 25, 2005 11:59
AW(I,J)=0.0
RATIO=(X(I)-XC(I))/DXMI(I+1)
IF(LW.EQ.13)RATIO=0.0
IF(N.EQ.2)U(I-1,J)=U(I,J)-RATIO*(U(I+1,J)-U(I,J))
IF(N.EQ.3)V(I-1,J)=V(I,J)-RATIO*(V(I+1,J)-V(I,J))
IF(N.EQ.4)W(I-1,J)=W(I,J)-RATIO*(W(I+1,J)-W(I,J))
IF(N.EQ.5)E(I-1,J)=E(I,J)-RATIO*(E(I+1,J)-E(I,J))
IF(N.EQ.6)D(I-1,J)=D(I,J)-RATIO*(D(I+1,J)-D(I,J))
IF(N.EQ.7)T(I-1,J)=T(I,J)-RATIO*(T(I+1,J)-T(I,J))
ENDIF
C WALL
IF(LW.EQ.14.OR.LW.EQ.16) THEN
AW(I,J)=0.0
DELTA=X(I)-XC(I)
AREA=R(J)*DYP(J)
UWAL=U(I-1,J)
VWAL=V(I-1,J)
WWAL=W(I-1,J)
VISWAL=VIS(I-1,J)
ANG=ATAN(W(I,J)/(V(I,J)+SMALL))
VT=(V(I,J)-VWAL)*COS(ANG)+(W(I,J)-WWAL)*SIN(ANG)
VTTAU=CMU**0.25*SQRT(ABS(E(I,J)))
YPLUS=VTTAU*DELTA*RHOP/VISWAL
EYPLUS=ELOG*YPLUS
TMULT=VISWAL*AREA/DELTA
TAUW(I-1,J)=VISWAL*VT/DELTA
IF(TURBUL)TAUW(I-1,J)=RHO(I-1,J)*VTTAU**2
IF(YPLUS.GT.11.6)TMULT=RHOP*CAPPA*VTTAU*AREA/ALOG(EYPLUS)
IF(N.EQ.2) THEN
SU(I,J)=AWNOW*UWAL+SU(I,J)
SP(I,J)=AWNOW+SP(I,J)
ELSE IF(N.EQ.3) THEN
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)