Page 318 - Introduction to Computational Fluid Dynamics

P1: ICD
0 521 85326 5
CB908/Date
0521853265appb
297
APPENDIX B. 1D CONDUCTION CODE
C COND(I) AND SPH(I) ARE DEFINED AT NODE P May 11, 2005 15:43
RR=R2-R1
DO 1 I=1,N
IF(X(I).LT. RR)COND(I)=200
IF(X(I).GT. RR)COND(I)=40
1 SPH(I)=SPHREF
RETURN
END
C *************************************************
SUBROUTINE SORCE
INCLUDE ’COM1D.FOR’
C *************************************************
C FORM PROBLEM DEPENDENT SOURCE TERM INCLUDING SU AND SP
DO 1 I=2,N-1
TERM=HPERI(I)*PERIM(I)*(XCF(I+1)-XCF(I))
SU(I)=SU(I)+TERM*TINF
SP(I)=SP(I)+TERM
1 CONTINUE
RETURN
END
C *************************************************
SUBROUTINE INTPRI
INCLUDE ’COM1D.FOR’
C *************************************************
RETURN
END
C *************************************************
C FUNCTION ROUTINES
C *************************************************
FUNCTION HPERI(II)
INCLUDE ’COM1D.FOR’
C H AT PERIMETER
I=II
HPERI=HPREF+X(I)*0.0+T(I)*0.0
RETURN
END
C --------------------------------------------
FUNCTION AREA(XX)
INCLUDE ’COM1D.FOR’
C AREA OF CROSS-SECTION
AREA=2*PI*(R1+XX)*THICK