• Nie Znaleziono Wyników

Symulacja cyfrowa do wyznaczania stałych szybkości reakcji w warunkach miareczkowania w warstwie dyfuzyjnej na wirującej elektrodzie dysk-pierścień. Cz. 2, Reakcje drugiego rzędu - Biblioteka UMCS

N/A
N/A
Protected

Academic year: 2021

Share "Symulacja cyfrowa do wyznaczania stałych szybkości reakcji w warunkach miareczkowania w warstwie dyfuzyjnej na wirującej elektrodzie dysk-pierścień. Cz. 2, Reakcje drugiego rzędu - Biblioteka UMCS"

Copied!
8
0
0

Pełen tekst

(1)

ANNALES UNIVERSITATIS MARIAE CURIE-SKLODOWSKA LUBLIN-POLONIA

VOŁ.XŁII/XŁIII,R SECTIO AA 1987/1988 Zakład Chemii Analitycznej i Analizy Instrumentalnej

Wydział Chemii UMCS

Joanna BASAK, Jarosław PENAR, Kazimierz SYKUT

Symulacja cyfrowa do wyznaczania stałych szybkości reakcji w warunkach miareczkowania w warstwie dyfuzyjnej na winiącej

elektrodzie dyak-plerścień. Częś<? II. Reakcje drugiego rzędu

Digital Simulationfor Determining Rate Constants in Diffusion Layer Titrationon the Rotating Ring-disc Electrode. PartII. Second Order Reactions

WSTĘP

W poprzedniej pracy [1) przedstawiono metodę symulacji cyfrowej do ba­

dania procesu EC I (reakcja chemiczna pierwszego rzędu poprzedzona procesem elektrodowym) na wirującej elektrodzie dysk-pierścień oraz omówiono zasadnicze części napisanego w Fortranie programu do symulacji tego procesu. Obecnie zaj- miemy się rozszerzeniem powyższego programu tak, aby był on przydatny do ba­

dania kinetyki reakcji homogenicznej drugiego rzędu w warunkach miareczkowania w warstwie dyfuzyjnej (proces EC II).

Od strony analitycznej problem opisu procesu EC II został w sposób przy­

bliżony rozwiązany przez Albery’ego i in. (2, 3, 4]. W pracy |4] podano sposób wysnacsania stałej szybkości reakcji drugiego rzędu z krzywych zależności %r od io dla różnych szybkości wirowania elektrody oraz dla różnych stężeń substan­

