Page 319 - Introduction to Computational Fluid Dynamics

P1: ICD
15:43
May 11, 2005
CB908/Date
0521853265appb
298
RETURN 0 521 85326 5 APPENDIX B. 1D CONDUCTION CODE
END
C --------------------------------------------
FUNCTION PERI(XX)
INCLUDE ’COM1D.FOR’
C PERIMETER
PERI=4*PI*(R1+XX)
RETURN
END
C *************************************************
SUBROUTINE RESULT
INCLUDE ’COM1D.FOR’
CHARACTER*20 HEADER
C *************************************************
HEADER=’ FINAL-TEMP ’
CALL PRINT(T,HEADER)
C EXTRACT PROBLEM DEPENDENT PARAMETERS IF ANY
QLOSS=ACF(2)*COND(1)*(T(1)-T(2))/(X(2)-X(1))
QMAX=2*PI*(R3**2-R1**2)*HPREF*(T(1)-TINF)
EFF=QLOSS/QMAX
WRITE(6,*)’ NUMERICAL SOLUTION ’
WRITE(6,*)’ QLOSS = ’,QLOSS,’ EFF = ’,EFF
C PLOT TEMP PROFILE
OPEN(12,FILE=’TEXT3.DAT’)
WRITE(12,*)’TITLE = ANNULAR FIN’
WRITE(12,*)’VARIABLES = XX TT ’
WRITE(12,*)’ZONE T = ZONE1, I = ’,N,’ ,F = POINT’
DO 51 J=1,N
51 WRITE(12,*)X(J),T(J)
CLOSE(12)
RETURN
END
C *************************************************
BLOCK DATA
INCLUDE ’COM1D.FOR’
C *************************************************
C LOGICAL DECLARATIONS
C *** DECLARE STEADY OR UNSTEADY AND SOLUTION METHOD
DATA STEADY,UNSTEADY,GAUSS,THOMAS/.TRUE.,.FALSE.,.FALSE.,.TRUE./
C --------------------------------------------
C CONTROL PARAMETERS
C FULLY IMPLICIT (PSI=1),FULLY EXPLICIT (PSI=0),SEMI IMPLICIT (0<PSI<1)