Matematické Fórum

Nevíte-li si rady s jakýmkoliv matematickým problémem, toto místo je pro vás jako dělané.

Nástěnka
22. 8. 2021 (L) Přecházíme zpět na doménu forum.matweb.cz!
04.11.2016 (Jel.) Čtete, prosím, před vložení dotazu, děkuji!
23.10.2013 (Jel.) Zkuste před zadáním dotazu použít některý z online-nástrojů, konzultovat použití můžete v sekci CAS.

Nejste přihlášen(a). Přihlásit

#1 10. 01. 2014 08:47

trakata
Zelenáč
Příspěvky: 6
Reputace:   
 

Výpočet čas.horizontu v jakém se dá očekávát 90%kontam. H2O v nadrzi

Ahoj členové matematického fóra, už rok řeším svůj projekt, kdy profesor mi ho pořád vrací zpět, a už jsem ve fázi zoufalosti, kdy opravdu už nevím jak dál pokračovat. Zadani projektu zni:

V přítoku do vodní nádrže se v poslední době objevil 2,2´PCB o koncentraci poloviční než je jeho rozpustnost (pokud je tato hodnota vyšší než 10-4 mol/l, uvažujte tuto koncentraci). Vypočtěte jak se bude tato látka distribuovat mezi vodu, vzduch a vodní sedimenty. Lze očekávat, že se bude látka významně bioakumulovat? Pokud ano, bude to spíše ve vodních organismech nebo v rostlinách v blízkosti nádrže?
Plocha nádrže 2 ha, průměrná hloubka 8m, průtok 0,5 m3/s, uvažujte různé scénáře distribuce do vzduchu. Kontaminace trvá dlouhodobě.

Rovnovaznou koncentraci vypocitanou mam -  Cout

A potrebuju jeste vypocitat: Výpočet časového horizontu v jakém se dá očekávát  90% kontaminaci vody v nádrží -

zkusila jsem to takto:
dM/dt = Qin * Cin – Qout * Cout - kCoutV
V: dC/dt = Q.Cin-Q.Cout/V-k.Cout
dc/c = Q.Cin-Q.Cout/V –k.d.t
=>
ln C/C0= Q.Cin-Q.Cout/V –k.t
ln0,9= Q.Cin-Q.Cout/V –k.t
t = ln 0,9/Q.Cin-Q.Cout/V –k


A ZPET MI PRISLA ODPOVED OD PANA PROF.
- profes. napoveda:

v časově závislé bilanci schází sedimenty, do kterých ta látka velmi
ochotně přechází. Bilanci si budete muset napsat pomocí fugacity,
protože vyhodnocovaný systém zahrnuje dvě složky (vodu a sedimenty v
nádrži). Na levé straně rovnice je změna množství v těchto složkách
(dM/dt=d(Cw.Vw+Cs.Vs)/dt=d(f.Zw.Vw+f.Zs.Vs)/dt) a napravo přítok -
odtok - reakce, odtok a reakce také zapsány pomocí fugacity. Pak tu
rovnici potřebujete zintegrovat

PROSIM O RADU JAK VYRESIT TUTO BILANCI.

Dekuji Monika

Offline

 

#2 10. 01. 2014 10:59

jelena
Jelena
Místo: Opava
Příspěvky: 30020
Škola: MITHT (abs. 1986)
Pozice: plním požadavky ostatních
Reputace:   100 
 

Re: Výpočet čas.horizontu v jakém se dá očekávát 90%kontam. H2O v nadrzi

Zdravím,

nezakládej, prosím, duplicity viz pravidla. Téma v sekci Chemie jsem smazala, 2. téma v matematice zamknu.

Je to takový specifický dotaz, jehož matematické dopracování asi nebude problémem, ale na sestavení rovnice je třeba si pořádně nakreslit bilanci a promyslet složky - pokusím se podívat přes víkend, ale bez záruky. Snad se také podívá i někdo z kolegů. Děkuji.

Offline

 

#3 10. 01. 2014 12:41

trakata
Zelenáč
Příspěvky: 6
Reputace:   
 

Re: Výpočet čas.horizontu v jakém se dá očekávát 90%kontam. H2O v nadrzi

↑ jelena: Ahoj, děkuji za info, vezmu to na vědomí. Co se týká bilance, budu ráda za každý postřeh.) Děkuji

