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

2. Modelo 1-D de flujo de agua y solutos en un surco de riego
2.3.5. Tratamiento de solutos
El flujo de solutos se representa a partir del modelo 1-D de advección-dispersión
propuesto por García-Navarro et al. (2000). La cantidad de sustancia en el flujo de agua se
especifica mediante su concentración volumétrica, C, y se admite que el transporte de esta
sustancia es mixto, no se desprecian los efectos difusivos en ningún caso.
En el contexto del modelo de agua propuesto anteriormente, C, la concentración
promediada en la sección transversal, es el objetivo principal y su evolución
espaciotemporal en el surco queda descrita por la ecuación (2.30a) o la (2.30b). Para el
modelado numérico se ha adoptado la forma expresada en (2.30a) que se reproduce aquí,
44
( AC)
∂(
QC)
∂
∂t
+
∂x
+ PzC =
∂ ⎛ ∂C
⎞
⎜ AE ⎟
∂x
⎝ ∂x
⎠
(2.69)
La resolución numérica de la ecuación anterior se lleva a cabo en tres fases: primero
se resuelve el término advectivo, en dos pasos sucesivos, y de forma acoplada con el
cálculo del modelo de circulación de agua, resolviéndose conjuntamente y satisfaciendo la
condición de estabilidad (ec. 2.78); a continuación, se resuelve el término dispersivo; por
último se resuelve el término debido a la infiltración.
• Fase 1: Advección
La parte advectiva de la ecuación (2.69) se recoge en la ecuación (2.70), y se acopla
con el modelo de circulación de agua y se resuelve como su término dinámico, es decir, se
le aplica el método de MacCormack.
Paso advectivo-predictor:
( AC)
∂(
QC)
∂
∂t
= −
∂x
(2.70)