! This component has been assigned Type Number 211. If that number conflicts with ! another user Type number, you will need to change it and recompile the appropriate ! dll. SUBROUTINE TYPE211(TIME,XIN,OUT,T,DTDT,PAR,INFO,ICNTRL,*) TOT00010 !DEC$ATTRIBUTES DLLEXPORT :: TYPE211 C***********************************************************************TOT00020 C STATIC BOILER TOT00030 C ref : "Boiler specification" IEA annex 10 AN10 8TOT00040 C***********************************************************************TOT00050 C M.DACHELET MAY 1987 MODIF MS1088:C3,C4 TOT00060 C***********************************************************************TOT00070 C MC BURNER FUEL MASS FLOW RATE TOT00080 C QU USEFUL POWER TOT00090 C QC CONSUMED POWER TOT00100 C CO2 CO2 CONTENT IN GAZ TOT00110 C EFF EFFICIENCY OF THE BOILER TOT00120 C ETAON NOMINAL EFFICIENCY TOT00130 C TETAON FRACTION OF CYCLE WHEN BURNER IS "on" TOT00140 C TETASB FRACTION OF CYCLE WHEN BURNER IS "on" DURING STAND-BY TOT00150 C EFF EFFICIENCY OF EQUIVALENT HEAT EXCHANGER TOT00160 C AU AU OF EQUIVALENT HEAT EXCHANGER TOT00170 C MW WATER MASS FLOW RATE TOT00180 C MW0 NOMINAL WATER MASS FLOW RATE TOT00190 C MC FUEL MASS FLOW RATE TOT00200 C MC0 NOMINAL FUEL MASS FLOW RATE TOT00210 C TWSU WATER SUPPLY TEMPERATURE TOT00220 C TWEX WATER EXHAUST TEMPERATURE TOT00230 C TAMB AMBIENT TEMPERATURE TOT00240 C TC FUEL TEMPERATURE TOT00250 C PC LOW CALORIFIC VALUE OF THE FUEL TOT00260 C NTU NUMBER OF TRANSFERT UNITS TOT00270 C TOT00280 C FOR MORE INFORMATION ON VARIABLES AND NOMENCLATURE : TOT00290 C "BOILER SPECIFICATIONS" (IEA - ANNEX 10) TOT00300 C***********************************************************************TOT00310 C NOTE : UNITS TOT00320 C TOT00330 C FOR CORRECT RUN , PLEASE USE COHERENT UNIT SYSTEM TOT00340 C TOT00350 C EXAMPLES: TOT00360 C TOT00370 C TIME s s h h TOT00380 C POWER W kW kJ/h kW TOT00390 C MASS kg kg kg kg TOT00400 C VOLUME m3 m3 m3 m3 TOT00410 C ENERGY J kJ kJ kWh TOT00420 C TOT00430 C***********************************************************************TOT00440 C INPUTS : 1 WATER SUPPLY TEMPERATURE TOT00450 C 2 WATER MASS FLOW RATE TOT00460 C 3 AMBIENT TEMPERATURE TOT00470 C 4 MEAN AQUASTAT TEMPERATURE TOT00480 C 5 GAMMA 0=OFF 1=ON TOT00490 C TOT00500 C OUTPUTS : 1 WATER EXHAUST TEMPERATURE TOT00510 C 2 WATER MASS FLOW RATE TOT00520 C 3 USEFUL POWER TOT00530 C 4 CONSUMED POWER TOT00540 C 5 MEAN FUEL CONSUMPTION TOT00550 C 6 TETAON TOT00560 C 7 TETASB TOT00570 C 8 EFFICIENCY TOT00580 C 9 EFFICIENCY 100% ON TOT00590 C 10 EFFICACITY OF THE EQ. HEAT EXCHANGER TOT00600 C 11 EQUIVALENT AU TOT00610 C TOT00620 C PARAM. : 1 BURNER FUEL MASS FLOW RATE TOT00630 C 2 CONCENTRATION OF CO2 IN GAS TOT00640 C 3,4,5 AU0, K1, K2 TOT00650 C 6,7,8 Yw, DYw, Kw TOT00660 C 9,10,11 MC0, MW0, CO20 TOT00670 C 12,13 C1,C2 TOT00680 C 14 CPW TOT00700 C 15,16 CPC,PC TOT00710 C TOT00720 C***********************************************************************TOT00730 C TOT00740 C PROPOSED VALUES : C1 3.294 kJ/K kg TOT00750 C C2 2.105 kJ/K kg TOT00760 C CPW 4.187 kJ/K kg TOT00790 C CPC 1.880 kJ/K kg TOT00800 C PC 42875.0 kJ/kg TOT00810 C TOT00820 C***********************************************************************TOT00830 DIMENSION XIN(5),PAR(16),OUT(11),INFO(15),C(4) TOT00850 REAL NTU,MW,MC,MW0,MC0,K1,K2,KB,KD,KE,KW TOT00860 INTEGER GAMMA TOT00870 DOUBLE PRECISION XIN,OUT TOT00880 ! Set the version information for TRNSYS IF (INFO(7).EQ.-2) THEN INFO(12) = 15 RETURN 1 ENDIF C TOT00890 C FIRST CALL TOT00900 C TOT00910 TOT00920 IF(INFO(7).EQ.-1) THEN TOT00930 INFO(6)=11 TOT00940 CALL TYPECK(1,INFO,5,16,0) TOT00950 ENDIF TOT00960 TOT00970 C TOT00980 C INITIALISATION OF INPUTS AND PARAMETERS TOT00990 C TOT01000 TOT01010 MC=PAR(1) TOT01020 CO2=PAR(2) TOT01030 AU0=PAR(3) TOT01040 K1=PAR(4) TOT01050 K2=PAR(5) TOT01060 YW=PAR(6) TOT01070 DYW=PAR(7) TOT01080 KW=PAR(8) TOT01090 MC0=PAR(9) TOT01100 MW0=PAR(10) TOT01110 CO20=PAR(11) TOT01120 C1=PAR(12) TOT01130 C2=PAR(13) TOT01140 CPW=PAR(14) TOT01170 CPC=PAR(15) TOT01180 PC =PAR(16) TOT01190 TOT01200 TWSU=XIN(1) TOT01210 MW=XIN(2) TOT01220 TAMB=XIN(3) TOT01230 TAQUA=XIN(4) TOT01240 GAMMA=INT(XIN(5)+0.1) TOT01250 TC=TAMB TOT01260 IF(MW.LT.0.) CALL TYPECK(3,INFO,0,0,0) IF(MW.LT.1.E-5) MW = 1.E-5 TOT01270 C TOT01280 C CHECK IF BOILER IS OFF TOT01290 C TOT01300 TOT01310 IF((GAMMA.EQ.0). OR .(TWSU.GT.TAQUA)) THEN TOT01320 CW=MW*CPW TOT01330 TWEX=(TWSU-(TWSU*YW*(1.-KW)/CW)+(YW*TAMB/CW))/(1.+(YW*KW/CW)) TOT01340 QU=CW*(TWEX-TWSU) TOT01350 QC=0. TOT01360 TETAON=0. TOT01370 TETASB=0. TOT01380 ETA=0. TOT01390 ETAON=0. TOT01400 EFF=0. TOT01410 AU=0. TOT01420 GOTO 20 TOT01430 ENDIF TOT01440 TOT01450 C TOT01460 C GENERAL CALCULATIONS TOT01470 C TOT01480 QB=MC*PC TOT01490 CPG=C1+C2/CO2 TOT01500 CPG0=C1+C2/CO20 TOT01510 CPA=CPG-CPC TOT01520 CG=MC*CPG TOT01530 CG0=MC0*CPG0 TOT01540 CW=MW*CPW TOT01550 CW0=MW0*CPW TOT01560 AU=AU0/(1.+K1*(CG0-CG)/CG+K2*(CW0-CW)/CW) TOT01570 NTU=AU/CG TOT01580 OMEGA=CG/CW TOT01590 UM0=1.-OMEGA TOT01600 EFF=(1.-EXP(-NTU*UM0))/(1.-OMEGA*EXP(-NTU*UM0)) TOT01610 KB=(CPG*TWSU-(CPA*TAMB+CPC*TC))/PC TOT01620 TWEX=TAQUA TOT01630 QU=CW*(TWEX-TWSU) TOT01640 TOT01650 TAV=KW*TWEX+(1.-KW)*TWSU TOT01660 KD=(DYW*(TAV-TAMB))/QB TOT01670 KE=(YW*(TAV-TAMB))/QB TOT01680 ETAON=EFF*(1.-KB)-KD-KE TOT01690 QUMAX=QB*ETAON TOT01700 TETAON=((QU/QB)+KE)/(EFF*(1.-KB)-KD) TOT01710 TETASB=KE/(EFF*(1.-KB)-KD) TOT01720 TOT01730 IF(QU.LE.QUMAX)GO TO 30 TOT01740 TOT01750 C TOT01760 C CALCULATION IN CASE OF FULL LOAD TOT01770 C TOT01780 TOT01790 CC=CW+((YW+DYW)*KW) TOT01800 CD=((CW*TWSU)+(QB*EFF*(1.-KB))-((YW+DYW)*(((1.-KW)*TWSU)-TAMB))) TOT01810 TWEX=CD/CC TOT01820 QU=CW*(TWEX-TWSU) TOT01830 TETAON=1. TOT01840 ETAON=QU/QB TOT01850 TOT01860 C TOT01870 30 CONTINUE TOT01880 C TOT01890 QC=QB*TETAON TOT01900 ETA=QU/QC TOT01910 TOT01920 C TOT01930 C OUTPUTS AND END TOT01940 C TOT01950 TOT01960 20 OUT(1)=TWEX TOT01970 OUT(2)=MW TOT01980 OUT(3)=QU TOT01990 OUT(4)=QC TOT02000 OUT(5)=QC/PC TOT02010 OUT(6)=TETAON TOT02020 OUT(7)=TETASB TOT02030 OUT(8)=ETA TOT02040 OUT(9)=ETAON TOT02050 OUT(10)=EFF TOT02060 OUT(11)=AU TOT02070 C TOT02080 RETURN 1 TOT02090 END TOT02100