Procesvoorontwerp Claus-proces

TECHNISCHE HOGESCHOOL DELFT.

LABORATORIU:.I VODR CHK.:lscrn:: TECHnOLOGIE.

PROCESVOORONT~ERP CLAUS-PROCES

In samenwerking ~et het Ingenieursbureau COMPRIhlO N.V. te Amsterdam.

DELFT, APRIL 1969.

G. van Os

Jac. van Beierenlaan 67, D~lft.

C.Vl.A.Schram

Oude Delft 29D, Delft.

P. \'[essel s

I SAMENVATTING

Voor een bestaande uitvoering van het straight-through claus-proces werd de stofbalans en ten dele

ook de warmtebalans over de verschillende onderdelen

van het proces en over het hele proces opgesteld.

Het claus-proces werd door middel van een

kompu-terprogramma berekend voor verschillende omstandighe-den.

-2-11 INHOUDSOPGAVE blz. 1 Samenvatting 1 11 Inhoudsopgave 2 111 Inleiding ₃

IV Nadere bespreking algol-programma 5

V Resultaten en konklusies 6

Bijlagen:

A Nadere beschouwing straight-through proces 8

B C E F Thermodynamika Thermodynamische gegevens

Aannamen bij de berekeningswijze Berekeningswijze brander eerste kondensor eerste line-burner eerste reaktor tweede kondensor

derde kondensor en demister incinerator totaal G Figuren en tabellen H Algol-programma JProcesberekening K Symbolenlijst L Literatuurlijst 11

14

16

11

11

21 22

25

21

28

29

32

33

42

52

51

59

III lHLElDlNG

De totale zwavelproduktie in de niet-kommunistische landen bedroeg in 1966 24,2 miljoen ton tegen een vraag van 25,0 miljoen ton (tab. 1 ,lito 1 ). Uit de schattin-gen voor 1970 volgt, dat er een zwaveloverschot dreigt te ontstaan. Een van de oorzaken hiervan is de sterk toene-mende zwavelwinning uit aardoliën.

Het grootste deel van de zwavel wo~dt gewonnen d.m.v. de ekonomisch gunstige frasch-processen. De winning van zwavel uit de bij de raffinage van aardolie vrijkomende zwavelwaterstof is niet in de eerste plaats een ekonom i-sche kwestie, maar is noodzakelijk vanwege de luchtver-ontreinigende en korrosieve eigenschappen van zwavelwater-stof. De oudere processen voor de winning van zwavel uit pyriet en gips zijn ekonomisch niet zo gunstig.

Tab. 2 en fig.

5

geven een indruk van de kost- en verkoopsprijzen van de zwavel bereid volgens de verschil-lende procedé's. Door het in de volgende jaren verwachte zwaveloverschot zal de prijs nadelig beïnvloed worden (lit. 1 ), waardoor de processen ter bereiding van zwavel

met een hoge kostprijs minder rendabel worden.

De zwavel wordt voor meer dan 80~ gebruikt voor de bereiding van zuavelzuur (fig.

6 ).

Voor de winning van zwavel uit zwavelwaterstof staan afhankelijk van de samenstelling van de voeding verschil-lende procedé's ter beschikking (lit. 2 ). Voor gas met

een zwavelwaterstofgehalte groter dan 15/0 en een gehalte aan lichte kooh/aterstoffen kleiner dan 2~; wordt het ver-beterde claus-proces toegepast. Het claus-proces kent twee hoofdvarianten, namelijk het split-flow proces voor

-4-gassen met een zwavelwaterstofgehalte kleiner dan 30% en het straight-through proces voor gassen met een zwavel

wa-terstofgehalte groter dan 50%.

Het processchema van het straight-through claus-pro-ces is afgebeeld in fig.11 •

De doelstelling van deze opdracht is het opstellen van een komputerprogramma van het in fig. 11 afgebeelde straight-through proces, teneinde dit proces voor ver-schillende omstandigheden snel te kunnen berekenen.

IV ALGOL-PHOGRAMilA

De berekeningen, uitgevoerd aan het straight-through claus-proces, zoals dit afgebeeld staat in fig.11 , werden samengebracht in een "ALGOL 60" komputerprogramma.

Dit progra~~a is als volgt samengesteld: 1 • Hoofdprogramma

2. Procedure 'burner'

3.

Procedure 'colire'

4.

Procedure 'incinerator'

5.

Getalband

In de verschillende onderdelen van het programma worden de volgende berekeningen uitgevoerd:

In de procedure 'burner' worden de konversie en de daardoor bepaalde samenstelling en temperatuur aan de uit-gang van de brander berekend.

In de procedure 'colire' wordt berekend: de zwavel-produktie in de kondensor, de in de line-burner benodigde hoeveelheden fuel gas en lucht, de samenstelling van het gas aan de uitgang van de line-burner, de konversie in de reaktor en de daarmee samenhangende samenstelling en tem-peratuur aan de uitgang van de reaktor. De procedure 'co-lire' wordt in de procesberekening tweemaal doorlopen.

In de procedure 'incinerator' wordt berekend: de zwa-velproduktie in de derde kondensor en demister, de in de incinerator benodigde hoeveelheden fuel gas en lucht voor de verbranding van het fuel gas, de zwavel en de zwavelwa-terstof, de samenstelling en zwaveldioxydegehalte aan de uitgang van de inGinerator, en de hoeveelheid lucht om te voldoen aan de specifikatie van het zwaveldioxydegehalte in het schoorsteengas.

Een voorbeeld van een volledige procesberekening staat afgebeeld in b~jl.J •

-6-V RESULTATEN EN KONKLUSIES

Met behulp van het komputerprogramma werd het claus-proces berekend voor een kapaciteit van circa 100 ton per dag. Bijlage J geeft een voorbeeld van een volledig uit-gevoerde procesberekening berekend met het komputerpro-gramma uit bijlage H.

Daarnaast werden de invloeden nagegaan van verschil-lende procesgrootheden, zoals temperatuur van de voeding, samenstelling van de voeding, en de temperatuur van de lucht.

Fig.10 geeft de invloed weer van de temperatuur van de voeding op de totale recovery. De recovery blijkt sterk afhankelijk te zijn van de voedingstemperatuur. Naarmate

,. -..

-de temperetuur hoger wordt, daalt de recovery sterker. Fig. 9 geeft de invloed weer van de luchttemperatuur op de recovery; een verhoogde luchttemperatuur geeft een daling van de recovery.

De invloed van het volumepercentage propaan en butaan in de voeding (fig.

3 )

blijkt vrij gering te zijn; een verhoging van het koolwaterstofgehalte geeft een geringe daling van de recovery.

Tenslotte werden globaal de invloeden nagegaan van het volumepercentage zwavelwaterstof in de voeding, de temperatuur aan de ingang van de eerste reaktor en de re-latieve vochtigheid van de lucht:

vol.

%

zwavelwaterstof 64,68

90,85

ingangstemp. eerste reaktor 510 OK 500 OK recovery 94,57

%

94,73

%

recovery 94,73

%

94,86

%

relatieve vochtigheid lucht (290 oK) 20

%

70

%

recovery

94,83

%

94,73

%

Verlaging van de temperatuur aan de ingang van de re-aktor geeft volgens de berekening een verhoogde recovery. De mogelijkheid bestaat echter, dat bij een te sterke tem-peratuurverlaging de kinetiek dermate verlaagd wordt, dat er in de reaktor geen evenwicht meer bereikt wordt, het-geen weer een daling van de recovery zou doen veroorzaken.

In het algemeen kan gekonkludeerd worden, dat in de

praktijk een recovery van

94,5

tot 95,0

%

verwacht mag

worden, afhankelijk van de verschillende grootheden. Slechts de verhoging van de temperatuur van de voeding kan een

-8-A NADERE B~SCHOUWING STRAIGHT THROUGH PROCES

De zwavel verbindingen in nafta en gasolie worden in een hydrotreater omgezet in zwavelwaterstof, hetgeen in een aminetreater verwijderd wordt. Het aldus vrijgemaak-te zwavelwavrijgemaak-terstof is de voeding voor het straight-through claus-proces.

Na het passeren van een knock-out-drum voor het vè~­

wijde re ri van water, bevat het gas 50 tot 90 vol.% zwavel-waterstof, maximaal 2 vol.% koolwaterstoffen (vnl. propaan en butaan) (lit. 2 ), is verzadigd met waterdamp, terwijl het resterende deel kooldioxyde is.

Dit gasmengsel wordt in een brander geleid, waarbij een derde deel van de zwavelwa_{terstof volledi b wordt}

om-gezet in zwaveldioxyde en water volgens de reaktie:

De voor de omzetting van de zwavelwaterstof en de verbranding van de koolwaterstoffen benodigde zuurstof wordt d.m.v. lucht in een stoechiometrische verhouding gedoseerd. Deze reakties zijn sterk exotherm en worden aflopend verondersteld. Fig.

14

geeft een beeld van de invloed van eventuele overmaat of ondermaat lucht (lit.3 ). De molaire verhouding van het ontstane zwaveldioxyde en de zwavelwaterstof dient precies twee te zijn. De gevorm-de zwaveldioxygevorm-de en gevorm-de zwavelwaterstof gatn in de brander de volgende evenwichtsreaktie aan:

Dit evenwicht is endotherm indien e=2 en exotherm indien

-

..."

-

_

.

r

\

~

e=6 of S, terwijl de

zwavelwaterstof/zwaveldioxyde-ver-=

"'houding konstant blijft.,

De zwaveldamp kan afhankelijk van de temperatuur als SS' S6 en S2 voorkomen (fig.

4 ,

lito

4 ).

Bij de in de brander optredende temperaturen (1200 tot 1400oC) komt zwaveldamp alleen in de vorm van S2 voor.

De stabiliteit van de vlam in de brander vereist een minimaal zwavelwaterstofgehal te van circa 40 vol.;~; bij lagere percentages maakt men daarom gebruik van het split-flow proces.

