Page 320 - Introduction to Computational Fluid Dynamics

P1: ICD
May 11, 2005
0 521 85326 5
CB908/Date
0521853265appb
299
APPENDIX B. 1D CONDUCTION CODE
DATA PSI,DELT,MXSTEP,ITERMX,RP,CC/1.0,5,100,500,1.0,1E-5/ 15:43
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/200.0,150.0,0.0,0.0,0.0,0.0/
DATA TINF,TINF1,TINFN,HPREF/25,0.0,0.0,20.0/
DATA CONDREF,RHOREF,SPHREF/1.0,1.0,1.0/
C --------------------------------------------
C GRID SPECIFICATION
DATA XCELL,XNODE/.TRUE.,.FALSE./
DATA N/8/
C PROBLEM DEPENDENT PARAMETERS (IF ANY)
DATA THICK/0.001/
DATA R1,R2,R3/0.0125,0.025,0.0375/
END

Library File LIB1D.FOR
C *************************************************
C THIS IS LIBRARY LIB1D.FOR - A. W. DATE
C *************************************************
SUBROUTINE MAIN
INCLUDE ’COM1D.FOR’
C *************************************************
WRITE(6,*)’ ***************************************’
IF(THOMAS)WRITE(6,*)’ SOLUTION BY TDMA’
IF(GAUSS)WRITE(6,*)’ SOLUTION BY GAUSS SIEDEL’
WRITE(6,*)’ ***************************************’
C*** CALCULATE CELL FACE COORDINATES, AREA AND VOLUME.
CALL GRID
C*** SPECIFY INITIAL TEMPERATURE DISTRIBUTION (USER FILE)
CALL INIT
ISTOP=0
IF(STEADY)PSI=1.0
IF(UNSTEADY)THEN
DO 101 I=1,N
101 TO(I)=T(I)
IF(PSI.EQ.0.0)ITERMX=0
ENDIF
TTIME=0.0
C*** BEGIN TIME STEP
TIMEMX=MXSTEP*DELT