SUBROUTINE TYPE75 (TIME,XIN,OUT,T,DTDT,PAR,INFO,ICNTRL,*) C* SUBROUTINE STINJ (Prop,P,MAir,TAirEnt,WAirEnt,MSteam,TSteam, C* & TAirLvg,WAirLvg,QTot,QSen,ErrStat) C*********************************************************************** C* Copyright ASHRAE. Toolkit for HVAC System Energy Calculations C* Revised for TRNSYS by N.Blair and R.Schwarz C*********************************************************************** C* SUBROUTINE: STINJ C* C* LANGUAGE: FORTRAN 77 C* C* PURPOSE: Calculate the performance of a dry steam C* humidifier. Outputs include leaving air C* temperature and humidity ratio. C*********************************************************************** C* INPUT VARIABLES SAMPLE VALUES C* XIN(1) MAiR Entering dry air mass flow rate (kg/s) 1.35 C* XIN(2) TAirEnt Entering air temperature (C) 15.5 C* XIN(3) WAirEnt Entering air humidity ratio (-) 0.005 C* XIN(4) MSteam Steam mass flow rate (kg/s) 0.01 C* XIN(5) TSteam Steam temperature (C) 150.0 C* C* OUTPUT VARIABLES C* OUT (1) TAirLvg Leaving air temperature (C) 17.2487 C* OUT (2) WAirLvg Leaving air humidity ratio (-) 0.0123243 C* OUT (3) QTot Total heat transfer rate (W) -27433.5 C* OUT (4) QSen Sensible heat transfer rate (W) -2396.16 C* OUT (5) ErrStat Error status indicator, 0 = ok, 1 = error (-) 0.0 C* C* PARAMETERS C* EffSat Saturation efficiency (-) 1.0 C* C* PROPERTIES C* CpAir Specific heat of air (J/kg C) C* CpVap Specific heat of water vapor (steam) (J/kg C) C* Patm Atmospheric pressure (Pa) C*********************************************************************** C MAJOR RESTRICTIONS: C C DEVELOPER: Shauna Gabel C Michael J. Brandemuehl, PhD, PE C University of Colorado at Boulder C C DATE: January 1, 1992 C C INCLUDE FILES: stinj.inc C SUBROUTINES CALLED: None C FUNCTIONS CALLED: HUMRATIO C SATPRESS C ENTHALPY C REVISION HISTORY: None C C REFERENCE: Clark, D.R. 1985. HVACSIM+ building C systems and equipment simulation program: C Reference Manual. NBSIR 84-2996, National C Institute of Standards and Technology, C Washingtion, D.C. C*********************************************************************** C INTERNAL VARIABLES: C cpMoist Specific heat of moist air (J/kg C) C wSatLvg Saturation humidity ratio at leaving air temp (-) C************************************************************************ DIMENSION XIN(5),OUT(5),PAR(1),INFO(15) DOUBLE PRECISION XIN OUT INTEGER PAR,INFO,IOPT,NI,NP,ND REAL PATM,CPAIR,CPWAT,CPVAP,CPLIQ,DVISCAIR,DVISCLIQ, &KAIR,KLIQ,RHOLIQ,HFG,RAIR,TKELMULT,TABSADD,PROP(16), &MAIR,SMALL CHARACTER*3 YCHECK(5),OCHECK(5) DATA small/1.E-9/ DATA YCHECK/'MF2','TE1','DM1','MF2','TE1'/ DATA OCHECK/'TE1','DM1','PW2','PW2','DM1'/ PARAMETER (PATM=1) PARAMETER (CPAIR=2) PARAMETER (CPWAT=3) PARAMETER (CPVAP=4) PARAMETER (CPLIQ=5) PARAMETER (DVISCAIR=6) PARAMETER (DVISCLIQ=7) PARAMETER (KAIR=8) PARAMETER (KLIQ=9) PARAMETER (RHOLIQ=10) PARAMETER (HFG=11) PARAMETER (RAIR=12) PARAMETER (TKELMULT=13) PARAMETER (TABSADD=14) PARAMETER (PAMULT=15) PARAMETER (PABSADD=16) PROP (PATM)=101325.0 PROP (CPAIR)=1006.0 PROP (CPWAT)=4186.0 PROP (CPVAP)=1805.0 PROP (CPLIQ)=4186.0 PROP (DVISCAIR)=.0000182 PROP (DVISCLIQ)=.00144 PROP (KAIR)=.026 PROP (KLIQ)=.604 PROP (RHOLIQ)=998.0 PROP (HFG)=2501000.0 PROP (RAIR)=287.055 PROP (TKELMULT)=1.0 PROP (TABSADD)=273.15 PROP (PAMULT)=1.0 PROP (PABSADD)=0.0 ErrStat = 0 IOPT = -1 NI = 5 !CORRECT NUMBER OF INPUTS NP = 1 !CORRECT NUMBER OF PARAMETERS ND = 0 !CORRECT NUMBER OF DERIVATIVES EFFSAT =PAR(1) MAIR =XIN(1) TAIRENT =XIN(2) WAIRENT =XIN(3) MSTEAM =XIN(4) TSTEAM =XIN(5) IF (INFO(7).EQ.-1) THEN CALL TYPECK(IOPT,INFO,NI,NP,ND) C CHECKS TO SEE IF USER'S INFO MATCHES CORRECT NUMBER CALL RCHECK(INFO,YCHECK,OCHECK) C CHECKS TO SEE IF INPUT AND OUTPUT UNITS MATCH) ENDIF C1*** Calculate the inlet air specific heat cpMoist = Prop(CpAir)+WAirEnt*Prop(CpVap) C1*** Calculate the leaving air temperature as a weighted average of inlet C1*** and steam temperatures TAirLvg = (MSteam*Prop(CpVap)*TSteam+MAir*cpMoist*TAirEnt)/ & (MSteam*Prop(CpVap)+MAir*cpMoist) C1*** Assuming all entering steam is absorbed by the air, calculate the C1*** outlet humidity ratio WAirLvg = WAirEnt + MSteam / MAX(MAir,small) C1*** Calculate the saturation humidity ratio at the leaving air C1*** temperature wSatLvg = HUMRATIO (Prop(Patm),SATPRESS(Prop,TAirLvg)) C1*** The leaving humidity ratio cannot exceed the saturation efficiency C1*** times the leaving humidity ratio at saturation WAirLvg = MIN (WAirLvg,wSatLvg*EffSat) C1*** Calculate heat transfer rates QTot = MAir*(ENTHALPY(Prop,TAirEnt,WAirEnt) & -ENTHALPY(Prop,TAirLvg,WAirLvg)) QSen = MAir*cpMoist*(TAirEnt-TAirLvg) OUT(1)=TAIRLVG OUT(2)=WAIRLVG OUT(3)=QTOT OUT(4)=QSEN OUT(5)=ERRSTAT RETURN 1 END REAL FUNCTION HUMRATIO (Patm,Pw) C*********************************************************************** C* Copyright ASHRAE. Toolkit for HVAC System Energy Calculations C*********************************************************************** C* FUNCTION: HUMRATIO C* C* LANGUAGE: FORTRAN 77 C* C* PURPOSE: Calculate the humidity ratio from water C* vapor pressure and atmospheric pressure C*********************************************************************** C* INPUT VARIABLES C* Patm Atmospheric pressure (Pa) C* Pw Partial water vapor pressure (Pa) C* C* OUTPUT VARIABLES C* HumRatio Humidity ratio (-) C*********************************************************************** C MAJOR RESRICTIONS: None C C DEVELOPER: Shauna Gabel C Michael J. Brandemuehl, PhD, PE C University of Colorado at Boulder C C DATE: January 1, 1992 C C INCLUDE FILES: None C SUBROUTINES CALLED: None C FUNCTIONS CALLED: None C C REVISION HISTORY: None C C REFERENCE: 1989 ASHRAE Handbook - Fundamentals C************************************************************************ C1*** Calculate the humidity ratio. HumRatio = 0.62198*Pw/(Patm-Pw) RETURN END REAL FUNCTION SATPRESS (Prop,T) C*********************************************************************** C* Copyright ASHRAE. Toolkit for HVAC System Energy Calculations C*********************************************************************** C* SUBROUTINE: SATPRESS C* C* LANGUAGE: FORTRAN 77 C* C* PURPOSE: Calculate saturation pressure of water C* vapor as a function of temperature C*********************************************************************** C* INPUT VARIABLES C* T Temperature (C) C* C* OUTPUT VARIABLES C* SatPress Saturation pressure (Pa) C* C* PROPERTIES C* TKelMult Multiplying factor to convert user T to Kelvin C* TAbsAdd Additive factor to convert user T to absolute T C* tKel = Prop(TKelMult) * (T + Prop(TAbsAdd)) C* PaMult Multiplying factor to convert user P to Pascals C* PAbsAdd Additive factor to convert user P to absolute P C* Pa = Prop(PaMult) * (P + Prop(PAbsAdd)) C*********************************************************************** C MAJOR RESTRICTIONS: 173.16 K <= Temp <= 473.15 K C C DEVELOPER: Shauna Gabel C Michael J. Brandemuehl, PhD, PE C University of Colorado at Boulder C C DATE: January 1, 1992 C C INCLUDE FILES: prop.inc C SUBROUTINES CALLED: None C FUNCTIONS CALLED: None C C REVISION HISTORY: None C C REFERENCE: 1989 ASHRAE Handbook - Fundamentals C C Hyland, R.W., and A. Wexler. 1983. C Formulations for the thermodynamic C properties of the saturated phases of H2O C from 173.15 K to 473.15 K. ASHRAE C Transactions, Vol. 89, No. 2A, pp. 500-519 C*********************************************************************** C INTERNAL VARIABLES: C tKel Temperature in Kelvin (K) C pascals Saturation pressure (Pa) C*********************************************************************** INTEGER PATM,CPAIR,CPWAT,CPLIQ,CPVAP,DVISCAIR, &DVISCLIQ,KAIR,KLIQ,RHOLIQ,HFG,RAIR, &TEKELMULT,TABSADD,PAMULT,PABSADD REAL PROP(16) PARAMETER (PATM=1) PARAMETER (CPAIR=2) PARAMETER (CPWAT=3) PARAMETER (CPVAP=4) PARAMETER (CPLIQ=5) PARAMETER (DVISCAIR=6) PARAMETER (DVISCLIQ=7) PARAMETER (KAIR=8) PARAMETER (KLIQ=9) PARAMETER (RHOLIQ=10) PARAMETER (HFG=11) PARAMETER (RAIR=12) PARAMETER (TKELMULT=13) PARAMETER (TABSADD=14) PARAMETER (PAMULT=15) PARAMETER (PABSADD=16) DATA C1/-5674.5359/,C2/6.3925247/,C3/-0.9677843E-2/ DATA C4/0.62215701E-6/,C5/0.20747825E-8/,C6/-0.9484024E-12/ DATA C7/4.1635019/,C8/-5800.2206/,C9/1.3914993/,C10/-0.048640239/ DATA C11/0.41764768E-4/,C12/-0.14452093E-7/,C13/6.5459673/ C1*** Convert temperature from user units to Kelvin. tKel = Prop(TKelMult)*(T+Prop(TAbsAdd)) C1*** If below freezing, calculate saturation pressure over ice. IF (tKel .LT. 273.15) THEN pascals = EXP(C1/tKel+C2+C3*tKel+C4*tKel**2+C5*tKel**3+C6* & tKel**4+C7*ALOG(tKel)) C1*** If above freezing, calculate saturation pressure over liquid water. ELSE IF (tKel .GE. 273.15) THEN pascals = EXP(C8/tKel+C9+C10*tKel+C11*tKel**2+C12*tKel**3+C13 & *ALOG(tKel)) ENDIF C1*** Convert pressure from Pascals to user units SatPress = pascals/Prop(PaMult) - Prop(PAbsAdd) RETURN END REAL FUNCTION ENTHALPY (Prop,TDB,W) C*********************************************************************** C* Copyright ASHRAE. Toolkit for HVAC System Energy Calculations C*********************************************************************** C* FUNCTION: ENTHALPY C* C* LANGUAGE: FORTRAN 77 C* C* PURPOSE: Calculate the enthalpy of moist air. C*********************************************************************** C* INPUT VARIABLES: C* TDB Dry bulb temperature (C) C* W Humidity ratio (-) C* C* OUTPUT VARIABLES: C* Enthalpy Enthalpy of moist air (J/kg) C* C* PROPERTIES: C* CpAir Specific heat of air (J/kg C) C* CpVap Specific heat of water vapor (J/kg C) C* Hfg Reference heat of vaporization of water (J/kg) C*********************************************************************** C MAJOR RESTRICTIONS Uses perfect gas relationships C Fit for enthalpy of saturated water vapor C C DEVELOPER: Shauna Gabel C Michael J. Brandemuehl, PhD, PE C University of Colorado at Boulder C C DATE: January 1, 1992 C C INCLUDE FILES: PROP.INC C SUBROUTINES CALLED: None C FUNCTIONS CALLED: None C C REVISION HISTORY: None C C REFERENCE: 1989 ASHRAE Handbook - Fundamentals C*********************************************************************** INTEGER PATM,CPAIR,CPWAT,CPLIQ,CPVAP,DVISCAIR, &DVISCLIQ,KAIR,KLIQ,RHOLIQ,HFG,RAIR, &TEKELMULT,TABSADD,PAMULT,PABSADD REAL PROP(16) PARAMETER (PATM=1) PARAMETER (CPAIR=2) PARAMETER (CPWAT=3) PARAMETER (CPVAP=4) PARAMETER (CPLIQ=5) PARAMETER (DVISCAIR=6) PARAMETER (DVISCLIQ=7) PARAMETER (KAIR=8) PARAMETER (KLIQ=9) PARAMETER (RHOLIQ=10) PARAMETER (HFG=11) PARAMETER (RAIR=12) PARAMETER (TKELMULT=13) PARAMETER (TABSADD=14) PARAMETER (PAMULT=15) PARAMETER (PABSADD=16) C1*** Calculate the enthalpy as a function of dry bulb temperature and C1*** humidity ratio. hDryAir = Prop(CpAir)*TDB hSatVap = Prop(Hfg) + Prop(CpVap)*TDB Enthalpy = hDryAir + W*hSatVap RETURN END