Offline

 

#4 10. 01. 2014 13:01

kaja.marik
Veterán
Příspěvky: 1915
Reputace:   57 
 

Re: Výpočet čas.horizontu v jakém se dá očekávát 90%kontam. H2O v nadrzi

No ale tam neni zadano, jak se to znecisteni dostava do vzduchu a jak do tech sedimentu. Rychlost je primo umerna koncentraci? Jaka je konstanta umernosti? Nebo jinak?

Offline

 

#5 10. 01. 2014 13:51

trakata
Zelenáč
Příspěvky: 6
Reputace:   
 

Re: Výpočet čas.horizontu v jakém se dá očekávát 90%kontam. H2O v nadrzi

↑ kaja.marik:

Zadáno bylo pouze toto, můžu sem dát co mám zatím vypočítané:

1,85/2= 0,925 mg/l = 0,000925 g/l
n=m/M 0,000925/223,098= 4,14*10-6mol/l =10-4  mol/l
                    4,14*10-3mol/m3       

Výpočet:
Potřebné informace, především chemické a fyzikálně chemické vlastnosti:

chemický název: 1,1'-Biphenyl, 2,2'-dichloro
poločas rozpadu: 9,3 let   
molekulová hmotnost:    M = 223.098 g/mol
rozpustnost ve vodě (25°C):    1,85 mg / l (25 ° C) 1,85→ 0, 00185 g/l , 1,85g/ m3
tenze nasycené páry (25°C):    0,841 Pa
teplota: 25°C→ T 25+273,15 = 298,15 K
log KOW:    log KOW = 4,73 
log KOC 3,65  = (0,74.log Kow + 0,15)

Popis sledované nádrže:               
plocha: 2 ha → 20 000m2
hloubka: 8m
vodní sedimenty dosahují do hloubky:  0,5 m
vzduch nad vodní plochou:  10 m
Přítok: 0,5 m3/s
Odtok: 0,5 m3/s

Výpočet vlastností ŽP:  objemy a hustoty složek ŽP:                     
    vzduch:             VA = 20000 x 10 = 200 000 m3
                               ρaA = 1,2 kg/m3

voda:                VW = 20000 x 8 = 160000 m3
                                ρW = 1 000 kg/m3

sedimenty:        VS = 20000*0,5 = 10000 m3
                         ρS = 2 400 kg/m3

Fugacitní kapacity do vzduchu, vody a sedimentu:
- vzduch:    ZA = 1 / RT = 1 / 8,314 * 298,15 = 4,03*10-4 mol/Pa.m3

- voda:    ZW = Cs / Ps = 8,29*10-3 / 0,841 = 9,85*10-3 mol/Pa.m3
CS = c / M = 1,85 g/m3 / 223,098 g/mol = 8,29*10-3 mol/m3

sedimenty:    ZS = Kp . ρS / H = Kp . ρS / pS / CS
    = (2,23*10-3* 2 400) / (0,841 / 8,29*10-3) = 0,05 mol/Pa.m3
       
Rovnovážná koncentrace ve vodě:
Bilance:     dM/dt = Qin * Cin – Qout * Cout - kCoutV
Ustálený stav:  0 = Q. Cin – Q . Cout -k . Cout .V
Cout = Q . Cin/Q+k.V
Cout = Q . Cin/Q+k.V
Cout = 1,585*108 *4,14*10-3/ 1,585*108 +0,0745*160000   
Cout =4,13*10-3mol/m3                                                                             

Výpočet fugacity:
Cout = 4,13*10-3 mol/m3
Zw = 9,85*10-3 mol/Pa.m3 (vypočteno u fugac. kapacit Za, Zw, Zs)   zdroj 5
f = Cout/Zw = 4,13*10-3/ 9,85*10-3 = 0,41Pa

Výpočet látkového množství C12H8Cl2 v jednotlivých složkách:
vzduch:      na= f.Za.Va =   0,41*4,03*10-4 * 200 000  = 33,046 mol         4%
voda:          nw = f.Zw.Vw = 0,41* 9,85*10-3* 160000 = 646,16 mol         73%
sedimenty:  ns =  f.Zs.Vs = 0,41* 0,05 * 10000  =  205  mol                      23%


A potrebuju jeste vypocitat: Výpočet časového horizontu v jakém se dá očekávát  90% kontaminaci vody v nádrží -
  ----jako napoveda mi prisla od prof. toto:

v časově závislé bilanci schází sedimenty, do kterých ta látka velmi
ochotně přechází. Bilanci si budete muset napsat pomocí fugacity,
protože vyhodnocovaný systém zahrnuje dvě složky (vodu a sedimenty v
nádrži). Na levé straně rovnice je změna množství v těchto složkách
(dM/dt=d(Cw.Vw+Cs.Vs)/dt=d(f.Zw.Vw+f.Zs.Vs)/dt) a napravo přítok -
odtok - reakce, odtok a reakce také zapsány pomocí fugacity. Pak tu
rovnici potřebujete zintegrovat

->coz me bohuzel me nic nerika:/

Offline

 

#6 11. 01. 2014 22:54

jelena
Jelena
Místo: Opava
Příspěvky: 30020
Škola: MITHT (abs. 1986)
Pozice: plním požadavky ostatních
Reputace:   100 
 

Re: Výpočet čas.horizontu v jakém se dá očekávát 90%kontam. H2O v nadrzi

Zdravím v tématu,

↑ trakata: děkuji za upřesnění, něco jsem snad vymyslila, ale něco mi není jasné.

Výpočet časového horizontu v jakém se dá očekávát  90% kontaminaci vody v nádrží

Návrhu pana učitele rozumím tak, že ve sledovaný časový okamžik odebereme vzorek vody a sedimentu (to je vyjádřeno odvození (dM/dt=d(Cw.Vw+Cs.Vs)/dt=d(f.Zw.Vw+f.Zs.Vs)/dt) - závěr tvoří levou část rovnice. Pravou část zřejmě chce popsat jako v materiálu (3. str. 2. vzorec)

To bys asi sestavila. Co nemám úplně jasno:

a) co znamená "90% kontaminaci vody v nádrží". Rozumím tomu správně, že 90 % vody obsahuje 2,2´PCB o stejné koncentraci jako v přítoku (nebo v koncentraci odpovídající rozpustnosti a přebytek jde do sedimentu?) Nádrž by měla být považována za dobře promíchávané prostředí (v modelu) a řeka za píst. Tedy jak interpretovat tuto část zadání?

Využije se tento údaj pro koncentraci odtokové vody? Že koncentrace odtokové je 0,9*koncentrace přítokové?

b) v bilanci napravo je složka "reakce" - uvažuje se tam nějak zadaný poločas rozpadu, nebo co bude v této složce použito?

Přispívá to nějak k řešení Tvého problému? Doufám, že i vážený kolega ↑ kaja.marik: přidá další poznámky.

No ale tam neni zadano, jak se to znecisteni dostava do vzduchu a jak do tech sedimentu.

podle zápisu pana učitele mám dojem, že v bilanci již vzduch neuvažuje. Ale nevím, že ho odmyslím, protože tak bylo v zadání, nebo ho ještě nějak předpokládá přepočíst předem: viz část zadání: "uvažujte různé scénáře distribuce do vzduchu". Vzduch ve výpočtu vyšel o řád méně, než voda. Nejvíce jsou sedimenty.

Děkuji za další komentář.

Offline

 

#7 12. 01. 2014 09:34

trakata
Zelenáč
Příspěvky: 6
Reputace:   
 

Re: Výpočet čas.horizontu v jakém se dá očekávát 90%kontam. H2O v nadrzi

↑ jelena:

rada bych poskytla vice informaci, ale rozumim tomu jako slepy k houslim..:-O .-)

A) Co se tyce bodu a, tak mi psal jeste jeden komentar a to je:
90%kontaminac rozumi:

90% kontaminace neznamená 90% koncentrace na vstupu, ale 90%
ustálené koncentrace
- musíte si prostě napsat časově závislou bilanci pro distribuci té
látky mezi vodu, vzduch a sedimenty. Bilanci pak řešit pro ustálený
stav (alespoň). Pro orientaci je dále užitečné řešit bilance jako
časově závislou a zjistit např. kdy koncentrace dosáhne 90% ustálené
hodnoty - jestli to jsou dny, měsíce, roky..


Využije se tento údaj pro koncentraci odtokové vody? Že koncentrace odtokové je 0,9*koncentrace přítokové?  - netuším:/


