2.2 Càlcul de la concentració a llacs, estuaris i mars

En el cas d'un llac, s'ha de treballar una mica més per poder obtenir la concentració de contaminant. Considerarem un llac com el de la Fig. 3,

Figura 3: Paràmetres per calcular la concentració d'un contaminant a un llac
\begin{figure}\leavevmode
\begin{center}
\leavevmode\epsfxsize =6cm
\epsfbox{llac.eps}
\end{center}
\protect\end{figure}
amb un abocament d'intensitat $T$, un volum $V$ i un cabal de sortida $R$. Si anomenem $Q$ a la quantitat de contaminant i $C$ a la seva concentració, està clar que la relació entre ambdues magnituds vindrà donada per
\begin{displaymath}
C = \frac{Q}{V}.
\end{displaymath} (4)

Pensem ara per quins motius varia la quantitat de contaminant $Q$. Per una banda tenim una variació deguda a l'abocament

\begin{displaymath}
\frac{dQ_a}{dt}=T,
\end{displaymath} (5)

una altra deguda a la quantitat de contaminant que es perd a través del flux de sortida
\begin{displaymath}
\frac{dQ_f}{dt} = - C R = - \frac{Q}{V} R
\end{displaymath} (6)

i una altra deguda a la descomposició del contaminant
\begin{displaymath}
\frac{dQ_d}{dt} = - Q \lambda.
\end{displaymath} (7)

Sumant aquestes equacions terme a terme obtenim la següent equació diferencial
\begin{displaymath}
\frac{dQ}{dt} = T - Q \left( \lambda + \frac{R}{V} \right),
\end{displaymath} (8)

la solució de la qual és
\begin{displaymath}
Q = \frac{T}{\lambda + R/V}
\end{displaymath} (9)

o aplicant l'Eq. 4, si es prefereix en termes de la concentració,
\begin{displaymath}
C = \frac{T}{V} \frac{1}{\lambda + R/V}.
\end{displaymath} (10)

Per fixar idees, fes les activitat 3, 4 i 5.

Taller de simulaciķ medi ambiental
2009-02-27