Page 367 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
0 521 85326 5
0521853265appc
346
APPENDIX C. 2D CARTESIAN CODE
TERM=(1-LN)/DYMI(J+1)+LN/(YC(J+1)-Y(J)) May 25, 2005 11:59
DN=(VISTN*PRTINV+VISN*PRINV)*RC(J+1)*DXP(I)*TERM
C *** CALCULATE CELL-PECLET NUMBERS
PECLW=CW/(DW+SMALL)
PECLE=CE/(DE+SMALL)
PECLS=CS/(DS+SMALL)
PECLN=CN/(DN+SMALL)
C *** CONVECTION SCHEMES
IF(UPWIND)THEN
AAW=1.0
AAE=1.0
AAS=1.0
AAN=1.0
ELSE IF(HYBRID)THEN
AAW=AMAX1(0.0,1.-0.5*ABS(PECLW))
AAE=AMAX1(0.0,1.-0.5*ABS(PECLE))
AAS=AMAX1(0.0,1.-0.5*ABS(PECLS))
AAN=AMAX1(0.0,1.-0.5*ABS(PECLN))
ELSE IF(POWER)THEN
AAW=AMAX1(0.0,(1.-0.1*ABS(PECLW))**5)
AAE=AMAX1(0.0,(1.-0.1*ABS(PECLE))**5)
AAS=AMAX1(0.0,(1.-0.1*ABS(PECLS))**5)
AAN=AMAX1(0.0,(1.-0.1*ABS(PECLN))**5)
ENDIF
C *** TOTAL COEFFICIENTS
AW(I,J)=DW*(AAW+AMAX1(PECLW,0.0))
AE(I,J)=DE*(AAE+AMAX1(-PECLE,0.0))
AS(I,J)=DS*(AAS+AMAX1(PECLS,0.0))
AN(I,J)=DN*(AAN+AMAX1(-PECLN,0.0))
1 CONTINUE
GO TO 2000
C COEFFICIENTS OF PRESSURE CORRECTION EQUATION
1000 DO 2 J=2,JNM
DO 2 I=2,INM
PP(I,J)=0.0
SP(I,J)=0.0
SU(I,J)=0.0
LW=NTAGW(I,J)/10
LE=NTAGE(I,J)/20
LS=NTAGS(I,J)/30
LN=NTAGN(I,J)/40
LB=1-NTAG(I,J)