• No results found

Berekening van de warmteproduktie en de soortelijke warmte van tuinbouwprodukten uit meetresultaten van de adiabatische calorimeter

N/A
N/A
Protected

Academic year: 2021

Share "Berekening van de warmteproduktie en de soortelijke warmte van tuinbouwprodukten uit meetresultaten van de adiabatische calorimeter"

Copied!
26
0
0

Bezig met laden.... (Bekijk nu de volledige tekst)

Hele tekst

(1)

S P R E N G E R I N S T I T U U T Haagsteeg 6, 6700 AA Wageningen

T e l . : 08370-19013

(Publikatie uitsluitend met

toestemming van de directeur)

RAPPORT NO. 2034

Ing. H. van der Krieke

BEREKENING VAN DE WARMTEPRODUKTIE EN DE

SOORTELIJKE WARMTE VAN TUINBOUWPRODUKTEN

UIT MEETRESULTATEN VAN DE ADIABATISCHE

CALORIMETER

Uitgebracht aan de directeur van het Sprenger Instituut.

Proj. no. 49

(2)

3. Invoer bij het computerprogramma 2

4. Berekening van de benodigde grootheden 5

4.1 Temperatuur - tijdkromme 5

4.2 Temperatuur - tijdafgeleide 5

4.3 Evenredigheidsconstante a 5

4.4 Massa bij de betreffende meettemperaturen 6

4.5 Warmteproduktie en soortelijke warmte 6

5. Uitvoer 7

6. Gebruik van ADIACAL in samenhang met

piotprogramma 7

7. Conclusies 8

8. Literatuur 8

9. Bijlagen 9

(3)

1

.

Samenvatting

Door gebruik te maken van het rekenprogramma "ADIACAL" is het mogelijk de warmteproduktie en de soortelijke warmte van tuinbouwprodukten te berekenen uit metingen met de adiabatische calorimeter en wel in het bijzonder met een correctie voor massaverlies en verlies van voelbare warmte als gevolg van het verdampen van vocht.

De warmteproduktie, niet gecorrigeerd voor vochtverlies, is eenvoudig en direct uit de meetgegevens af te leiden. Men wil echter over de gecorrigeerde warmteproduktie beschikken omdat alleen dit gegeven kan worden gebruikt voor het vaststellen van het te installeren koel vermogen bij gekoelde bewaarplaatsen. Er kunnen belangrijke verschillen optreden tussen de ni et-gecorrigeerde en de gecorrigeerde grootheden (zie bijlagen). Het aantal meetgegevens dat per meting ingevoerd moet worden is bij

deze versie van het rekenprogramma aanmerkelijk geringer dan bii de versie vermeld in lit. 1. Over het temperatuurtraject 5°C t/m 40°C behoeven slechts 56 temperaturen te worden geponst. Bij het programma met het

oude rekensysteem werd tijdens het doorlopen van het temperatuurtraject ieder uur de temperatuur afgelezen en dienden alle temperaturen te

worden geponst. Het huidige programma levert dus een belangrijke tijdsbe-sparing op.

(4)

?.. Inleiding

Voor de wijze waarop de warmteproduktie wordt gemeten

m e

t behulp van

de adiabatische calorimeter wordt verwezen naar lit. 1 en 2.

Het rekenprogramma "ADIACAL" berekent de warmteproduktie en eventueel

de soortelijke warmte bij de temperaturen 5, 10, 15, 20, 25, 30, 35, 40°C

en schept de mogelijkheid een correctie aan te brengen op de

warmte-produktie in verband met massaverlies en verlies van voelbare warmte door

verdamping. De verdampingswarmte van water is hoog en bedraagt ca. 2500 kJ/kq.

De voor de verdamping benodigde warmte moet geleverd worden uit de door

het produkt geproduceerde warmte. Daarom heeft verdamping van een kleine

hoeveelheid vocht een grote invloed op het meetresultaat. Tijdens de

op-slag van produkt, waarbij warmteproduktie gepaard naat met verdampinq,

resulteert een effectieve warmteproduktie, die sterk kan verschillen

van de "absolute", slechts door de temperatuur en niet door vochtverlies

beheerste warmteproduktie (zie lit. 1, 2 ) .

Het rekenprogramma verricht een aantal op elkaar volgende handelingen:

-invoer van de gegevens,

-berekening van de temperatuur-tijdkromme,

-berekening van de temperatuurstijging per tijdseenheid bij de genoemde

meettemperaturen,

-berekening van de dampdrukdeficiten nodig voor de verdeling van het

massaverlies over het temperatuurtraject,

-berekening van het over de meetperiode gesommeerde dampdrukdeficit en

van de evenredigheidsconstante a, waarmee het momentane verlies aan

vocht uit het op dat tijdstip heersende dampdrukdeficit kan worden afgeleid,

-berekening van de massa bij de verschillende temperaturen,

-berekening van de warmteproduktie en de soortelijke warmte uit het

voorgaande,

-uitvoer.

S. Invoer

