Vista/Abrir - RiuNet
Vista/Abrir - RiuNet
Vista/Abrir - RiuNet
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
C<br />
t<br />
C<br />
d<br />
C<br />
v C D C Q f x<br />
C C C<br />
(2.71)<br />
Donde Q ( x, t) v( x, t) C( x, t)<br />
C<br />
es el flujo dispersivo condicional. De forma<br />
similar se obtienen las expresiones para las condiciones de contorno:<br />
C( x,0) C ( x )<br />
x <br />
C<br />
0<br />
C( x,0) C ( x , t)<br />
x <br />
C<br />
D<br />
D C( x, t) n( x) W ( x , t)<br />
x<br />
d<br />
C<br />
<br />
v( x, t) C( x, t) Dd C( x, t) QC( x, t) n(x) P( x , t)<br />
x<br />
C C C <br />
El flujo dispersivo Q C viene dado por la expresión implícita:<br />
2<br />
1<br />
C<br />
3<br />
(2.72)<br />
Q ( x, t) ( x, t; s, ) Q ( x, t; s, ) dsd<br />
C<br />
0 <br />
C y C<br />
t<br />
0 <br />
( x, t; s, ) C( s, ) dsd<br />
- ( x, t; s, ) C( s, ) dsd<br />
0<br />
t<br />
0 3<br />
T<br />
- ( x, t; s, ) Q ( x, t; s, ) n( s)<br />
ds<br />
d<br />
+<br />
t<br />
t<br />
<br />
<br />
t<br />
0<br />
<br />
<br />
3<br />
C<br />
C y C<br />
C<br />
C<br />
( x, t; s, ) C( s, ) n( s)<br />
dsd<br />
C<br />
C<br />
C<br />
(2.73)<br />
Donde:<br />
( x, t; s, ) G( x, t;<br />
s, ) v( x, t)<br />
C<br />
T<br />
( x, t; s, ) G( x, t; s, ) v( x, t)<br />
v<br />
( s, )<br />
C<br />
( x, t; s, ) G( x, t; s, ) v( x, t)<br />
f( s, )<br />
C<br />
C<br />
C<br />
C<br />
Esta función aleatoria de Green, G(x,t;s,τ), satisface la ecuación estocástica de<br />
convección-dispersión sujeta un cero inicial y unas condiciones de contorno dadas<br />
(Morales-Casique y Neuman, 2006).<br />
Para poder reducir este sistema acoplado a una única ecuación de transporte, hay<br />
que imponer una serie de restricciones físicas. Si el campo de velocidades v es<br />
estacionario, los momentos de la función de Green dependen sólo de los incrementos en<br />
el espacio y el tiempo (Dentz y Tartakovsky, 2008). Es decir, de G( xs,<br />
t<br />
.<br />
Igualmente, las integrales que contienen ( x, t; s, )<br />
son de un orden inferior y pueden<br />
despreciarse en (2.73). Del mismo modo, para que el campo de velocidades sea<br />
estacionario, el dominio debe ser infinito, por lo que Ω = Ω ∞ y así las integrales<br />
C<br />
43