Dyfuzja w fazie gazowej wieloskładnikowa.doc

Wykłady z inżynierii chemicznej – Dyfuzja

Maksymilian Pająk i Wojciech Skrzypiński

Procesy dyfuzyjne

              Wiele procesów inżynierii chemicznej, w których zachodzi międzyfazowa wymiana substancji, najczęściej określa się jako procesy wymiany masy. Należą do nich absorpcja, desorpcja, adsorpcja, ekstrakcja, nawilżanie gazów, destylacja i rektyfikacja oraz wiele podobnych. Teoretyczny opis tych procesów czy inaczej mówiąc operacji jest trudniejszy od opisu procesów wymiany ciepła. Już sama ich systematyka jest bogatsza i można ją sformułować według różnych kryteriów. Można mówić o procesach wymiany masy, które przebiegają: jednokierunkowo, przeciwkierunkowo, z wymianą jednego składnika, z wymiana wielu składników, z równoczesną reakcja chemiczną, z równoczesną wymianą ciepła, izotermicznie, z efektem cieplnym, a to i tak nie wypełnia całej ich różnorodności. W różnych operacjach przemysłowych można przecież spotkać przypadki o charakterze mieszanym. Oprócz różnorodności układów dwufazowych, jak na przykład faza gazowa i ciekła, faza ciekła i faza stała czy dwie różne fazy ciekłe, może występować także różnorodność substancji, które podlegają wymianie masy. Jako przykład można wymienić dwa procesy, w których z tego samego surowca jakim jest gaz koksowniczy raz można usuwać benzen, a innym razem amoniak. Jeśli dodatkowo uwzględni się fakt, że w procesach wymiany ciepła najczęstszym elementem aparatury jest rura, a w procesach wymiany masy, gdzie chodzi o zapewnienie jak największego pola powierzchni międzyfazowej i jak najdłuższego czasu kontaktu faz, stosuje się aparaty o różnorodnych wypełnieniach, o czym szczegółowo będzie mowa w dalszej części, to twierdzenie o bardziej skomplikowanym opisie procesów wymiany masy staje się w pełni uzasadnione.

              Procesy wymiany masy najczęściej wiążą się z dyfuzyjnym ruchem molekuł, stąd teoretyczny opis tych operacji opiera się na prawach dyfuzyjnego ruchu masy, a procesy zyskały nazwę procesów dyfuzyjnych. Matematyczny opis wszystkich operacji dyfuzyjnych pojawiających się w inżynierii chemicznej i procesowej ma zatem wspólne podstawy teoretyczne.

Mechanizmy dyfuzyjnego ruchu masy

              Mechanizmy ruchu masy wykazują pewne podobieństwa do mechanizmów ruchu ciepła. Przy wymianie ciepła można wyróżnić przewodzenie, konwekcję i promieniowanie. W przypadku wymiany masy rozróżnia się tylko dwa mechanizmy, tj. dyfuzję i konwekcję.

              Aby w procesach wymiany masy zachodził transport substancji (przepływ określonego rodzaju molekuł) konieczny jest pewien gradient stężenia transportowanego składnika wzdłuż drogi poruszania się tych molekuł, czyli ruch masy odbywa się zawsze od ośrodka (miejsca) o stężeniu wyższym do ośrodka (miejsca) o stężeniu niższym.

              Dyfuzja jest zjawiskiem molekularnym i podobnie jak przewodzenie ciepła przebiega w warstwach nieruchomych lub poruszających się ruchem laminarnym. Dyfuzja najłatwiej zachodzi w fazie gazowej, trudniej w fazie ciekłej, a najtrudniej w fazie stałej. Szybkość dyfuzji rośnie ze wzrostem temperatury, maleje ze spadkiem ciśnienia i jest mniejsza dla molekuł o większej masie molowej.

              Dyfuzyjny ruch masy może występować, jeśli w dwóch różnych punktach przestrzeni występuje różnica stężeń, różnica temperatur lub różnica ciśnień. Mówimy wówczas o dyfuzji stężeniowej, termodyfuzji lub dyfuzji ciśnieniowej. Znany jest także dyfuzyjny ruch masy pod wpływem różnicy natężenia pola elektrycznego, taki przypadek w przyrodzie obserwuje się w komórkach organizmów żywych i wówczas mówi się o elektrodyfuzyjnym ruchu masy.

              Konwekcja lub ruch masy przez konwekcję odbywa się w ten sposób, że dana faza znajdująca się w ruchu burzliwym przenosi pewne porcje płynu z miejsc o stężeniu wyższym do miejsc o stężeniu niższym, zupełnie podobnie jak w procesach wymiany ciepła. Jednakże w procesach międzyfazowej wymiany masy tak opisana konwekcja nie zachodzi. W obecności drugiej fazy na ich granicy następuje spowolnienie ruchu płynu aż do ruchu laminarnego. Zatem wymiana masy do drugiej fazy następuje poprzez laminarne warstwy graniczne istniejące pomiędzy dwiema fazami. W tych warstewkach granicznych zanika konwekcja, a pojawia się dyfuzja. Zatem transport składnika z głębi burzliwej fazy do granicy międzyfazowej, to najpierw konwekcja a potem dyfuzja. Taki ruch składnika odbywający się w obrębie jednej fazy nazywa się wnikaniem masy. Jeśli połączy się dwa procesy wnikania masy w sąsiadujących fazach, to mamy do czynienia z przenikaniem masy. Zatem przenikanie masy odbywa się następująco najpierw składnik z głębi (z rdzenia) jednej fazy wskutek wnikania dociera do granicy międzyfazowej, przechodzi do drugiej fazy na granicy i w końcu wnika do głębi drugiej fazy. Schematyczny opis tego zjawiska można przedstawić graficznie za pomocą wykresu stężenia (np. ułamków molowych) od drogi.