De gegevens nodig voor de berekening, die achtereenvolgens gevraagd worden,

zijn:

-tijdsinterval in uren, dat wil zeggen de tijd tussen twee opeenvolgende

temperatuurwaarnemingen.

(5)

-aantal uren na laatste waarde, dat wil zeggen de tijd in uren, die ligt

tussen de laatst ingevoerde temperatuurwaarneming en het einde van de meting.

-het temperatuurverschil in °C, dit is het constante temperatuurverschil dat wordt onderhouden tussen het door het meetvat geleide en met waterdamp verzadigde gasmengsel en het produkt. Bijvoorbeeld een temperatuurverschil van 2°C bij een relatieve vochtigheid van + 90%. -massa van vat 1 begin, eind (kg). Het begin- en eindgewicht van het

pro-dukt in meetvat 1 wordt gevraagd in kg.

-de begintemperatuur in °C, hiermee wordt de produkttemperatuur bij de start van de meting bedoeld.

-laatste temperatuur, waarbij opgave van de warmteproduktie wordt gewenst. Dit is het laatste temperatuurniveau in °C (een veelvoud van 5°C), waarbij de warmteproduktie berekend moet worden.

-produkt, hier dient de produktnaam te worden gegeven. -pointer C (is de soortelijke warmte bekend of niet bekend)

0 - wel bekend 1 - niet bekend

Is de soortelijke warmte bekend (0), dan wordt deze opgevraagd (in J/kgK). Daarna start het inlezen van de ponsband met daarop de temperatuur en tijdgegevens.

De gegevens zijn als volgt op de ponsband vastgelegd:

Rond het meetpunt (dit is van de reeks waargenomen temperaturen de temperatuur die zich het dichtst bij één der gewenste temperaturen 5, 10, 15 40°C bevindt) worden drie lagere en drie hogere temperaturen geponst. De

interval periode tussen twee geponste temperaturen is via de invoer met het toetsenbord gegeven. Na de rij van 7 meetpunten wordt 99999 geponst en de rij wordt voorafgegaan door de tijd in uren, die verstreken is

tussen de laatst waargenomen temperatuur bij het voorgaande meetpunt en de eerst waargenomen temperatuur van het actuele meetpunt.

Voorbeeld: 99999

60 aantal uren van begintemperatuur tot 4,7°C 4,7 \

4,8

V

4,9

5 / serie van 7 meetpunten 5,1

5,2 5,3

(6)

99999 80 9,4 9,6 9,8 10 10,2 10,4 10,6

aantal uren van 5,3°C tot 9,4°C

serie van 7 meetpunten

99999

enz.

De ponsband wordt afgesloten met X-off.

Ieder getal (temperatuur) op de ponsband dient te worden afgesloten met een CR/LF.

Indien de soortelijke warmte berekend moet worden uit de meting (1 via toetsenbord), is nog een aantal extra gegevens nodig, nl.: -massa van vat 2, begin, eind: het begin- en eindgewicht van het produkt

in meetvat 2 in kg.

-stookvermogen in Watt, het via de verwarmingsdraad toegevoerde stook-vermogen in meetvat 2.

Nu wordt na de eerste ponsband een tweede ponsband of tweede deel van de ponsband ingelezen, waarop de gegevens, die nodig zijn voor de berekening van soortelijke warmte zijn vastgelegd. Dit is een reeks temperaturen, waargenomen bij het tweede meetvat met het reeds bekende tijdsinterval.

Deze reeks beslaat het gehele temperatuurtraject dat door de calorimeter wordt doorlopen. Het aantal temperaturen, dat geponst moet worden in dit geval is gering, omdat in het gestookte meetvat de eindtemperatuur veel eerder wordt bereikt.

(7)

4. Berekening van de benodigde grootheden

liL?§!2§!S§OlQ9_vaG_de_temgeratuur-tiidkrgiT)nie

De temperatuur-tijdkromme van de calorimeter (T in °C en t in uren), die vaak een S-vorm vertoont, wordt benaderd door een vijfdegraads polynoom: T = a + bt + et2 + dt3 + et4 + ft5. Door middel van matrixberekening worden de coëfficiënten a, b, c, d, e en f bepaald.

^i?_Berekemng_yan_de_temgeratuur:tijdafge^eide —rp- in - r — De afgeleide van de temperatuur naar de tijd bij een temperatuur T4 wordt uit zeven temperatuurgegevens als volgt berekend:

(~af)

T

= 60 At (" T l + 9 T2 " 4 5 T3 + 4 5 T5 - 9 T6 + T? )

4

Hierin is T^ een temperatuur van de reeks 5°C, 10°C ...40°C,

At een tijdsinterval tussen twee temperaturen (zie invoer). Het programma berekent de warmteproduktie met afgeleiden van de gemeten temperatuur-tijdcurve, dus temperaturen die verkregen zijn uit ingevoerde meetpunten. Berekening van de afgeleiden uit de berekende temperatuur-tijd-kromme (zie 4.1) blijkt vooral bij de lage temperaturen grote afwijkingen te geven. In feite is de benadering met een vijfdegraads polynoom in dit gebied te goed en slingert de kromme zich enigszins rond de ingevoerde meetpunten.

