El sistema de ecuaciones de Saint-Venant y Richards del riego por gravedad: 2. Acoplamiento numérico para la fase de avance en el riego por melgas

Se presenta una solución numérica monótona para el acoplamiento interno de las ecuaciones de Saint-Venant y Richards en el riego por melgas. La solución de las ecuaciones de Saint-Venant se aproxima utilizando un esquema lagrangiano en diferencias finitas; mientras que para la ecuación de Richards se utilizan elementos finitos para la integración en el espacio y diferencias finitas implícitas para la integración en el tiempo. Se muestra que una discretización de la ecuación de momentum, basada en aproximaciones de las derivadas y los coeficientes mediante ponderaciones en tiempo y espacio, puede generar cambios de monotonía en la vecindad de la cabecera de la melga; asimismo, se muestra que dicha situación no se presenta si se utiliza una discretización basada en: a) las derivadas en el espacio y la pendiente de fricción en una celda de cálculo se aproximan adelante en el tiempo; b) las derivadas en el tiempo se aproximan utilizando una forma ponderada; y c) los coeficientes son calculados en el tiempo anterior.