Page 358 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
0 521 85326 5
0521853265appc
337
APPENDIX C. 2D CARTESIAN CODE
IF(N.EQ.2)SP(I,J)=AENOW+SP(I,J) May 25, 2005 11:59
AE(I,J)=0.0
IF(N.EQ.2)U(I+1,J)=0.0
IF(N.EQ.3)V(I+1,J)=V(I,J)
IF(N.EQ.4)W(I+1,J)=W(I,J)
IF(N.EQ.5)E(I+1,J)=E(I,J)
IF(N.EQ.6)D(I+1,J)=D(I,J)
IF(N.EQ.7)T(I+1,J)=T(I,J)
ENDIF
C EXIT
IF(LE.EQ.23.OR.LE.EQ.25) THEN
AE(I,J)=0.0
RATIO=(XC(I+1)-X(I))/DXMI(I)
IF(LE.EQ.23)RATIO=0.0
IF(N.EQ.2)U(I+1,J)=U(I,J)+RATIO*(U(I,J)-U(I-1,J))
IF(N.EQ.3)V(I+1,J)=V(I,J)+RATIO*(V(I,J)-V(I-1,J))
IF(N.EQ.4)W(I+1,J)=W(I,J)+RATIO*(W(I,J)-W(I-1,J))
IF(N.EQ.5)E(I+1,J)=E(I,J)+RATIO*(E(I,J)-E(I-1,J))
IF(N.EQ.6)D(I+1,J)=D(I,J)+RATIO*(D(I,J)-D(I-1,J))
IF(N.EQ.7)T(I+1,J)=T(I,J)+RATIO*(T(I,J)-T(I-1,J))
ENDIF
C WALL
IF(LE.EQ.24.OR.LE.EQ.26) THEN
AE(I,J)=0.0
DELTA=XC(I+1)-X(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)=AENOW*UWAL+SU(I,J)
SP(I,J)=AENOW+SP(I,J)