li3_Berekemng_van_de_eyenredighei^dsc

~ "" ~ ~— "" ——— — — — — — — — — — — — — — — —— — — — — — — — — — — ______.._.. _ ..__,.—

De afgeleide van de massa naar de tijd -gr- wordt evenredig verondersteld met het dampdrukdeficit Ap.

De nodige verzadigingsdampdrukken worden berekend met de Magnus-formule: pverz = 1 02'7 8 5 7 + 2 3 7 ^ I T (T in °C en P in Pa)

Het dampdrukdeficit AP is: Ap = Pyerzj " pverzxi

Hierin is Tl de temperatuur waarbij de in het meetvat gevoerde lucht is verzadigd met waterdamp. Deze lucht wordt daarna opgewarmd tot de temperatuur van het meetvat T. Tl volg T als functie van de tijd met een

voor iedere meting vast gekozen temperatuurverschil, dat afhangt van de gewenste relatieve vochtigheid.

(8)

kan nu berekend worden uit: a

Am

J Ap dt

waarin: Am = totale massaverlies, dat gedurende de meting is opgetreden.

T = tijdsduur van de meting

Ap is een functie van de tijd, verkregen door de Magnus-formule toe

te passen op de temperaturen T en Tl. Beide functies van de tijd.

De functie voor T is berekend volgens paranraaf 4.1.

4i!_§Êr§!S§QiQ9_v§Q_de_massa_bij_d^-^

De massa b i j de meettemperatuur T wordt berekend met

"T * \tart • - f

i M t

waarin t = tijdstip waarop de meettemperatuur T is bereikt.

^§_Wari!)teBrodukti

p

e_en_soortelijke_wa

De ongecorrigeerde formules

warmteproduktie bij temperatuur T:

W

T

=

1

V [ C T - ^ |

T

en soortelijke warmte bij temperatuur T:

P

C

T

=

l

m

2tat " - a t ; J T

(mW/kg)

(J/kgK)

rat ~ar

gaan met gebruik van de voorgaande afleidingen over in:

WT

=

C

T

^ [ ' T

"öT

dT

1.

-%-

°

I A P I

] T

1

"3SÜÜ"

W>2 P + Qv (a\apiTl\g-a^ap2ml )

"h m

2

j - g ^

(-air" - a r y

(mW/kg)

(J/kgK)

Hierin is:

Wj = warmteproduktie bij temperatuur T

C-j- = soortelijke warmte bij temperatuur T

T = Temperatuur T

x

= temperatuur in vat 1

T

2

= temperatuur in vat 2

t = tijd

QvT= verdampingswarmte van water bij temperatuur T

m = massa index 1 of 2 betekent vat 1 of vat 2

* = gecorrigeerd voor vochtverlies op het

moment dat temperatuur T is bereikt

a = evenredigheidsconstante (index 1 of 2; vat 1 of 2)

mW/kg

J/kg K

°C

oc

h

J/kg

k

9water

/ P a h

(9)

-7-5. Uitvoer

De uitvoer van het programma bestaat uit:

a. massa en afgeleide van de temperatuur-tijdcurve vs. temperaturen (tabel)

b. warmteproduktie en soortelijke warmte vs. temperaturen (tabel)

c. evenredigheidsconstante a

d. grafieken van temperatuur vs. tijd en warmteproduktie vs. temperatuur.

De grafieken worden gemaakt met het programma "PLOTMUS" (cassette 29, blok 2 ) .

6. Gebruik van ADIACAL in samenhang met het piot-programma

Het gebruik van het programma "ADIACAL" is als volgt:

1. het inlezen van het programma,

2. het plaatsen van een lege datacassette voor opvang van de output

voor zover nodig bij het programma PLOTMUS,

3. inzetten van de ponsband,

4. runnen van ADIACAL (voor de input zie pag. 2)

Voor meerdere produkten evenzoveel cassettes gebruiken

5. clear het programma ADIACAL,

6. lees het programma PLOTMUS is,

7. runnen van PLOTMUS

Nu wordt gevraagd welke variabelen men als x en y variabelen wil gebruiken.

Voor elk produkt moet er nu worden ingevoerd:

(grafiek no.) (x-var.) (y-var.) (serie no.)

warmteproduktie

(gecorrigeerd)-temperatuur

warmteproduktie (niet

gecor-rigeerd) - temperatuur

temperatuur als functie

van tijd

Als bijlagen zijn een listing van het programma en de uitvoer voor enige

produkten opgenomen.

1 2 3 1 1 4 2 3 1 1 1 2

(10)

7. Conolusies

-Dit programma berekent de warmteproduktie en soortelijke warmte op

twee manieren, nl. niet gecorrigeerd voor massaverlies (= vochtverlies)

tijdens de meting, maar ook met vochtcorrectie. Dat deze gecorrigeerde

warmteproduktie aanzienlijke verschillen met de niet gecorrigeerde

warmteproduktie kan hebben, blijkt uit de voorbeelden in de bijlage 2.

-Het feit dat voor de gecorrigeerde warmteproduktie niet alle, maar

ongeveer 50 temperaturen geponst moeten worden, betekent een grote

tijdsbesparing in vergelijking met het rekenprogramma, dat opgesteld

werd voor het vorige computersysteem.

8. Literatuur

1. Rudolphij J.W., W. Verbeek, F.H. Fockens.

Measuring heat production of respiring produce under normal and

CA-storage conditions withan adiabatic calorimeter.

Lebensmittel-Wissenschaft und Technolog]ie, 10, 153-158 (1977)

2. Verbeek W., J.W. Rudolphij.

De bepaling van de warmteproduktie van tuinbouwprodukten met behulp

van een adiabatische calorimeter.

Koeltechniek 70 (1977) nr. 11 (november)

Wageningen, 12-10-1978

HvdK/LB

(11)

L i s t i n g van het programma ADIACAL b i j l a o e l 20 REM (ADIACAL)

30 DIM A(70),D(40)

y

E<3),M<40),Ml<9),C<8),P<40),P 1<8),X<9),P2