De reaktiegassen worden in een waste-heat-boiler afgekoeld tot een temperatuur, die afhankelijk is van de druk van de op te '..;ekken stoom. De zwavel wordt in de eerste kondensor gekondenseerd (circa 60 vol.% van de zwavel in de voeding) en naar een zwaveltank afge-voerd. Aangezien de viskositeit van vloeibare zwavel binnen het temperatuurtrajekt 170 250 oe een maximum vertoont (fig. 13,lit.

5 ),

dient de uitgangstempera-tuur niet binnen dit trajekt te vallen. Daar bij tempe-raturen boven 250 oe de dampspan~ing van de vloeibare zwavel groter is dan bij 160 oe, wordt een uitgangstem-peratuur van circa 160 oe genomen.

De zwavelmist, die uit deze kondensor komt, wordt niet door een demister afgevangen, aangezien dit een sterke korrosie aan de gaasbedden van de demister geeft.

In de kondensor kan gekoeld worden met water of af-gewerkte stoom.

In een line-burner wordt fuel-gas met lucht verbrand, waarna de verbrandingsgassen met de gassen uit de

konden-sor worden gemengd. Voor andere herverhittingsmethoden, I ~

zie li t. 2,3,6,11. De temperatuurverhoging is nodig om de j..

~~1.k~~~~L--, " IJ?!

minimale aansla~temperatuur van de katalytische reaktie

Sq

-te bereiken (lit~1,S). Als katalysator wordt gesinterd bauxiet-gruis gebruikt. In de reaktor treedt de volgen-de evenwichtsreaktie op (fig. 12 ):

-10-waarin ex de dissociatiegraad van S8 is (fig. 4 ). De

uitgangstemper~tuur van de reaktor is 250 - 300 °C.

De gevormde zwavel wordt ~~~ln- tweede kondensor afgescheiden (circa 28 ~. ï~ van de zwavel in de voeding), waarna de rea~tiegassen in een twe8de line-burner wor-den opgewarmd om vervolgens in een tweede reaktor geleid te worden. De in de tHeede reaktor gevormde zwa.vel wordt in een derde kondensor afgescheiden ( circa.

5

.~.% van de zwavel in de voeding) .

De drie kondensors zijn in één apparaat ondergebracht. De zwavelmist wordt in tegenstelling tot na de eer-ste en tweede kondensor nu wel in een demieer-ster afgeschei-den (efficiency 99~).

Bij een temperatuur van ongeveer 600 °C, die gerea-liseerd wordt door verbranding van fuel-gas met lucht, wordt in de incinerator de nog in het gas aanwezige zwa-vel en zwazwa-velwaterstof met 100~ overmaat lucht omgezet in zwaveldioxyde.

Na verdunning van het gasmengsel met lucht ter af-koeling en ter verla.ging van het zwaveldioxydegehalte, wordt het gas via de schoorsteen afgevoerd.

B THl:;H,iWDYN;'.lUKA

A. De reaktiewarmten als funktie van de temperatuur wer-den met de volgende vergelijkingen berekend (lit.

9 ):

T.

t, HO =t, HO + jt,CPdT

T s

T s

B.

De evenwichtsliggingen als ~unktie van de temperatuur werden berekend met (lit~10):

dIn K (T) = P 2 = a. + b.T + e.T l. l. 1. In

K

_{= -nT} .6

H~

+ ~.ln la T + . 6 b . 6 c 2 2R· T + 6R. T + I

waarbij de integratiekonstante I werd berekend d.m.v.:

.6G

0

=

-RT

.ln K

(T )

S 8 P S

t,.n Het verband tussen K en K luidt: K =K .p

P x P x

C. Voor de koolwaterstoffen propaan en butaan, waarvoor in de literatuur geen soortelijke warmten als funktie van de tem~eratuur werden gevonden, werden deze funk-ties berekend uit gegevens (lit. 12) met een door Janz beschreven methode (lit.14 ):

-12-2 e a + bT + cT p de d2e ~ _{b + 2eT} ----.:Q 2e d'r dT2

Hieruit zijn nu a, b en c te berekenen.

D. Voor de berekening van de konversie werd gebruik

ge-maakt van de gekoppelde warmte- en stofbalans volgens

Kramers en Westerterp (lit. 9 ):

Cp

=rv.e .

i ~ p~

o

E. Voor de berekening van de hoeveelheden water in het

voedingsgas en in de lucht afhankelijk van de

verza-digde dampspanning werd gebruik gemaakt van:

voor een ideaal gas

I

r

zodat

I

P

=

P

-14-C THERHODYNArnSCHE GEGEVENS (lit. 9,12,13)

I. Soortelijke warmte : c _p

₌

a + bT + cT 2 + diT (caltfmOl)

a b c d zwavel (S2) 3,58 0,00030 zwavel (S6) 19,25 0,Q02664 zwavel (S8) 25,OQ Q,003552 zwavelwaterstof 7,20 0,0036Q zwaveldioxyde 7,70 0,00530 8,3.10

-7

water (damp) 8,22 0,Q0015 _{1 ,34.1}

°

-6 kooldioxyde 1 0,34 0,OQ274 195500 zuurstof 8,27 0,000~58 18770 stikstof 6,50 Q,00100 4,76 0,052 -6 propaan 14,7.10 butaan 7,53 0,0654 18,7.10 -6

11. Vormingsenthalpie:

6H~

(25

°c,

1 atm) (kcal/mol) zwavel (S2) 31 ,02 _ zwavel (S8) 27,090 zwavel (S6) 27,78 zwavelwaterstof

-4,77

zwaveldioxyde -70,94 water -57,1979 kooldioxyde -94,05 1 .;I zuurstof

°

stikstof

°

p!,opaan -24,820 n-butaan -29,812

111. Vrije vormingsenthalpie: l'lGo(250C,1atm) s zwavel (S2) 19,36 zwavel (Se) 12,770 zwavel (S6) 13,97 zwavelwaterstof -7,85 zv/aveldioxyde -7 1 ,68 water -54,635 propaan -5,614 butaan -3,754

IV.

Verbrandingswarmte: waterstof 57,7979 methaan 191,759 etheen 316,195 ethaan 341,261 prope.en 460,4 28 propaan 488,527 butaan 635,384 z"Tavel (S6) 453,42 zwavel (S ) 8 594,61 zwavelwaterstof 66,17 (kcal/mol) (kcal/mol)

-16-E AANNMdEN BIJ DE BEREK~NINGSWIJZE

Brander:

Bij de berekening van de konversie en de temperatuur

aan de uitgang van de brander werd verondersteld, dat het

thermodynamische evenwicht bereikt werd. Deze

veronderstel-ling wordt gerechtvaardigd door gegevens uit de praktijk.

Mogelijk vindt dit zijn oorzaak in de hoge temperaturen,

die in de brander optreden. In de literatuur werden geen

gegevens betreffende de kinetiek gevonden.

l,-,-<~

t,JS

Nevenreakties (S, 11, 12,13 ) werden, ter