cji C reagującej z substancją В generowaną na dysku. W skrócie metoda polega na wyznaczeniu tsw. kinetycznego prądu dysku »£>,*, tj. wartości prądu dysku w momencie gdy powierzchnia reakcji substancji В i C znajduje się na wewnętrznej krawędzi pierścienia [3,4] (powierzchnią reakcji określa się obszar, gdzie [B| = |C|).

Przy takim prądzie dysku przy stałej prędkości wirowania elektrody w (dla da­

nego stężenia substancji C) prąd pierścienia Ir jest uzależniony od ilości substancji В przenikającej w obszar dominacji substancji C, a więc w tych warunkach ir

(2)

44 Joanna Вмак, Jarosław Penar,К aa i mi er«Sykut

zależy od szybkości reakcji pomiędzy В i C. Mając wyznaczone wartości iD<k dla różnych w i dla różnych stężeń substancji C, określa się wartości kinetycznego prądu pierścienia 1я,к z krzywych ip vs. Îr. Mając te wartości można określić zależność Nk od w/Cc, gdzie Nk = |*K,fc/*D,fc| a Cp jest stężeniem substancji C. Zgodnie s teorią Albery’ego i in. [4] Nk wyraża się wzorem:

Nk = _ F(e)]wfc2co (1)

gdzie ri i r3 oznaczają promień dysku i promień dysku z przerwą izolacyjną, D jest współczynnikiem dyfuzji, u - lepkością kinematyczną roztworu a F(a) jest funkcją wymiarów elektrody [5]. Widać, że wykres zależności Nk od ш/Ср powinien być linią prostą przechodzącą przez początek układu. Jest to prawdą jedynie dla bardzo małych wartości ш/Ср. Dla większych wartości w/C wartość Nk nie rośnie nieograniczenie (jak to wynika ze wzoru) lecz asymptotycznie dąży do wartości No [4,5] (efektywności zbierania elektrody). Wyznaczenie z zależności Nk vs. ш/Ср polega na określeniu nachylenia prostoliniowego odcinka tej zależności, które jest proporcjonalne do l/«3.

Prater i Bard przedstawili metodę wyznaczania stałej szybkości fc2 opartą na symulacji cyfrowej [6]. Metoda polega na symulacyjnym wyznaczeniu krzywych roboczych Nk od x2m (gdzie x2 = к2С0льГ1и113О-11\0Ы)-'*13, m = С£/С°, Сд oznacza początkowe stężenie substancji elektroaktywnej) przy skoku potencjału na dysku do wartości prądu granicznego. Otrzymuje się w ten sposób (dla danej elektrody) rodzinę krzywych - dla każdej wartości parametru m jedna krzywa.

Ustalając w eksperymencie wartość ip na wartość graniczną ipj (graniczny prąd dyfuzyjno-konwekcyjny) należy rejestrować zmiany Nk (= |»hjd,j|) przy zmianie w i Cp. Wartość k2 można otrzymać przez wyznaczenie z krzywych roboczych wartości x2fn dla każdego Nk i wykreśleniu zależności x2fn °d ш/Ср. Nachylenie prostej jest wprost proporcjonalne do l/k2.

Naszym celem jest zaprezentowanie innego sposobu wykorzystania symulacji cyfrowej do wyznaczania k2. W symulacji przez nas proponowanej na dysk nakłada się prąd o wartości takiej, aby powierzchnia reakcji znajdowała się na wewnętrznej krawędzi pierścienia. Dla otrzymanych z analizy doświadczalnych krzywych Îr vs. ip wartości Nk [4] znajdujemy odpowiadające im wartości x2fn na roboczej krzywej symulacyjnej Nk od X2m. Następnie wykreślamy zależność x2"» od Cc/w (zależność prostoliniowa). Nachylenie prostej jest wprost proporcjonalne do l/k2.

Metoda ta ma kilka zalet w stosunku do metody Albery’ego i in. Między innymi charakteryzuje się większą dokładnością i jest użyteczna bez nakładania ograniczeń na wartości ш/С^. Zaletą jej jest także to, że dla danej elektrody mamy jedną krzywą roboczą Nk vs. x2m bez względu na wartość parametru m.

U Pratera i Barda dla każdego parametru m należy wykonać osobne obliczenia.

Prócz tego, jak wskazuje Albery i in. [4] dla wyznaczania k2 najkorzystniej jest яustalić” pozycję powierzchni reakcji na wewnętrznej krawędzi pierścienia, gdyż wtedy pierścień nsbiera” jedynie tę część substancji B, która nie zdążyła przerea- gować z substancją C w rejonie jej dominacji. Tak więc opisywana metoda łączy ze

(3)

Symulacja cyfrowa... II Reakcje drugiego rcedu 45

sobą lale tę metody Albery’ego i in. (praca w sakresie prądów kinetycinych) oral saletę metod symulacyjnych (brak ogranicseń sakresu stosowalności).

SYMULACJA CYFROWA PROCESU EC II ZE SKOKIEM PRĄDU NA DYSKU

Symulowany proces moina opisać równaniami:

A ± ne —» В dysk (2)

В + С —♦ D + E warstwa dyfusyjna (3)

B ± ne —» A pierścień (4)

gdsie sakłada się, ie substancje В, C i D są elektronieaktywne. W stosunku do procesu ЕС I powyisiy proces róini się jedynie reakcją (3). Do symulacji opisa­

nej popnednio naleiy kilka imian swiąsanych ■ П-nędowością reakcji w warstwie dyfusyjnej, skokiem prądu na dysku oras dyfusją i konwekcją substancji C. Wpro- wadsamy parametr m, określony laleinością [6]:

gdsie i Сд osnacsają pocsątkowe stęienia substancji C i A (stęienia substancji A i C w głębi rostworu). Pocsątkowe stęienie substancji C w kaśdej komórce siatki priestrsennej będsie równe [1]:

FC1( J, K) = m

Dyfuiję i konwekcję substancji C uwzględniamy analogicinie jak dyfuiję i konwekcję substancji A i В [1, 6]. Zmianie musi ulec parametr swiąsany s siybkością reakcji homogenicsnej. Dla reakcji II rsędu mamy:

- dCc/di = —dCB/dt = k2CB CG

Wyraiając to prsy pomocy parametrów symulacyjnych (tsn. Са/Сд = ЕС,

= ЕВ) otrzymamy [1, в]:

-ДГС = -A.FB = k2FB СсЫ

Ponieważ (1, 7] Де = tk/L = vV’p-V’w-^-^O.Sl)-3/3 а Сс = FC С°А, to —&FC = —&FB = kitbC^FB FC/L. Bezwymiarowy parametr k2tkC^

oinacimy XKTC [6]. Mamy:

XKTC = k2tkC^ = k2C^u~1i/1f3DA1,9(0.5i)-2/3 A więc ostatecsnie

FB1(J, K) = FB1(J, K)(l - XKTC FC1(J, K)/L)

(4)

46 Joanns Валак, Jarosław Psnar, KazimieraSykut

FC1(J, K) = FC1(J, K)(i - XKTC FB1(J, K)/L)

Ostatnia i mian a dotyczy warunków brzegowych przy skoku prądu na dysku. Na dysku mamy:

FA1(1,1) = FA1(1,1) - FLUX - DMA(FA1(\, 1) - FA1(2,1))

FB1(1,1) = FB1(1,1) + FLUX - DMB(FB1(1,1) - FBI(2,1)) /471(1,1) = /471(1,1) - DMC(FC1(1,1) - FC1(2,1))

W przerwie izolacyjnej i na pierścieniu warunki brzegowe są takie jak po­

przednio [1) (oczywiście należy uwzględnić dyfuzję substancji C do elektrody).

Parametr F LU X określa ilość substancji A jaka uległa reakcji elektrodowej na dysku w czasie At prsy założeniu, że na dysk nałożony jest prąd .kinetyczny”

Ponieważ ogólnie [7]

FLUX = DMA(FA1(2,1) - FA1(1,1)) = ZD(DMA/L)1'2

a [1, 7]

ZD = io/tO.Siy^nFAC^D^u1/2,,-1/0 (ZD jest prądem dysku wyrażonym przez wielkości bezwymiarowe) orai

»„ = .вгпГАС^Р2/3«1/2!/-1/«

(*»» jest prądem dysku w stanie stacjonarnym) to możemy napisać

»d/».. = (0.51)1/3/0.e2F££7X(£/DMA)1/2 czyli

FLUX = (*Ъ/»..)0.в2/(0.51)1/3(рА/А/£)1/з

Ponieważ t£> musi być równy prądowi .kinetycznemu” dysku, który spełnia zależność [2, 5]:

|.p,*| = M/(l-/*(«))

gdzie M = nFAD2/3(O.51)1/3w1/2iz-1/eC,0/1.288 a /’(«) jest funkcją zależną od wymiarów elektrody, przy czym a = (rj/ri)3 — 1, to ostatecznie możemy zapisać :

FLUX = [m/(l - F(a)]0.776(DA/A/£)1/2

(5)

Symulacja cyfrowa ... IIReakcjedrugiego rzędu 47

Ryc. 1. Wykree zależności Nk od XKTCm dla elektrody o wymiarach n = 0.3 cm, r2 = 0.33 cm, rs = 0.42 cm. o m = 0.00871, m = 0.004354, Д m = 0.006531, x

m = 0.001088.

WYNIKI I DYSKUSJA

Dla danej elektrody wartość ЛГ* jest taka sama dla układów s tą samą wartością x2fn> РГ1У czym prąd dysku musi być tak dobrany, aby powierschnia re­

akcji anajdowała się w tej samej pozycji dla każdego z tych układów. Innymi słowy układy > tą samą wartością x2m i tym samym położeniem powierzchni reakcji dają te same wartości Nk. Dzięki temu w symulacji procesu ЕС II prsy skoku prądu na dysku do wartości, pny której powierschnia reakcji osiąga wewnętrzną krawędś pierścienia, otrzymujemy jedną krsywą robocsą zależności Nk od ХКТС m, bes wsglądu na wartość m. Krsywą taką prsedstawia ryc.l .

Opisaną tutaj metodą uśyto do wyznaczenia stałej szybkości reakcji roda- nowania floroglucyny. Na ryc.2 przedstawiono zależność ХКТС m od C^/w, prsy czym wartości ХКТС m odpowiadają doświadczalnie wyznaczonym wartościom Nk.

Rycina 3 przedstawia eksperymentalną krsywą N od w/Cg. Zgodnie z zależnością Albery’ego, początkowa prostoliniowa część tej krzywej odpowiada

(6)

48 Jonin» Baiłfc, Jarosław Psnar,Kaaimisr» 8ykut

Ryc. 2. Wykres saleinofci XKTC m od С$,/ы [mol s/dm*rad] |—] Cg = 1.3 10“*

[mol/dm*], (o) Cg = 1.7 10"* [mol/dm*], (Д) Cg = 8.7 10^[mol/dm8]

■aleinoéci (1). W pracy [8] porównano teorię Albery’ego i in. ■ metodę symula- cji cyfrowej, ingerując imianę współcsynnika .339 we wiorie (1) na .424, dsięki ciemu wyniki symulacji są bardsiej i godne i teorią. Analiiując ryc. 2 i 3 widać więksiy lakres oru więkssą precyiją metody symulacji cyfrowej .

(7)

Symulacja cyfrowa ... IIReakcje drugiego rt«du 49

Ryc. 3. Wykres saleinoéci Nyodw/Ca [rad dm3/mol b|,[—j C" = 1.3 10 3 [mol/din3|, (A) Cg = 1.7 10~3 [mol/dm3), (o) Cg = 8.710~ł |mol/dm3|

LITERATURA

(1) PenarJ., Sykut K., Annale«Unir. Manat Cwie-SUodowtka, 43, 44, (1988/1989) [2] AlberyW.J., В rucken« tein S,Johnson D.C.; TVune Faraday Soc., 62, 1938, (1966) [3] Albery W.J., Brucken«tein S.;tbd.,62,2584,(1966)

[4] Albery W.J., HitchmanM.Ł.,Ulstrup .;ibid, 65, 1101, (1969).

[5] AlberyW.J., HitchmanM.L.; Finy dite elictrodtf, Claredon Pre»«,Oxford 1971.

[6] Prater K.B., Bard A.J.; J. Bectroclum. Soc., 117, (3),335, (1970).

[7] Prater K.B., BardA.J.; «h<f.ll7. (2), 207,(1970).

[8] Albery W.J., Drury J.S;J.Cham.Soc Faraday Trant.,68,(3),(1972).

SUMMARY

Thia work present« a new way of usage the digital simulation for second order kinetics study on rotating ring-discelectrode.

Złocone w Redakcji23II 1990

(8)

Cytaty

Powiązane dokumenty

Reakcję bromowania przeprowadzono ogrzewając molowe ilości 1,2,4-triazolu lub jego pochodnych z dwukrotnym nadmiarem bromu w lodowatym kwasie octowym.. (llb) i

nieniem suchą pozostałość zadano wodą i ekstrahowano bromową pochodną benzenem. Otrzymywanie

Oznaczając x = aktywność atomów wodoru w metalu, y = aktywność elektronów na powierzchni metalu, z = aktywność nie nasyconych drobin po przyłączeniu elektronów, ale

Gdy produkt generowany na dysku ulega w roztworze reakcji homogenicznej I rzędu, to prąd graniczny pierścienia zmienia się proporcjonalnie do zmian prądu dysku, przy

Celem niniejszej pracy jest szczegółowe opisanie metody Pratera i Barda oraz prezentacja zasadniczych części napisanego w języku Fortran programu do symulacji procesu, w

Badania chronaksymetryczne są cennym uzupełnieniem metody odruchowo-warunkowej i mogłyby znaleźć zastosowanie w warunkach klinicznych, szczególnie przy ustalaniu tła

wanie sferolitów w analizowanym przypadku może wiązać się z faktem, że tworząca się w reakcji substancja bezpostaciowa Cr2Mn3O9 posiada silne skłonności do adsorpcji

21 Giddens zauważa, że ta dysproporcja dopiero od niedawna jest na korzyść kobiet, do XIX wie- ku znacznie dłuższym życiem cieszyli się mężczyźni (2006: 181)... radzą sobie