home *** CD-ROM | disk | FTP | other *** search
/ The World of Computer Software / World_Of_Computer_Software-02-387-Vol-3of3.iso / h / house_ii.zip / FOR / FUEL.FOR < prev    next >
Text File  |  1992-04-08  |  7KB  |  182 lines
  1.       SUBROUTINE FUEL 
  2. C *** CALCULATES PROPERTIES OF COMBUSTION FOR GAS FURNACE PROGRAM 
  3.       IMPLICIT REAL(A-H,M-Z)
  4.       DIMENSION MW(14),HF(14),A(14),B(14),C(14),D(14) 
  5.      +,NOC(13)
  6.       PARAMETER (CPAIR=.24) 
  7. CMDK BLKCM1
  8. CMDK BLKCM2
  9. CMDK FUEL1
  10. CMDK PRT1
  11.       PVF(Z)=EXP(16.72048-9400.7858/(Z+459.7))
  12.       DATA A/8*0.,  -.168288E-6,  -.438419E-7,  -.509447E-7,
  13.      +       0.,  -.536138E-7,  0./ 
  14.       DATA B/6.09436E-3,  10.1670E-3,  14.1216E-3,  16.7486E-3, 
  15.      +       17.4328E-3,  21.4615E-3,  21.0380E-3,  40.2128E-3, 
  16.      +       .159751E-2,  .518974E-3,  .555976E-3,  .1488E-3, 
  17.      +       .999914E-3,  -3.018458E-3/ 
  18.       DATA C/6.96554,  10.8508,  15.722, 21.4222, 21.0855,
  19.      +       26.0134,  26.1784,  29.1591,  9.46994,  6.65257, 
  20.      +       7.05349,  6.99778,  7.52025,  -9.69913/
  21.       DATA D/4152.99,  8577.87,  12169.23,  14986.09,  15514.83,
  22.      +       17899.39,  18739.14,  21833.76,  3209.73,  3226.06,
  23.      +       3150.49,  3103.78,  3703.66,  19683.4/ 
  24.       DATA HF/-32200.,  -36425.,   -44676.,  -57870.,  -54270., 
  25.      +        -66456.,  -63000.,  -71928.,  -169293.,  3*0.,
  26.      +        -104036.,  -112971./
  27.       DATA MW/16.042,  30.068,  44.094,  2*58.120,  2*72.146, 
  28.      +        86.172,  44.010,  28.016,  32.,  2.016,2*18.016/
  29. C *** WATER IN FUEL 
  30.       PH2O=PVF(TAIN)
  31.       NH2OP= PH2O/BARPSI
  32.       NH2OF= NH2OP/(1.-NH2OP) 
  33. C *** VOLUME OF DRY FUEL PER HOUR 
  34.       VFDRY= (1.-NH2OP)*FUELCM*60.
  35. C *** MOLECULAR WEIGHT OF DRY FUEL
  36.       MWF=0.
  37.       DO 1 I=1,12 
  38. 1     MWF= MWF+XF(I)*MW(I)
  39.       SGFUEL=MWF/28.85
  40. C *** MASS & MOLES OF DRY FUEL PER HOUR AND WET FUEL
  41.       MF= SGFUEL*(.07664)*VFDRY 
  42.       WNF= MF/MWF 
  43.       WFUEL= MF+18.016*NH2OF*WNF
  44. C *** AIR HUMDITY 
  45.       WSAT= (.90029E-5*TWBI-.68134E-3)*TWBI+.01939
  46.       W= ((1093.-.556*TWBI)*WSAT-CPAIR*(TCOMB-TWBI))/(1093.-.444*TWBI)
  47.       IF ( IQEXT .EQ. 0) W= WSAT
  48. C *** MOLES OF PRODUCTS OF COMBUSTION LESS THOSE DUE TO AIR 
  49.       FNO(9)=XF(1)+2.*XF(2)+3.*XF(3)+4.*(XF(4)+XF(5))+5.*(XF(6)+XF(7))+ 
  50.      + 6.*XF(8)+XF(9) 
  51.       FNO(10)= XF(10) 
  52. C  FNO(11) IS NEG. OF NO. OF MOLES OF O2/MOLE OF FUEL(=-21 A/79)
  53.       FNO(11)= -2.*XF(1)-3.5*XF(2)-5.*XF(3)-6.5*(XF(4)+XF(5))-8.*(XF(6) 
  54.      + +XF(7))-9.5*XF(8)-.5*XF(12)+XF(11) 
  55.       FNO(12)=0.
  56.       FNO(13)= NH2OF+2.*XF(1)+3.*XF(2)+4.*XF(3)+5.*(XF(4)+XF(5))+6.*
  57.      + (XF(6)+XF(7))+7.*XF(8)+XF(12)
  58. C *** STOICHIOMETRIC AIR
  59.       WAIR= -FNO(11)/.21*WNF*28.85
  60.       WSTO= WAIR*(1.+W) 
  61.       WSA=(-1.)*FNO(11)/.21 
  62. C        WMSA IS NO. OF MOLES OF DRY AIR/MOLE FUEL
  63.       WMSA=WSA*(1.+EXCESS)
  64. C *** FUEL ENERGY IN
  65.       ENINF= 0. 
  66.       XF(13)=NH2OF
  67.       DO 2 I=1,13 
  68. 2     ENINF= ENINF+ (H(I,TAIN)-H(I,77.)+HF(I))*XF(I)
  69. C *** ENTHALPY PER MOLE OF DRY AIR
  70.       ENAIR=0.
  71.       FN(10)=.79
  72.       FN(11)= .21 
  73.       FN(12)=0. 
  74.       FN(13)= 7.622*.21*W 
  75.       DO 3 I=10,13
  76. 3     ENAIR= ENAIR+(H(I,TCOMB)-H(I,77.)+HF(I))*FN(I)
  77. C *** CONSTANT PORTION OF FLUE ENTHALPIES 
  78.       AS=0. 
  79.       BS=0. 
  80.       CS=0. 
  81.       DS=0. 
  82.       DO 4 I=9,13 
  83.       AS= AS+A(I)*FNO(I)
  84.       BS= BS+B(I)*FNO(I)
  85.       CS= CS+C(I)*FNO(I)
  86. 4     DS= DS+(D(I)-H(I,77.)+HF(I))*FNO(I) 
  87.       AAIR= A(10)*.79+A(11)*.21+A(13)*FN(13)
  88.       BAIR= B(10)*.79+ B(11)*.21+B(13)*FN(13) 
  89.       CAIR= C(10)*.79+C(11)*.21+C(13)*FN(13)
  90.       DAIR= (D(10)-H(10,77.)+HF(10))*.79+(D(11)-H(11,77.)+HF(11))*.2
  91.      +1+(D(13) -H(13,77.)+HF(13))*FN(13)
  92. C *** TOTAL FLUE ENTHALPY PER MOLE OF DRY FUEL IS OF THE FORM: H=((AS*T+
  93. C *** BS)*T+CS)*T +DS+(((AAIR*T+BAIR)*T+CAIR)*T+DAIR)*NAIR: H2O IS
  94. C *** IGNORED.FOR TOTAL HEAT INPUT EVALUATE H AT TCOMB BY USING ABOVE AN
  95. C *** H2O CONDENSED 
  96.       HFACT= ((AS*TCOMB+BS)*TCOMB+CS)*TCOMB+DS
  97.       HFACTA=((AAIR*TCOMB+BAIR)*TCOMB+CAIR)*TCOMB+DAIR
  98.       HFG= (B(14)*TCOMB+C(14))*TCOMB+D(14)
  99.       PSAT=PVF(TCOMB) 
  100.       FA= PSAT/BARPSI 
  101.       NOT=FNO(9)+FNO(10)+FNO(11)+FNO(13)
  102.       FB=(FA*(NOT+WMSA*(1.+FN(13)))-FNO(13)-WMSA* 
  103.      +FN(13))/(FA-1.)*HFG 
  104.       IF(ISAVE.EQ.0)RETURN
  105. C 
  106. C  SAVE VALUES FOR FLUE CALCULATIONS
  107.       WNFC=WNF
  108.       WFLUE=WSTO*(1.+EXCESS)+WFUEL
  109.       WNFWG=WNF/WFLUE 
  110.       WMSAC=WMSA
  111.       NOC(9)=FNO(9) 
  112.       NOC(10)=FNO(10) 
  113.       NOC(11)=FNO(11) 
  114.       NOC(12)=FNO(12) 
  115.       NOC(13)=FNO(13) 
  116.       NO13=FNO(13)
  117.       N13=FN(13)
  118.       NOTC=NOT
  119.       NOTTC=NOT+WMSA
  120.       ASC=WMSA*AAIR 
  121.       BSC=WMSA*BAIR 
  122.       CSC=WMSA*CAIR 
  123.       DSC=WMSA*DAIR 
  124.       DO 40 I=9,12
  125.       ASC=ASC+A(I)*NOC(I) 
  126.       BSC=BSC+B(I)*NOC(I) 
  127.       CSC=CSC+C(I)*NOC(I) 
  128.    40 DSC=DSC+(D(I)-H(I,77.)+HF(I))*NOC(I)
  129.       RETURN
  130. C 
  131. C * * * * * * * * * * * 
  132. C 
  133.       ENTRY FUELA 
  134. C 
  135.       PSAT=PVF(TFLUE) 
  136.       PH2O=NO13*BARPSI/NOTTC
  137.       IF(PSAT.GT.PH2O)THEN
  138. C      NO CONDENSATION IN FLUE GAS
  139.           NH2OL=0.
  140.           NOC(13)=NO13
  141.       ELSE
  142. C       CONDENSATION IN FLUE GAS
  143. C          ASSUME FLUE GAS IS SATURATED THROUGHOUT
  144.           FA=PSAT/BARPSI
  145.           NH2OL=(FA*(NOTC+WMSAC*(1.+N13))-NO13-WMSAC*N13)/(FA-1.) 
  146.           NOC(13)=NO13-NH2OL
  147.       IF(NOC(13).LT.0.)NOC(13)=0. 
  148.           ENDIF 
  149.       ASV=ASC+A(13)*NOC(13) 
  150.       BSV=BSC+B(13)*NOC(13) 
  151.       CSV=CSC+C(13)*NOC(13) 
  152.       DSV=DSC+(D(13)-H(13,77.)+HF(13))*NOC(13)
  153. C  QVSLOS IS SENSIBLE HEAT LOSS AT FLUE (INCLUDES COMP DUE TO WAT VAP)
  154. C  QVLLOS IS LATENT HEAT LOSS OF UNCONDENSED WAT VAP) 
  155.       HCMB=(((ASV*TCOMB+BSV)*TCOMB+CSV)*TCOMB+DSV)*WNFWG
  156.       HTFLU=(((ASV*TFLUE+BSV)*TFLUE+CSV)*TFLUE+DSV)*WNFWG 
  157. C      CPFLU=(HTFLU-HCMB)/(TFLUE-TCOMB+.0000001) 
  158.       QVSLOS=(HTFLU-HCMB)*WFLUE 
  159. C   AS MORE AND MORE WAT VAP IS CONDENSED FROM FLUE GAS, NOC(13) WILL 
  160. C     APPROACH NH2OL
  161.       QVLLOS=NOC(13)*HFG*WNFC 
  162.       IF(QVLLOS.LT.0.)QVLLOS=0. 
  163.       RETURN
  164.       END 
  165.       FUNCTION H(I,T) 
  166.       DIMENSION A(14),B(14),C(14),D(14) 
  167.       DATA A/8*0.,  -.168288E-6,  -.438419E-7,  -.509447E-7,
  168.      +       0.,  -.536138E-7,  0./ 
  169.       DATA B/6.09436E-3,  10.1670E-3,  14.1216E-3,  16.7486E-3, 
  170.      +       17.4328E-3,  21.4615E-3,  21.0380E-3,  40.2128E-3, 
  171.      +       .159751E-2,  .518974E-3,  .555976E-3,  .1488E-3, 
  172.      +       .999914E-3,  -3.018458E-3/ 
  173.       DATA C/6.96554,  10.8508,  15.722, 21.4222, 21.0855,
  174.      +       26.0134,  26.1784,  29.1591,  9.46994,  6.65257, 
  175.      +       7.05349,  6.99778,  7.52025,  -9.69913/
  176.       DATA D/4152.99,  8577.87,  12169.23,  14986.09,  15514.83,
  177.      +       17899.39,  18739.14,  21833.76,  3209.73,  3226.06,
  178.      +       3150.49,  3103.78,  3703.66,  19683.4/ 
  179.       H=((A(I)*T+B(I))*T+C(I))*T+D(I) 
  180.       RETURN
  181.       END 
  182.