Obliczenia technologiczne w przemyśle chemicznym - Krzysztof Schmidt-Szałowski, Krzysztof Krawczyk, Jan Petryk, Jan Sentek

Kup ebooka

104.00 zł
83.20 zł (67,60 zł najniższa cena z 30 dni)

-
Proszę czekać

1.4. Procesy chemiczne w reaktorach przemysłowych

Do podstawowych aparatów stosowanych w instalacjach chemicznych należą reaktory. Przebieg i właściwości procesów zachodzących w reaktorach chemicznych można przedstawić na przykładzie kilku prostych reaktorów modelowych, do których należą:

- zbiornikowy reaktor zamknięty o zupełnym wymieszaniu reagentów, - zbiornikowy reaktor przepływowy o zupełnym wymieszaniu reagentów, - reaktor rurowy o przepływie tłokowym.

Modelowe reaktory są odpowiednikami obiektów stosowanych w praktyce, mającymi ich podstawowe właściwości, jednak z pominięciem szczegółów ich konstrukcji (rys. 1.2).

Rys. 1.2. Podstawowe modele reaktorów: a) zbiornikowy reaktor zamknięty, b) zbiornikowy reaktor przepływowy, c) reaktor rurowy

1.4.1. Zbiornikowy reaktor zamknięty działający przy zupełnym wymieszaniu reagentów

Jest to model reaktora (rys. 1.2a) działającego w trybie okresowym, który spełnia następujące założenia:

Proces prowadzi się w układzie zamkniętym, tak że całkowita masa reagentów w reaktorze jest stała. Objętość czynna reaktora jest stała Zawartość reaktora jest skutecznie mieszana tak, że proces zachodzi w układzie jednorodnym, w którym stężenia reagentów, temperatura i inne wielkości opisujące stan układu mają jednakowe wartości w całej objętości.

Oprócz wymienionych założeń modelowych przyjęto następujące założenia dodatkowe:

a) W reaktorze zachodzi reakcja o równaniu

1A = 1B

b) Początkowy stopień przemiany substratu A x0 = 0, czyli początkowa ilość produktu n0[B] = 0. c) Temperatura reagentów jest stała w ciągu całego czasu trwania procesu.

Przebieg procesu charakteryzują następujące wielkości:

? - czas bieżący (zegarowy) [s];

? - czas przebywania reagentów w stanie reakcji [s];

- początkowa, aktualna i końcowa liczba moli substratu A w reaktorze;

, - stężenie molowe substratu A [mol/m3], w chwili początkowej i w toku procesu;

, - równowagowy i końcowy stopień przemiany A w B [mol/mol];

- początkowa i średnia szybkość reakcji [mol/(m3-s)].

W toku procesu stopień przemiany wzrasta, jego wartość przybliża się do wartości równowagowej a szybkość reakcji odpowiednio maleje, zbliżając się do zera. W praktyce proces przerywa się po upływie czasu ?, zwanego czasem przebywania reagentów w stanie reakcji (rys. 1.3). Średnią szybkość reakcji od początku do zakończenia procesu oblicza się następująco:

(1.12)

Rzeczywistą (chwilową) szybkość reakcji, zmieniającą się w toku procesu, określa zależność:

(1.13)

Znak minus w tym równaniu wynika z przyjętego założenia, że szybkość reakcji ma wartość dodatnią, natomiast ubytek substratu ma wartość

Rys. 1.3. a) Zmiany stopnia przemiany x w toku procesu w reaktorze zamkniętym; b) zmiany szybkości reakcji r w toku procesu w reaktorze zamkniętym;? - czas bieżący, ? - całkowity czas przebywania reagentów w stanie reakcji,x* - równowagowy stopień przemiany, r - średnia szybkość reakcji w okresachod 0 do ? oraz od 0 do ?

ujemną. W układzie o stałej objętości rzeczywistą szybkość reakcji można też określić następująco:

(1.14)

Zależność (1.14) można przekształcić, wprowadzając stopień przemiany x (rys. 1.3):

(1.15)

Korzystając z zależności (1.15), można obliczyć czas ?, w którym uzyskuje się określony stopień przemiany x. Po przekształceniu otrzymuje się równanie

(1.16)

Przeprowadzając całkowanie tego równania w granicach: ? od 0 do ? orazx od 0 do xk

(1.17)

otrzymuje się zależność:

(1.18)

lub

(1.19)

W przypadku procesu izotermicznego czas przebywania ?, w którym uzyskuje się założony stopień przemiany można obliczyć przez wprowadzenie do wyrażenia podcałkowego w równaniu (1.19) zależności szybkości reakcji od stopnia przemiany określonej przez równanie kinetyczne dotyczące założonej temperatury.

1.4.2. Zbiornikowy reaktor przepływowy działający przy zupełnym wymieszaniu reagentów

Za pomocą modelu zbiornikowego reaktora przepływowego (rys. 1.2b) można przedstawić przebieg procesu prowadzonego w trybie ciągłym w warunkach, które są zgodne z założeniami 2 i 3 poprzedniego modelu, tj. przy stałej objętości układu i jednorodnym składzie mieszaniny reagentów. Zakłada się, że proces przebiega w stanie stacjonarnym, co oznacza, że wielkości opisujące stan mieszaniny reagentów są niezmienne w czasie trwania procesu. Przez reaktor przepływa ciągły strumień reagentów o natężeniu WE na wejściu oraz WY na wyjściu. Ubytek substratów zużywanych w toku reakcji jest uzupełniany przez strumień E wprowadzany do reaktora, a wytworzony produkt i nieprzereagowany substrat są odprowadzane w strumieniu wychodzącym Y. W stanie stacjonarnym liczba moli substratu A w reaktorze nie zmienia się W takich warunkach szybkość reakcji r można obliczyć na podstawie wielkości strumieni substratu doprowadzanego WE[A] i odprowadzanego WY[A], z uwzględnieniem objętości czynnej reaktora

