transporte de solutos en el flujo de agua en riego por surcos - Helvia ...

Paso advectivo-corrector:
2. Modelo 1-D de flujo de agua y solutos en un surco de riego
( ) ( ) ( ) m
p
m
AC AC − ∆t
QC
i
( AC)
c
i
= (2.71)
i
p
i
1
i+
2
= ( AC)
− ∆t(
QC)
(2.72)
La variación de (AC) por el término advectivo es el promedio del resultado
obtenido en ambos pasos:
• Fase 2: Dispersión
( AC)
mac
i
p
1
i−
2
p
c
( AC)
i + ( AC)
i
= (2.73)
2
Para resolver numéricamente el término de dispersión (ec. 2.74) se aplica un
esquema estándar de diferencias centradas, en el que se utilizarán los valores calculados en
el paso anterior:
∆t
⎡
⎢
∆x
⎢
⎢⎣
⎣
disp
mac
( AC)
− ( AC)
= ⎡(
AE)
i
i
∂
( AC )
∂t
⎤
⎥⎦
=
∂ ⎛ ∂C
⎞
⎜ AE ⎟
∂x
⎝ ∂x
⎠
mac ( C)
− ( C)
∆x
− ⎡
⎢⎣
⎤
⎥⎦
mac ( C)
− ( C)
∆x
⎤
⎥
⎦
(2.74)
m+
1
mac
m+
1
mac
1
i−1
i
i
i+
1
i−
( AE)
1
(2.75)
i+
2
2
La evaluación de las cantidades ⎡(
AE ) 1 ⎤
i−
⎢⎣ 2 ⎥⎦
forma lineal:
y ⎡(
AE) 1 ⎤
i+
⎢⎣ 2 ⎥⎦
se realiza de una
⎡
⎢⎣
( )
m+
1
( AE)
( AE)
m+
1
m+
1
m+
1
1 ⎤ ⎡ i + i−1
⎤
AE i = ⎢
(2.76)
2
− ⎥⎦
2 ⎥
⎣
⎦
45