Page 366 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
0 521 85326 5
0521853265appc
345
APPENDIX C. 2D CARTESIAN CODE
PRTINV=1.0/PRTN May 25, 2005 11:59
DO 1 J=2,JNM
DO 1 I=2,INM
SU(I,J)=0.0
SP(I,J)=0.0
C *** DIFFUSION COEFFICIENTS AND INTERPOLATED VALUES
LW=NTAGW(I,J)/10
LE=NTAGE(I,J)/20
LS=NTAGS(I,J)/30
LN=NTAGN(I,J)/40
C **** LAMINAR VISCOSITY
TERMW=(X(I)-XC(I))/VIS(I,J)+(XC(I)-X(I-1))/VIS(I-1,J)
VISW=DXMI(I)/TERMW*(1-LW)+LW*VIS(I-1,J)
TERME=(XC(I+1)-X(I))/VIS(I,J)+(X(I+1)-XC(I+1))/VIS(I+1,J)
VISE=DXMI(I+1)/TERME*(1-LE)+LE*VIS(I+1,J)
TERMS=(Y(J)-YC(J))/VIS(I,J)+(YC(J)-Y(J-1))/VIS(I,J-1)
VISS=DYMI(J)/(TERMS+SMALL)*(1-LS)+LS*VIS(I,J-1)
TERMN=(YC(J+1)-Y(J))/VIS(I,J)+(Y(J+1)-YC(J+1))/VIS(I,J+1)
VISN=DYMI(J+1)/TERMN*(1-LN)+LN*VIS(I,J+1)
C **** TURBULENT VISCOSITY
IF(TURBUL)THEN
TERMW=(X(I)-XC(I))/VIST(I,J)+(XC(I)-X(I-1))/VIST(I-1,J)
VISTW=DXMI(I)/TERMW*(1-LW)+LW*VIST(I-1,J)
TERME=(XC(I+1)-X(I))/VIST(I,J)+(X(I+1)-XC(I+1))/VIST(I+1,J)
VISTE=DXMI(I+1)/TERME*(1-LE)+LE*VIST(I+1,J)
TERMS=(Y(J)-YC(J))/VIST(I,J)+(YC(J)-Y(J-1))/VIST(I,J-1)
VISTS=DYMI(J)/TERMS*(1-LS)+LS*VIST(I,J-1)
TERMN=(YC(J+1)-Y(J))/VIST(I,J)+(Y(J+1)-YC(J+1))/VIST(I,J+1)
VISTN=DYMI(J+1)/TERMN*(1-LN)+LN*VIST(I,J+1)
ENDIF
C *** CONVECTION COEFFICIENTS
CW=FINTW(RHO,I,J)*FINTW(U,I,J)*R(J)*DYP(J)
CE=FINTE(RHO,I,J)*FINTE(U,I,J)*R(J)*DYP(J)
CS=FINTS(RHO,I,J)*FINTS(V,I,J)*RC(J)*DXP(I)
CN=FINTN(RHO,I,J)*FINTN(V,I,J)*RC(J+1)*DXP(I)
C **** DIFFUSION COEFFICIENTS (ALLOWANCE FOR BLOCKED REGIONS )
TERM=(1-LW)/DXMI(I)+LW/(X(I)-XC(I))
DW=(VISTW*PRTINV+VISW*PRINV)*R(J)*DYP(J)*TERM
TERM=(1-LE)/DXMI(I+1)+LE/(XC(I+1)-X(I))
DE=(VISTE*PRTINV+VISE*PRINV)*R(J)*DYP(J)*TERM
TERM=(1-LS)/DYMI(J)+LS/(Y(J)-YC(J))
DS=(VISTS*PRTINV+VISS*PRINV)*RC(J)*DXP(I)*TERM