(1.20)

Rozpatrując przebieg procesu w tym modelu reaktora, przyjęto następujące założenia dodatkowe:

w reaktorze zachodzi reakcja o równaniu 1A = 1B; strumień E wprowadzany do reaktora zawiera tylko substrat A, czyli a początkowy stopień przemiany

Reaktor znajduje się w stanie stacjonarnym, dlatego wszystkie wielkości określające warunki procesu (masa i stężenie reagentów, natężenie przepływu, temperatura i in.) mają stałe wartości. Mieszanina reagentów ma skład jednorodny, zatem stopień przemiany x substratu A w całej objętości reaktora ma jednakową wartość, taką samą jak w strumieniu odprowadzanym z reaktora

(1.21)

W takich warunkach rzeczywista szybkość reakcji r jest w całej objętości układu jednakowa i równa średniej szybkości reakcji :

(1.22)

Na podstawie kinetycznej charakterystyki reakcji mającej postać , można obliczyć objętość czynną reaktora, w którym przy określonej wielkości strumienia substratu WE[A], w określonej temperaturze uzyskuje się stopień przemiany

(1.23)

1.4.3. Rurowy reaktor przepływowy o przepływie tłokowym

Określenie "przepływ tłokowy" dotyczy teoretycznego modelu ruchu płynów (cieczy lub gazów), który w rzeczywistych warunkach nie występuje, lecz bywa stosowany do przybliżonego opisu działania niektórych aparatów przepływowych. Przepływ płynu w takim modelowym aparacie (rys. 1.2c) spełnia następujące warunki:

1) w całym przekroju poprzecznym strumienia prędkość liniowa płynu jest jednakowa, a więc strumień porusza się tak, jak sztywny tłok;

2) w kierunku przepływu, czyli wzdłuż strumienia, płyn nie ulega mieszaniu.

Taki reaktor ma kształt rury o długości L i objętości (rys. 1.4). Pole przekroju poprzecznego strumienia S jest w całym reaktorze jednakowe.

Przyjęto ponadto następujące założenia:

a) w reaktorze zachodzi reakcja o równaniu 1A = 1B; b) strumień wprowadzany do reaktora zawiera tylko substrat A, czyli a początkowy stopień przemiany c) reaktor znajduje się w stanie stacjonarnym, czyli wartości parametrów (stężenia, natężenia przepływu i in.) są stałe w czasie prowadzenia procesu, ale w różnych miejscach w reaktorze ich wartości mogą być różne; d) proces zachodzi w jednakowej temperaturze.

W rurowym reaktorze przepływowym, do którego w sposób ciągły doprowadza się strumień substratu A i odprowadza się strumień zawierający reagenty A i B, stopień przemiany wzrasta w przepływającym wzdłuż reaktora strumieniu reagentów. Stopień przemiany substratu x w dowolnym miejscu wzdłuż reaktora może być wyznaczony na podstawie miejscowego natężenia przepływu substratu A:

(1.24)

(1.25)

gdzie: - natężenie przepływu substratu A w strumieniu doprowadzanym do reaktora (wyrażone odpowiednio w [mol/s] lub w [g/s]); - natężenie przepływu substratu A w odległości l od miejsca, w którym strumień substratu E dopływa do przestrzeni czynnej reaktora (wyrażone odpowiednio w [mol/s] lub w [g/s]).

Rys. 1.4. Model reaktora rurowego o przepływie tłokowym; W[A] i W[B] - natężenie przepływu substratu A i produktu B; pozostałe oznaczenia w tekście

