Testovac´ı u´lohy pro ovˇeˇren´ı funkˇcnosti chemick´eho modulu programu Flow123d ˇ Luk´aˇs Zedek, Jan Sembera 20. prosinec 2010
Abstrakt Pˇredkl´ adan´ a zpr´ ava pˇredstavuje pˇrehled funkcionalit a v´ ysledky proveden´ ych test˚ u propojen´ı simulaˇcn´ıch softwar˚ u Flow123D pro transport kontaminant˚ u podzemn´ı vodou a Semchem pro simulaci chemick´ ych proces˚ u.
Obsah 1 Testy chemie
3
1.1
Poˇc´ ateˇcn´ı a okrajov´e podm´ınky, materi´alov´e vlastnosti, . . . . . . .
3
1.2
Modely s mobiln´ı porozitou . . . . . . . . . . . . . . . . . . . . .
4
1.2.1
Vypnut´ a chemie, referenˇcn´ı model . . . . . . . . . . . . .
4
1.2.2
Kinetika s nulov´ ymi kinetick´ ymi mocninami . . . . . . . .
5
1.2.3
Kinetika - implicitn´ı schema - dechlorace trichlorethenu .
6
1.2.4
Kinetika - explicitn´ı schema - dechlorace trichlorethenu .
7
1.2.5
Rovnov´ aha - protol´ yza vody . . . . . . . . . . . . . . . . .
8
1.2.6
Simulace obou typ˚ u reakc´ı . . . . . . . . . . . . . . . . . .
9
Modely s imobiln´ı porozitou . . . . . . . . . . . . . . . . . . . . .
10
1.3.1
Vypnut´ a chemie, referenˇcn´ı model . . . . . . . . . . . . .
10
1.3.2
Kinetika - im-/explicitn´ı schema - dechlorace trichlorethenu 11
1.3.3
Rovnov´ aha - protol´ yza vody . . . . . . . . . . . . . . . . .
12
1.3.4
Simulace obou typ˚ u reakc´ı . . . . . . . . . . . . . . . . . .
13
1.3
2 Testy transportu ovlivnˇ en´ eho chemi´ı
14
2.1
Poˇc´ ateˇcn´ı a okrajov´e podm´ınky, materi´alov´e vlastnosti . . . . . .
15
2.2
Model s mobiln´ı porozitou . . . . . . . . . . . . . . . . . . . . . .
16
2.2.1
Transport bez vlivu chemie . . . . . . . . . . . . . . . . .
16
2.2.2
Vliv kinetiky s nulov´ ymi kinetick´ ymi mocninami . . . . .
17
2.2.3
Vliv kinetiky - implicitn´ı schema - dechlorace trichlorethenu 18
2.2.4
Vliv kinetiky - explicitn´ı schema - dechlorace trichlorethenu 19
1
2.3
2.2.5
Vliv rovnov´ ahy - protol´ yza vody . . . . . . . . . . . . . .
20
2.2.6
Simulace obou typ˚ u reakc´ı . . . . . . . . . . . . . . . . . .
21
Model s dvoj´ı porozitou . . . . . . . . . . . . . . . . . . . . . . .
22
2.3.1
Transport bez vlivu chemie . . . . . . . . . . . . . . . . .
22
2.3.2
Vliv kinetiky - implicitn´ı schema - dechlorace trichlorethenu 23
2.3.3
Vliv kinetiky - explicitn´ı schema - dechlorace trichlorethenu 24
2.3.4
Vliv rovnov´ ahy - protol´ yza vody . . . . . . . . . . . . . .
25
2.3.5
Simulace obou typ˚ u reakc´ı . . . . . . . . . . . . . . . . . .
26
Pro poˇc´ıt´ an´ı chemie mus´ı b´ yt zapnut´ y transport.
2
Kapitola 1
Testy chemie 1.1
Poˇ c´ ateˇ cn´ı a okrajov´ e podm´ınky, materi´ alov´ e vlastnosti, . . .
Obr´ azek 1.1: Geometrie modelu pro simulaci reakc´ı bez transportu Pro testov´ an´ı funkˇcnosti chemick´eho modulu programu Flow123D byla vytvoˇrena sada testovac´ıch u ´loh s geometri´ı obsahuj´ıc´ı 3 tetraedrick´e elementy (viz obr´azek 1.1). Na tyto elementy byly zad´any takov´e okrajov´e podm´ınky proudˇen´ı, aby nedoch´ azelo k transportu (tj. na cel´e geometrii byl nastaven nulov´ y sp´ad piezomet-
3
rick´e v´ yˇsky). Spoleˇcnˇe pro vˇsechny testy chemie byly pouˇzity vstupn´ı soubory sit.mtr, sit.msh, sit.fbc, sit.ngh, sit chem.bcd, sit chem.tic. Testovac´ı u ´lohy chemie zahrnovaly simulaci: • kinetick´e reakce reprezuntuj´ıc´ı dechloraci trichlorethenu, 2M nO4− + C2 Cl3 H → 3Cl− + 2CO2 + 2M nO2 + H + s nulov´ ymi kinetick´ ymi mocninami reaguj´ıc´ıch l´atek, s nenulov´ ymi kinetick´ ymi mocninami, implicitn´ı schema, s nenulov´ ymi kinetick´ ymi mocninami, explicitn´ı schema, • rovnov´ aˇzn´e reakce pˇredstavuj´ıc´ı protol´ yzu vody, H + + OH − H2 O, • obou typ˚ u reakc´ı najednou.
1.2
Modely s mobiln´ı porozitou
Modely pracuj´ıc´ı s geometriemi, u kter´ ych byla uvaˇzov´ana pouze mobiln´ı porozita, byly poˇc´ıt´ any s reviz´ı 783 programu Flow123d.
1.2.1
Vypnut´ a chemie, referenˇ cn´ı model
Aby bylo v´ ysledky simulace chemick´ ych reakc´ı s ˇc´ım porovnat, byl vytvoˇren referenˇcn´ı model transportu kontaminant˚ u bez chemick´ ych reakc´ı. Jednotliv´e chemick´e l´ atky lze v grafu identifikovat podle poˇc´ateˇcn´ı podm´ınky. V poˇrad´ı od nejmenˇs´ı poˇc´ ateˇcn´ı podm´ınky k nejvyˇsˇs´ı pˇr´ısluˇs´ı jednotliv´e kˇrivky l´atk´am C2 Cl4 , M nO4− , C2 Cl3 H, H + , OH − , Cl− . Podle oˇcek´av´an´ı z˚ ustaly koncentrace sledovan´ ych l´ atek po celou dobu simulace rovn´e hodnotˇe dan´e pr´avˇe poˇc´ateˇcn´ı podm´ınkou (viz 1.2). cesta: . . ./Test ulohy/simulace/samotna chemie/mobil por/flow M CH off/ vstupn´ı soubory: flow M CH off.ini v´ ystupn´ı soubory: chemie M off.pos, transport M off.pos
4
Obr´ azek 1.2: Geometrie modelu pro simulaci reakc´ı bez transportu
1.2.2
Kinetika s nulov´ ymi kinetick´ ymi mocninami
Prvn´ı simulovanou chemickou rekc´ı byla dechlorace trichlorethenu manganistanem. Tato reakce m´ a kinetick´ y pr˚ ubˇeh. Pro jednoduchost byly mocniny v modelu kinetick´e reakce nejprve nastaveny na hodnotu 0. Tomu odpov´ıd´a oˇcek´avan´ y line´ arn´ı pokles (resp. n´ ar˚ ust) koncentrac´ı sledovan´ ych l´atek, kter´e se reakce u ´ˇcastn´ı. Spr´ avnost pˇredpokladu potvrzuje graf 1.3. Je v nˇem patrn´ y pokles koncentrace manganistanu a trichlorethenu, zat´ımco koncentrace Cl− roste. cesta: . . ./Test ulohy/simulace/samotna chemie/mobil por/flow M CH kin zero/ vstupn´ı soubory: flow M CH kin zero.ini v´ ystupn´ı soubory: chemie M kin zero.pos, transport M kin zero.pos
5
Obr´ azek 1.3: Geometrie modelu pro simulaci reakc´ı bez transportu
1.2.3
Kinetika - implicitn´ı schema - dechlorace trichlorethenu
Dechlorace trichlorethenu byla simulov´ana tak´e s nenulov´ ymi mocninami v modelu kinetick´e reakce. Pokles mnoˇzstv´ı reaktant˚ u a n´ar˚ ust koncentrace produktu se projevil i v tomto pˇr´ıpadˇe, ale pr˚ ubˇeh zmˇeny koncentrac´ı si nezachoval line´arn´ı charakter (viz 1.4). V´ yvoj koncentrac´ı byl v tomto pˇr´ıpadˇe poˇc´ıt´an za pouˇzit´ı implicitn´ıho v´ ypoˇcetn´ıho schematu. cesta: . . ./Test ulohy/simulace/samotna chemie/mobil por/flow M CH kin impl/ vstupn´ı soubory: flow M CH kin impl.ini v´ ystupn´ı soubory: chemie M kin impl.pos, transport M kin impl.pos
6
Obr´ azek 1.4: Geometrie modelu pro simulaci reakc´ı bez transportu
1.2.4
Kinetika - explicitn´ı schema - dechlorace trichlorethenu
Pr˚ ubˇeh v´ yvoje koncentrac´ı z´ıskan´ y simulac´ı dechlorace tetrachloethenu pomoc´ı explicitn´ıho schematu ukazuje obr´azek 1.5. V´ ysledky dosaˇzen´e s pouˇzit´ım implicitn´ıho a explicitn´ıho schematu se shoduj´ı. cesta: . . ./Test ulohy/simulace/samotna chemie/mobil por/flow M CH kin expl vstupn´ı soubory: flow M CH kin expl.ini v´ ystupn´ı soubory: chemie M kin expl.pos, transport M kin expl.pos
Obr´ azek 1.5: Geometrie modelu pro simulaci reakc´ı bez transportu
7
1.2.5
Rovnov´ aha - protol´ yza vody
Otestov´ an´ı simulace rovnov´ aˇzn´e reakce bylo provedeno na modelu protol´ yzy vody. Vzhledem k u ´vaze, kter´ a pˇredpokl´ad´a rychlost rovnov´aˇzn´e reakce za nekoneˇcnˇe velkou vzhledm k rychlosti transportu, se d´a oˇcek´avat skokov´a zmˇena koncentrac´ı reaktant˚ u a produk˚ u v prvn´ım kroku simulace. Graf 1.6 tento pˇredpoklad potvrzuje. Modr´e kˇrivky odpov´ıdaj´ıc´ı koncentarc´ım reaktent˚ u (H + , OH − ) se hned v prvn´ım kroku simulace ust´al´ı na hodnotˇe r˚ uzn´e od poˇc´ateˇcn´ı podm´ınky. cesta: . . ./Test ulohy/simulace/samotna chemie/mobil por/flow M CH kin eq vstupn´ı soubory: flow M CH eq.ini v´ ystupn´ı soubory: chemie M eq.pos, transport M eq.pos
Obr´ azek 1.6: Geometrie modelu pro simulaci reakc´ı bez transportu
8
1.2.6
Simulace obou typ˚ u reakc´ı
Obr´ azek 1.7 ukazuje v´ yvoj koncentrac´ı ovlivnˇen´ y obˇema typy dˇr´ıve samostatnˇe testovan´ ych chemick´ ych reakc´ı (dechlorac´ı a protol´ yzou). cesta: . . ./Test ulohy/simulace/samotna chemie/mobil por/flow M CH both vstupn´ı soubory: flow M both.ini v´ ystupn´ı soubory: chemie M both.pos, transport M both.pos
Obr´ azek 1.7: Geometrie modelu pro simulaci reakc´ı bez transportu
9
1.3
Modely s imobiln´ı porozitou
Modely pracuj´ıc´ı s geometriemi, u kter´ ych byla uvaˇzov´ana du´aln´ı porozita, byly poˇc´ıt´ any s reviz´ı 431 programu Flow123d.
1.3.1
Vypnut´ a chemie, referenˇ cn´ı model
Jako referenˇcn´ı model byla opˇet nejprvˇe provedena simulace bez chemick´ ych reakc´ı. Obr´ azek 1.8 ukazuje srovn´an´ı s v´ ysledky modelu bez chemie, spoˇc´ıtan´eho s geometri´ı s mobiln´ı porozitou. Postupn´ y pokles koncentrac´ı v mobiln´ıch p´orech aktu´ alnˇe uvaˇzovan´eho modelu je vyvol´an mechanismem difuze ˇr´ızen´ ym procesem v´ ymˇeny koncentrace s imobiln´ımi p´ory. cesta:. . ./Test ulohy/simulace/samotna chemie/imobil por/flow I CH off/ vstupn´ı soubory: flow I CH off.ini v´ ystupn´ı soubory: chemie I off.pos, transport I off.pos
Obr´ azek 1.8: Geometrie modelu pro simulaci reakc´ı bez transportu
10
1.3.2
Kinetika - im-/explicitn´ı schema - dechlorace trichlorethenu
V´ ysledky simulace dechlorace trichlorethenu poˇc´ıtan´e jak za pomoci implicitn´ıho, tak explicitn´ıho schematu a porovnan´e s koncentracemi modelu bez chemick´ ych reakc´ı zachycuje obr´azek 1.9. cesta: . . ./Test ulohy/simulace/samotna chemie/imobil por/flow I CH kin impl/ vstupn´ı soubory: flow I CH kin impl.ini v´ ystupn´ı soubory: chemie I kin impl.pos, transport I kin impl.pos
cesta: . . ./Test ulohy/simulace/samotna chemie/imobil por/flow I CH kin expl/ vstupn´ı soubory: flow I CH kin expl.ini v´ ystupn´ı soubory: chemie I kin expl.pos, transport I kin expl.pos
Obr´ azek 1.9: Geometrie modelu pro simulaci reakc´ı bez transportu
11
1.3.3
Rovnov´ aha - protol´ yza vody
Simulace protol´ yzy vody na geometrii s dvoj´ı porozitou ukazuje zˇretelnˇe vyˇcerp´an´ı iont˚ u OH − a pokles mnoˇzstv´ı iont˚ u H + v prvn´ım kroku simulace (viz 1.10). cesta: . . ./Test ulohy/simulace/samotna chemie/imobil por/flow I CH eq/ vstupn´ı soubory: flow I CH eq.ini v´ ystupn´ı soubory: chemie I eq.pos, transport I eq.pos
Obr´ azek 1.10: Geometrie modelu pro simulaci reakc´ı bez transportu
12
1.3.4
Simulace obou typ˚ u reakc´ı
V´ yvoj koncentrac´ı obou typ˚ u reakc´ı l´atek ovlivnˇen´ ych protol´ yzou a dechlorac´ı, bez transportu, ukazuje graf 1.11. cesta: . . ./Test ulohy/simulace/samotna chemie/imobil por/flow I CH both/ vstupn´ı soubory: flow I CH both.ini v´ ystupn´ı soubory: chemie I both.pos, transport I both.pos
Obr´ azek 1.11: Geometrie modelu pro simulaci reakc´ı bez transportu
13
Kapitola 2
Testy transportu ovlivnˇ en´ eho chemi´ı V r´ amci propojov´ an´ı softwar˚ u Flow123d a Semchem byla navrˇzena sada testovac´ıch u ´loh na provˇeˇren´ı funkˇcnosti propojen´ı a na zjiˇstˇen´ı ˇcasov´e n´aroˇcnosti v´ ypoˇct˚ u model˚ u reakˇcn´ıho transportu s jednotliv´ ymi druhy reakc´ı. Simulace transportu ovlivnˇen´eho chemick´ ymi reakcemi byla odzkouˇsena na modelech kolonov´eho experimentu s uvaˇzovanou jak pouze mobiln´ı, tak i dvoj´ı porozitou. Geometrie modelu obsahovala 45 tetreaedrick´ ych element˚ u (viz 2.1. Spoleˇcnˇe pro vˇsechny testy reakˇcn´ıho-transportu byly pouˇzity vstupn´ı soubory sit trans.mtr, sit transc.msh, sit trans.fbc, sit trans.ngh, sit trans.bcd, sit trans.tic. V´ ysledky jsou zaznamen´ any vˇzdy v podobˇe tˇr´ı graf˚ u k jednomu pˇr´ıkladu. Grafy zachycuj´ı v´ yvoje reakcemi a transportem ovlivnˇen´ ych koncentrac´ı v elementu na zaˇc´atku kolony (kontaminovan´ y u ´plnˇe), v elementu veprostˇred (kontaminovan´ y ˇc´asteˇcnˇe) a v elementu na konci kolony (vliv poˇc´ateˇcn´ı kontaminace a konvekce).
14
2.1
Poˇ c´ ateˇ cn´ı a okrajov´ e podm´ınky, materi´ alov´ e vlastnosti
Sp´ ad piezometrick´e v´ yˇsky na obou konc´ıch kolony ˇcinil 100 m a propustnost horniny v kolonˇe byla nastavena na 0.01 m · s−1 . D´elka kolony je 45 m. D´elka ˇcasu simulace je 200 s. Poˇc´ ateˇcn´ı podm´ınky byly nastaveny shodnˇe s modelem samostatn´ ych reakc´ı (viz 1.1).
Obr´ azek 2.1: Geometrie modelu pro simulaci reakcemi ovlivnˇen´eho transportu
15
2.2
Model s mobiln´ı porozitou
Modely pracuj´ıc´ı s geometriemi, u kter´ ych byla uvaˇzov´ana pouze mobiln´ı porozita, byly poˇc´ıt´ any s reviz´ı 783 programu Flow123D.
2.2.1
Transport bez vlivu chemie
V´ yvoj koncentrac´ı l´ atek ve zvolen´ ych elementech bez zahrnut´ı vlivu chemie je pro srovn´ an´ı zanesen v kaˇzd´em grafu. cesta: . . ./Test ulohy/simulace/transport/mobil por/flow M Tr off/ vstupn´ı soubory: flow M Tr off.ini v´ ystupn´ı soubory: sit M Tr off.pos, transport M Tr off.pos ˇ cas v´ ypoˇ ctu: 0:17.74
16
2.2.2
Vliv kinetiky s nulov´ ymi kinetick´ ymi mocninami
V´ ysledky modelu transportu kontaminant˚ u se zahrnut´ım kinetick´e reakce (dechlorace trichloethenu) s nulov´ ymi mocninami v matematick´em popisu ukazuj´ı grafy 2.2, 2.3 a 2.4. Je zˇrejm´e, ˇze na zaˇc´atku kolony m´a chemie v´ yrazn´ y vliv, zat´ımco v elementech bl´ızk´ ych konci kolony je pˇrevaˇzuj´ıc´ı vliv transportu. cesta: . . ./Test ulohy/simulace/transport/mobil por/flow M Tr kin zero/ vstupn´ı soubory: flow M Tr kin zero.ini v´ ystupn´ı soubory: sit M kin expl.pos, transport M Tr kin zero.pos ˇ cas v´ ypoˇ ctu: 2:46.82
Obr´ azek 2.2: Koncentrace v elementu na zaˇc´atku kolony
Obr´ azek 2.3: Koncentrace v elementu uprostˇred kolony
17
Obr´ azek 2.4: Koncentrace v elementu na konci kolony
2.2.3
Vliv kinetiky - implicitn´ı schema - dechlorace trichlorethenu
Zmˇeny koncentrac´ı l´ atek ve druh´em, dvac´at´em a ˇctyˇric´at´em elementu, ovlivnˇenn´e dechlorac´ı a poˇc´ıtan´e podle implicitn´ıho schematu jsou zaznamen´any v grafech 2.5, 2.6 a 2.7. cesta: . . ./Test ulohy/simulace/transport/mobil por/flow M Tr kin impl/ vstupn´ı soubory: flow M Tr kin impl.ini v´ ystupn´ı soubory: sit M kin impl.pos, transport M Tr kin impl.pos ˇ cas v´ ypoˇ ctu: 10:45.75
Obr´ azek 2.5: Koncentrace v elementu na zaˇc´atku kolony
18
Obr´ azek 2.6: Koncentrace v elementu uprostˇred kolony
Obr´ azek 2.7: Koncentrace v elementu na konci kolony
2.2.4
Vliv kinetiky - explicitn´ı schema - dechlorace trichlorethenu
Model reakˇcn´ıho transportu s kinetikou poˇc´ıtanou podle explicitn´ıho schematu poskytl v´ ysledky, kter´e jsou zaneseny do graf˚ u 2.8, 2.9 a 2.10. Je zˇrejm´e, ˇze v´ ypoˇcet podle explicitn´ıho schematu probˇehl rychleji neˇz v´ ypoˇcet stejn´e u ´lohy podle schematu implicitn´ıho. cesta: . . ./Test ulohy/simulace/transport/mobil por/flow M Tr kin expl/ vstupn´ı soubory: flow M Tr kin expl.ini v´ ystupn´ı soubory: sit M Tr kin expl.pos, transport M Tr kin expl.pos ˇ cas v´ ypoˇ ctu: 2:36.91
19
Obr´ azek 2.8: Koncentrace v elementu na zaˇc´atku kolony
Obr´ azek 2.9: Koncentrace v elementu uprostˇred kolony
2.2.5
Vliv rovnov´ ahy - protol´ yza vody
Model transportu zohledˇ nuj´ıc´ıho protol´ yzu vody mˇel v´ ysledky zobrazen´e v grafech 2.11, 2.12 a 2.13. cesta: . . ./Test ulohy/simulace/transport/mobil por/flow M Tr eq vstupn´ı soubory: flow M Tr eq.ini v´ ystupn´ı soubory: sit M Tr eq.pos, transport M Tr eq.pos ˇ cas v´ ypoˇ ctu: 10:50.68
20
Obr´ azek 2.10: Koncentrace v elementu na konci kolony
Obr´ azek 2.11: Koncentrace v elementu na zaˇc´atku kolony
2.2.6
Simulace obou typ˚ u reakc´ı
Vliv zohlednˇen´ı rovnov´ aˇzn´e i kinetick´e reakce v modelu transportu kontaminant˚ u ukazuj´ı grafy 2.14, 2.15 a 2.16. cesta: . . ./Test ulohy/simulace/transport/mobil por/flow M Tr both/ vstupn´ı soubory: flow M Tr both.ini v´ ystupn´ı soubory: sit M Tr both.pos, transport M Tr both.pos ˇ cas v´ ypoˇ ctu: 43:31.95
21
Obr´ azek 2.12: Koncentrace v elementu uprostˇred kolony
Obr´ azek 2.13: Koncentrace v elementu na konci kolony
2.3
Model s dvoj´ı porozitou
Modely pracuj´ıc´ı s geometriemi, u kter´ ych byla uvaˇzov´ana du´aln´ı porozita, byly poˇc´ıt´ any s reviz´ı 431 programu Flow123d.
2.3.1
Transport bez vlivu chemie
V´ ystup referenˇcn´ıho modelu transportu kontaminant˚ u bez simulace chemick´ ych reakc´ı je souˇc´ ast´ı kaˇzd´eho grafu. cesta: . . ./Test ulohy/simulace/transport/mobil por/flow I Tr off/ vstupn´ı soubory: flow I Tr off.ini v´ ystupn´ı soubory: sit I Tr off.pos, transport I Tr off.pos ˇ cas v´ ypoˇ ctu: . . .
22
Obr´ azek 2.14: Koncentrace v elementu na zaˇc´atku kolony
Obr´ azek 2.15: Koncentrace v elementu uprostˇred kolony
2.3.2
Vliv kinetiky - implicitn´ı schema - dechlorace trichlorethenu
V´ ysledky modelu transportu kontaminant˚ u se zahrnut´ım kinetick´e reakce, imobiln´ı porozity a poˇc´ıtan´eho podle implicitn´ıho numerick´eho schematu, ukazuj´ı grafy 2.17, 2.18 a 2.19. cesta: . . ./Test ulohy/simulace/transport/mobil por/flow I Tr kin impl/ vstupn´ı soubory: flow I Tr kin impl.ini v´ ystupn´ı soubory: sit I kin impl.pos, transport I Tr kin impl.pos ˇ cas v´ ypoˇ ctu: . . .
23
Obr´ azek 2.16: Koncentrace v elementu na konci kolony
Obr´ azek 2.17: Koncentrace v elementu na zaˇc´atku kolony
2.3.3
Vliv kinetiky - explicitn´ı schema - dechlorace trichlorethenu
Pouˇzit´ım explicitn´ıho numerick´eho schematu pro v´ ypoˇcet zmˇeny koncentrac´ı dechlorac´ı v modelu transportu s imobiln´ı porozitou jsme dospˇeli k v´ ysledk˚ um zobrazen´ ym v grafech 2.20, 2.21 a 2.22. cesta: . . ./Test ulohy/simulace/transport/mobil por/flow I Tr kin expl/ vstupn´ı soubory: flow I Tr kin expl.ini v´ ystupn´ı soubory: sit I Tr kin expl.pos, transport I Tr kin expl.pos ˇ cas v´ ypoˇ ctu: . . .
24
Obr´ azek 2.18: Koncentrace v elementu uprostˇred kolony
Obr´ azek 2.19: Koncentrace v elementu na konci kolony
2.3.4
Vliv rovnov´ ahy - protol´ yza vody
V´ ysledky modelu transportu kontaminant˚ u se zahrnut´ım rovnov´aˇzn´e reakce a mobiln´ı porozity jsou zobrazeny v grafech 2.23, 2.24 a 2.25. cesta: . . ./Test ulohy/simulace/transport/mobil por/flow I Tr eq vstupn´ı soubory: flow I Tr eq.ini v´ ystupn´ı soubory: sit I Tr eq.pos, transport I Tr eq.pos ˇ cas v´ ypoˇ ctu: . . .
25
Obr´ azek 2.20: Koncentrace v elementu na zaˇc´atku kolony
Obr´ azek 2.21: Koncentrace v elementu uprostˇred kolony
2.3.5
Simulace obou typ˚ u reakc´ı
V´ ystup modelu transportu kontaminant˚ u rozˇs´ıˇren´eho o simulaci rovnov´aˇzn´ ych i kinetick´ ych reakc´ı ukazuj´ı grafy 2.26, 2.27 a 2.28. cesta: . . ./Test ulohy/simulace/transport/mobil por/flow I Tr both/ vstupn´ı soubory: flow I Tr both.ini v´ ystupn´ı soubory: sit I Tr both.pos, transport I Tr both.pos ˇ cas v´ ypoˇ ctu: . . .
26
Obr´ azek 2.22: Koncentrace v elementu na konci kolony
Obr´ azek 2.23: Koncentrace v elementu na zaˇc´atku kolony
Obr´ azek 2.24: Koncentrace v elementu uprostˇred kolony
Obr´ azek 2.25: Koncentrace v elementu na konci kolony
27
Obr´ azek 2.26: Koncentrace v elementu na zaˇc´atku kolony
Obr´ azek 2.27: Koncentrace v elementu uprostˇred kolony
Obr´ azek 2.28: Koncentrace v elementu na konci kolony
28