vereenvou-'Ç'(\"\\ C~t di ging van de berekening, verondersteld niet op te treden. De brander werd verondersteld adiabatisch te werken.

Line-burner en incinerator:

Bij de berekening van de line-burners en de

incine-rator werden de verbrandingswarmten van de verschillende

komponenten van het fuel-gas onafhankelijk van de

ratuur verondersteld, aangezien de invloed van de

tempe-ratuur op de verbrandingswarmten gering is en dit de

bere-kening aanzienlijk vereenvoudigd.

Reaktoren:

Bij de berekening van de reaktoren werd de

tempera-tuur aan de ingang van de reaktor aan de hand vnn

waar-den uit de literatuur (lit.S,11) en de praktijk aangenomen. Om analoge redenen als bij de brander is ook hier de

be-rekende konversie de evenwichtskonversie

De nevenreaktie zijn ook hier verwaarloosd, alhoewel

bij hogere kooldioxydegehalten in de voeding de invloed

van de nevenreakties toeneemt, hetgeen zich in de reaktor

manifesteert door een hogere temperatuur aan de uitgang,

F BEfillKENING::i\<fIJZE

Brander:

In de brander treden de volgende reakties op (lit. ):

H 28 +

i-0

2 «8°2 + H20 (1 ) -2H 28 + °2 2H2O + S 2 (5) 8 2 + 2°2 28°2 (6) 2H 28 + 8°2 - - - - 2H2O +

3/

282 (7) -CO 2 + H28 ==== C08 + H20 (a) CO₂ + 2H₂8 ==== C8 2 + 2H2O (11 ) À.C08 + _8°2 ==== 2.C0₂ +311..8 2 (12 ) C8 2 + 8°2 ---- CO2 +

3/

282 (13) C

₃

Ha + 5°2 3C0₂ +

4

H2O (9) -C

4

H10 + 6-?zo2

4

C02 + 5H2O (10)

-Daar voor de berekening van de konversie en de tem-peratuur aan de uitgang van de brander slechts begin- en eindtoestand beschouwd werden, zijn alleen de reakties (1,

7, 9,

10) verondersteld op te treden. De reakties (a,11,12,13 ) werden verwaarloosd. Het reaktieverloop in de brander is dan als volgt schematisch voo te stellen:

--'-l

~ . : . '," ,:--,,~''':.!.~.~. ~:':. :'.:.! ;' ,;",,: ~ •

. .

.

~ . \ . . . ..

f

'

••

De

t

tbi tb1

t

voeding lucht

I

I

3H 8 ₂ + _{3/2°2 8°2} + H 20 C

₃

H

a

+ 5°2

3

C02 +

4

H2O C

4

H10 + 6t02

4

C02 + 5H2O exotherm (x

=

1)

I

tb1

t

2H₂S + S02 + H 20 ====

3

H2O + endotherm

I

tb2

T

voeding bestaat uit:

zwavelwaterstof b [1] kmol/uur

propaan b [2] kmol/uur

n-butaan b

[3]

kmo1/uur

kool dioxyde b [4] kmol/uur

water

deze werd berekend uit de verzadigde bij de temperatuur van de voeding en

+ 2H 2S

3/

2S 2 dampspanning bedraagt: :.;" 1

'{ b

I. ( 0-~~ tJg ~ .1

•

-19-De benodigde zuurstof kmol/h

Stikstof in toegevoerde lucht: 3, 768(Jb [1) +5b [2] +6·~b[3]) kraol/h

Water in toegevoerde lucht :4,768(}b[1]+5b[2]+61b[3])/

(1/betha.O,01.exp(-4882,7(1/tbl-1/373,2))-1) kmol/h

Bij de berekening werden alle hoeveelheden betrokken op drie molen zwavelwaterstof: zwavelwaterstof 3b[1]/b[1] ₃ propaan 3b[2]/b[1] pb n-butaan 3b [3]/b[1]

₌

qb kool dioxyde 3b [4]/b [1J

₌

mb zuurstof 3/2+5pb+6~qb rb stikstof 3,768rb lb water (totaal) kb ~ (3+pb+qb+mb)/(p/exp(-48 82 ,7(1/tbi-1/373,2))-1)+ +4,768rb/(1/(betha.O,01.exp(-4882 ,7(1/tbl-1/373,2))-1)

De temperatuur (tb1) werd berekend met een gekoppelde

stof-en warmtebalans (lit.

9 ):

tb1

I

tbi (3c +rb.c +kb.c +lb.c +mb.c +pb.c + PH2S P02 PH20 PN2 PC02 _PC3H8

~H(1)

=

-112104-3,685.tb1+7,32.10-4.tb12+1,7.10-7.tb13+ -28155/tb1 cal/mol

~H(9) =

-493646+17,79.tb1-0,02216.tb12+6,71.10-6.tb1 3+ +493650/tb1 cal/mol

~H(10)=

-641671+21,17.tb1-0,0277.tb12+8,47.10-6.tb1 3+ +660000/tb1 cal/mol

De .hoeveelheden ges, betrokken op drie molen zwavelwater-stof in de voeding, worden na de verbranding:

zwavelwaterstof zwaveldioxyde water stikstof kooldioxyde 2 1 kb+4pb+5qb+1 lb mb+3pb+4qb

De konversie (x) en de temperatuur (tb2) aan de uitgang van de brander werden berekend met een gekoppelde stof-en warmtebalans (lit. 9 ): tb2

I

(2c +c +(kb+4pb+5qb+1).c +lb.c + PH2S PS02 . _{PH20 P1T2} tb1 . +(mb+3pb+4qb) . c ) _odT + x ol:I H(2 ) tb2 = 0 PC02 x-( ) 6 2 ' -6 3 / t.H,\2 tb2=13583,2+5,9 5.tb2-0,00544.tb2 +1,17.10 . • tb2 cal mol

De samenstelling van het gas uit de brander is:

zwavelwaterstof h2s=2b{1] .(1-x)/3 kmol/uur

zwaveldioxyde so2=b[1] .(1-x)/3 kmol/uur

water h2o=(kb+4pb+5qb+2x+1 ). b [1] /3 kmol/uur kooldioxyde co2= (mb+ 3pb+4qb) • b [1] /3 kmol/uur

stikstof n2 =lb.b[1]/3 kmol/uur

-21-Eerste kondensor:

De temperatuur aan de uitgang van de kondensor werd

niet berekend, doch er werd een waarde voor aangenomen. De zwavelopbrengst in de eerste kondensor bedraagt:

op1 = b [1] -h2s-so2-perc1 • b [1] kmol/uur

hierin is perc1 de uit de kondensor komende heid zwavelmist en zwaveldamp als fraktie van de

hoeveel-heid zwavel in de voeding. Voor de zwavelmist zijn alleen ~

praktijkwaarden bekend.

De zwaveldamp in het gas aan de uitgang van de konden-sor bevat bij die temperatuur (tuc1) zowel S8 als S6' De

afhankelijkheid v_an de dissociatiegraad van de temperatuur

volgt uit (lit.10 ) :

In

K

P 41n(4 a )/3+ln(p)/3-ln(3-3 a )-ln(3+ a )/3

=

-4907,213/t+O,336In(t)+4,654

Voor t

=

uitgangstemperatuur kondensor

hieruit

a

= alfix.

tuc1 volgt

liet gas aan de uitgang van de kondensor heeft wat be-treft zwavelwaterstof, zwaveldioxyde, waterdamp, koolstof-dioxyde en stikstof dezelfde samenstelling als het gas aan

de uitgang van de brander; de hoeveelheid zwaveldamp

be-draagt nu:

zwavel (S8) s8 = (1-alfix).perc1.b[1]/S kmol/uur

Eerste line-burner:

De temperatuur aan de ingang van de line-burner werd

gelijkgesteld aan de temperatuur (tuc1) aan de uitgang

van de eerste kondensor.

Het fuel-gas voor de line-burner heeft de volgende

samenstell ing:

waterstof H

2 fg [1] kmol/1 kmol totaal

methaan CH

4 fg [2] kmol/1 kmol totaal

etheen C

2H

₄

fg (3) kmol/1 kmol totaal

ethaan C

2H

6

fg [4] kmol/1 kmol totaal

propeen C

3H6 fg [5] kmol/1 kmol totaal

propaan _C

3HS fg [6] kmol/1 kmol totaal

n-butaan C

4

H10: fg [7] kmol/1 kmol totaal

stikstof N₂ fg [8] kmol/1 kmol totaal

kooldioxyde CO₂ fg[9] kmol/1 kmol totaal

Voor de berekening werden de volgende grootheden ingevoerd:

tir1 tir1 v [1] =

J

c dT v[6J =

/

c dT PH2S PS8 tuc1 tuc1 tir1 tir1 v [2] =

/

c dT v [7]

J

c dT PS02 PS6 tuc1 tuc1 tir1 tuc1 v[3]

/

c dT v[s] =

J

c dT PH20 PC02 tuc1 tli tir1 tuc1 v

[4]

/

c d'r v[9)

J

c dT PN2 PH20 tuc1 tli tir1 tuc1 v [5] =

J

c dT v [10] =

J

c dT PC02 PN2 tuc1 tli

-23-De hoeveelheid zuurstof, nodig voor de verbranding van 1 kmol fuel-gas:

De hoeveelheid water, die vrijkomt bij de verbranding van 1 kmol fuel-gas:

w I2J = fg [1] +2fg [2] +2fg [3] + 3fg [4J +3fg[5] +4fg t6] +5fg [7]

De hoeveelheid water in 1 kmol gas met een relatieve

voch-tigheid betha en een temperatuur tli:

w(3]

=

1/(betha.o,01.exp(-4882 ,1(1/tli-1/373,2)))-1

De hoeveelheid koolstofdioxyde, die vrijkomt bij de

ver-branding van 1 kmol fuel-gas:

Bij de verbranding van 1 kmol fuel-gas komt aan warmte vrij:

.. 1 [5 J

=

1 03 • ( 51 , 1979 • f g [ 1] + 1 91 , 159 • f g [2] + 31 6, 1 95. f g [3] + + 34 1 , 261 • f g [4] +460 , 420 • f g [5J +

4

8 8 , 527 • f g

r

6J +, 635 , 384 • f g [7

J )

Bij de berekening va.n de line-burners werden de verbrandings-warmten van het fuel-g&s bij 25 oe en 1 atm genomen.

De warmtebalans over de line-burner, bij volledige verbran-ding van het fuel-gas met een stoechiometrische hoeveelheid lucht, wordt nu:

h2 s • v [1] +s 02. v [2] +h2 o. v [3] +n2 • v [4] + c02. v [5] + ( ( 1-al fix) • perc 1 • b [ 1J /8) • v [6J +al fix. perc 1 • b [1] • v [7J / 6+n. ( v [5J +v [8J ) • w [4]+n. w [2J • (v [3) +v [9] ) +4,168. n. w [1] • (v [3] +v [9] )

/w

[3J + +3, 7G8.n.w[1]. (v[4] +v [10] )+n.ftda]. (v[4]+v [10] )+n.Vl[5] =0

Hieruit is de benodigde hoeveelheid fuel-gas (n kmol/uur) te berekenen.

Na de verbranding van n = na kmol/uur fuel-gas worden de hoeveelheden van de verschillende komponenten in kmol/uur:

z'tlavel wa terstof a [1] h2s

zwaveldioxyde a [2] s02

water a [3]

=

h2o+n.w[2] +4, 76e.n.w(1] /w[3]

stikstof a [4] n2+3, 76S.n.w (1) +n.fg[S]

koolstofdioxyde a [5] co2+n. w [4]

zwavel (Se) a [6] (1-alfix). pere1 • b [1] /S zwavel (S6)

e[7]

₌

alfix. pere1 • b [1] /6

-25-Reaktor:

In de reaktor kunnen de volgende evenwichtsreakties optreden:

2H S 2 + SO 2 2H20 + 3(1- a)/eSe + 3 a/6S6 ( 3) 2COS + SO ₂ 2C0 2 + 3(1-a)/eS e + 3 a/6s6

,

çM1

CS

2 + S02 ==== 2C02 + 3(1-a)/eS e +

3

a/6S

6

~

CS 2 + 2H2O CO2 + 2H2S (>6) CO 2 + H2S ---- COS + H20 /

(Sr

se ----

4/3

3₆ ( 17)

De evenwichtsreakties

(e,

;14, 15,16) werden bij de bere-kening verwaarloosd.

Bij de berekening van de konversie en de temperatuur aan de uitgang van de reaktor werden de hoeveelheden betrokken op

1 kmol/uur S02 : zwavelwaterstof a (1] / a[2] 2 zwavel dioxyde a (2) / a [2] 1 water a

[3J /

a [2] krl stikstof a [4] / a [2] lrl kooldioxyde a [5)/ a (2) mrl zwavel (Se) a [6]/ a [2] nrl

zwavel (S6) a[7]/ a[2] orl

De gekoppelde warmte- en stofbalans over de reaktor luidt (lit. 9 ): tr1

J

(2c +c krl.c +lrl.c +mrl.c +nrl.c + . PH2S PS02 PH20 PN2 PC02 _Pse tl.r1 orl.c )_OdT+x.H(3)tl=0 (1) Ps6 x- r

De evenwichtsligging van de reaktie ( 3 ) als funktie van

de temperatuur (lit. 4 ):

In K

P = 2In(krl+2x)+3(1- a ).ln(1- a )/8+3(1- a ).ln(3x/8+

+perc1. b(1] .a[2]/8)/8+~ a In(4 a/3)+~ a .lnC3x/8+

+perc1.b[1] ,a[2]/8)+(S/8-1/8 a ).ln(krl+lrl+mrl+ +3-5x/B- 0:.x/8+(1+ 0:/3)(perc1.b[1]/(B.a[2] )))+ -ln(4)-3In(1-x)-(5/B- 1/B).ln(p) = 12886,62/trl+1,87 2 .1n(trl)-O,00274.trl+2,943. 10- 7 • • tr12-24,S34+ a.(-1836,90/trl+O,11S3.ln(trl)+1,779) (2)

De evenwichtsligging van de reaktie (17 ) als funktie van

de temperatuur (lit.10,4):

In K

=

4ln(4 a)/3+ln(p)/3-ln(3-3 0:)-ln(3+ 0:)/3 =

=

-4907,213/tr1 + O,336.ln(tr1) + 4,654 (3)

Met behulp van de vergelijkingen (1), (2), en (3) is het

mogelijk de konversie, de temperatuur en de dissociatie-graad aan de uitgang van de reaktor te berekenen.

De samenstelling van het gas aan de uitgang van de reaktor:

zwavelwaterstof: h2s=2(1-x).a[2] kmol/h

zwaveldioxyde so2=(1-x).a[2) kmol/h

water h2o=a[3] +2. x. a(2) kmol/h

kooldioxyde co2=a[5] kmol/h

stikstof n2 =a[4) kmol/h

zwavel (Ss) sS =(1- a)(3x.a[2] /S+perc1. b (1] /S) kmol/h

zwavel (S6) s6 =4 a.(3 x • a [2J/S+perc1.b[1]/B)/3 knol/h

-2i{-Tweede kondensor:

Evenals bij de eerste kondensor werd de temperatuur

aan de uitgang niet berekend, d~ch werde hiervoor een waa

r-de aangenomen, die gelijkgesteld werd aan de temperatuur

van de ingang van de tweede line-burner.

De zwavelopbrengst in de tweede kondensor bedraagt:

op2 b[1]-op1 -h2s-so2-perc1. b [1]

Hierin zal perc1 nu een andere waarde hebben als bij de

eerste kondensor.

De gassamenstelling aan de uitgang van de tweede

kondensor is : zwavelwaterstof zwaveldioxyde water koolstofdioxyde stikstof zwavel (SS) zwavel (8

6

)

h2s kmol/uur s02 kmol/uur h20 kmol/uur c02 kmol/uur n2 kmol/uur (1-a ). perc1 • b [1] /Skmol/uur

a .

perc1 • b [1]

/6

kmol/uur De waarden voor h2s, s02, h20, c02, en n2 worden nu g

e-substitueerd in a 1 t/m a

5 ,

waarna de berekeningen van

de tweede line-burner en de tweede reaktor volkomen analoog

Derde kondensor en deoister:

Voor de temperatuur tuc1 aan de uitgang van de derde kondensor werd, evenals voor de eerste en tweede kondensor,

een waarde aangenomen.

De uit de derde kondensor afkomstige mist wordt in de demister nagenoeg geheel afgevangen. De hoeveelheid zwa-veldamp, die de demister passeert, wordt berekend uit de dampspanning van de zwavel als funktie van de temperatuur (fig.

7

,lito

5 )

en bedraagt perc .b1 • De zwavel opbrengst in de derde kondensor en demister bedraagt:

s-perc.b1

De gassamenstelling aan de uitgang van de demister is:

zwavelwaterstof h2s kmel/uur

zwaveldiox.fde s02 krool/uur

'vla ter h20 kmol/uur

koolstofdioxyde c02 kmol/uur

stikstof n2 kmel/uur

zwavel (s ) ₈ s8 = p~rc.b1.(1-alfix)/8 kmol/uur zwavel (S6) s6 perc.b1.alfix/6 kmol/uur

-29-Incinerator:

Het gasmengsel wordt in de incinerator opgewarmd tot

een temperatuur tiu door verbranding van fuel-gas van

de-zelfde samenstelling als bij de line-burners gebruikt wordt.

De benodigde lucht en fuel-gas dient opgewarmd te

wor-den vanaf de temperatuur tbl. De nog a~nwezige_ zwavel en

zwavelwaterstof wordt met 100

%

overmaat lucht volgens (

6 )

en ( 1 ) in zwaveldioxyde omgezet. De hierbij vrijkomende warmte dient tevens voor de opwarming.

Voor de berekening van de benodigde hoeveelheid lucht

en fuel-gas werden de volgende grootheden ingevoerd:

tiu tiu y[1] _{. / c dT} y[2]

J

c :dT PS02 PH20 tuc1 tuc1 tiu tiu y(3)

I

_{c dT} y[4] ::: / cp dT P , tuc1 N2 tuc1 02 tiu tuc1 y[5] _{/ c dT} y [6J :::

J

c dT PC02 ,_PH20 tuc1 tbl tuc1 tuc1 y[7] =

J

c dT y[8]

J

c dT PC02 PN2 tbl tbl tuc1 y[9]

/'

c dT P02 tbl

De hoeveelheid zuurstof, nodig voor de verbranding van 1 kmol fuel-gas bedraagt:

De hoeveelheid zuurstof, nodig voor de verbranding van

zwavel en zwavelwaterstof, bedraagt:

t [2]

=

1~h2s+8s8+6s6

De hoeveelheid water in 1 kmol gas met een relatieve

voch-tigheid betha en een temperatuur tbl:

t[3]

=

1/betha.O,01.exp(-4882 ,7(1/tbl-1/373,2))-1

De hoeveelheid water, die vrijkomt bij de verbranding van 1 kmol fuel-gas:

t

[4]

= fg [1] +2fg [2] +2fS [3] +3fg [4J +3fg [5J +4fg [6] +5fg [7]

De hoeveelheid koolstofdioxyde, die vrijkomt bij de

ver-branding van 1 kmol fuel-gas:

De hoeveelheid warmte, die vrijkomt bij de verbranding

van 1 kmol fuel-gas:

t[6]

=

103 • (fg [1J. 57, 7979+fg[2] .191, 759+fg[3]. 316, 195+

+fg[4].341,261+fg[5].460,420+fg[6].488,527+fg[7J.635,384)

De totale warmtebalans over de incinerator wordt:

n.

«

t [4J +4, 76t

[1J

/t [3] )(y[2] +y [6] )+( 3, 76t [1] +fg[8] ) •

• (y [3J +y[8] )+t [5] • (y [5J +y [7] ) ) +y [1] • (so2+h2s+8s8+6s6) +

+y[2] .(h2o+9,53.t[2J/t[3J+h2s)+y[6J.9,53.t[2l/t[3]+ +y[3J • (n2+7, 53. t [2J )+yJs] .7,53 .,t[2] +(y [41+y [9]). t [2] +

+y[5] .c02

=

n.t[6]+s8.L1H +s6.L1H +h2s.L1H

c

-31-Hieruit is de benodigde hoeveelheid fuel-gas ( nc

=

n in kmol/uur ) te berekenen.

Het gas aan de uitgang van de incinerator heeft de vol-gende samenstelling in kmol/uur:

zwaveldioxyde 902 s02+h2s+8.s8+6.s6

water h20

₌

h20+4, 76(n.t C1J +2t[2] )jt[3J+n.t f4J stikstof n2 n2+3, 76(n. t [1] +2t [2] )+n. fg[S]

kool dioxyde c02 co2+n.t(5]

zuurstof 02 1}h2s+8s8+6s6

De koncentratie aan zwaveldioxyde bedraagt in p.p.m.

ppm

=

106.S02/(h20+so2+o2+c02+n2)

De quench-air, nodig om aan de specifikatie spec (in p.p.m.) van'het gehalte aan zHaveldioxyde in het afgas te voldoen, bedraagt in kg/uur :

Totae.l:

De totale recovery (rec) bedraagt in procenten:

ree

=

100.(b1-so2)/b1

De totale zwavelopbrengst bedraagt in kmol/uur:

b1-so2

en is gelijk aan :

op1+op2+s-pere.b1

De totale benodigde hoeveelheid fuel-gas bedraagt:

nd := na+nb+ne kmol/uur

is: nd.(2.fg[1J+16.fg[2]+28.fg[3]+30.fg[4]+4 2 .fg[5]+44.

G FIGUREN EN TABELL~"'N

1966 1967 1968 1969 1970 (est.) (est.) (est.) World Summary

Demand 25.0 26.0 27.8 29.4 31.0

Supply 24.2 25.7 27.7 30.2 35.1

Surplus (deficit) (0.8) (0.3) (0.1) 0.8 4.1

Supply Demand Balance in North America U. S. Demand (4.75%/yr, 1967·70) 9.2 9.3 10.0 10.6 11.0 Production 9.3 9.3 10.1 11.5 14.1 Surplus (deficit) 0.1 0.1 0.9 3.1 Canada Demand 1.2 1.2 1.3 1.4 1.5 Production 2.2 2.7 3.3 3.7 4.2 Surplus (deficit) 1.0 1.5 2.0 2.3 2.7 Mexico Demand 0.2 0.2 0.3 0.4 0.5 Production 1.6 1.8 1.9 2.0 2.2 Surplus (deficit) 1.4 1.6 1.6 1.6 1.7

Surplus (deficit) For North America 2.5 3.1 3.6 4.8 7.5

Source: Arthur D. LiWe, Inc. estimates and various publications.

tab. 1 (lit. 1

)

Mini.

Mfg. FOB

1970 1975 eost ROI~' SG&A ~~

(million Ion€; tons - ( $/Iong ton sulfur sulfur equi.) equivalent)

Frasch Low 3.0 3.0 7 3 4 . 14

Medium 5.2 5.8 11 4 4 19

High 2.8 2.8 15 8 4 27

Sour gas Natural gas 0.8 1.1 8 7 2 17

Refinery 0.8 1.0 15 5 2 22

Smelter gas 0.5 0.7 6 12 2 20

Pyrites 0.4 0.5 20 15 3 38

Gypsum 0.35 0.35 25 10 3 38

Other native 0.1 .1 25 10 3 38

Utility stack gas 0.1 .3 25 18 3 46

Total Production 14.1 15.8 Total Demand';~' 11.0 13.5

* To obtain 15% pre·tax return on fixed investment.

C<O 4.7?%LY~r:9wth 1967·1975.

tab. 2 (lit. 1 )

-33-.30 .20 .JO

9~·OO~~~~~~--~~~~~~~~~~--~--1.0

1.5

2.0

O.S

- - VOl.UME PROCENT

[%] - -..

~ fig. 3 I . O r - : : : : - - - : : : - - - - , ~ 0.9 t= ~ 0.8 Cl: I&. lal 0.7 .... ~ 0.6 Cl: o 0.5 lIJ Cl: ~ 0.4 Cl) lIJ

g:

O.~

....

! 0.2

...

Cl: f 0.1 BOILING POINT OF SULFUR

1

I

I

I . 0.Oo.-...:;;:b:---:2;;-;!OO=---~:----:O-40bO~--;;5~OO:::--~~--~:-O:-";;"""="80!:-::O,.-l TEMPERATURE - oe fig.

4

(lit.

4 )

-34-c ~ a: E ~ E 'ë ~ c .~ .0 0 g t i. cri ei l! -" ~ a: ~r---ï J2 28 24 20 16 12 -lew Cort FrHCh 0 0 2 10 9.3 9

--8

--3 2 1 0 Hi9h Drm.nd

Unc:luding Net Ext'Qru}

lew Otm.nJ - Utility StJICk Gu

$ o U R G A $ Medium Con F"sçh S M--E R l E T F E I R N G E R Y A G S A $ 10 Unclvding Net Exporn) Ht;lt"lCon FUKh 11 12

Millions or long Tons of Sulfur Av.iliblt. by Souroe

fig.

5

(lit. 1 )

13

NON·ACID USES 0.9 - SUlFUR IN RUBBER AND OTHER USES;·

--- S02.nd S03 USES OTHER CHEMICAL

-

I-USES INClUOING

EXPLOSIVES 1.7

~ _{r-CELl. ANO RAYON 0.3} HF 0.3

..--_r-I RON ANO STE E l 0.3

-

r-Ti0 2 0.5

--

-PETROLEUM REFINING 0.8 SUlFURIC

ACID F _ _{-AMMaN. SUlFATE 0.7}

H~4 [ f -8.4 R T I L 4.5 I _{.-I-PHOSPHATES 3.8} Z E R TOTAL fig.

6

(lit. 1 )

-35-p G Y P S U M y & R I N T A E T S I V E 14 Ptfcent of Tout 10 19 3 3 3 5 9 ~ 100 PTOfit 15'""Pu·llx ~e'urn on Fill:t<1 'm.Mtrr-.rnt Cort. Mfg. Cost, DeDr .• 5.a11!\ 'nd Admirl. Ex~SoeS

(

.

! f r I

I

[. f

r

t

"

r

r

~

i

, , , I I:II ~ ~~ , , , I ' 9 •• ,c .. · .• ·_".·.j~ c'.,·i"l-· .cc:>;: _. >o,i··~c··L". j. • . ~'";.'. " C ' ". c-··:,jC_êc c~c ,:.'.

:~

~0~~S~5~i~~'~~

"'S

l~~1;~'~?~~~~~!2~

:: .-'~' ... .. ·oL· - :_~ - : __ 5 •• ~ .~

=

~~ ~-; __ ~~~~~.C~ __ =-~~-~~ . -3.

l

2_"

i' i '

1

'

1

1

'

'

1'11

'

1

'

1

'

I

,

i~i

'

ililiiiiiiii'i

'

i

I

r

~

200 -' . 210 220 230 240

J

L~

___

~~_~.

~~~~

~em:::.u.:: :~~

~

.

~~".~.~~

.

,J

Fig. 10. Vapor pressure of monoclinie sulphur. (From K. Neumann, Zei+' Phys. Chem .• A 171, 416-420 (1934).)

fig.

1

(lit.

₅

)

36

-\ /

L.

_RERD

-

...

/i

,

r

""

+

1

...

..

"'"

..

_COURt:

"

_~

""

~

BUR.NEI<..

....

HODFDPROGRRHJfR ti,. r

_...

r INelA/ERf/TOR. ~ ~ ~

....

SCHEMA RLGOl-6Q

PEOG/ÇAHI1R

fig.8

-31-f

,30 ,20 ,10 91.{.00 L L _ & _ ' _ L L L . . . . L L

-2BO

290

300

310

9

5.00

f

9~'00 ,80

,60

, 'iO .20

- - tbl

[OK]

•

fig.

9

DE INVLOED WW OF rEHPI:RflTVUR Vluv DE

VO~D/NG OP DE RECOVERY.

93.

00 ~-..L-_'"""'---L_-'----t._--'---l ___ ~ _ _ _ . 280 300 320

340

360

- - fhi

[OK] - -...

~ fig 10

-38-I ~ \0 I ~ 1 2,. '--

.,..--.

V///..I: 1 j~ j

!>~

I TEE/:) I 2 .. ,

r)~··

.

L1?oC

,

t. I

R

[X

F.G

0

...

-3 lf

L

5 6

:

I

,

i

c)(

,

~

/

/

;J/,J

r

c,:' ( ;..

"

.,,:f{\io' ij/,

I

r

-

-

/

I

13,)" I

to

oe.

..

_..

I

~X'i,i'"

',' ~

'I

..

_/1'1

lt

"

''I' .. ( ~ "T J

"

'i(,y ! (~

...

t- i -8 9 20 " - - U

,

_S

'I. '

S

"'.

'rS

9°

,'1 \ 'r$ \."

~

J1.~l;

.

('ow \...·..,.~.l· .,;- .,t .J

.•

/ .

ti

fig. 11

-

_.

~,. L. 'F.G

-.... --.

_{L -}

₇

,

.'l

--

~ ,': . . f. "" --,. -'\ •

c..

'F.G . I

-t

KNoeI( OUT Dt<UH.

2

8UR.Ne:R.

3

Wr:lSTF HEAT 801LéR. l.f l/tIE 8UR.NER..

5

Ie Rf3~kTOR.

6

:J.E' R,éF/KrOf! .

t

LJNF 8URNER..

8

IF CONDéIVSDR. .

9

2~CDNDEN~OR.. JO .31:' CDNDENSOR..

11

DEN/STeR. /2 INCIIVERfJToR.

tn

~ilht

"\ I('e~~'" ~\ .~t:." '. L. ,Ir 1.,.

I

.

ST 12 _{" Gil}

10 z 0 Vi _! " ~!~, o I u i-. 70 65 60 55 5!,}",OO CATAL'I'TIC RfGION 5CO TEMFERATIJR[ _.;:

fig. 5. Theoretica! equi!iorium conversion cf hydrogen sul!,d& 10 vopor sullur by s&lective·

oxidation with Ihe stoichiometrie oir according to Ihe over-all equatian 2H,S -i-0, ~ 2H,.o

+

2/e5,_ 100 Rt;COVE RY [%]

9

5 I I I ,

"

fig. 12

,

,

,

,

,

, ...

-

....

-.

(lit.

4 )

t -10 0 +10 EXCE"SS AIR [%] 1 ReSOLUTE STDICJ..IIOHETR.Y, UPPl?R LlNIT OF CONVER.$ION.

2 nCTUR~.

fig. 14 (lit. 3 )

-4

0-r'- ----.. ~ ... ~ .... -.' . I · ' rr- ... ·~ .. . . __ .--... ...... , " ' _.; ....

'-

.

- .

'

-p ..

_'..

,

, ,... v ~ ... ,... S ... ~ .§ cs l:( ;; l' ~. I I I I : 1 , r1~~~--~,~.~ I~~~~--~, 1 I I T~',,: --;-~ 1-1--+; -+, -t,~·-"--..L.-H-++--0;-- -;-;_

-t

+

-++--:--: -::---1-+:"""', -;--~ ~~--+ __ 4 __ ..L.~ __ ~~-+~..L.4 ___ ~1_ I-I....;...:...:...jl-~+-H--+...:....+-~ ~....J.~_ 1 "

r+

.:;: 1 , _, ' I-t-

~.---14

,

-1:....l.

1

---'---

.

-+Ll...!...l.J-+; --';-+--

,~

0. 00501--1--'

I'

-~

,:-

.:

_--:I--

I-~

-_-

..:.~

'

_""""'-_';:'..:.""

'

~t=f.-../-::"":.+_-..:

,,:,

,

::j;:

.:j

I---~_+-L~ , I ~ ~ r ~~~~,-+.~.-~---,--~'_I~-+--~++-~":'l 1--4-~---1I--4-~~~~_I~- 1---t1~1~~~+-..:. . ..:.I~+--+~'-+ __ ,_ I I , .1 I ,of or i I I ' I 0.0000 L.l."":"".l--2':"'5 O-:"-...L.., O.LO...L...~---'-:--,.L5 0,-'---'1'--... 1 I I • j I I Tcmpcrature. OF. , I i I i , ; I I I , i I ' . ..!... I I I I I 70~~~~~!~;~-+~~--~~~4---~-+---~--~1-~~ i I I 60 50 40 30 20 10 I o [7. 300 ( i I I I I 1 , _i I I I , I \ . I , I \ I , , _! I , I i\: I , , , I ; . \ I _I I I , , _I I I I I /400

:in'(

I ·T ! , I I I , , , I 1 I I I 1 I , I , I I I

,

.

I I .", I , , -...::::' I I 450 500 550 Tcmperature. 'F. I ; i , 1 --..; 600 : ! , i I I I , ! I I , I , , I , , , I I I , I , I , I i 1 , i 1 I 650 I I , , I I • I i , I ' I 1 I I I I I i I I i 1 ~ ________ . , __ ~ ____ ... __ " .. , ... ~ .. --1 Fig. 5. Viscosity of liquid sulphur. High temperature range.

(From R. F. Bacon and R. Far:elli, J. Am. Chem. Soc., 65,

639-648. (1943).)

fig.13 (lit. 5 )

-41-'

.

...

" I" 1 H ALGOL-PROGRM~ ==:11 ,

'hngin'

'procedure'

burmr;

'code';

• 0 roe e du re' coL

i

re; , c 0

de'

j

'procedure'

incinerator;

'code';

'real' p, tbi, tbL, betha,

pr,

tiPI, tUCl, tLi,perCI, hl,

-

h 2S,

soz,h zo, n2,coz,S,

x,

tb23. trl3, aLfa,na, nb, ne, nd,

OPl,

opz.

OP3,

n.perc.spcc,rec,02,pp~,pi,tiu;

, ar r

Cly • -

b [ 1 : 4

J.

f

9 [

1:

9 ] ;

rdam:

read{p,tbi,tbL,betha,b,fg);

orint ( • ,

cLaus-oLant") ;

nLcr(s);

print(··

print("

v

as k

0 ( 2, 1,

b eth a) j

print("

proc

es

b

erekenin

9

input

data"); nLcr(3);

reL. vochtigheid Lucht:

tt);

space(l

0);

sp ac

e (

3);

wr

i

te ( , ,

p

r

oc en

t f f ) ;

n

Lc

r

(1 ) ;

te-nperatuur Lucht

: " ) ;

-

space(lo);

v

as

k 0 ( 3, 1 , t b L ) ;

print(t.

sp ac e ( 2 ) ;

tenperatuur voeding

'

wr i

te ( , ,

9

rad en

:tI)jSpace(lo);

k " ) ;

n L

c

r ( 1 ) j

v

as

k 0 ( 3, 1 , t

bi) ; space( 2) j

wri

tee'

t

graden

~\")

j nLcr( 2);

print("

steLLi n

r

l

print("

wri te(' ,

samonsteLLing voeding

krooL/uur

samen

fu eL

gas

moL fraktie") ; nLcr(

1 ) ;

-

h2s:

"~);

VaskO(3,4,b[1]);

hz : ti);

v as

k 0 ( 1 , 4 , f fJ [ 1 ]) ; p

r

in

t ( , ,

C3 : ' , ) j S P

ac

e (

5 ) ;

v

as

k

0 ( 3,

-4,

b [ 2 ]); s

p

ac e ( 2

0 )

j

wr

i

t 0 ( ,

'c

h -4 : ' , ) ;

s p ac e (

5) j

v

as

k 0 ( 1, 4, f 9 [ 2

J)

j

n

L

c r (

1 ) ;

Sp ac e ( 15 );

wr

i te ( , ,

C 4 : • , ) ;

sp ac e ( 5 );

v

as

K

0 (

3, 4,

b [ 3

J) ;

space(zo);

write("c2h4

: " ) ; soace(5); vasl<ü\1,4,fg[,]);

nLcr(l); space(15); write(ttc02 : " ) ; space(s);

vasko(3,~,b[4]);

space(2o);

write("czh6

: " ) ; space(s);

v as

k

0 ( 1.

4,

f

Cl [ 4 J);

n

L

c r (

1 ) ;

sp

ac

e ( 5 5 )

j

wr

i

te ( • 'e 3h

6 : • , ) ;

sp

ac e (

5) ;

vasko(1~4,fg[5));

nLcr(I); space(55); write(

f t

_c3hs

_:It);

sp ac e ( 5 );

v

as

k

0

(1 4 ,

f

9

[6

J);

n

L c

I" ( 1) ; .

s pace ( 5 5 ) ;

wr

i

t

e ( • , C 4 h 1 0 : t , ) ;

sp

ac e ( 5) ; v

as

k 0 ( 1, 4, f 9 [ 7

J) ;

-

n L

cr (

1 ) ;

space(5s); write("n2

: " ) ; space(5); vasko(l,4,fg[s]);

nLcr(l)j space(ss)j

write("co2

: " ) ; space(s);

vasko( 1,4, fg[ 9J); npag;

111j~.~

hu rn er (p, tb i, tbL, b, b eth a, h 2S, S 02, h 20, n 2, c 02, x. tb 2

3) ;

Pf

nLcr(3);

-

-bI:

=

b

[1];

tL

i:

=

t bL ;

-r

P.

ad (

n

r, tir 1

f

tu Cl,

P

e rc

1) i

oOl:=b[1 J-h2S-S02-percl.b[I]; space(lO);

'end' ;

-43-write(" eerste

kond~

nsor"); nler( 3); spaee(1 0);

writ!:(t

·

'zwavelop

b

rengst

: " ) ; spaee(lo); vasko(3,3,Opl);

sp

ac

e ( 1);

wr

i t e ( , • I<m

0

L

l

u

u r' t ) ; n L c r (

2) ;

sp

ac

e ( 10) ;

'

wri

te (t 'u i t

~an

gs

te

:np

eratuu r :

t, );

sp ae e (

1 0) ;

vasko(3,1,tuel); spaee(3);

write("

graden kit);

npag;

nler(3)

'

; sp

a

ee(lo);

v1

rite("eerste reaktor"); nLer(3);

spaee(l~);

w

riteCtleerste Lineburner

l f );

nLer(z);

spaee(lo); './rite("in

ga

n

g

stanperatuur : " ) ; spaee(10);

vasko(3,1,tuCl); spaee(3);

vl

rite("graden

l<tI);

nLer(l)

'

; space(lo); 'Nr

i

te(·'uitgangstenpero.tuur : " ) ;

soa e e ( I

0 ) ;

v

as

k

0 (

3, 1, tir

1 );

S

P

ae

e (

:3); ~'I

rit e ( , ,

9

rad en

k");

nLer(2);

~~eoLire(pr,

tirl, tue

1,

f

g

,

~Li,

betha,. h

zs!~o:,

h

zo,

•

nz,eoz,s,perel.bl, trI 3,

x,

alfa. n). na.-n, nler( 4),

spaee(lO);

write("tweede

kondensor"); nLep();

read(pr', tirl,tucl,perel);

op

2 :

=

b rIJ

-op 1

-h

2S

-S

02

-p

ere

1 • b [ 1 ] ;

spaee(lO); wr'iteC"zVlaveLoobren

g

st

: " ) ; spaee(lo);

vasko(3,3,OPZ); spaee(l); write ("kr.lOL/uur");nLerC

·

l);

sp ac e (

1 0 ) ;

wr i

t e ( , • u i t

~J

an gs

t

e

m

p

er a tu u r :

I , ) ;

s p ae e (

1 :J ) ;

vaskO(3,I,tuCl); space(3);

write("graden

kt');

npag; nLer(3); spaee(lo);

v

/rite("tweede reaktor");

nLer(s); space(lo);

v/riteCtttweede

Lineburner''');

n L e r ( 3 ); sp ae e (

1 0 ) ;

wr i t e ( , , in

9

an gs

t em p er

a

tu

u

r

:

I , ) ;

space(lo); v

a

sko(3,l,tuel); sp

a

ce(3);

v

/rite("graden

kl

' ) ;

nLer(!)

·

; space(lo);

w

rite("uit

ga

ngstanperatuur : " ) ;

spa c e (

1 0 ) ; v

ö.S

k

0 ( 3 ti, tip

1 ) ;

S

P

ae e ( 3) ;

wr

i t

e

(

t t 9

r

Cl d

en

1<' t ) ;

n

L

cr

( z ) ;

eOLire(pr .. tirl, tUC1, f

g

, tLi, betl1a,I1zs,SOZ,

h

20,

n 2, CO 2.

S, P

ere 1 , bi. trI :3~

x,

aL fa,

n); nb: =n ;

rea d (

0

i, t

i

u,

p ere,

tu Cl.

sp

ee );

n

L C

r ( 4 ) ;

ineinerator(tuel,pi.tiu,perc.bl,fg,hzs,betha,

tb

L • s 0 Z,

11

zo. n. rec. 0 2,

P

ere

1 ,

n

2, C

02, S ) ;

ne:=n; nLer(3); pom:=sozo,o6/(hzo+02+S0Z+

C02+n2); soaee(lO); 'ölritc("quench air

space(lo) ;dri

k

o(

3.

(32+3.

76tdS)/4.

76.S02-, 6.

(1

/spee-l

lop

m

»;

write(" kg/uur");

nLcr(

1);

• t t ) •

.

.

nd:=na+nb+nc; space(lo);

v

/riteCttfueL

g

as totaaL

: " ) ;

spaee( 10); vas:

(

o(

3t

2, nd~(

zaf

g

[ 1 ]+16 .f

g

[ 2 J+2s.f

g

[ 3]+ 30.

f 9 [ 4

J+

4 Z. f ][ 5 ]+ 4 4. f g[ 6 ]+ S 8a f 9 [ 7 ]+ Z a. f 9 [ 8

J

+ 4 4. f g[ 9]) ) ;

==a3,

'orocedure' burner(p,tbi,tbL,b,betha,hzs,soz,hzo,nz, C02,x,tbz3);

'real' p,tbi,tbL,betha,hls,soz,hzo,l1Z,COl,x,tbz3i 'array' bi

'boqin"reaL' ab,bb,cb, db, eb, fb, gb,hb, kb, Lb,nb,pb,qb, Zbl,zbz, Zb3, tb 1 , t bn, s, zb 4 1 , Zb 4 2, zb ~ 3, tb Z 1, tb z 2, tb z 3 ;

'orocedure' deLta (tbz. zb4);

'real' tbz, zb 4; -

-'beçJin' -x:=(eb*(tbz-tbl)+fb/z*(tbz 'powepl - l-tbl -'power' Z)+gb/3. ( tb l i p 0 \'/ e r t - 3 -tb 1 - 'p 0 VI er ' 3 ) -h b .. (l / t b l -1 / t bI) ) / ( -1 3 5 8. z-5.965*tb2+o.00544*tb2*tb2-1.17~)-6*tb2 'power' 3); 'end' j ' i f ' x 'qreater-'O.999'then' X:=O.999; • i f' x 'L ess r o. (101 • th en t -X: = O. 0 {) 1 ; zb z : =-5 077 /tb 2 + 3. 00 2* Ln (tb 2) -0. I) C 273 * tb 2+ Z.9 "10 -7" tb z * tb 2-1 1 .3 4 ; test( zbz); Zb3:=z*Ln(kb+l+4*pb+5*qb+z*X)+3/zotn(1.s.x)+o.50Ln(p)-Ln(4)-3*Ln(1-x)-O. s*Ln(<I+o. 50 x+kb+Lb+lilb+7*pb+9*qb);

tes t ( zb 3) ; zb 4: = zb 3 - Zb z ; t es t ( zb 4) ; pb:=3*b[2]/b[l]; qb:=3 0 b[3]/b[1]; - . kb:· = 3 0 ( b [ 1 ] + b [ 2 ] + b [ 3 ] + b [ 4 ] ) / ( p

I

ex p ( - 4 6 6 l . ·7 a ( 1 / t b i-I / 3 7 3» -1 ) / b[lJ+4.7660(3/2+Sopb+6.Soqb)/(1/(bethaaO.Ol.exp(-~88Z.7a ( I/tb L -1 /37 3) ) ) -1 ) ; -Lb:=3. 76 SO(3/2+50 pb+6. 50 qb); mb:=b[ 4 ]/b[ 1 ]03;

a

b:

=

3 4. 0 0 + k

tJ.

8. 2 2 + L b - 6. S 0 +11 b

*

1 0 • 3 4

+

p b a 46 .11 +q b .. 6 1 .Z8;· -bb: = O. 0 1 1 2 +kb. O. 00 0 1 5 + L b * O. 00 1 +m b* O. 0 0 Z 7 4+ 0 IJ s 3* pb + o. 06 7 oq b ; Cb:=l. 34,0-6t1>kb-l. 47 8"..,-s*pb-18. 7,o-6 a qb; db : = - 2 8 1 5 5 - 1 9 5 S 0 0 ":Tl b - 9 3 8 -S 0 0 pb -1 2 zoo s 0 q b ; eb:= 3 O. 3 z+!<b. 8. 2 z+Lb-6. 5 C fmb*1 O. 3 4+pb*6 3.9 o+qb* 8 2. 46 ; f b :

=

O. 0 1 2 6 5 + I, b.o. 0 0 0 1 5 + Lb. o. 0

a

1 ~1 bo.

a •

0 0 Z 7 4 + pb. o. 0 0 8 8 2 +q b 00.01171; gb:=5. 1'0-7+1. 34,o-6*I<b+pb.5. 36 to-6+qb.6 )/0-6; hb:=-(195 50o*mb+s86 sOoopb+7 3 20oooqb) j - . .

zbl :=ab.tbi+bbotbi.tbi/2+cb.tbi 'po'./er' 3/,-db/tbi+ll 21 a 4.6

+ 4 9 36 46 .p b+6 4 1671 oq b ;

print(" hoofdbrandert f _{) ; nLcr(3);}

p rin t ( • • d ru kin b r an der· : ' , ); s p ac e ( 1 0 ) ;

vasko(l.Z,p); space(6); rtrite("atat . ) ; nLcr(l);

rrint(" t8lilperatuur voeding : " ) ; sn8ce(lo); vasko( 3, I, tbi); space( 5); write(ft graden ~~"); nLcr( 2);

cl

s

t -:t : " , ' f • 1"\.-',' .:. I , I I '. I . . . . \ 'I ,, __ 1 . _ • • '. c,' r ..., , -' , c , .' ,." ,l ... ~.' .. 4... '-' \ t t

-45-!<; OL/LI') t' 112 C 3 (', """ '1"' '.;' 1) 1,..1 ... '. J ' \ 1..' J :-':;:}CC;( 1 :;) ,t ., .. -_ - ", ... ~

..

-/ ...

43;

, )

.

9l' .LI~ ) ' \ , , , _ • -.J ; ;; Lc i) ( =~

,

"

.~~\,' :;!J t";:.j '~~~ ::.~~(;

,

) ~ 2 0 ;, Le r (1 ) ;

.

'",", ::.Iv .t_ "'" l '-, ,~! ~" • I v • .l..!J , CG2 r_~r:L~:t(t' :""l'J'Jsr' c , -,.,~ ' l (1 ' , . J ,.J '-..4'.J ,- '.' I \ /

..

P C02 _{; •} S _..y _{,,) G} nrC_~ I _{, u} \--'1 _{• '-}-'-" 1 _{3 .. /} ,. _{... ,.} "L'~" _{(...i. J C.I..} j"_{I) ;} \ I ,. e Ld.' ç; t', t ~ r~ 1 • 'LU Cl . 't L~ .• J, 2S. S 02. f'l 2 O. n 2. CO Z. S. p

c

r'

c

1 , bIt t r I 3. x. aL f'l, ~Jc:trlLi, ;) ; t art'~ty ' '!-iJ;

, rea L' él r I ,

r;

r 1 , crI .. rJ t" 1 • n rIt 0'" 1 • 1\ rIt L rIt ~,~ rIt aLf J.;-). aL -;-i x t

él L fin. z r' ) • Z t~ 2 t Z f' 3, Z '"' .:" ZIt S 6. S 6 • t t"" 1 1 • t rI 2 t Z r 4 1. Z r 4 2 t Z r 4 3 ; , ê! r r' i:lj' ;::, [ 1 : 7

J

t ',/ [ 1 : 5

J

t V [ 1 : 1 û

J

;

' nrocerlur'e' rJcLtar('ct·l, zr4);

, r

e

Cl L t t r1 • Z t' 4 ;

'

beqi

n'

arnho'D

:

Zrl:=-49G7. 213ltrl+O. 336oLn(t rl)+4.654j

a

Lf3.:=o.5;

t

es

t

(

a L'f (1) ; aLfan:=ilLfa-(4/3.Ln(4.aLfa)+1/3.Ln( pr)-Ln(3-3.~LfQ)-1/3§ fLOn: deLft: L n ( 3 + Cl L f <1 ) - n' 1 ) / ( 4/ ( 3" ::J. L fa) + 1 / ( 1 - él V: a ) -) / ( 9 + 3 ~ él L 'f a) ) ; ' i f' äbs(LlL-;-a- aLfan) fL oss ' 10-4 'tl1e~lr ' go

t

of

'fLop;

al

f a:

=

él!,-;' U;l ;

l i f t

a

lf

o.

tL

e

ss'

0

I

tflen

t élLfa.:=O. l ; t :~o

to'

arn

;le:l

j

x:==(-ar'll) (t t'l - t::'r'l)-td"'l/Z$.(t t·l '(io'.'!et~t 2-t5,r1 f:)O',!é)r'

2)-C r' 1 / 3" ( t t" 1 t~) 0 ij G '" , 3 - 't Ü' 1 t DO' I eH'" , 3) + cl r 1 lil ( 1 /'~ r 1 -1 / t i, r 1 ) ) / ( - 2 5 6 0 4 + 36 5 G (1 Cl. L f an + (). 2 3 5 (1 Cl. L -;-31'1" t r I - o. 0 0 5 4S" t '" 1 • iJ 0 './ e t~ , 2 + 1 • 17_K) -6 (1

t

j' 1 '~, 0 \ler f 3); ' i f '

X

' w<:at er t 0.999 ' tIHm' )(:=0.999; lift

>

:

cLess' 0.001 ' th

en

t X:=O.O Ol ; Z r 2 :

=

2 .. L n ( k r' 1 + 2 .. X ) + 3" ( 1 - éJ. L f an) / SI) L n ( 1 - Cl L f 2.n ) + 3 <> ( 1 - aLf a:1 ) /

8

I!f L n (

3 /

a. x + iJ

G

rc

1 .1 /

3

"

fJ

I/a [ 2 ] )

+ û •

so

a

L -;-

an

i!i L

n ( 4 /3

§ aLf

a

n )

+0. 5oéJ.L'fûl~.Ln( 3/8.)<+perC1~1/8(~b 1/a[ 2 ])+( 5/S-1/s.aLfan) • L n ( f< r 1

+

L r 1 +~:;, t' 1 + 3 - 5 /

30::

-1 / éj"

aL

f

;1"1 e X + ( 1 + 1 / 3 I) éÜ f Q, 1 ) ~ ( O. 1 / 8 \'f b I /a [ 2

J

) ) -

L n ( 4 ) - 3 .. L n ( 1 -x ) - ( 5 / Ó -1 / 8 (; u

L-

f

~,li1 ) I) L i 1 ( P (I) ;

t

8S t ( zr 2) ; zr 3 :

=

1 2 8 B 6 .. 6 2

/t

r~ 1 + 1. 87 2(1., L 1'1 ( t rI) - 0 • 0 ij 27 4. t r 1 + 2 .. 9 -4 3 I() -7 I) tr'l tpv\ïGt'" 2-24, ~34+éJ.LféJ.nc:.(-la36. ~ o/trl+O J183~Ll (trl)+ 1.719); t os

t

(

zr' 3 ) ; zr 4: =Zj"3-Zr 2; , en d t ; Z 1 ~ ~ - 4 9 fJ 7 • 2 I 3/ t i r 1

+

O. 3 3 6 ~ l11 (

t

ir' 1 ) + 4 .6 5 .j ; aLTJ.n:=-=o. 5; tos~(~Lf~.~~) ;

á L f i~~

:

=:l L': :~:~ - ( 4 / 3 ~. U-: ( 4~· û L f in) + 1 / 39-L i ( IJ r ) - V~ ( 3 - 3 ~ :} L'f l11 )

-1 / 3 (I; Lil ( 3 + ~j

L":-

':;'71 ) - 21 ) / ( 4 / ( 3 (ot Zi L -; ~~ i 1 ) + 1 ;' ( 1 - 2. L f ~:'1 ) -1 / ( 9 + 3 n

aVf-::":.) ) ;

t ." .(: • . ~.,

co

'

.

L.!O ~., - ') L.t: ~ " ) ' I 1 ' ) ( " " t _ . t -!o', 'y.~ t ' .:1"'" 7 Q. t ',: L :; ,:: ','''' J'

....;.. 1 ...11; v\. ,'.J. • • _11 U I _~" v .... · 0 0 JC- ~ \..1' · .... 11 - J V ' " _ • _ _

él V~ i Ti :

=

~ L'~' ~L:: ;

t ~ f ' cl V'

i :

-

,

" V}

s

s·

û t t:~

'

'}

Ti • û. Vi"

i;

1 :

= ;;

.

Cl ; t ij J t :) t cl.:; t,-j--t j

\Ir

1 J::::7 .2 :~ (-':ii·l-:LiCl )+:l .001 3::~(tirl · pc','I':;r' 2-tUCl

-41-v [

2

J : ::::

7. 7 0 Ir> (t i r' 1 - tu Cl) + O. 0 0 2 6 S. (t ir 1 •

n

°

VI 8

r'

2 -tu C 1 t p 0\'18

r •

2)-8.3,()-7/3*Ctirl 'po',/ert _{, )-tUCl 'po\'1cr' 3);}

v[ 3]:=8. 22"" (tir'l-tUCI )+0 J)OOC7 s-(tirl 'pO\'/erl _2-tUCl

'oower' 2)+1.34₁₀-6 /3.(tirl 'po\'/Cr' )-tUCl 'power');

~~: J:=6. set,! (tirl-tucl )+0.000 5.(tirl 'power' Z-tUCl 'power'

v( 5 J:=1 O. 3 H(tirl-tucl )+0.00 137.(tirl 'power' 2-tUCl

'por/er' 2)+195500fJ (1/tirl-I/-CUCl);

V[6]:=25.00*(tirl-tuCl)+O.001776.(tirl 'power' 2-tUCl

• P 0 1,'1 er' 2);

V[7]:=19.Z5*(tirl-tuCl)+O.OOI33z.(tirl 'power' z-tUCl , po \'1er ' 2);

v[e]:= 10. 34*(tUCl-tLi)+o.00137.(tUCl 'power' z-tLi 'power'z

) + 1 95 S 0 Co ( I/tU Cl -I /t t., i) ;

v [

9 ] :

=

8. 2 2 & ( tu C 1 - t

LU

+

o.

0 0 0 0 7 5. (

tu

Cl' P 0 \'I

e r

t

z

-tL

i '

p 0

we

r

I

z)

+1. 34_,0-6.(tucI'power' 3-tLit'poi'/er' 3)/3; ,

V rIO]: = 6. 50 Cl (tuc 1 -t Li) + 0.000 50 (tue 1 • power' z-t L i' power' 2) ;

w[ 1 ]:=3 .5":'7)[ 1 J+2~~j[ 2

J+

_{3 ..} "fJr 3 J+1 é>,:--;':Jf -i =-H. 5.fèi[ s]+ sofg[ 6 J+6. 50 fg(7 J;

wr

2 ]:=fg[ 1 J+2 o fg[ 2 J+zofg[ 3 ]+3"-fg[ 4]+ 3.fg[ 5 J+ 4.fg[ 6 ]+ s.fD[7]; , , -w[ 3]:=I/(betl,ao O.Ol*exp(-488Z*(1/tLi-1/37 3»)-1; w[ 4 ] : = f 9 [ 2

J

+ 2* f g[ 3 ] + 2* f g[ 4 ]+ 3*' f g[ 5 J+ 3. f g[ 6 ]+ 4. f g[ 7 ] +f 9 [ 9 ] ;

wr

5 J:=IO 3o (fg[ 1 ].57.797 9+fg[ Z J.l 91.7 5 ~+ffj~3 ].316.1 95+-fg[ 4

J.

3 4 1. 2 6 1 + f 0 [ 5

J ....

46 O. 4 2 e + f 9 [6

J.

4 8 8. 5 27 +T 9 [ 7 J

*

6 3 5. 3 8 4) ; n:

=

C h 2 S ..

v [

1 ] + S 0 Z •

v [

2 ] +h zoo v [ 3 ] + n 2. V [ 4 ) + C 0 20 V [ 5 ]+ V [6

J.

(1-aLfix)*oerClol/aobl+4/30 aLfixoperClol/8.bl*V[7])/(

wr 5 J-v[ 3 J*v/[ 2 ]-4. 16 o\l[ 1 ]/i/[ 3 ]~vr 3 ]-3.76 O\'v'[ 1 ].V[

4]-fÇl[

8 Jov[" J-v[ 5

J*w[

4

J-

v

/[

4 ].v[ 8 ]-tw[

z

J+4.7 6 aow[ 1 ]/I'/[ 3 ])ov[ 9 J

- ( f (j [ 8 ] + 3. 7 6 8 - w[ 1 ] ) •

v [

1 0 ]) ; a(1]:=h2S; a[2J:=s02;

ar

3 ] : =11 Z û+ n. w[ 2 J+ 4. 16 8. w[ 1 ]

/

'i

JL

3 ]. n ; a[ 4 ):=n2+3. 76 aonow[ 1 ]+nofg[ 8]; a[ 5 ]: =eo 2+no1,'/[ 4 ]; a[6]:=(1-aLfix)operCl.1/a.bl; a[7]:=4/30 aLfix-perClel/8 o bl; test ( a) ;

nLcr(I); space(10); v/riteC"benodigde Lucht

Sp ae e C 1 0 ) ;

• • I ) •

.

,

v

ns l< 0 ( 3. 2. (3 Z + 3. 76 • 28 ). n. \'I[ 1 J); sp ac e( 2) ;

write(l'k~J/uur"); nLcr(l); space(lo); write("benodigde

f LJ

e

L g

as

:

I , ) ;

s

p

ac e (

1 0 ) ;

v

as k

° (

3. 3, ( ze. f g [ 1 ]+ 16" f ~[ 2] + z 8 -f (j[ 3 ]+ 30" f g[ 4 J+ 42.

fg[ 5]+ 4 uf0[ 6 J+5 8*-fg[ 7 ]+2s o fg[ a J+4 ".fg[ 9 ]).n) ;

space(l); .. ,rite('tl<g/uur"); nLer(4);

space( 1 (l); write(" reaktor"); nLer( 3);

sp ac e ( 1 0 ); ~'1r i te ( t , dru k : ' t ) ; sp ae e ( 1 0 ) ;

v

as

k 0 (1 2 t p

r) ;

spac:e(4); write(ttatat l _{) ; nLcr(l);}

space(lo); i/rite("ingangste-nperatuur : " ) ; space(lo)i

fLup: 'begin' 'end' ; fLap:

'end' ;

nLcr(z); Vlrtte("k.îl0L/uur hzs S02 hzo C02

nz S8 S6 stot"); nLcr(2);

wr i te ( , , i nv 0 er ' • ); V llS 1\0 ( 3, 3 a[ 1 J); v as k 0 ( 3, 3 a[ 2 ]) ;

vas!<o( 3, 3, a[ 3]); vasko( 3, 3, a[ 5]); vasko( 3,3, a[

4]);

vasko( 3. 3,

ar

6 J); vaskO( 3, 3, a[ 7 ] ) ; vasko( 3,3,6 .a[ 7 ]+80 a[ 6

J) ;

nrl:=a[6 ]/él[

d;

Orl :=a[ 7 ]/a[ 2];

krl :=a[ 3

J

/a[

2]; Lrl :=a[ 'I J/a[ 2]; :nrl:=a[5]/a[Z]i ar 1 : = Z 2. 1 C

+

3. Z Z -!( rl +6. S 0 ~ Lr 1

+

1

o.

3 4

-m

r 1 + Z 5. 0 con r 1 + 19.zS-0rl; b rl : : 0 o. () 1 2 5+ o. 000 IS" k r 1 +0 • Ol. Lr 1 + o. 0 0 21 ot -171 r 1 + o. 0 0 3 5 5. nrl +0. (JO 266 4.0rl; crI :=-8 .31()-7+1. 3 4_fQ-6.krl; . drl:=-1955000:nrl; trll:=tirl; c1eLtar(trll,Zr41); tr12:=6s0; deLtar(tr12,zr42); tr13:=(tr'11-trl z)/(zr4 2-zr41 )oZr41+trll jdeLtar(tr13. zr43); ' i f ' abs(zr43)tLess'o.oo5'then' 'goto' fLap;

'if' zr43 'Less' 0 'then"

trlz:=tr13; zr42:=zr43j

'go

tot fLup; trll:=trl)iZr41:=zr43; 'go tot fLup;

h2S:=2.(I-:~)"',lr z]; S02:~(1-x)-a[

z];

hzo:=a[ 3]+z.x-a[ z]; c02:=a[s]; n2:=a[4J; . s 8: = ( 1 -q L fan) 0 ( 3/8" x. a[ 2 ]+ p erc 1-*1 /8" b 1- ) ; S 6 :

=

4 / 3" aLf an iJ-( 3 / 8 .. X Co a [ 2 ] + P G r'c 1 • 1 / 80 b 1 ) ; S:=8~S8+6·S6 ; p rin

t ( "

u

i t,/ 0 e r " );

v

as k 0 ( 3. :3, h 2S ); Ij

as

1< 0 ( 3. 3. S 0 2 ) ; v ?S k 0 ( 3. 3, f, 2 {; ~; v a~, k 0 ( 3, 3, C 0 ~ ) ; ~~ us 1-( 0 ( 3. 3 ~ 1: 2 ) ;

véiskO(3. 3,sa), vaSr-o(3,3.S6). -Va""kO(3. 3,S)~

nLcr(2}; rJr;int(" uitgangste-:Jperatuur :'1);

sp ac e( 10) ;

vasl·W(3.2.tr13); SpaCG(2); write("graden kt l ) ; .

print(" - conversie :11); space(lO);

v

as k 0 (2 2. x. 1 co) ;