A (40),XI(56,6)

f

Yl(56),G< 6,56>,I < ó,6),Bl< 6),S < 6),B < 40)

40 INPUT "TIJD IN UREN ",H

GO INPUT "AANTAL UREN NA LAATSTE WAARDE",Hl

f 60 INPUT "TEMP.VERSCHIL IN 'C",T1

70 INPUT "MASSAVAT1 BEGIN,EIND",M<1),M

ÖOINPUT "BE GINTE MPER A TUU R",TG

£ 90INPUT "EINDTEMPERATUUR WAARBY WARMTEPRODUKTIE GEWENST",T2

.1.00 INPUT "PRODUKT",U$

110 REM C ONBEKEND 'Z = l

f 120 REM C BEKEND Z=0

160 INPUT "BESLISSINGSVARIABELE C",Z

161 IF Z = l THEN 170

A 162 INPUT "SOORTELIJKE WARMTE",Cl

170 MAT A=ZER

e

ISODATA LOAD /618„A<)

190 IF Z = 0 THEN 240

200 INPUT "MASSA VAT 2 BEGIN,EIND",Ml<1),Ml

210 INPUT "STOOKVERMOGEN IN WATT ",P

A 220 MAT B=ZER

230 DATA LOAD /61S,B( )

240 FOR 1=1 TO 70

£ 250 IF A(I)=0 THEN 260:NEXT I

260 N=I-1

270 1=70:NEXT I

m

230 FOR 1=1 TO 255

290 IF B<I)=0 THEN 300:NEXT I

300 G=I-1

Q 310 1=255:NEXT I

320 X=0

330 FOR 1 = 1 TO N •-- • ••

Q 340 IF A(I)<>99999 THEN 360

350 X=X+7+A(I+l>

360 NEXT I

Q 370 Y=X+H1

330 X<1)=0

390 R=0

0 400 FOR 1=1 TO N

410 I F A( 1 ) 0 9 9 9 9 ? THEN 480 420 R=R + 1 0 4301F R O I THEN 460 440 A2=3

450 GOTO 470

0 460 A2=6

470 X(R+1)=X(R)+A(I+1)+A2*H

480 NEXT I

£ 490 -R=0

500 A1=0

510 11=1

o

I"'-)

0

530 B=T5t2

540 FOR 1=1 TO N

550 IF A(I)=99999 THEN 570

560 GOTO 600

570 1=1)1

580 A1--A1+A(I)

590 GOTO 760

600 11=11+1

610 XI(1,1)=1

620 X1<1,2)=1

630 X1<1,3)=1

640 Xl(l,4)=l

641 X1<1,5)=1

• ' « * — ^ • "

(12)

I 6 8 0 X l < I lf3 ) = A l t 2 6 9 0 X 1 < I 1 , 4 ) = A 1 * 3 7 0 0 X 1 ( I 1 , 5 ) = A 1 * 4 j 7 1 0 X H I l , 6 > = A l + 5 7 2 0 Y i < I l ) * A ( î ) 7 3 0 S = S + Y K I 1 ) 7 4 0 B = B + Y H I l ) t 2

7S0 A1=A1+H

760 NEXT I

770

MAT G= TRN(X1>

780MAT I=G*X1

790 MAT I=INM<I)

300MAT S=G*Y1

810MAT B1=I*S

820 N2=I1

830 11=0

840 FOR 11=1 TO N2

8 5 0 R i = R l + < Y l < I l > ~ B l < l ) - - ß l < 2 > * X l < I l , 2 > ~ B l < 3 > * X i a i , 3 > ~ B l < 4 > * X l < I Q 1 , 4 > - - B l ( 5 ) * X l < I l , 5 > - B l < 6 > * X l < I l , 6 > > * 2 8 6 0 NEXT I I