Rozkład stopnia przemiany x i szybkości reakcji r wzdłuż reaktora przedstawiono na rys. 1.5. Odległość od miejsca, w którym strumień substratu E dopływa do przestrzeni czynnej reaktora zaznaczono na osi l równoległej do osi reaktora. Całkowitą długość przestrzeni czynnej reaktora oznaczono L. Do obliczenia średniej szybkości reakcji w tym reaktorze korzysta się z takich samych zależności, jak dla przepływowego reaktora zbiornikowego, ponieważ dotyczą one każdego z reaktorów przepływowych w stanie stacjonarnym [wynik (1.22)]. Rzeczywistą (miejscową) szybkość reakcji w miejscu określonym przez wartość l można obliczyć, gdy znana jest charakterystyka kinetyczna reakcji. W tym celu korzysta się z bilansu substratu A wewnątrz odcinka reaktora (rys. 1.4) o objętości d? i długości dl (przy czym Dopływa do niego strumień substratu W[A], a odpływa strumień A pomniejszony o ubytek spowodowany przebiegiem reakcji w elemencie objętości d?, czyli . Szybkość reakcji w tym elemencie objętości można określić następująco:

(1.26)

Wprowadzając stopień przemiany substratu A (x), otrzymuje się następujące zależności:

(1.27)

(1.28)

Odpowiednie podstawienie prowadzi do wyniku, który przedstawia szybkość reakcji w elemencie objętości d?:

(1.29)

Rys. 1.5. Rozkład wartości stopnia przemiany x i szybkości reakcji r wzdłuż reaktora rurowego o przepływie tłokowym

Po przekształceniu otrzymuje się równanie przedstawiające przyrost stopnia przemiany dx w objętości d?:

(1.30)

Ponieważ szybkość reakcji można wyrazić następująco:

(1.31)

Z równania (1.29) wynika zależność

(1.32)

Korzystając z tego równania, po przeprowadzeniu całkowania w granicach:? od 0 do ?r oraz x od 0 do otrzymuje się zależność, z której można obliczyć objętość czynną reaktora przy założonej wartości końcowego stopnia przemiany

(1.33)

W celu obliczenia czynnej objętości reaktora zależność (1.33) trzeba uzupełnić przez wprowadzenie charakterystyki kinetycznej reagującego układu w postaci równania

Warto zauważyć, że wiele ze stosowanych reaktorów przemysłowych działa w sposób zbliżony do jednego z przedstawionych modeli. Na przykład w reaktorach ze złożem ziarnistego katalizatora, przez które przepływa strumień gazowych reagentów, proces przebiega podobnie, jak w modelu reaktora rurowego o przepływie tłokowym.

Stosowane oznaczenia

a - udział molowy składnika w mieszaninie [mol/mol]

b - ciśnienie atmosferyczne [MPa]

c - stężenie molowe [mol/dm3]

c - ciepło właściwe [kJ/(kg-K)]

G - masowe natężenie przepływu; produktywność instalacji [kg/h]

H - entalpia ogrzania [J], [kJ]

?H - strumień entalpii ogrzania [J/s], [kJ/h]

J - pojemność cieplna strumienia masy [kJ/(h-°C)]

k - stała szybkości reakcji (jednostka zależy od rzędu reakcji)

l - długość drogi strumienia wzdłuż reaktora rurowego [m]

L - długość przestrzeni czynnej reaktora rurowego [m]

n - liczba moli [mol]

p - ciśnienie w aparaturze [MPa]

qr - molowe ciepło rozcieńczania [kJ/kmol]

QL - strumień entalpii parowania [kJ/h]

Qr - strumień entalpii reakcji R [J/s]; [kJ/h]

?Qi - strumień entalpii tworzenia składnika i [J/s]; [kJ/h]

r - szybkość reakcji [mol/(m3-s)]

r - średnia szybkość reakcji [mol/(m3-s)]

S - wielkość powierzchni poprzecznego przekroju aparatu [m2]

T - temperatura [°C]

T0 - temperatura odniesienia [°C]

? - czynna objętość aparatu [m3]

V - objętościowe natężenie przepływu strumienia [m3/s]

u - liniowa prędkość strumienia [m/s]

W - molowe natężenie przepływu [mol/s]; [kmol/h]

x - stopień przemiany substratu [mol/mol]

x* - równowagowy stopień przemiany substratu [mol/mol]

? - czas przebywania reagentów w stanie reakcji [s]; [h]

? - średni czas przebywania reagentów w stanie reakcji [s]; [h]

?[A] - ubytek substratu A [mol]

?x - przyrost stopnia przemiany substratu [mol/mol]

? - gęstość strumienia [mol/(m2-s)]

? - wydajność przetwarzania substratu w produkt [mol/mol]

?i - selektywność wytwarzania produktu i [mol/mol]

? - czas [s]; [h]

? - selektywność wytwarzania produktu w określonym miejscu lub w określonym czasie [mol/mol]

1.6. Bilanse technologiczne

Bilans technologiczny procesu produkcyjnego składa się z równań bilansowych, w których występują wielkości strumieni masy lub energii (entalpii), doprowadzanych i odprowadzanych z instalacji lub jej wyodrębnionej części, np. pojedynczego aparatu. Równania te zapisuje się w postaci: przychód = rozchód Dla każdego z występujących w procesie pierwiastków można ułożyć pierwiastkowe równanie bilansowe w którym uwzględnia się obecność tego pierwiastka w każdym ze strumieni, niezależnie od tego, w jakiej postaci występuje. Wielkości strumieni występujące w równaniach bilansowych odnoszą się do określonej podstawy bilansu, którą może być wielkość dobowej (godzinowej, rocznej) produkcji albo ilość przetwarzanego surowca, wytwarzanego produktu itp. Bilans technologiczny może dotyczyć całej instalacji produkcyjnej lub jej części. Zaznacza się to przez wskazanie określonego obszaru bilansowania (OB). W równaniach bilansowych uwzględnia się te strumienie, które przecinają granice obszaru bilansowania (rys. 1.6).

Rys. 1.6. Reaktor przepływowy; znaczenie symboli w tekście

1.6.1. Bilans procesu prowadzonego w trybie ciągłym

Gdy proces jest prowadzony w instalacji działającej w trybie ciągłym, pierwiastkowe równanie bilansowe dotyczące pierwiastka A zawartego w strumieniu E wprowadzanym do układu oraz w strumieniu odprowadzanym Y mają postać:

(1.34)

(1.35)

Symbole WE[A] i WY[A] oznaczają natężenie strumieni pierwiastka A przepływających przez granice obszaru bilansowanego OB. W równaniu (1.34) natężenie strumienia wprowadzanego E i odprowadzanego Y zostało wyrażone w jednostkach molowych, np. [mol/s]. Symbole GE[A] i GY[A] w równaniu (1.35) oznaczają natężenie tych samych strumieni wyrażone w jednostkach masy, np. [g/s].

Pozostałe symbole oznaczają:

?? [s] - okres, którego dotyczy bilans (godzina, doba itd.); ?n[A] [mol] - przyrost liczby moli pierwiastka A (dodatni lub ujemny) w obszarze OB w okresie ??; ?m[A] [g] - przyrost masy pierwiastka A (dodatni lub ujemny) w obszarze OB w okresie ??.

Gdy proces przebiega w warunkach stacjonarnych, wielkość zasobu każdego ze składników znajdujących się wewnątrz obszaru bilansowanego nie zmienia się z upływem czasu, czyli:

1.7. Bilans strumieni masy w instalacji produkcyjnej

Przykład 1. Produkcja kwasu siarkowego(VI) z siarki

Proces wytwarzania kwasu siarkowego(VI) z siarki przebiega zgodnie ze schematem instalacji, który w uproszczonej postaci przedstawiono na rys. 1.7. Proces ten składa się z następujących operacji:

Utlenianie (spalanie) siarki w strumieniu powietrza w piecu PK (oddział I):

(R1) S + O2 = SO2

Utlenianie SO2 w reaktorze katalitycznym AK (oddział II):

(R2) SO2 + 0,5 O2 = SO3

Absorpcja SO3 w wodzie z wytworzeniem kwasu siarkowego(VI) w absorberze WA (oddział III):

(R3) SO3 + H2O = H2SO4

Rys. 1.7. Schemat instalacji stosowanej do wytwarzania kwasu siarkowego(VI) z siarki. I - oddział utleniania siarki (reakcja R1), II - oddział utleniania SO2 (reakcja R2), III - oddział absorpcji SO3 (reakcja R3); znaczenie symboli w tekście

Stosowane oznaczenia

W[S], W[SO2], W[SO3], W[H2SO4], W[O2], W[N2], W[pow] i W[H2O] - molowe natężenie przepływu strumieni, odpowiednio, siarki, SO2, SO3, H2SO4, O2, N2, powietrza i wody [kmol/h] G[H2SO4], G[H2O] i G[KS] - masowe natężenie przepływu strumieni H2SO4, wody i roztworu kwasu siarkowego(VI) w wodzie [kg/h]

Miejsce strumieni w instalacji zaznaczono na schemacie (rys. 1.7).

Założenia

Z-1.

Proces wytwarzania kwasu siarkowego(VI) z siarki prowadzi się w instalacji, której uproszczony schemat przedstawiono na rys. 1.7. W instalacji prowadzonej w trybie ciągłym proces przebiega w warunkach stacjonarnych. Produktem jest roztwór kwasu siarkowego(VI) o stężeniu 96% mas. H2SO4 wytwarzany w ilości 12 000 kg/h (strumień KS).

Z-2.

W piecu PK prowadzi się reakcję R1, w której strumień doprowadzanej siarki S ulega całkowitemu przetworzeniu w ditlenek siarki. Do pieca PK doprowadza się strumień osuszonego powietrza P1 o składzie [mol/mol]: O2 - 0,21; N2 - 0,79. Z pieca PK doprowadza się strumień gazów G1 o składzie [mol/mol]: SO2 - 0,12; O2 - 0,09; N2 - 0,79.

Z-3.

Strumień gazów G1 kieruje się do reaktora katalitycznego AK, w którym prowadzi się reakcję R2, uzyskując końcowy stopień przemiany SO2 w SO3

Z-4.

W celu regulacji temperatury w złożu katalizatora, do reaktora AK doprowadza się strumień osuszonego powietrza P2, którego wielkość

Z-5.

Strumień gazów G2 odprowadzany z reaktora AK zawierający SO2, SO3, O2 i N2 kieruje się do absorbera WA, w którym następuje całkowite przetworzenie SO3 w H2SO4 (reakcja R3). SO2 pozostaje w gazach kominowych GK odprowadzanych z instalacji.

Potrzebne dane

Masa molowa M [kg/kmol]

Cel zadania

Należy obliczyć wielkość i skład strumieni w instalacji wytwarzającej kwas siarkowy(VI) z siarki zgodnie ze schematem (rys. 1.7) i przedstawić ogólny bilans masy tego procesu.

Podstawa bilansu

Natężenie strumienia KS zawierającego roztwór kwasu siarkowego(VI) o stężeniu 96% mas. H2SO4, który odprowadza się z instalacji

GKS = 12 000 kg/h

Rozwiązanie

Całkowity bilans strumieni masy w procesie wytwarzania kwasu siarkowego(VI) obejmuje pierwiastkowe bilanse siarki, tlenu, wodoru i azotu w instalacji produkcyjnej (obszar bilansowania OB1, rys. 1.7).

Wytwarzany strumień produktu KS składa się z kwasu siarkowego(VI) i wody. Wielkość i skład strumienia KS wynika z podstawy bilansu i założenia Z-1.

Strumień kwasu siarkowego(VI)

GKS [H2SO4] = 12 000 - 0,96 = 11 520 kg/h, czyli

WKS [H2SO4] = 117,6 kmol/h (1.36)

Strumień wody

GKS [H2O] = 12 000 - 0,04 = 480 kg/h, czyli

WKS [H2O] = 26,7 kmol/h (1.37)

Pierwiastkowy bilans siarki [S]

Siarka doprowadzana do instalacji w strumieniu S (obszar OB1) utlenia się w piecu PK do ditlenku siarki (założenie Z-2), który następnie utlenia się do tritlenku siarki w reaktorze AK (obszar OB2). W reaktorze AK uzyskuje się stopień przemiany SO2 w SO3 x = 0,98 (założenie Z-3). Strumień tritlenku siarki wytworzony w reaktorze AK jest w całości przetwarzany w kwas siarkowy(VI) w absorberze WA (obszar OB3). Ditlenek siarki znajdujący się w strumieniu gazów G2 opuszczających reaktor AK odprowadza się z instalacji w gazach kominowych GK, a zatem wielkość strumienia SO2 w gazach kominowych jest taka sama, jak w strumieniu G2 odprowadzanym z reaktora AK. Można ją obliczyć na podstawie założonego stopnia przemiany SO2 w SO3 w tym reaktorze. Tritlenek siarki zawarty w strumieniu G2 jest w całości przetwarzany w kwas siarkowy(VI) (założenie Z-5), dlatego korzystając z wyniku (1.36), można obliczyć natężenie strumienia WG2[SO3]:

WG2 [SO3] = WKS[H2SO4] = 117,6 kmol/h

Z założonej wielkości stopnia przemiany SO2 w reaktorze AK (x = 0,98) wynika natężenie strumienia SO2

czyli

(1.38)

Siarka doprowadzana do pieca PK zostaje w całości przetworzona i odprowadzana w postaci SO2 (założenie Z-2), dlatego strumienie WS[S] i WG1[SO2] mają jednakową wielkość:

(1.39)

Wielkość strumienia SO2 w gazach kominowych wynika z różnicy:

(1.40)

Na podstawie wyniku (1.39) i założonego składu strumienia gazów G1 odprowadzanych z pieca PK (założenie Z-2) można obliczyć wielkość składowych strumienia G1:

(1.41)

Skład strumienia G1:

(obliczony poprzednio) (1.42)

(1.43)

(1.44)

Tabela 1.1. Pierwiastkowy bilans siarki [S] w obszarze OB1

PRZYCHÓD

ROZCHÓD

Strumień [S]

kmol/h

kg/h

Strumień [S]

kmol/h

kg/h

S [siarka]

120

3 840

KS [H2SO4]

117,6

3 763

GK [SO2]

2,4

77

Suma przychodu

120

3 840

Suma rozchodu

120,0

3 840

Pierwiastkowy bilans wodoru [H]

W odprowadzanym z instalacji strumieniu KS znajduje się wodór będący składnikiem kwasu siarkowego(VI) i wody [wyniki (1.36) i (1.37)].

Strumień wodoru [H] zawartego w kwasie siarkowym(VI):

(1.45)

Strumień wodoru [H] zawartego w wodzie rozcieńczającej:

(1.46)

Całkowita ilość wodoru [H] zawarta w strumieniu KS:

(1.47)

Doprowadzany do instalacji wodór jest zawarty w strumieniu wody (W) kierowanym do absorbera WA. Wielkość strumienia wody wynika z bilansu wodoru:

(1.48)

Tabela 1.2. Pierwiastkowy bilans wodoru [H] w obszarze OB1

PRZYCHÓD

ROZCHÓD

Strumień wodoru [H]

kmol/h

kg/h

Strumień wodoru [H]

kmol/h

kg/h

W [H2O]

288,6

288,6

KS [H2SO4]

235,2

235,2

KS [H2O]

53,4

53,4

Suma przychodu

288,6

288,6

Suma rozchodu

288,6

288,6

Pierwiastkowy bilans tlenu [O]

Na podstawie składu strumienia G1 kierowanego z pieca PK do reaktora AK i wielkości strumienia siarki S doprowadzanej do pieca PK [wynik (1.39)] można obliczyć wielkość strumieni składników zawierających tlen. Do reaktora AK wprowadza się strumień gazów G1 o natężeniu 1000 kmol/h [wynik (1.41)] i o składzie wynikającym z założenia Z-2: SO2 - 0,12; O2 - 0,09; N2 - 0,79 [mol/mol]. Na tej podstawie można obliczyć wielkości poszczególnych składników strumienia G1 [wyniki (1.38) oraz (1.42)-(1.44)].

Do pieca PK doprowadza się strumień powietrza P1 zawierający tlen i azot. Skład tego strumienia wynika ze składu strumienia G1 [wyniki (1.42)-(1.44)]:

(1.49)

[wynik (1.44)]

W celu regulacji temperatury reagentów do reaktora AK doprowadza się dodatkowy strumień powietrza P2 o natężeniu WP2 = 0,5 - WP1 = 500 kmol/h.

Strumień ten zawiera:

(1.50)

(1.51)

Do absorbera WA doprowadza się strumień wody W, którego wielkość wynika z bilansu wodoru [wynik (1.47)]:

(1.52)

(1.53)

Do instalacji (obszar OB1) doprowadza się następujące strumienie zawierające tlen:

- strumień P1 doprowadzany do pieca PK: [wynik (1.49)] - strumień P2 doprowadzany do reaktora AK: [wynik (1.50)] - strumień W doprowadzany do absorbera WA: WW[H2O] = 144,3 kmol/h [wynik (1.53)]

Całkowita ilość tlenu [O] doprowadzanego do instalacji:

(1.54)

Strumienie odprowadzane z instalacji (obszar OB1), w których znajduje się tlen:

- strumień tlenu [O] zawarty w produkowanym kwasie siarkowym(VI) [wynik (1.36)]:

(1.55)

- strumień tlenu [O] zawarty w wodzie rozcieńczającej [wynik (1.37)]: - w gazach kominowych GK odprowadza się SO2 w ilości w którym strumień tlenu [O] jest równy 4,8 kmol/h.

Nieprzetworzony w toku procesu tlen zostaje odprowadzony z instalacji w strumieniu gazów kominowych GK. Strumień tlenu WGK[O2] można obliczyć na podstawie pierwiastkowego bilansu tlenu [O] przez porównanie strony przychodu i rozchodu

Z tego równania wynika wielkość strumienia tlenu zawartego w gazach kominowych

(1.56)

Tabela 1.3. Pierwiastkowy bilans tlenu [O] w obszarze OB1

PRZYCHÓD

ROZCHÓD

Tlen [O]

kmol/h

kg/h

Tlen [O]

kmol/h

kg/h

P1 [powietrze]

420,0

6 720

KS [H2SO4]

470,4

7 526

P2 [powietrze]

210,0

3 360

KS [H2O]

26,7

427

W [woda]

144,3

2 310

GK [SO2]

4,8

77

GK [O2]

272,4

4 360

Suma przychodu

774,3

12 390

Suma rozchodu

774,3

12 390

Pierwiastkowy bilans azotu [N2] w obszarze OB1

We wszystkich strumieniach w instalacji azot występuje w postaci niezwiązanej. Jest wprowadzany w strumieniach powietrza P1 i P2. Odprowadza się go w strumieniu gazów kominowych GK. Równanie bilansowe azotu ma postać:

[wynik (1.44)]

[wynik (1.51)]

Suma przychodu

Suma rozchodu

Ogólny bilans instalacji (obszar OB1)

Przychód

Siarka

czyli (tab. 1.1)

Powietrze (P1 + P2)

w tym:

Tlen: (tab. 1.3)

Azot: (z bilansu azotu)

Suma tlenu i azotu (P1 + P2): 43 260 kg/h

Woda

(tab. 1.3)

Suma przychodu:

Rozchód

Kwas siarkowy(VI) (96% H2SO4) 12 000 kg/h (podstawa bilansu)

Gazy kominowe 37 700 kg/h

Składniki gazów kominowych:

Azot

Tlen

SO2

Suma rozchodu:

Ogólny bilans instalacji został przedstawiony na wykresie strumieniowym (rys. 1.8). Rozbieżność między całkowitymi wielkościami przychodu i rozchodu wynika ze stosowanej dokładności obliczeń strumieni masy (1 kg/h). Wielkości przedstawione na wykresie strumieniowym podano z dokładnością do 10 kg/h.

Rys. 1.8. Wykres strumieniowy przedstawiający ogólny bilans instalacji stosowanej do wytwarzania kwasu siarkowego(VI) z siarki. Wielkość strumieni podano z dokładnością do 10 kg/h

1.9. Bilans strumieni entalpii w reaktorze przepływowym

Przykład 2. Utlenianie amoniaku do tlenków azotu

W procesie wytwarzania kwasu azotowego(V) z amoniaku pierwszym etapem jest wysokotemperaturowy, katalityczny proces, w którym w wyniku utleniania amoniaku w mieszaninie z powietrzem otrzymuje się tlenek azotu NO. Proces ten prowadzi się w reaktorze przepływowym RK (rys. 1.9) pod zwiększonym ciśnieniem. Stosowane w tym procesie katalizatory mają postać pakietu siatek wykonanych ze stopów platyny. Z wielu reakcji zachodzących w toku procesu, w tym przykładzie uwzględniono dwie, które mają największy wpływ na jego efekt cieplny:

(R1) 2 NH3 + 2,5 O2 = 2 NO + 3 H2O

(R2) 2 NH3 + 1,5 O2 = N2 + 3 H2O

Rys. 1.9. Reaktor RK do utleniania amoniaku; znaczenie symboli w tekście

Stosowane oznaczenia
W[NH3], W[NO], W[O2], W[N2], W[H2O] - natężenie przepływu strumieni, odpowiednio NH3, NO, O2, N2, H2O [kmol/h] - stopień przemiany amoniaku w NO w reaktorze RK (reakcja R1) [mol/mol] - strumień entalpii ogrzania [kJ/h] - strumień entalpii tworzenia [kJ/h] - strumień entalpii reakcji [kJ/h]

Potrzebne dane

Entalpia tworzenia w temp. 25 °C [kJ/kmol]

NH3 -46-103

NO 90-103

H2Og (para) -240-103

Średnie ciepło molowe [kJ/(kmol-K)]

NH3 36

O2 28

N2 28

NO 30

H2Og (para) 38

Założenia

Z-1. Proces utleniania amoniaku w reaktorze RK jest prowadzony w trybie ciągłym w warunkach stacjonarnych.

Z-2. Do reaktora RK doprowadza się strumień E zawierający mieszaninę amoniaku z powietrzem o składzie: 60 kmol/h NH3 oraz 500 kmol/h powietrza, w którym znajduje się 100 kmol/h O2 i 400 kmol/h N2.

Z-3. Cały strumień amoniaku doprowadzany do reaktora RK ulega przetworzeniu w produkty reakcji R1 i R2.

Z-4. Stopień przemiany amoniaku w NO a pozostała część amoniaku zostaje przetworzona w reakcji (R2). Z reaktora RK odprowadza się strumień produktów Y o składzie:

NO 57 kmol/h; N2 401,5 kmol/h; H2O 90 kmol/h; O2 26,5 kmol/h.

Z-5. Temperatura strumienia E doprowadzanego do reaktora RK:

Z-6. W reaktorze RK proces przebiega w warunkach adiabatycznych (bez przenikania ciepła przez ściany reaktora).

Cel zadania

Przyjmując temperaturę odniesienia należy obliczyć:

bilans strumieni entalpii w procesie utleniania amoniaku; temperaturę strumienia Y produktów odprowadzanych z reaktora RK.

Podstawa bilansu: Natężenie strumienia amoniaku doprowadzanego do reaktora RK

Rozwiązanie

1. Bilans strumieni entalpii

Przychód [kJ/h]

Entalpia ogrzania składników strumienia E

Łącznie (1.61)

Entalpia tworzenia składników strumienia E

(1.62)

Rozchód [kJ/h]

Entalpia ogrzania składników strumienia Y

Łącznie (1.63)

Entalpia tworzenia składników strumienia Y

Łącznie (1.64)

Wielkość strumienia entalpii reakcji wynika z bilansu procesu w reaktorze RK

Zestawienie pozycji bilansu entalpii w reaktorze RK przedstawiono w tab. 1.4.

Tabela 1.4. Bilans strumieni entalpii w procesie utleniania amoniaku przy podstawie bilansu: WE[NH3] = 60 kmol/h i temperaturze odniesienia T0 = 25 °C

Przychód

kJ/h

Rozchód

kJ/h

Strumień E

Strumień Y

Entalpia ogrzania składników strumienia E

Entalpia ogrzania składników strumienia Y

NH3

108 - 103

NO

1 450 - 103

O2

140 - 103

O2

620 - 103

N2

560 - 103

N2

9 550 - 103

H2O

2 900 - 103

Razem:

Razem:

808 - 103

14 520 - 103

Entalpia tworzenia amoniaku

Entalpia tworzenia NO i H2O

5130 - 103

-2 760 - 103

-21 600 - 103

Razem:

-16 470 - 103

Suma przychodu:

Suma rozchodu:

Entalpia strumienia E

Entalpia strumienia Y

-1 952 - 103

-1 950 - 103

(Rozbieżność między przychodem i rozchodem wynika z zaokrąglenia niektórych wielkości występujących w bilansie)

2. Temperatura strumienia Y produktów odprowadzanych z reaktora RK

Ogólne równanie bilansu entalpii w reaktorze RK ma postać:

(1.65)

Po podstawieniu odpowiednich wartości strumieni entalpii otrzymuje się równanie bilansowe

Z tego równania wynika następująca temperatura strumienia Y zawierającego produkty reakcji odprowadzane z reaktora RK

(1.66)

1.8. Bilans entalpii

Racjonalna gospodarka energią należy do podstawowych problemów produkcji przemysłowej. W przemyśle chemicznym sposób organizacji procesów produkcyjnych zależy w znacznym stopniu od efektu energetycznego prowadzonych reakcji. W przypadku reakcji endotermicznych potrzebną energię doprowadza się ze źródeł zewnętrznych. Instalacje, w których prowadzi się procesy egzotermiczne, są samowystarczalne w zakresie bilansu energii (mają charakter autotermiczny), a uzyskany znaczny nadmiar ciepła można wykorzystać do różnych celów. Ekonomiczna wartość zasobu entalpii zawartej w materiale zależy od temperatury, w której może być przekazywana. Entalpię zawartą w materiałach o wysokiej temperaturze można wykorzystać do wielu celów. Entalpii materiałów o niskiej temperaturze nie wykorzystuje się praktycznie. Gdy trzeba ją odprowadzić z instalacji, stosuje się czynniki chłodzące, np. powietrze lub wodę. Dokumentem przedstawiającym sposób gospodarowania energią w instalacji produkcyjnej jest bilans strumieni energii lub entalpii.

Jeżeli w instalacji produkcyjnej proces przebiega pod stałym ciśnieniem, to jako miarę energii zawartej w przetwarzanych materiałach można stosować ich entalpię, która występuje w dwóch postaciach. Jedną z nich jest entalpia ogrzania materiału, czyli ilość ciepła potrzebnego do jego ogrzania do określonej temperatury. Drugą postacią jest entalpia tworzenia substancji zawartych w materiale, wynikająca z ich struktury molekularnej i stanu skupienia. Entalpią tworzenia nazywa się entalpię reakcji syntezy danej substancji z substratów będących substancjami prostymi i występujących w postaci trwałej w warunkach syntezy. Na podstawie entalpii tworzenia reagentów można obliczyć efekt cieplny zachodzących między nimi reakcji. Zgodnie z przyjętą konwencją, entalpia tworzenia reakcji egzotermicznych jest wyrażana liczbą ujemną.

1.8.1. Entalpia ogrzania zasobu i strumienia materiału

Entalpię ogrzania ?H zasobu materiału o masie m i temperaturze T określa się następująco:

(1.57)

gdzie: m - masa zasobu materiału [g]; c - ciepło właściwe pod stałym ciśnieniem [J/(g-°C)]; T - temperatura [°C]; T0 - umownie przyjęta temperatura odniesienia [°C], w której entalpia ogrzania substancji jest równa zeru.

Wielkość zasobu materiału może być wyrażona w jednostkach masy, np. [g, kg], lub w jednostkach molowych [mol, kmol]. Iloczyn [J/(°C)] jest pojemnością cieplną zasobu materiału.

Natężenie strumienia materiału może być wyrażane w jednostkach masy, np. [kg/s] lub jednostkach molowych, np. [kmol/s]. Wielkość strumienia entalpii ogrzania przenoszonego przez strumień materiału o natężeniu W [kmol/s] lub G [kg/s] w temperaturze T można wyrazić za pomocą wzorów:

(1.58)

(1.59)

gdzie: C - ciepło molowe pod stałym ciśnieniem [J/(mol - °C)]; c - ciepło właściwe pod stałym ciśnieniem [J/(g - °C)].

Wielkości [J/(°C - s)] oraz [J/(°C - s)] wyrażają pojemność cieplną strumienia materiału.

Ciepło właściwe i ciepło molowe zależą od temperatury. W obliczeniach, które nie wymagają dużej dokładności, w celu uproszczenia rachunku można stosować ich wartości średnie w zakresie od temperatury odniesienia T0 do aktualnej temperatury zasobu lub strumienia materiału.

1.8.2. Entalpia reakcji

W obliczeniach technologicznych efekt cieplny reakcji jest traktowany inaczej niż w obliczeniach fizykochemicznych. W wyniku reakcji egzotermicznych następuje wzrost entalpii ogrzania i wzrost temperatury reagentów, dlatego w bilansie cieplnym procesu zostaje to zapisane po stronie przychodu. Natomiast efekt cieplny reakcji endotermicznej umieszcza się po stronie rozchodu. Przy obliczaniu bilansu cieplnego złożonego procesu, w którym zachodzi wiele reakcji, całkowity efekt cieplny procesu dotyczący wszystkich zachodzących reakcji, może być przedstawiony w postaci różnicy między sumami entalpii tworzenia substratów i produktów:

(1.60)

gdzie: - liczba moli reagenta i (substratu lub produktu); - entalpia tworzenia składnika i [J/mol].

Sumowanie dotyczy wszystkich reagentów biorących udział w reakcji. W przypadku procesu egzotermicznego ma wartość dodatnią.

Entalpia tworzenia substancji chemicznych zależy od temperatury, jednak przedstawiony sposób obliczania efektu cieplnego reakcji wymaga, żeby w równaniach bilansowych były stosowane wartości entalpii tworzenia reagentów określone w takiej temperaturze, jaką obrano za temperaturę odniesienia.

1.3. Charakterystyka procesów chemicznych

Stosowane w tym podręczniku określenie "reakcja chemiczna" dotyczy samej przemiany substratu w produkt, natomiast określenie "proces chemiczny" obejmuje oprócz reakcji także towarzyszące jej zjawiska przenoszenia masy i energii. Do opisu procesów chemicznych stosuje się wielkości określające warunki, w których proces przebiega, jego postęp i dynamikę oraz końcowy wynik.

Wielkościami, które określają warunki procesu, są m.in. temperatura, ciśnienie, stężenia reagentów, natężenie przepływu strumieni masy lub energii.

Miarą postępu procesu jest stopień przemiany określający, jaka część substratu została przetworzona w produkt.

Dynamikę procesu określa szybkość reakcji chemicznych i przemian fizycznych zachodzących w aparaturze.

Wynik procesu charakteryzuje jego wydajność produkcyjna, zwana produktywnością, czyli ilość produktu wytwarzanego w jednostce czasu lub natężenie odprowadzanego strumienia produktu, a także wydajność surowcowa, sprawność energetyczna i inne.

1.3.1. Stan układu reagentów

Układ zawierający mieszaninę reagentów może znajdować się w jednym ze stanów, które określa się jako:

stan równowagi termodynamicznej, będący także stanem równowagi chemicznej, w którym skład mieszaniny reagentów nie ulega zmianie, gdyż wypadkowa szybkość zachodzących w nim reakcji jest równa zeru; stan reakcji, w którym skład układu zmienia się wskutek zachodzących reakcji chemicznych. Gdy proces zachodzi samorzutnie w reaktorze zamkniętym, odizolowanym od otoczenia, to stan reagującego układu przybliża się do stanu równowagi chemicznej; stan zamrożenia (zahamowania) reakcji chemicznych, w którym równowaga nie występuje, lecz reakcje nie zachodzą, gdyż w danych warunkach (np. z powodu niskiej temperatury) ich szybkość jest równa zeru.

1.3.2. Stopień przemiany, wydajność surowcowa

Miarą postępu procesu chemicznego jest stopień przetworzenia substratu, zwany stopniem przemiany, który wskazuje, jaka część początkowej ilości substratu została przetworzona w produkt. Jeżeli reakcję o równaniu

1A ? 1B

prowadzi się w zamkniętym reaktorze zbiornikowym, stopień przemiany substratu A w produkt B, czyli można obliczyć na podstawie następujących definicji:

(1.1)

(1.2)

gdzie: n0[A], m0[A] - początkowa liczba moli lub początkowa masa substratu A; n[A], m[A] - liczba moli lub masa substratu A, który pozostał nieprzetworzony w toku reakcji.

Wynik zakończonego procesu przedstawia się za pomocą wydajności surowcowej ?, czyli wydajności przetworzenia substratu A w produkt B, którą w przypadku procesu prowadzonego w reaktorze zamkniętym można określić, podobnie jak stopień przemiany, na podstawie ubytku substratu A (w jednostkach molowych lub w jednostkach masy):

(1.3)

(1.4)

gdzie: - liczba moli lub masa substratu A po zakończeniu procesu.

W procesie prowadzonym w reaktorze przepływowym w warunkach stacjonarnych, końcowy stopień przemiany określa się na podstawie wielkości strumienia substratu A doprowadzanego do reaktora oraz strumienia nieprzereagowanego substratu A odprowadzanego z reaktora:

(1.5)

lub: (1.6)

WE[A] i WY[A] oznaczają wielkość strumienia substratu A, odpowiednio: wprowadzanego do reaktora oraz odprowadzanego, wyrażone w jednostkach molowych (np. mol/s), a GE[A] i GY[A] oznaczają wielkości tych strumieni wyrażone w jednostkach masy (np. g/s).

1.3.3. Termodynamiczna charakterystyka układu reagentów

W stanie równowagi termodynamicznej zależności między stężeniami (ściślej: między aktywnościami) reagentów przedstawia się za pomocą stałej równowagi reakcji K. Na przykład stała równowagi reakcji odwracalnej o równaniu stechiometrycznym

aA + bB = cC + dD

ma postać:

(1.7)

Symbole , , i  oznaczają aktywności reagentów (w postaci bezwymiarowej) w stanie równowagi. Stała K jest termodynamiczną stałą równowagi, której wartość zależy tylko od temperatury. W obliczeniach technologicznych, jeśli nie powoduje to istotnego błędu, bywa stosowana stała równowagi w postaci:

(1.8)

gdzie: , , i  są stężeniami molowymi reagentów w stanie równowagi.

W przypadku, gdy reagenty występują w fazie gazowej, stała równowagi może mieć postać:

(1.9)

gdzie: , , i  są ciśnieniami cząstkowymi reagentów w stanie równowagi.

Błąd, który wynika z zastosowania Kp zamiast termodynamicznej stałej równowagi K, jest na ogół niewielki i można go pominąć, jeżeli ciśnienie w układzie jest zbliżone do ciśnienia normalnego (0,1 MPa). Na podstawie stałej równowagi K można wyznaczyć równowagowe wartości stężeń molowych lub ciśnień cząstkowych reagentów.

1.3.4. Kinetyczna charakterystyka układu reagentów w stanie reakcji

W obliczeniach technicznych za miarę szybkości reakcji przyjmuje się ilość przetworzonego substratu lub ilość wytworzonego produktu w ciągu jednostki czasu, w jednostce objętości układu. Wielkość tę wyraża się w jednostkach:

; ; lub analogicznych.

Jako szybkość reakcji można przyjąć szybkość przetwarzania dowolnego substratu lub wytwarzania dowolnego produktu. Przyjmując za podstawę ilość składnika A, który jest substratem przykładowej reakcji 1A = 1B, jej szybkość można zdefiniować następująco:

(1.10)

gdzie: oznacza ilość substratu A, który został przetworzony w czasie , a  jest czynną objętością reaktora. Symbol oznacza średnią szybkość reakcji odnoszącą się do całej objętości czynnej reaktora i do czasu trwania procesu Szybkość reakcji może zmieniać się w toku procesu, dlatego korzysta się z wielkości będącej chwilową (rzeczywistą) szybkością reakcji:

(1.11)

gdzie: oznacza liczbę moli dowolnego substratu lub produktu reakcji. Jeśli dotyczy ona substratu po prawej stronie równania (1.11) umieszcza się znak "-", a jeśli produktu to znak "+".