JESTE MUZU PRIDAT, CO JSEM PANU PROF. TAKE POSILALA, A TO MI TAKE NEUZNAL (pomoc je od kolegyne, ja jsem to nesestavovala)

dM/dt= Qin*Cin-Qout*Cout-Qair*Cair-Qsed*Csed                                           M=c*V
dM=dc*V
Objem systému zůstává konstantní, tudíž se mění pouze koncentrace.

dC* V/dt= Qin*Cin - Qout*Cout- Qair*Cair - Qsed*Csed

Vše vydělím objemem systému. Nemůžu uvažovat, že vzduch je součást systému, jelikož ho je nekonečně mnoho - počítám pouze objem nádrže a objem sedimentu.

dC*V/dt*V = (Qin*Cin/V - Qout*Cout /V - Qair*Cair /V - Qsed*Csed/ V)       
dC/dt=(Qin*Cin/V - Qout*Cout /V - Qair*Cair /V - Qsed*Csed/ V)       
dC= [(Qin/V)*Cin]-[(Qout/V)*Cout]-[(Qair/V)*Cair]-[(Qsed/V)*Csed]*dt
dC/[(Qin/V)*Cin]-[(Qout/V)*Cout]-[(Qair/V)*Cair]-[(Qsed/V)*Csed]= dt

Provedu substituci:
dC/p =dt

Derivací se zbavím konstant:
dp=d*(Qin/V*Cin)-d*( Qout/V*Cout)-d*( Qair/V*Cair)-d*( Qsed/V*Csed)
        kons.=0                   

Koncentrace ve vzduchu je hodně malá vzhledem k tomu, že se vzduch stále vyměňuje - proto tento součinitel můžu zanedbat.

                    Uvažuji, že se za určitý čas bude koncentrace v sedimentu konstantní -                         sediment se nasytí =0

dp= -d * (Qout/V) *Cout)
Uvažuji, že koncentrace na výstupu je 90% konc. na vstupu.

dp=d Coutx(-Qout/V)= dCoutx-(Qout/V)  : dp=d Coutx -(Qout/V)
dCout= - dp/Qout/V
-dp/(Qout/V)xp = dt
-dp/p = dt * Qout/V

Integruji:
p kon. ∫ dp/p = t t∫ -( Qout/V)dt
p poč.

ln p/pc= -Qout/V*t
Dosadím subst. p

(Q/V):Cin  - (Q/V*Cout) / (Q/V):Cin  = Qout/V *t     ((Q/V):Cin   a (Q/V):Cin  -skrtnuto)
poč.stav

Předpokládám, že Cout je 90% konc. na vstupu
ln Cin - Cout / Cin =( - Qout /V) t
                                           Uvažuji koncentraci 10-4 mol/l = 0,0001mol/l = 0,1 mol/m3
Odlogaritmuji:
Cin-Cout/ Cin = e-(Qout /V ) t
0,1-0,09 /0,1 = e -0,5/160000*t
0,1= e -0,000003125* t
t= 736827, 2298 h
736827, 2298 / 24 = 30701,13458 dní
30701,13458 / 365 = 84, 11269 let.

b) v bilanci napravo je složka "reakce" - uvažuje se tam nějak zadaný poločas rozpadu, nebo co bude v této složce použito? - netuším:/

c)"uvažujte různé scénáře distribuce do vzduchu"   - tzn. napr:
Scénář:
Kolik by přešlo látky do vzduchu, kdyby se ustálila rovnováha s vodou?

Offline

 

#8 12. 01. 2014 14:39

jelena
Jelena
Místo: Opava
Příspěvky: 30020
Škola: MITHT (abs. 1986)
Pozice: plním požadavky ostatních
Reputace:   100 
 

Re: Výpočet čas.horizontu v jakém se dá očekávát 90%kontam. H2O v nadrzi

↑ trakata:

děkuji za další upřesnění. Bohužel, mám dojem, že záležitost není na internetovou diskusi - pokládám stejné otázky, co jste již řešili, Také je vidět, že jste už probrali různé varianty modelu.

Snad se podařilo odkazem upřesnit vzorec bilance s využitím fugacity, ale už více bych do tématu nepřinesla. Lepší bude poprosit učitele o osobní konzultaci, než mailem.

Ať se podaří.

Offline

 

Zápatí

Powered by PunBB
© Copyright 2002–2005 Rickard Andersson