870 B=B-St2/Il

Q 880 K1=SQR<<B-R1)/B>

890 SELECT PRINT 215(112)

900PRINT HEX(OE);"ADIACAL"

Q 910 PRINT HEX(OE);"PRODUKT=";U*

920 PRINT

930 PRINT "CORR. C0EFF.=";K1

940 P1=0

950 FOR T3=l TO Y STEP H

96Q Y3=B1 < 1 ) +B1 <2)«T3+B14-3)*T3f2+B1 <4) *T3t3+Bl < 5)*T3t4+Bl < 6) *T3t

5

980 Pi=Pl + U O t < 0 . 6609+7. S*Y3/<237.3+Y3)>-10t<0. 6609+7. S*.<Y3

-Tl>/<237.3+Y3-Tl))>«H

990 NEXT T3

1000 K=0

1010 R=0

1020 T=5

1030 FOR 1=1 TO N

1040 IF A(I+1>=99999 THEN 1200

1050 IF ABS(A(I)--T)>ABS(A(I + 1)-T) THEN 1210

1060 K=K+1

1070 P(K+1)=0

1080 D<K) = <A(I~2)-8*A(I-l)+3*A<I + l)--A(I+2))/(12«H>

1090 R=R+1

1100 P2(K)s<i0t(0.6609+7.5*T/(237.3+T))-10t<0.6609+7.5*<T--T1 >/<2

37 .3+T--Tl)))*H

1 1 1 0 FOR T 3 = l TO X(R+1> STEP H 1 1 2 0 Y 4 = B l < i > + B l < 2 > * T 3 + B i < 3 > * T 3 t 2 + B i < 4 > * T 3 t 3 + B i < S > * T 3 t 4 + B K 6 ) * T 3 * 5

1130 P(K+1)=P(K+i)+<10t(0.6609+7.S*Y4/(237.3+Y4)>-10t(0.6609+7

.5*(Y4-T1)/< 237.3+Y4-T1)))«H

1140 NEXT T3

1150 F=<M<1)-M>/<P1>

1160 M<K+1>=M<K)~F*(P<K+1)-P<K)>

1170 T=T+S

1180 IF T O T 2 THEN 1210: I=N: NEXT I : GOTO 1220

1190 GOTO 1210

1200 1=1+2:GOTO 1210

1210 NEXT I

1220 IF Z--=0 THEN 1410

1230 L=l

1240 P(1)=0

H

o

! -n

H

(13)

#

j

. \ V X )

1250 T=ü' '" " "-• -' • -••••• g 1260 P2 = 0 1270 FOR 1=1 TO G 1280 P2=P2+(10t(0.6609+7.5*B(I)/(237.3+B(I) >>-10t<0.6609+7.5*<B< g I )-Tl)/(237.3 + B(I)-Tl)))*H 1290 NEXT I 1300 FOR 1=1 TO G £ 1310 P2(L + l)=P2 < L +1) + (lOt < 0.6609+7.5*B(G)/(237.3+B < G >)>-10* < 0.66 09+7.5*(B(G)-T l)/(237.3+B(G)-Tl>))*H 1320 IF ABS(B(I)-T)>ABS(B(I+1)~T> THEN 1390 0 1330 E(L) = (B(I-2)-8*B<I-1)+8*B(I+1)~B(I+2))/(12*H) 1340 PI<L)=<lOt(0.6609+7.5*T/<237.3+T))~10t(0.6609+7.5*(T-Tl)/(2 3 7.3+T-Tl>))*H 0 1350 F1=<M1(1)-M1)/P2 1360 M1(L+1)^M1(L)"F1*(P2(L+1)-P2(L)) 1370 T=T + S 0 1330 L=L + 1

1390 IF T<T2 THEN 1400 : I=G:NEXT I:GOTO 1410 1400 NEXT I 0 1410 DATA 8 1420 DATA 2438.6E3,2476.9E3,2465.1E3,2453.4E3,2441.7E3,2429.9E3, 2418.0E3,2406.2E3 0 1430 READ N 1440 L=0 1470PRINTU3ING 2310 0 1480 T,L = 0 1490 FOR K=l TO N 1500 L = L + 1 1510 T = T + 5 1520 PRINTUSING 2280,TjD(K);E<L);M(K+1);Ml(L+l) 1530 NEXT K Q 1540 PRINT • • -1550 PRINT 1560PRINTUSING 2320 0 1570 T,L=0

1531 DATA SAVE OPEN "PLOT.DAT" 1532 DATA SAVE 4,2

0 1533 DATA SAVE 8

1584 DATA SAVE "TEMP.('C)","W.PROD.<MW/KG)","W.PROD.*(MW/KG)", "TYD(UREN)"

0 1S35DATA SAVE U*

1536 FOR K=l TO N

1590 IF M(K+1)>0 THEN 1600:NEXT KîGOTO 1721

