Page 316 - Introduction to Computational Fluid Dynamics

P1: ICD
CB908/Date
0521853265appb
295
APPENDIX B. 1D CONDUCTION CODE
C EXACT SOLUTION 0 521 85326 5 May 11, 2005 15:43
AM=SQRT(HPREF*PERIM(2)/CONDREF/ACF(2))
QLOSS=SQRT(HPREF*PERIM(2)*CONDREF*ACF(2))*(T1-TINF)*TANH(AM*AL)
EFF=TANH(AM*AL)/(AM*AL)
WRITE(6,*)’ EXACT SOLUTION ’
WRITE(6,*)’ QLOSS = ’,QLOSS,’ EFF = ’,EFF
C NUMERICAL SOLUTION
QLOSS=ACF(2)*CONDREF*(T(1)-T(2))/(X(2)-X(1))
QMAX=2*AL*BREADTH*HPREF*(T(1)-TINF)
EFF=QLOSS/QMAX
WRITE(6,*)’ NUMERICAL SOLUTION ’
WRITE(6,*)’ QLOSS = ’,QLOSS,’ EFF = ’,EFF
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.,.TRUE.,.FALSE./
C --------------------------------------------
C CONTROL PARAMETERS
C FULLY IMPLICIT(PSI=1),FULLY EXPLICIT(PSI=0),SEMI IMPLICIT (0<PSI<1)
DATA PSI,DELT,MXSTEP,ITERMX,RP,CC/1.0,5,100,500,1.0,1E-5/
C --------------------------------------------
C BOUNDARY SPECIFICATION
DATA T1SPEC,Q1SPEC,H1SPEC/.TRUE.,2*.FALSE./
DATA TNSPEC,QNSPEC,HNSPEC/.FALSE.,.TRUE.,.FALSE./
DATA T1,TN,QB1,QBN,HB1,HBN/225.0,205.0,0.0,0.0,0.0,0.0/
DATA TINF,TINF1,TINFN,HPREF/25,0.0,0.0,15.0/
DATA CONDREF,RHOREF,SPHREF/45.0,1.0,1.0/
C --------------------------------------------
C GRID SPECIFICATION
DATA XCELL,XNODE/.TRUE.,.FALSE./
DATA N,AL/7,0.02/
DATA XCF/0.0,0.0,0.2,0.4,0.6,0.8,1.0,43*1.0/
C PROBLEM DEPENDENT PARAMETERS (IF ANY)
DATA BREADTH,THICK/0.2,0.002/
END