APPENDIX APROGRAM LISTINGThe computational program EQUIL is listed in this appendix in Control DataSTAR-100 FORTRAN language ].4 (an extension of ANSI FORTRAN for the ControlData STAR-100 computer). This program consists of one main program and sixsubroutines.PROGRAM EQUIL(INPUT,OIJTPlIT,TAPE5=INPIJT,TAPF6=nlJTPIJT)DIMENSION CIS{500,24)COMMON/E1/T(5OO),P(5OO),WMIX(5OO),CL(5OO,6),SH(500)COMMON/E2/HI(500,24),MOLFF(500,24),CPI(500,24)CnMMON/E4/MW(24),SYMB(24),AA(24,6),MWEL(6),AAA(500,24,_)COMMON/ER/NE,NSCOMMON/EI4/VIS(500),CNND(500),SIG(500),CP(500)COMMON/EI5/RR,CRIT,RCRIT,MNS,MMM,MNE,MNFI,MXCOMMON/EI6/ENT(500,24),CPI1(500,24)FOtJIVALENCE (CIS(1,1),ENT(1,1))REAL MW,MWEL,MOI.EFRR=1.987CRIT=I.E-6RCRIT=.I*CRITMI=500C T IS IN DEGREE K AND P IS IN ATM.T(1:100)=4000.T(I01;100)=6000.T(201;300)=12570.P(1:MI)=6.3549C THE ELEMENTS ARE IN THE ORDER OF E-,H,HE,C,O AND N.CL 1,1;M1}=1.E-10CL 1,2;MI)=.4229649999CL 1,3;MI)=.1008CL 1,4;M1)=°4518CL 1,5:M1)=.02443CL 1,6:M1)=5.E-6CALL READ(MI)MNS=MImNSMMM=NE+IMNE:MI:_NEMNEI:MMM_MIMX=MMMmMMMmMICALL THERMO(MI)KTEST=OCALL CHEO(KTEST,MI)CALL TP(MI)DO 10 I=I,NSC CIS IS THE MASS FRACTION OF SPFCIE I. HI IS IN J/KG.CW=MW(1)CW1=4184./CWCIS(1,1:MI)=MOLEF(I,I:MI)mCW/WMIX(I:M1)10 HI(1,I:MI}=HI(I,I;M1)x:CWIC HERE SH IS IN J/KG, CP IS IN J/KG K, VlS IS IN N SEC/M2, COND ISC IN W/M K.SH(I;M1)=SH(1;M1)$4184.8
[1]
Sanford Gordon,et al.
Computer program for calculation of complex chemical equilibrium compositions
,
1972
.
[2]
C. Wilke.
A Viscosity Equation for Gas Mixtures
,
1950
.
[3]
Selmer M. Johnson,et al.
Chemical Equilibrium in Complex Mixtures
,
1958
.
[4]
W. E. Nicolet.
User's manual for RAD/EQUIL/1973: A general purpose radiation transport program
,
1973
.
[5]
C. Stroud,et al.
Chemical equilibrium of ablation materials including condensed species
,
1969
.
[6]
R. Pike,et al.
Thermodynamic properties in polynomial form for carbon, hydrogen, nitrogen, and oxygen systems from 300 to 15000 K
,
1970
.