0 1600 T=T+S 1610 READ V 1620 L=L+1 0 1630 IF Z=0 THEN 1670 1640 C(K) = < M(K + l)*P + V*(F*P2< K)*M1< L + l)/3600-Fl*Pl(L)*M <K + l)/3600 ))/(M(K +1)*M1(L + 1 ) * <E (L)-D(K))) 0 1650 Ul=(V*F*P2<K)/(M<l^l)) + ( <M(K + l>*C<K)+4200)/(M(K + l))) *D(K >>/3.6 1660 GOTO 1700 0 1670 Wl=((<M(K+l)*Cl+4200)/(M<K+l>)) *D(K))/3.6 1630 W=(V*F*P2(K)/(M(K + l)) + ((M(K + l)*C1 + 4200)/M(K + l)) *D< K ))/3.6 0 1690 Z1=W1*24*3600/<1000*4.19002) 1700 Z3=W*24*3600/(1000*4.19002) 1701 DATA SAVE T,W,W1,0 0 1710PRI.NTUSING 2290, T; C<K) J Ml; U; Zl; Z3 1720 NEXT K

1721DATA SAVE END

0 172ÓDATA SAVE U$

1727 FOR T3=1T0 Y STEP. 5

1728Y3-B1<1)+B1<2)«T3+B1(3)*T3t2+Bl(4)*T3t3 + Bl(5)*T3t4+Bl< 6)*T3t dft 5

(14)

1745 PRINTUSING 2350, "ALPHA 1==" ; F, "ALPHA 2=";F1

1748 DATA SAVE END

1749 REWIND

1750 END

2280 %«««.««

2290 %«««.««

«ttaaa.tt

2300

**&«8HM8 88. 888888

2310 % T, 'C DT/DT1,'C/H

2, KG

888.888

8888.8 8888.8

888.888

8888.8

D T / D T 2 , ' C / H ttfltl.tttiti

88888.8

ft

2320 %T,'C

C,J/KG'C

W*,MW/KG

W,MW/KG

Ml,KG M

W*,KCAL/TON

24H

2350 %

W,KCAL/TGN24H

8888888 88.888888

0

; -

T

T

(15)

-^^•"MaMMWWiii i.4*^.*u-m, eg 00 «3 O O LU e;

Û

Z

<!

T,

or

2

0

Ü

H

h iL

Ü

Ü3

H O

<EÙ.

o o ca o Ü » o o o o o o o C N - T O C O O O O O ü C ^ C s C N C N p O O O n CM f-4 CN N O O O O X N U » O »•O C-4 O H • Q O N O O O o o o o o o o o o o o o o o o o o o o Cl o O X ro N o u KI s rs ». Cs \~ Os Q Û ii - * Lu Lu u o o O O O T H O O O O rtMMNOOOO « » » U '. * 9 : 9 a O O Q O O O O O CC cc o u u O O O O O O Ö O o o o o o o o o H-.L0 o m O U) O 10 o •n TH r-4 r-4 ro ro 'T

x

•T f-4 O H V , M N < r H rs o c-i <N CJ Kl LO -o co i : Ki -o -ri r-s X sT CN O J -N 'T CO Kl -N —J < • • • <C Os CO TH O U CN Os CO »O ü so <• M O • N i n o ' O * C3 ü s o <T «O 5T 3 Kl tO Ki 'T X. sO tn so o ».T-H Ki L O CO 3 C3 \ 3 uo sQ uo CO TC » • « « ».N sQ CO CO * r-4 so es fs, 3 1-1 C-4 sT rs Ü O O O O v. o o o o

u

o o o o o o o o o • » > • s i f l O ü l O : •> t-t »-I f -4 co LU n L ü ^ - , ^ - , H- • :• X «o <o -l-> +-> te- « io o * ? 5 en O) o o o o o o Ü: o o LU Os LU rs. ca o M o o: o ££ « • O O O CJ LU P f« C-4 f - <c <x LU X X M Cu Cu 2 - J _J II <c <c *

(16)

9.75 r

9.00

8.25 ..

7.50 .

6.75 •

6.00 ..

5.25 .

4.50 .

3.75 ..

3.00

2.25 ..

1.50 .

J-0.20 0.40 0.60 0.80 1.00 1.20 1.40 1.60 1.80 2.00 2.20 2.40 2.60

LEGENDE GRAFIEK 1

X -AS (schaalfactor 10 ) - TEMP.C'C)

Y -AS (schaalfactor 100 ) - W.PROD.(MW/KG)

(17)

Ni et-gecorrigeerde warmteproduktie - temperatuurcurve van

handgeoogste champignons

9.00 8.25 • 7.50 -. 6.75 .. 6.00 . 5.25 .. 4.50 ,. 3.75 • 3.00 • 2.25 •• 1.50 .. 0.75 .. ) 1 ! T t ! ' ' ' ' ' ' ' 0.20 0.40 0.60 0.80 1.00 1.20 1.40 1.60 1.80 2.00 2.20 2.40 2.60 LEGENDE GRAFIEK 2

