ČVUT v Praze Fakulta stavební Studentská vědecká a odborná činnost Akademický rok 2005/2006
Transport a akumulace solí v porézních stavebních materiálech
Jméno a příjmení studenta :
Jan Bažil
Ročník, obor :
4., Konstrukce a dopravní stavby
Vedoucí práce :
Ing. Z. Pavlík, Ph.D.
Katedra :
Stavební mechaniky
1. Úvod Návrh stavebních konstrukcí, ať jsou to stavby pozemní, podzemní či vodohospodářské, vyžaduje spolupráci mnoha technických a průmyslových oborů. Prvotním cílem při návrhu každé konstrukce je zajištění její funkčnosti a bezpečnosti během celé doby plánované životnosti. Ukazuje se, že právě nízká životnost bývá největším problémem mnoha staveb. Například mosty pozemních komunikací jsou navrhovány na období sto let. Velmi často se však stává, že mosty, které nejsou ani v polovině své životnosti jsou ve stavu, který v lepším případě vyžaduje rozsáhlou rekonstrukci, která se vynaloženými prostředky téměř vyrovná stavbě mostu nového, v horším případě je nutné přistoupit k demolici takovéhoto objektu. Častou příčinou snížení životnosti konstrukcí je působení roztoků solí ve struktuře jakéhokoliv typu stavebních materiálů. Typickým příkladem solí, které ohrožují stavební konstrukce jsou chloridy (sodný, vápenatý), dusičnany, sírany, hydrogenuhličitany a ionty amonné, sodné, vápenaté a hořečnaté. Ostatní ionty se vyskytují pouze vyjímečně. Roztoky solí se mohou dostat do struktury materiálu například kapilárním vzlínáním podzemních vod, solením vozovek a chodníků v zimním období, působením kyselých dešťů a živých organismů (zvířecí exkrementy, mikroorganismy) atd. Soli jsou do vnitřní struktury materiálu transportovány pouze rozpuštěné ve vodě. V případě, že jsou soli ve struktuře materiálu ve formě vodného roztoku, nijak materiál neohrožují. Nebezpečí pro strukturu materiálu nastává v okamžiku, kdy se voda z roztoku odpaří a sůl začne krystalizovat. Pokud sůl krystalizuje na povrchu materiálu, dochází sice k zhoršení estetického působení, ale samotná konstrukce není bezprostředně ohrožena. Problém nastává až v případě, kdy soli krystalizují uvnitř pórového systému materiálu, protože krystalizační tlaky, v závislosti na velikosti a uspořádání porézní struktury materiálu, mohou přesáhnout pevnost materiálu, a tím způsobit jeho destrukci. Pokud se zaměříme na betonové konstrukce, pak jejich životnost nejvíce ohrožují chloridy a to hned dvojím způsobem. Jednak degradují beton samotný, jednak způsobují korozi ocelové výztuže. Koroze ocelové výztuže v betonu je závislá na koncentraci iontů, např. hydroxidových, chloridových, v pórovém roztoku, a také na difúzi kyslíku z okolní atmosféry a vlhkosti krycí vrstvou. Ocelová výztuž v betonu na svém povrchu vytváří v prostředí s vysokou koncentrací hydroxidových iontů pasivační vrstvu z oxidů železa, které tvoří ochrannou bariéru proti dalšímu pronikání agresivních iontů z roztoku, čímž dojde k zastavení koroze výztuže.
Obr. 1 - Degradace mostní římsy vlivem působení chloridů. 1
Transport iontů solí porézní strukturou stavebního materiálu je většinou autorů považován za difúzní proces. Je-li koncentrace iontů na povrchu materiálu vyšší než v jeho struktuře, pak koncentrační gradient vyvolává tok iontů směrem do porézního materiálu. Na základě stanovení difúzních koeficientů je možno předpovědět transport solí porézní strukturou stavebních materiálů a tím předejít jejich případné degradaci a prodloužit životnost stavebních konstrukcí.
2. Základní modely transportu solí v porézních materiálech Většina používaných modelů popisujících difúzi chloridových iontů je jednoduchá a je založena na použití Fickovy difúzní rovnice pro koncentraci chloridů s konstantním difúzním součinitelem a s konstantními počátečními a okrajovými podmínkami ∂C = D ⋅ (div grad (C )), ∂t C (0, t ) = C 0 ,
(1)
C (∞ , t ) = 0, C (0, x ) = 0,
kde
(
)
D je součinitel difúze m 2 / s ,
C je koncentrace roztoku (kg / m3 ), C0 je koncentrace na exponovaném povrchu (kg / m 3 ), x je vzdálenost od exponovaného povrchu (m ).
V případě, že uvažujeme pouze jednorozměrný transport, přejde rovnice (1) do tvaru
∂ 2C ∂C ∂ ⎛ ∂C ⎞ = D⋅ ⎜ ⎟ = D⋅ 2 . ∂x ∂t ∂x ⎝ ∂t ⎠
(2 )
Tato rovnice má jednoduché řešení ve tvaru (např. [1])
⎡ ⎛ x C ( x , t ) = C 0 ⎢1 − erf ⎜ ⎝ 2 Dt ⎣
⎞⎤ ⎟⎥, ⎠⎦
(3)
což je důvodem častého použití tohoto typu modelů. Největším nedostatkem těchto modelů je předpoklad konstantního součinitele difúze soli a zanedbání vlivu transportu vody na transport chemických sloučenin. Z naměřených koncentračních profilů nemůže být téměř nikdy získána jediná hodnota difúzního součinitele, zvláště pokud jsou tato měření prováděna v delším časovém úseku, kdy se difúzní součinitel poté jeví jako funkce času. Přes uvedené nedostatky našly tyto modely uplatnění v mnoha praktických aplikacích, neboť lze vypočítané hodnoty difúzního součinitele porovnat pro různé typy materiálů i pro vzorky materiálů vystavené rozdílným podmínkám okolního prostředí s odlišnou charakteristikou.
2
Za opačný extrém popisu transportu solí porézním materiálem lze považovat komplexní modely zahrnující při popisu transportu solí všechny sekundární vlivy, které ovlivňují kombinovaný transport tepla, vlhkosti a solí. Mezi tyto propracované modely patří např. model Grunewaldův [2]. Ač by se na první pohled mohlo zdát, že takto komplexní modely představují ideální řešení, není tomu tak. Je třeba experimentálně stanovit celou řadu vstupních parametrů, jejichž měření je časově náročné a pro některé materiály téměř neuskutečnitelné, protože je třeba naměřit tyto parametry v závislosti na teplotě, vlhkosti a koncentraci solí. Tyto modely je tedy rozumné použít v případě, že potřebné vstupní parametry máme již k dispozici. Jistý kompromis mezi těmito dvěma postupy popisu transportu solí představuje model založený na nelineární difúzní rovnici ve tvaru ∂C = div (D (C ) ⋅ div grad (C )), ∂t
(4 )
s počátečními a okrajovými podmínkami tvaru C (0, t ) = C1 ,
C (∞, t ) = C2 , C (0, x ) = C2 .
Za předpokladu, že transport solí probíhá na principu difúze, lze provést výpočet součinitele difúze v závislosti na koncentraci solí z naměřených koncentračních profilů obdobnými inverzními metodami jako při určení součinitele vlhkostní vodivosti v závislosti na vlhkosti nebo součinitele tepelné vodivosti v závislosti na teplotě [3]. Mezi nejjednodušší metody řešení tohoto problému patři Boltzamann-Matanův [4] princip založený na aplikaci Boltzmannovy transformace ve tvaru
η=
x
, 2 t ω (η ) = C ( x, t ).
(5) (6)
Za podmínky konstantních počátečních podmínek a Dirichletových okrajových podmínek na obou koncích vzorku a předpokladu jednorozměrného transportu solného roztoku je poté parabolická difúzní rovnice transformována do tvaru ∂C ∂C ∂η dC 1 = ⋅ = ⋅ , ∂x ∂η ∂x dη 2 t ∂C ∂C ∂η dC ⎛ x = ⋅ = ⋅⎜− ∂t ∂η ∂t dη ⎝ 4t t
(7 ) x ⎞ dC ⎛ ⋅⎜− ⎟= ⎠ dη ⎝ 2 t
⎞ 1 dC ⎛ η ⎞ ⋅ ⎜ − ⎟. ⎟⋅ = ⎠ 2t dη ⎝ 2t ⎠
(8)
Úpravou vztahu (8) získáme obyčejnou diferenciální rovnici ve tvaru d dη
⎛ dω ⎞ dω ⋅ ⎜⎜ D = 0, ⎟⎟ + 2η ⋅ dη ⎝ dη ⎠
(9)
s podmínkami 3
ω (0) = C1 , ω (∞ ) = C2 . Předpokládáme-li, že v daném čase t 0 = konst. známe rozložení koncentrace solí v materiálu C ( x,t 0 ) , můžeme použít druhou transformaci ve tvaru
z = 2η ⋅ t 0 ,
(10)
C (z ) = ω (η ).
(11)
Dosazením této transformace do rovnice (7) získáme rovnici
d ⎛ dC ⎞ z dC ⋅⎜D ⋅ =0 ⎟+ dz ⎝ dz ⎠ 2t 0 dz
(13 )
s podmínkami C (0) = C1 ,
C (∞ ) = C2 . Řešením rovnice (13) získáme vztah pro výpočet součinitele difúze solí v závislosti na koncentraci solí. D (C 0 ) =
∞
1 ⎛ dC ⋅ ∫⎜z⋅ ⎛ dC (z = z 0 )⎞⎟ z0 ⎝ dz 2t 0 ⋅ ⎜ ⎝ dz ⎠
⎞ ⎟ dz . ⎠
(14 )
Dalším krokem k přesnějšímu popisu transportu solí porézní strukturou materiálů je uvažování vlivu transportu vody na transport solí a vlivu vázání iontů solí na stěny porézního prostoru. Model zahrnující tyto vlivy je označován jako difúzně– advektivní [5]. Transport solného roztoku je popsán systémem dvou parabolických rovnic, a to bilancí hmotnosti solí vyjádřenou jako ∂ (wC ) = div (wD ⋅ grad (C f ∂t
)) − div (C
f
)
⋅v −
∂C b , ∂t
(15 )
a bilancí hmotnosti vlhkosti ve tvaru ∂w = div (κ ⋅ grad (w )), ∂t
kde
(16 )
(
)
D je součinitel difúze solí m 2 / s ,
(
)
C f je koncentrace volných iontů ve vodě kg / m3 ,
(
)
Cb je koncentrace vázaných iontů na stěnách pórů matriálu kg / m 3 , 4
w objemový obsah vlhkosti (m3 / m3 ),
κ je součinitel vlhkostní vodivosti (m 2 / s ), v je Darcyho rychlost (m / s ).
Darcyho rychlost kapalné fáze je dána vztahem
v = −κ ⋅ grad (w).
(17 )
Vstupní parametry C f , Cb , w v rovnicích (15) a (16) určíme experimentálně, zbylé dvě materiálové charakteristiky určíme pomocí inverzní analýzy. Problém opět zjednodušíme na případ jednorozměrného transportu solného roztoku. Rovnice (15) a (16) poté převedeme na tvar ∂C f ∂ (wC ) ∂ ⎛ ⎜⎜ wD ⋅ = ∂t ∂x ⎝ ∂x ∂w ∂ ⎛ ∂w ⎞ = ⎜κ ⎟. ∂x ∂x ⎝ ∂x ⎠
(
)
⎞ ∂ ∂C b ⎟⎟ − Cf ⋅v − , ∂t ⎠ ∂x
(18 ) (19 )
Pro řešení soustavy dvou parabolických rovnic použijeme opět Boltzmann–Matanovy metody s předpokladem konstantních počátečních podmínek a Dirichletových okrajových podmínek. Po aplikaci Boltzmannovy transformace na rovnice (18) a (19) dostáváme
2
dC d ⎛ dw ⎞ d ⎛ ⎜⎜ D ⋅ w f ⎟⎟ + 2 ⎜⎜ C f ⋅ κ dη ⎝ dη ⎠ dη ⎝ dη
2
d ⎛ dw ⎞ dw ⎜⎜ κ ⎟⎟ + η = 0. dη ⎝ dη ⎠ dη
d (w ⋅ C f ) ⎞ dCb dC f ⎟⎟ + η +η = 0, dη dC f dη ⎠
(20) (21)
Druhou Boltzmannovu transformaci ve tvaru z = η t0 provedeme
(22) za
podmínky,
že
v čase
t = t0
známe
w( x, t 0 ), C f ( x, t 0 ), Cb ( x, t 0 ) .
Dostáváme dw ⎞ d ⎛ ⎜ w ⋅ C f ⋅κ ⎟+2 dz ⎠ dz ⎝
dC f ⎛ ⎜⎜ D ⋅ w dz ⎝
2
d dz
2
dw ⎞ z dw d ⎛ = 0. ⎜ κ (z 0 ) ⎟+ dz ⎠ t 0 dz dz ⎝
⎞ z d (w ⋅ C f ) z dC b dC f ⎟⎟ + ⋅ + ⋅ = 0, dz t 0 dC f dz ⎠ t0
(23 )
(24 )
Z řešení rovnice (24) určíme součinitel vlhkostní vodivosti v závislosti na vlhkosti
5
∞
1 ⎛ dw ⎞ ⋅ ∫⎜z κ (z0 ) = ⎟dz , ⎛ dw ⎞ z0 ⎝ dz ⎠ ( z 0 )⎟ 2t 0 ⎜ ⎝ dz ⎠
(25)
kde κ (z 0 ) = κ (w0 , C f 0 ), w0 = w( z0 , t0 ), C f 0 = C f ( z0 , t0 ) . Dosazením vztahu (25) do rovnice (23) a řešením této rovnice získáváme vztah pro výpočet součinitele difúze solí v závislosti na koncentraci solí.
D ( z0 ) = −
C f ( z0 ) ⋅ κ ( z0 ) w( z0 )
dC f dz
dw (z0 ) dz +
( z0 )
∞ ⎛ d (w ⋅ C f ) dCb dC f ⎞ ⎟dz , ⋅ ∫⎜ + ⎜ dz ⎟ dC f dC dC f b ⎠ ( z0 ) z0 ⎝ 2t0 ⋅ w( z0 ). dz
1
D ( z 0 ) = D (w0 , C f 0 ) .
(26)
(27)
3. Popis zkoumaného materiálu V této práci jsem se zabýval studiem transportních a akumulačních parametrů solí pro materiál na bázi kalcium silikátu. Jedná se o materiál, který se běžně používá jako vnitřní tepelná izolace. Vyrábí se chemickou reakcí oxidu vápenatého, oxidu křemičitého, příměsí a vody v autoklávu při vysokých tlacích. Tento výrobní proces umožňuje vyrobit materiál s velmi jemnou porézní strukturou zcela homogenního výsledného složení: hydrát kalcium silikátu 75 – 90%, minerální plnivo 0 – 20%, celulóza 3 – 6%. Pro přesný popis studovaného materiálu jsem provedl stanovení základních materiálových charakteristik.
4. Popis experimentu K měření základních materiálových parametrů byly připraveny vzorky o rozměrech 50x50x50 mm, které byly vysušeny ve vakuové sušárně při 65 °C, aby nedošlo k rozkladu hydrátu kalcium silikátu. Vzorky byly opakovaně váženy na analytické váze s přesností 0,001 g. Po ustavení hmotnosti vzorku na ± 0,01g, tzv. počáteční hmotnost vzorku, byly vzorky umístěny do nádoby s destilovanou vodou tak, aby byly zcela ponořené, a ze sestrojené aparatury byl pomocí vakuové pumpy odčerpán vzduch. Po třech hodinách bylo vakuum odpojeno a vzorky byly ponechány dalších 24 hodin zcela ponořené ve vodě. Zvážením vzorků byla zjištěna tzv. nasycená hmotnost vzorku, a jejich následným ponořením na závěsném zařízení do vody a zvážení byla stanovena tzv. Archimedova hmotnost. Ze zjištěných hmotností byl nejprve vypočítán objem vzorku dle rovnice V=
msat − ma
ρH O
(m3),
(28)
2
6
kde ρ H 2O je hustota vody (kg/m3), msat je hmotnost zcela nasyceného vzorku (kg), ma je hmotnost nasyceného vzorku ponořeného ve vodě (kg). S použitím vypočítaného objemu (28) byla pro každý vzorek stanovena objemová hmotnost, hustota matrice a otevřená pórovitost dle následujících vztahů
ρ=
m0 (kg/m3), V
ρ mat = ψ=
(29)
m0 (kg/m3), V (1 −ψ )
(30)
msat − m0 (-), Vρ H 2O
(31)
kde m0 je hmotnost vysušeného vzorku (kg). Pro aplikaci difúzně-advektivního modelu Beara a Bachmata [5] bylo nezbytné experimentálně stanovit profily obsahu vlhkosti, koncentrace solí a vazebnou izotermu solí. Nejprve byly nařezány vzorky o rozměru 40x40x160 mm, které byly na podélných stranách parotěsně a vodotěsně odizolovány epoxidovou pryskyřicí. Po vytvrdnutí izolace byly vzorky umístěny do vakuové sušárny a vysušeny při 65 °C. Vysušené vzorky se nechaly volně vychladnout na vzduchu. Poté byly upevněny do stojanu a neizolovanou stranou byly 1 – 2 mm ponořeny do penetrujícího média. Jako penetrující médium byla zvolena destilovaná voda a 1M roztok chloridu sodného, viz. Obr. 2. V daných časových intervalech (30, 60 a 90 minut) byly vzorky po 20 mm rozřezány, epoxidová pryskyřice odstraněna a takto upravené vzorky byly ve váženkách zváženy a následně vysušeny při teplotě 110°C do ustálené hmotnosti. Ze zjištěných hmotností byl vypočítán obsah vlhkosti pro každý vzorek dle rovnic u=
msat − m0 (kg/kg), m0
(32)
w=
msat − m0 ρ (m3/m3), m0 ρ H 2o
(33)
kde
m0 je hmotnost vysušeného vzorku (kg), msat je hmotnost nasyceného vzorku (kg), ρ je objemová hmotnost materiálu (kg/m3), ρ H 2o je hustota vody (kg/m3).
Vzorky, které byly nasyceny roztokem chloridu sodného byly po vysušení a zvážení ponořeny do váženky s 200 ml vroucí destilované vody, jak ukazuje Obr.3. Váženky byly poté parotěsně utěsněny páskou a ponechány po dobu dvou měsíců
7
v laboratoři při teplotě 25 ± 2°C. Poté byl ve výluhu stanoven obsah chloridů pomocí iontově selektivní elektrody přístrojem inoLab pH/ION 740, Obr. 4. Vazebná izoterma je důležitým vstupem do matematických modelů popisujících transport solných roztoků strukturou stavebních materiálů. Pro určení vazebné isotermy iontů, byla využita modifikovaná metoda Tanga a Nilssona [6, 7], kteří navrhli postup založený na adsorpci iontů z roztoku o známé koncentraci. Vysušené vzorky byly umístěny do váženek s 200 ml roztoku chloridu sodného o známé koncentraci v rozsahu 0,01M – 1M. Váženky byly parotěsně utěsněny páskou a ponechány v laboratoři při 25 ± 2°C. Po ustavení rovnováhy byl stanoven obsah chloridů ve výluhu pomocí iontově selektivní elektrody přístrojem inoLab pH/ION 740, který odpovídá obsahu volných chloridů c f (mg/l). Po té byl vypočítán obsah vázaných chloridů dle rovnice cb =
kde
V (c0 − c1 ) ρ (kg/m3 vzorku), m0
(28)
V je objem roztoku soli (m3), c0 , c1 je počáteční a rovnovážná koncentrace roztoku chloridů (kg/m3) m0 je hmotnost vysušeného vzorku (kg), ρ je měrná hmotnost materiálu (kg/m3). Měřením roztoků s různými počátečními koncentracemi solí c0 jsou získány
( )
body funkce cb = cb c f , tzv. vazebná izoterma iontů.
Obr. 2 Nasákání vzorků.
8
Obr. 3 Vážení a vyluhování vzorků.
Obr. 4 Přístroj inoLab pH/ION 740 s iontově selektivní elektrodou.
9
5. Výsledky a diskuse Výsledky měření základních materiálových parametrů jsou prezentovány v Tab. 1. Ze získaných výsledků je zřejmé, že se jedná o vysoce porézní materiál, u kterého je možné očekávat vysokou kapilární aktivitu z pohledu transportu vlhkosti a solného roztoku. Zjištěná objemová hmotnost koresponduje s hodnotou otevřené pórovitosti. Tab. 1 Základní materiálové parametry Testovaný materiál
Objemová hmotnost
Hustota matrice (kg/m3)
230,6
Kalcium silikát
Otevřená pórovitost (-) 0,87
1986,3
Experimentální výsledky stanovení vlhkostních profilů a celkových koncentračních profilů solí demonstrují Obr. 5-9. Obr. 10 prezentuje experimentálně stanovenou vazebnou izotermu chloridů pro kalcium silikát. 3.50 0.5 hod 1 hod 1.5 hod
Obsah vlhkosti [kg/kg]
3.00 2.50 2.00 1.50 1.00 0.50 0.00 0
0.02
0.04
0.06
0.08
0.1
0.12
0.14
Vzdálenost [m]
Obr. 5 Profil obsah vlhkosti pro penetraci destilované vody.
10
0.16
3.5 0.5 hod 1 hod 1.5 hod
Obsah vlhkosti [kg/kg]
3 2.5 2 1.5 1 0.5 0 0
0.02
0.04
0.06
0.08
0.1
0.12
0.14
0.16
Vzdálenost [m]
Obr. 6 Profil obsah vlhkosti pro penetraci 1 M roztoku NaCl.
3
3
Objemová vlhkost [m /m ]
0.8
0.6
0.4
0.2
0 0
0.0005
0.001
0.0015
0.002
0.0025
0.003
0.0035
0.004
x/√t [m/s1/2]
Obr. 7 Vlhkostní profil v Boltzmannově formě pro nasákání destilovanou vodou.
11
3
3
Objemová vlhkost [m /m ]
0.8
0.6
0.4
0.2
0 0
0.0005
0.001
0.0015
0.002
0.0025
0.003
0.0035
0.004
x/√t [m/s1/2]
Obr. 8 Vlhkostní profil v Boltzmannově formě pro nasákání 1M roztoku NaCl.
30.0
3
Obsah soli [kg/m ]
25.0 20.0 15.0 10.0 5.0 0.0 0
0.0005
0.001
0.0015
0.002
0.0025
0.003
0.0035
0.004
x/√t [m/s1/2]
Obr. 9 Profil celkového obsahu solí vyjádřený v Boltzmannově formě.
12
140
Cb [kg/m
3 (vzorku)]
120 100 80 60 40 20 0 0
50
100
150
200
Cf [kg/m3(roztoku)]
Obr. 10 Vazebná izoterma chloridů pro kalcium silikát. Experimentálně stanovené profily vlhkosti a koncentrace chloridů a vazebná izoterma chloridů byly použity pro výpočet součinitele vlhkostní vodivosti a součinitele difúze solí pomocí výše uvedené inverzní analýzy. Vstupní data upravená pomocí metody lineární filtrace jsou uvedena na Obr. 11 a 12. Obr. 13 ukazuje závislost vlhkostní vodivosti na obsahu vlhkosti stanovené pomocí výše popsané inverzní analýzy. Z průběhu je patrné, že získané výsledky odpovídají hodnotám pro vysoce kapilárně aktivní materiály, což je v souladu se změřenou celkovou otevřenou pórovitostí materiálu. Obr. 14 prezentuje součinitel difúze chloridů v závislosti na obsahu volných chloridů vypočítaný pomocí funkce κ (w) na obr. 13. Ze získaných výsledků je evidentní, že vypočítaný součinitel difúze chloridů je o několik řádů vyšší než součinitel difúze většiny iontů ve vodě. Tento fakt si lze vysvětlit pouze tak, že koncentrační rozdíl chloridů není jedinou hnací silou difúzního mechanismu v roztoku solí, ale významnou roli zde hrají patrně i další jevy, například elektromigrace iontů či povrchová difúze iontů soli na stěnách porézního prostoru studovaného materiálu.
13
0.7 0.6
0.4
3
3
w [m /m ]
0.5
0.3 0.2 0.1 0 0
0.0005
0.001
0.0015
0.002
0.0025
0.003
0.0035
x/√t [m/s1/2]
Obr. 11 Vlhkostní profil v Boltzmannově formě pro nasákání 1 M roztoku NaCl. 18 16
12 10
3
C [kg/m ] vzorku
14
8 6 4 2 0 0
0.0005
0.001
0.0015
0.002
0.0025
0.003
1/2
x/√t [m/s ]
Obr. 12 Koncentrační profil chloridů v Boltzmannově formě.
14
0.0035
1.0E-05
2
κ [m /s]
1.0E-06
1.0E-07
1.0E-08
1.0E-09 0.10
0.20
0.30
0.40 3
0.50
0.60
0.70
3
w [m /m ]
Obr. 13 Součinitel vlhkostní vodivosti kalcium silikátu v závislosti na obsahu vlhkosti.
2
D [m /s]
1.0E-06
1.0E-07
1.0E-08 21.50
22.50
23.50
24.50
25.50
C [kg/m3(roztoku)]
Obr. 14 Součinitel difúze solí jako funkce koncentrace solí.
15
26.50
6. Závěr Transport vlhkosti a solí porézní strukturou stavebních materiálů představuje dosud nepříliš probádaný problém. Z tohoto důvodu jsou teoretické a experimentální studie zabývající se kombinovaným transportem vody a solí aktuálními tématy při popisu transportních jevů ve stavebních materiálech. Analýza provedená v této práci ukázala, že popis transportu solí v materiálech musí být vždy proveden v kombinaci s transportem vody. Zanedbání vlivu transportu vody a adsorpce iontů solí na stěnách porézního prostoru na transport iontů solí, s čímž se běžně setkáváme v jednoduchých difúzních metodách založených na aplikaci Fickovy difúzní rovnice, může vést k nejednoznačným výsledkům při inverzní analýze koncentračních profilů chloridů, a následně k odchýlení od reality. Použití takto zjednodušených modelů by mělo být omezeno pouze na referenční účely pro základní porovnání transportních vlastností různých typů materiálů a vždy by mělo být bráno v úvahu, že při jejich aplikaci dochází k zanedbání některých velmi důležitých faktorů. Tato práce vznikla na základě finanční podpory grantů GAČR č. 103/04/P085 a 106/04/P255. Seznam použité literatury [1]
Carslaw H. S., Jaeger J. C., 1959, Conduction of Heat in Solids, Oxford – Clarendom Press.
[2]
Grunewald, J., DELPHIN 4.1 - documentation, theoretical fundamentals, TU Dresden, Dresden, 2000.
[3]
Pavlík Z, Jiřičková M, Černý R, Rovnaníková P 2004 Computational Analysis of the Parameters of Coupled Water and Chloride Transport in Cement Mortar on the Basis of Experiemental Water and Chloride Profiles, Cement Combinations for Durable Concrete, Thomas Telford Publications, London, pp. 611-620.
[4]
Matano C. On the relation between the diffusion coefficient and concentration of solid metals. Jap. J. Phys., Vol. 8, 1933, pp. 109-115.
[5]
Bear, J. – Bachmat, Y., 1990, Introduction to Modelling of Transport Phenomena in Porous Media, Vol 4. Kluwer, Dordrecht.
[6]
Tang, L.- Nilsson, L.O., 1993, Chloride binding capacity and binding isotherms of OPC pastes and mortars. Cement and Concrete Research, 23, pp. 247253.
[7]
Jiřičková M., Černý R., 2004, Chloride binding isotherms of different types of building materials. Proceedings of CIB 40 Meeting 2004, Glasgow Caledonian University, Vol. 1, pp. 8.1-8.10.
16