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 / HUMDF.FOR < prev    next >
Text File  |  1992-04-07  |  4KB  |  127 lines
  1.       SUBROUTINE HUMDF(TI)
  2. C 
  3. C  HUMIDIFIER MODEL 
  4. C  REF. BULLOCK,C.E.,RESIDENTIAL CENTRAL HEATING SYSTEM HUMIDIFIER
  5. C       COMPUTER SIMULATION AND APPLICATION GUIDELINES, ASHRAE PAPER
  6. C       NO. AT-78-8, NO.3, 1978
  7. C  NOT USED BY BCL....
  8. CMDK CPAIR
  9. CMDK NZN
  10. CMDK HUMIDC
  11. CMDK WIBLK
  12. C 
  13.       REAL LE,MWOUT,MWIN
  14.       DATA AREA,THK/144.,2.25/
  15. C ***   HUMIDIFIER READ - IN
  16. C ***    MWI - SUPPLY WATER FLOW RATE, LB/HR
  17. C ***    TWI - SUPPLY WATER TEMP,DEG F
  18. C ***    CFMBP - BYPASS AIR FLOW RATE, CFM
  19. C ***    HSTC - HUMIDISTAT TIME CONSTANT, MIN 
  20. C ***    SETHS - HUMIDISTAT SET POINT, %
  21. C ***    DIFFHS - HUMIDISTAT DIFFERENTIAL , % 
  22.       READ(5,*) MWI,TWI,CFMBP,HSTC,SETHS,DIFFHS 
  23.       IF ( CFMBP .LE. 0.) GO TO 147 
  24.       WRITE(60,2250) MWI,TWI,CFMBP,HSTC,SETHS,DIFFHS
  25.       LE = 0.9
  26.       VOL = AREA*THK/1728.
  27.       VEL = CFMBP/AREA*144. 
  28.       HDAV = 8.217*VEL+2135.
  29.       HSTC = HSTC/60. 
  30.       WA1=4.5*CFMBP 
  31.       F1=HDAV*VOL/WA1 
  32.       F1LE= F1*LE 
  33.       VOLWA1=VOL/WA1
  34.       HDAVLE= HDAV*LE 
  35.       DMDTW=.1
  36.       F1B= F1/(1.+F1/2.)
  37.       F1C= (1.-F1/2.)/(1.+F1/2.)
  38.       F2= F1LE*CPAIR
  39.       CPA2=1./(CPAIR+F2/2.) 
  40.       HDAVLC= HDAVLE*CPAIR
  41.       QSPP= (TWI-32.)*MWI 
  42.   147 CONTINUE
  43.            RETURN
  44. C 
  45.       ENTRY HUMDF1(TI)
  46. C ***   SET HUMIDIFIER STATES 
  47.       TWOUT = TWI 
  48.       QSP=0.
  49.       QWA=0.
  50.       QHUM=0. 
  51.       WBULK=WI(1) 
  52.       HIN=CPAIR*TI+1061.*WI(1)
  53.       TBULK=TI
  54.       RETURN
  55. C 
  56.       ENTRY HUMDF3(TI)
  57. C 
  58. C ***      $$$$$$$$$  HUMIDIFIER CALCULATIONS  $$$$$$$$$
  59.       IF (CFMBP .LE.0.0) GO TO 630
  60.       IF(IHUM.EQ.1) GO TO 580 
  61. C ***   IF HUMIDIFIER IS OFF, INLET=OUTLET CONDITIONS 
  62.       QSP=0.
  63.       QHUM=0. 
  64.       QWA=0.
  65.       TWOUT = TWI 
  66.       WBULK=WAC 
  67.       TBULK=TB
  68.       GO TO 630 
  69. C ***   START OF HUMIDIFER FINE TIME STEPS
  70. 580   LTW=0 
  71.       HIN=CPAIR*TB+1061.*WAC
  72.       ITW=1 
  73. C ***   START OF HUMIDIFIER LUMP CALCULATIONS 
  74. 582   TM = TWOUT
  75.       IF(TM.LT.32.) TM = 32.
  76.       F = (((1.14204D-9*TM-5.73544D-8)*TM-8.2068D-5)*TM+4.57702D-2)*TM
  77.      +    -6.95562
  78.       WSM = EXP(F)
  79.       WBULK = F1B*WSM+F1C*WAC 
  80.       WBAR = (WAC+WBULK)*.5 
  81.       HGM= 1061.+.45*TM 
  82.       R1 = HDAV*(WSM-WBAR)
  83.       F3 = R1*HGM*VOLWA1
  84.       TBULK= (F2*(TM-TB*.5)+HIN-1061.*WBULK+F3)*CPA2
  85.       IF(TBULK.LT.TWOUT) TBULK=TWOUT
  86.       TBAR = (TB+TBULK)*.5
  87.       R2 = HDAVLC*(TM-TBAR)+R1*HGM
  88.       MWOUT = MWI - R1*VOL
  89.       AH = R2*VOL 
  90.       BH = MWI*(TWOUT-TWI)
  91.       CH = R1*VOL*(TWOUT-32.) 
  92.       DF3DTM= (R1*VOL*.45+HGM*DMDTW)/WA1
  93.       DWBULK=F1B*WSM*(((4.56816D-9*TM-17.20632D-8)*TM-16.4136D-5)*TM+ 
  94.      * 4.57702D-2)
  95. C **** MWIN NOT SET IN ORIGINAL PROGRAM 
  96.       MWIN=0. 
  97.       DTBULK= (F2-1061.*DWBULK+DF3DTM)*CPA2 
  98.       ABC= (AH+BH-CH)/((HDAVLC*(1.-DTBULK*.5)-.55*R1)*VOL+MWIN+ 
  99.      * (1093.-.55*TM)*DMDTW)
  100.       IF (LINES .LT. ITER) WRITE(60,2330) ABC,TM,DMDTW,TB 
  101.       IF(ABS(ABC).LT.0.01)GO TO 620 
  102.       IF (ITW .EQ. 2) GO TO 584 
  103.       ITW=2 
  104.       ABC1=ABC
  105.       DM=MWOUT
  106.       TWOUT=TWOUT-ABC 
  107.       GO TO 582 
  108. 584   TWOUT=TWOUT+ABC1*ABC/(ABC-ABC1) 
  109.       DMDTW= (MWOUT-DM)/ABC1
  110.       LTW=LTW+1 
  111.       ITW=1 
  112.       IF( LTW .LT. 10) GO TO 582
  113.       WRITE(60,3030)
  114.       STOP
  115. 620   QHUM = MWOUT*(TWOUT-TWI)+(MWI-MWOUT)*(1093.+.45*TBULK-TWI)
  116.       QSP =QSPP 
  117.       QWA = MWOUT*(TWOUT-32.) 
  118. 630   CONTINUE
  119.       RETURN
  120. 2250  FORMAT (/' HUMIDISTAT PARAMETERS:'/ '   MWI       TWI     CFMBP 
  121.      $    HSTC     SETHS     DIFFHS   '/ F8.2,F10.2,F10.1,F10.2,
  122.      $ 2F10.4)
  123.  2330 FORMAT('  ABC=',F8.3,'  TWOUT=',F8.2,'  DMDTW=',F8.3, 
  124.      +'  TB=',F8.2) 
  125.  3030 FORMAT('  NO CONVERGENCE IN HUMIDIFIER WATER TEMP LVG') 
  126.       END 
  127.