Page 343 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
May 25, 2005
11:59
CB908/Date
0521853265appc
322
C 0 521 85326 5 APPENDIX C. 2D CARTESIAN CODE
STEP=2.0*R(JN)
XATTCH=0.0
COND=VISCOS*SPHEAT/PR(7)
DO 10 I=2,INM
XX=X(I)-XC(1)
IF(U(I,JNM).GT.0.0.AND.U(I-1,JNM).LT.0.0)THEN
BB=(U(I,JNM)-U(I-1,JNM))/DXMI(I)
AA=U(I,JNM)-BB*X(I)
XATTCH=-AA/BB-XC(1)
ENDIF
FLUX=QW(I,JN)
SHEAR=ABS(TAUW(I,JN))
UPLUS=U(I,JNM)/SQRT(SHEAR/DENSIT)
YPLUS=(Y(JN)-Y(JNM))*SQRT(SHEAR/DENSIT)/VISCOS
TW=T(I,JN)
ANUM=0.0
DEN=0.0
DEN1=0.0
DO 11 J=2,JNM
ANUM=ANUM+T(I,J)*ABS(U(I,J))*R(J)*DYP(J)
DEN=DEN+ABS(U(I,J))*R(J)*DYP(J)
DEN1=DEN1+U(I,J)*R(J)*DYP(J)
11 CONTINUE
TB=ANUM/DEN
UBAR=DEN1*2/R(JN)**2
ANU=FLUX/(TW-TB+SMALL)*2.0*R(JN)/COND
XX=XX/STEP
REY=UBAR*DENSIT*2*R(JN)/VISCOS
ANUTH=0.0123*REY**0.874*PR(7)**0.4
ANUR=ANU/ANUTH
WRITE(6,*)XX,ANUR,UPLUS,YPLUS
10 CONTINUE
WRITE(6,*)’ ANUTH = ’, ANUTH,’ REY = ’,REY

XATTCH=XATTCH/STEP
WRITE(6,*)’ ATTACHMENT X = ’,XATTCH
IBEG=2
IEND=35
ITOT=IEND-IBEG+1
OPEN(12,FILE=’EXPN.DAT’)
WRITE(12,*)’TITLE = EXPANSION - KRALL’