Question

cryogenic engineering subject research question.

Develop a mathematical model of Claude system driving pulse tube refrigerator using R507 refrigerant Nd discuss its impact on system.....

Answer 1

C *************************************************

PROGRAM ONED

INCLUDE ’COM1D.FOR’

C *************************************************

OPEN(6,FILE=’OO’)

WRITE(6,*)’ **********************************************’

WRITE(6,*)’ ADHESION OF PLASTIC SHEETS - PROB1-CHAPTER2’

WRITE(6,*)’ **********************************************’

DO 1 I=1,N

1 XCF(I)=XCF(I)*AL

CALL MAIN

CALL RESULT

STOP

END

C *************************************************

SUBROUTINE INIT

INCLUDE ’COM1D.FOR’

C *************************************************

C GIVE INITIAL GUESS AT TIME=0.0 OR AT ITER=0 FOR STEADY STATE

TIN=30

DO 1 I=1,N

T(I)=30

IF(I.EQ.1.OR.I.EQ.N)T(I)=250

1 CONTINUE

RETURN

END

C *************************************************

SUBROUTINE NEWVAL

INCLUDE ’COM1D.FOR’

C SET NEW VALUES OF HB1,HBN,QB1,QBN,TINF1,TINFN OR SOURCES

RETURN

END

C *************************************************

SUBROUTINE PROPS

INCLUDE ’COM1D.FOR’

C *************************************************

C COND(I) AND SPH(I) ARE DEFINED AT NODE P

DO 1 I=1,N

RHO(I)=RHOREF

COND(I)=CONDREF

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

SU(I)=SU(I)+0.0

1 CONTINUE

RETURN

END

C *************************************************

SUBROUTINE INTPRI

INCLUDE ’COM1D.FOR’

CHARACTER*20 HEADER

C *************************************************

WRITE(6,*)’ TIMESTEP = ’,NTIME,’ TOTAL TIME = ’,TTIME

C PRINT TEMPERATURES AT THE CURRENT STEP

HEADER=’ TEMP ’

CALL PRINT(T,HEADER)

C STORE MID-POINT TEMPERATURE

DUM1(NTIME)=T(4)

C GIVE STOP CONDITION

IMID=4

IF(T(IMID).GT.140)ISTOP=1

RETURN

END

C *************************************************

C FUNCTION ROUTINES

C *************************************************

FUNCTION HPERI(II)

INCLUDE ’COM1D.FOR’

C H AT PERIMETER

I=II

HPERI=HPREF*0.0+X(I)*0.0+T(I)*0.0

RETURN

END

C --------------------------------------------

FUNCTION AREA(XX)

INCLUDE ’COM1D.FOR’

C AREA OF CROSS-SECTION

AREA=1.0+0.0*XX

RETURN

END

C --------------------------------------------

FUNCTION PERI(XX)

INCLUDE ’COM1D.FOR’

C PERIMETER

PERI=0*XX

RETURN

END

C *************************************************

SUBROUTINE RESULT

INCLUDE ’COM1D.FOR’

CHARACTER*20 HEADER

C *************************************************

HEADER=’ FINAL-TEMP ’

CALL PRINT(T,HEADER)

HEADER=’ X(I) ’

CALL PRINT(X,HEADER)

HEADER=’ XCF(I) ’

CALL PRINT(XCF,HEADER)

C EXTRACT PROBLEM DEPENDENT PARAMETERS IF ANY

WRITE(6,*)’ PRINT MID-POINT TEMPERATURE’

DO 1 I=1,NTIME

TT=FLOAT(I)*DELT

1 WRITE(6,*)TT,DUM1(I)

TNOW=DUM1(NTIME)

TOLD=DUM1(NTIME-1)

TT=FLOAT(NTIME-1)*DELT

TIME=(140-TOLD)/(TNOW-TOLD)*DELT+TT

WRITE(6,*)’ TIME FOR ADHESION = ’,TIME

RETURN

END

C *************************************************

BLOCK DATA

INCLUDE ’COM1D.FOR’

C *************************************************

C LOGICAL DECLARATIONS

DATA STEADY,UNSTEADY,GAUSS,THOMAS/.FALSE.,.TRUE.,.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/0.0,10,10000,500,1.0,1E-5/

C --------------------------------------------

C BOUNDARY SPECIFICATION

DATA T1SPEC,Q1SPEC,H1SPEC/.TRUE.,2*.FALSE./

DATA TNSPEC,QNSPEC,HNSPEC/.TRUE.,2*.FALSE./

DATA T1,TN,QB1,QBN,HB1,HBN/250.0,250.0,0.0,0.0,0.0,0.0/

C DATA TINF,TINF1,TINFN,HPREF/25,150,250,12.0/

DATA CONDREF,RHOREF,SPHREF/0.25,1300,2000.0/

C --------------------------------------------

C GRID SPECIFICATION

DATA XCELL,XNODE/.TRUE.,.FALSE./

DATA N,AL/7,0.01/

DATA XCF/0.0,0.0,0.2,0.4,0.6,0.8,1.0,43*1.0/

END