transporte de solutos en el flujo de agua en riego por surcos - Helvia ...
transporte de solutos en el flujo de agua en riego por surcos - Helvia ...
transporte de solutos en el flujo de agua en riego por surcos - Helvia ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
2. Mo<strong>de</strong>lo 1-D <strong>de</strong> <strong>flujo</strong> <strong>de</strong> <strong>agua</strong> y <strong>solutos</strong> <strong>en</strong> un surco <strong>de</strong> <strong>riego</strong><br />
38<br />
∂U<br />
=<br />
∂t<br />
f<br />
don<strong>de</strong> f f agrupa todos los términos <strong>de</strong>l mo<strong>de</strong>lo dinámico sin fricción ni infiltración,<br />
f i conti<strong>en</strong>e los términos <strong>de</strong> infiltración,<br />
y f r los <strong>de</strong> fricción,<br />
f f<br />
f<br />
⎛ 0 ⎛ Q<br />
⎞<br />
⎞<br />
∂ ⎜ 2<br />
= −<br />
⎟<br />
⎜<br />
⎟<br />
∂ ⎜ Q 1<br />
+ ⎟<br />
⎝ gAS ⎠ x gAh<br />
0<br />
⎝ A 2 ⎠<br />
⎛− Pz⎞<br />
= ⎜ ⎟<br />
⎝ 0 ⎠<br />
+<br />
f<br />
i<br />
+<br />
f<br />
r<br />
(2.42)<br />
(2.43)<br />
f i (2.44)<br />
f<br />
r<br />
⎛ 0<br />
= ⎜<br />
⎝−<br />
gAS<br />
f<br />
⎞<br />
⎟<br />
⎠<br />
(2.45)<br />
Una vez que se han clasificado los <strong>flujo</strong>s y términos fu<strong>en</strong>te <strong>de</strong> las ecuaciones <strong>en</strong><br />
estos tres términos, se resolverán <strong>de</strong> forma in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>te y <strong>por</strong> or<strong>de</strong>n, los valores <strong>de</strong> las<br />
variables obt<strong>en</strong>idos con <strong>el</strong> término dinámico se utilizarán para calcular la infiltración, y una<br />
vez que se actualic<strong>en</strong> las variables (t<strong>en</strong>i<strong>en</strong>do <strong>en</strong> cu<strong>en</strong>ta la infiltración) se proce<strong>de</strong>rá al cálculo<br />
<strong>de</strong> la fricción, <strong>en</strong> estos pasos se sigue <strong>el</strong> algoritmo:<br />
con,<br />
m<br />
i<br />
f<br />
i<br />
( ) m<br />
f<br />
∆ U = ∆t<br />
(2.46)<br />
inf<br />
i<br />
f<br />
i<br />
( ) f<br />
f<br />
∆ U = ∆t<br />
inf<br />
(2.47)<br />
r<br />
i<br />
( ) inf<br />
f<br />
r<br />
i<br />
i<br />
∆ U = ∆t<br />
(2.48)<br />
∆U<br />
= ∆U<br />
+ ∆U<br />
+ ∆U<br />
f<br />
i<br />
inf<br />
i<br />
r<br />
i<br />
(2.49)