X -AS (schaalfactor 10 ) - TEMP.('C)

V -AS (schaalfactor 100 ) - W.PROD.*(MW/KG)

(18)

2.50 4.

2.25 f

2.00 f

1.75 l

1.50 f

1.25 f

1 . 0 0 . 0 . 7 5 0 . 5 0 0 . 2 5 '. i

0.00

' m • » -• • f t . . • • • » • • _ j • • r,_ • • 1 t , . . ..t..,

-0.10 0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 1.10

LEGENDE GRAFIEK 3

X -AS ( s c h a a l f a c t o r 100 ) - TYD(UREN)

Y -AS ( s c h a a l f a c t o r 10 ) - TEMP.('C)

•••«»CHAMPIGNONS HANDßEOQGST

(19)

^ "•iltriti'iriVhfcV H- O- 04 r O ^ O M O "s, • • t « » • _ l LO f-4 CS CO CO CO <r Lo i s rn Lo I N co CJ -«o rs r-4 LO co c-4 iC H H H f J O o o o o o o o o o o o o o o o o o o o o o o o » . o o o o o o o o - r r 3 K i r n 5J3 i a o o CO LO C--1 CS '-0 »H O O «r <r - ? r o r-o K Ï o o 1~< 1 - i 1-1 T H O O Ct o r-\ fS _ l » < r <r-CJ O v ; u-j SN * 3 Cs IS -T i-i TH < I C>- O <• <• O• TH CO o-, TH <r SJ i-I T-i i - l

Ld

CO

u

HO

M» L i CJ o o o o o o o o o o o o o o o o o o o o o o o o o o o s Q xxt

0

S.T3 >is J Î. y

ui

i; II

Jh

w »J J T A r i co IS CO CJ C-* !> !>• CN 3 II LL u. UJ o s CJ -• <r rs ». r-j CM i-t o o i - » • Û o o H Q LO r s -T u-j O O i a o o rs o o o o co o o o o o o a i v i o o o o o o o o o o o o CJ o o o o o o o o :c ». LO o LQ o 3 i— TH T-i r-j 0 O LO O •r-i i--i C-4 K i 1*0 <r CD s O <f r l Ifl O &• o 3 i-t :e r-o 9 « 3 CJ X.-* N . 3 <r j » ».•er * r-4 3 O 's CJ o ü » s o *-i • t u rs o--ro LO LO e s a 1 r s LO r-4 -sr o o » * o o LO rs i - i 3 CO LO o a o o o o o o CJ » o o * » s LO O LO ». H •tH i H o * o r-4 T H 0s K I a Cu <i O o o o a LO F-4 o 1-1 1-1 «o * 1-t co o 3 O o o a o !"0 CO LU n _ I CC U I

>ST^

!— . . - u « m CJ CL CL o > > s i i -- ^ fl) 8) 4-> 4-> û i «o re o X X ÏT C7> O) C i CO O LU rs o LU TH O CJ o o w O O L i O O ÛC « » o o o CJ LU CD i n r-4 H - <E <C LU x ir H L L Ù . •z. - 1 _ l ii <r <r *

(20)

1.40

f

1 . 3 0 y. < 1.20 -1.10 '. 1.00 . 0.9 0 '•• 0.80 '.

0.70

'.. 0.60

i-0.50

'••

0.40

'..

0.30

'.. J.

JL

J-

J_

J.

X

JL

J- J

0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 2.25 2.50 2.75 3.00 3.25

LEGENDE GRAFIEK 1 X -AS ( s c h a a l f a c t o r 1 0 ) - TEMP.(»C) Y -AS ( s c h a a l f a c t o r 1 0 0 ) - W.PROD.(MW/KG) •••-SCHORSENEREN

(21)

Ni et-gecorrigeerde warmteproduktie - temperatuur curve van schorseneren

1.05 0.97 .. 0.90 .. 0.82 ".. 0.75 .. 0.67 •• 0.60 -0.52 .. 0.45 . 0.37 .. 0.30 •• 0.22 ..

J L.

0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 2.25 2.50 2.75 3.00 3.25

LEOENDE GRAFIEK 2

X -AS (schaalfactor 10 ) - TEMP.C'C)

Y -AS (schaalfactor 100 ) - W.PROD.*(MW/KG)

•.»«SCHORSENEREN

(22)

4.00 .

3.60 •

3.20

•-2.80 •

2.40 .

2.00 ..

1.60 .

1.20

-0.80 •

0.40 •

0.00 •

-L.

_1 L

-L- L

-L

-L

-0.10 0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 1.10

LEGENDE GRAFIEK 3

X -AS (schaalfactor 1000 ) - TYD(UREN)