Rys. 1. Rozkład stężeń przy przenikaniu masy z fazy y do fazy x

Dwie fazy są rozdzielone granicą międzyfazową „i”. W pobliżu granicy każda z faz wytwarza warstewkę laminarną o grubości sy lub sx. Daleko od granicy międzyfazowej stężenia w każdej z faz są stałe, a w warstewkach granicznych zmieniają się wskutek dyfuzji.

Dyfuzja w fazie gazowej

Teoria dyfuzji gazów rozwinęła się wraz z teorią kinetyczną gazów. Wymienić tu należy takich badaczy jak Maxwell, Stefan, Mayer, Sutherland i inni).

Dyfuzja może być procesem ustalonym bądź nieustalonym. W pierwszym przypadku stężenie transportowanego składnika jest wyłącznie funkcją miejsca:

,

a w przypadku dyfuzji nieustalonej jest dodatkowo funkcją czasu:

,

Weźmy pod uwagę warstwę gazu znajdującą się w spoczynku, ograniczoną dwiema płaszczyznami I i II. Przez warstwę tę dyfunduje składnik mieszaniny gazowej A. Ruch składnika występuje na skutek istniejącej różnicy stężeń. Do wyrażenia stężeń można wykorzystać koncentrację (stężenie) oznaczaną symbolem cA [kmol/m3]. Rozpatrzmy dyfuzję ustaloną. Zatem stężenia w przekrojach I i II są niezmienne w czasie, tzn.:
cAI = const, cAII = const, a ponadto cAI > cAII.

Przemieszczanie się cząsteczek A w kierunku dyfuzji (w kierunku malejącego stężenia) wynika z różnicy stężeń składnika A w przekrojach I i II, można zatem mówić o gradiencie stężenia składnika A wzdłuż drogi. Stężenie składnika można wyrazić jako jego prężność pA, a więc można mówić o gradiencie prężności dpA/ds, który ze względu na spadek prężności wzdłuż drogi ma wartość mniejszą od zera. Gradient prężności składnika A należy traktować jako miarę oporu na jaki napotykają cząsteczki substancji A podczas dyfuzji wśród innych cząsteczek. Dla mieszaniny dwuskładnikowej Maxwell i Stefan opór dyfuzji w kierunku s opisali zależnością:

                            (1)

gdzie:

dpA               - strata prężności składnika A w kierunku zgodnym z drogą s na drodze o długości ds, Pa

              cA, cB              - stężenie składnika A lub B, kmol/m3,

              uA, uB              - prędkość (ruchliwość) cząsteczek A lub B, m/s,

lAB              - współczynnik proporcjonalności charakterystyczny dla układu A i B, [kg×m3/(kmol2×s)]

Natomiast dla przypadku, w którym cząsteczki A dyfundują przez mieszaninę składników B, C, ..., N, odpowiednie równanie opisujące gradient prężności składnika A ma postać:

                            (2)

Oba powyższe równania wynikają z teorii kinetycznej gazów.

              Aby za pomocą równań opisać ile cząsteczek substancji A przemieszcza się wzdłuż drogi s w jednostce czasu dokonajmy przekształceń powyższego równania.

