Page 369 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
0 521 85326 5
0521853265appc
348
APPENDIX C. 2D CARTESIAN CODE
AN(I,J)=RHON*(RC(J+1)*DXP(I))**2/SUMN*2.0*LB May 25, 2005 11:59
ENDIF
2 CONTINUE
2000 CONTINUE
RETURN
END
C *******************************************
SUBROUTINE SORCE(NNV)
INCLUDE ’COM2D.FOR’
C *******************************************
DIMENSION PROD(IT,JT)
N=NNV
GO TO (10,20,30,40,50,60,70),N
C *** FOR PRESSURE CORRECTION
10 DO 11 J=2,JNM
DO 11 I=2,INM
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)
SM=CE-CW+CN-CS
IF(UNSTDY)SM=SM+(RHO(I,J)-RHOO(I,J))/DELT*VOL(I,J)
SU(I,J)=SU(I,J)-SM*(1-NTAG(I,J))
11 CONTINUE
GO TO 1000
C *** FOR U-VELOCITY
20 DO 21 J=2,JNM
DO 21 I=2,INM
DPDX=(FINTE(P,I,J)-FINTW(P,I,J))/DXP(I)
SU(I,J)=SU(I,J)-DPDX*VOL(I,J)*(1-NTAG(I,J))
21 CONTINUE
GO TO 1000
C *** FOR V-VELOCITY
30 DO 31 J=2,JNM
DO 31 I=2,INM
DPDY=(FINTN(P,I,J)-FINTS(P,I,J))/DYP(J)
SU(I,J)=SU(I,J)-DPDY*VOL(I,J)*(1-NTAG(I,J))
VISP=VIS(I,J)+VIST(I,J)
IF(AXISYMM)SP(I,J)=SP(I,J)+VISP/R(J)**2*VOL(I,J)
31 CONTINUE
GO TO 1000
C *** FOR W-VELOCITY