Y -AS (schaalfactor 10 )= TEMP.('C)

(23)

-T C-i •z. O ï - -O o CO N . , 3 a - j o <r co <z <t O £>• U N - < J -Si Cr- Cl K ! a • a LO to Ci I N PS rH CO 1-1 M ü CN n «s- 10 FN co o o o o o o o o o o o o o o o o c a o o o o o o © o 1 ^ s a i a s a s - . . » . o © o o o o o o x: > r c o T - i < r r N o - o o CD TH O O *>• CO fN o O v ; LO 10 lO 'T <- *r O O s s a a a a a a a a T-) »O < i -O < l -O »O o o o Η N. O !N O O- O fN —I a a a a a a <î H O -O N H Ci Ci <1 Ci O . H H O x -T c-i o Ci ro a »• r-4 f-'j <r io ~o tN * 3

U

J

Lfe

Q

€ v

^33 3B3 U J

€û

HO

0 iV and su« a - o o o o o o o o ». o o o o o o o o C-i o o o o o o o o ^— a a a a a a a » O O O O O O O O O N. t -Û co 0-- N, i v Ci-, o i-0 o <r o o o . ». CN ra <i -M to o- o o

r-i o i-i a~< c-i CN ei o o

Q O O O N. O O O O Lü O O O o o o o o o o •î CJ o o o o o o o o Jl* *• 3 3 3 3 3 3 3 3 :;:' ». to o i-O o io o LO O

D r TH T-I a c-i ro r-o <r

CD ü O- fN O CN N , a a i • 3 c-i <r n -r z: r-o fN r-j co CO a fN r H a r*) -T o ». -rt TH CN c-i n 3 CD - v N 3 r-o <r <i co .Tü 3 3 3 3 ». CN »0 CO CN *: T-i LO 0~' »O 3 H H H n O *. CD O O O O ü a a a a N. O O O O »« CJ o o o o o o o o C O 3 B 3 a N LO O u".i O o 3 -o • * " " " : p} o t o o o a LO C-i <r K i « o <r ro o 3 O O o a O n co u i—i _ j LU

???

X (0 (8 O O- O .

°>>

z> *-- ^ 0) CU 4-> 4-> Ü (0 (O o s s o - v D i r a - * -V - * »*_-^^»—^ Q cs: co o u <r o LU i-l o CD O O M O O CC O O CC a a o o o CD Lü CD r-ï CN H <r <r Lü X X H Û . L I z: _J _J *

(24)

5.20 ..

4.80 ..

4.40 .

4.00 ..

3.60

-3.20 ,.

2.80

2.40

-2.00 .

1.60

-1.20

JL JL -L

J-0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 2.25 2.50 2.75 3.00 3.25

LEGENDE GRAFIEK 1

X -AS (schaalfactor 10 ) - TEMP.CC)

Y -Af) (schaalfactor 100 )» W.PROD. (MW/KG)

(25)

Ni et-gecorrigeerde warmteproduktie - temperatuur curve van witlofwortelen

3.75 r 3.50 .. 3.25 3.00 2.75 . 2.50 ,. 2.25 .. 2.00 . 1.75 • 1.50 •. 1.25 • 1.00 .. J. Jl L 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 2.25 2.50 2.75 3.00 3.25 LEGENDR GRAFIEK 2

X -AS (schaalfactor 10 )= TEMP.('C)

Y -AS (schaalfactor 100 )= W.PROD.*(MW/KG) •••=VlTLOFWORTELEN

(26)

4.00 3.60 -3.20 2.80 .. 2.40 2.00 1.60 1.20 -0.80 • 0.40 .. 0.00 .. X X X X -0.25 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 2.25 2.50 2.75 LEGENDE GRAFIEK 3 X -AS ( s c h a a l f a c t o r 1 0 0 ) = TYD(UREN) Y -AS ( s c h a a l f a c t o r 1 0 ) = TEMP.C'C) •WITLOFWORTELEN

Referenties

GERELATEERDE DOCUMENTEN

fulvum, en dat deze effec- tor in staat is een aantal cysteïne-proteasen te remmen die vereist zijn voor de basale afweer van tomaat tegen verscheidene schimmels.. Ook voor de

Ook de financiering met 1,5 miljoen euro door ZonMw van het zogenaamde WIDE-project waarin het Antoni van Leeuwenhoek, het UMC Utrecht en de Hartwig Medical Foundation

Cardiovascular risk reduction in high- risk pediatric patients: a scientific statement from the American Heart Association Expert Panel on Population and Prevention Science;

In het zuiden waar (geheel tegen het dominante beeld in de historiografie in) kern- gezinnen domineerden, waren er betrekkelijk weinig dienstboden in de bevolking, terwijl in

An additional sign that Australian manufacturing was indeed in a more advanced state of development in the early twentieth century is the fact that it exported harvest machinery

Inmiddels zijn in het materiaal van 1961 enkele wijzigingen aangebracht door het toepassen van de nauwkeuriger formule voor het bepalen van de gebouwenkosten. Bovendien is het

Avis de la Commission paritaire relatif à l’avis de la Chambre de chiropraxie concernant le niveau de formation pour atteindre le profil requis en

These life expectancy tables, survival charts and also the available prognostic scores at the start of dialysis (REIN score [France], aREIN score [NBVN]) can help to inform