Zdefiniujmy:

Gęstość strumienia składnika A:               [kmol/(m2×s)]

Gęstość strumienia składnika B:               [kmol/(m2×s)]

Gęstość strumienia składnika N:               [kmol/(m2×s)]

Prężność składnika A:              

              oraz             

                            (3)

Kinematyczny współczynnik dyfuzji:

[m2/s]

, ...,

              stężenie ogólne mieszaniny A, B, ..., N [kmol/m3]

Dynamiczny współczynnik dyfuzji:

[kmol/(m×s)]

, ...,

                            (4)

Oznaczmy stosunki strumieni (stosunki gęstości strumieni) odniesione do gęstości strumienia IA:

, , ...,

              Jeśli kierunek dyfuzji wybranego składnika będzie zgodny z kierunkiem dyfuzji składnika A, to strumień dyfuzji tego składnika jest dodatni, a w przypadku dyfuzji w kierunku odwrotnym do składnika A – znak ujemny. Stosunki strumieni oczywiście też będą podlegały tej regule. Zatem po przekształceniach uzyskuje się zależność:

                            (5)

Stąd:

                            (6)

              Wstawmy:

, , ...,

Zatem po wprowadzeniu ułamków molowych uzyskuje się równanie:

                            (*7)

Analizując powyższe równanie widać, że strumień dyfundującego składnika A jest proporcjonalny do jego gradientu, a ponadto zależy od stężeń pozostałych składników obecnych w mieszaninie, od ich współczynników dyfuzji, ich kierunku dyfuzji oraz od stosunków strumieni. Gdyby współczynniki dyfuzji wszystkich składników mieszaniny gazowej były jednakowe i równe , to równanie przekształca się do postaci:

                            (8)

Zdefiniujmy sumę względnych strumieni wszystkich składników:

                            (9)

Wyrażenie w nawiasie ulega zmianie na skutek zmienności ułamka .Gdyby wyrażenie w nawiasie pomnożyć przez , a jednocześnie gradient ułamka podzielić przez ten sam wyraz, to uzyskalibyśmy zależność:

                            (10)

Powyższa zależność, co prawda oparta na dość dużym uproszczeniu (równości współczynników dyfuzji), pozwala zaproponować jak uzyskać stałość pierwszego członu w równaniu oznaczonym (*) w przypadku, gdy nie czynimy już uproszczeń, tj dla nierównych współczynników dyfuzji. Na tej podstawie równanie (*) można zapisać w postaci:

                            (11)

Przez analogię do wcześniejszych stwierdzeń można przyjąć, że pierwszy ułamkowy człon równania jest wielkością prawie stałą i nazywa się go zastępczym współczynnikiem dyfuzji, zaś drugi ułamek wyraża „siłę” napędową procesu, w której wielkość znajdująca się w mianowniku wskazuje, że strumień składnika jest tym mniejszy im większy jest mianownik yf lub inaczej wyrażenie w mianowniku działa hamująco na proces dyfuzji składnika A. Ostatecznie otrzymamy zależność na gęstość strumienia dyfuzji składnika A przez nieruchomą warstwę gazu w postaci:

                            (12)

Jeśli zachodzi potrzeba wyznaczenia gęstości strumienia dyfuzji składnika A pomiędzy przekrojami I i II, to powyższe równanie należy scałkować. W procesie ustalonym , zatem strumień pomiędzy przekrojami I i II odległymi o drogę s wynosi:

                            (13)

              dla:             

                            (14)

                            (15)

Średnia logarytmiczna definiowana jest następująco:

                            (16)

Ostatecznie otrzymamy zależność:

                            (17)

Możliwych jest wiele przypadków dyfuzji, które zostaną teraz omówione szczegółowo. Rozpatrzmy kolejno wybrane mechanizmy procesu dyfuzji i wyprowadźmy ostateczne wzory, z których można obliczyć gęstość strumienia dyfuzji składnika.

Przypadek I:

dyfuzja składnika A przez składnik inertny B

,              ,                ,    

,
 

                            (18)

Otrzymuje się kilka równoważnych zależności dla różnych sposobów wyrażania stężeń składnika A:

              (19)

              (20)

              (21)

              (22)

Określmy równanie, według którego zmienia się stężenie w funkcji drogi:

              Przekształćmy równanie:

                            (23)

                            (24)

                            (25)

                            (26)

Zmiana stężenia w procesie jednokierunkowej dyfuzji składnika A przez inert B odbywa się po krzywej exponencjalnej.

Przykład 1: