101
CONGRESO IMTA 2013
Habida cuenta de los resultados de Sauce-
do et al. (2006), sobre el grado de aproxi-
mación de la hipótesis del tiempo de con-
tacto en el riego por melgas, es posible
simplificar más la aproximación emplean-
do la ecuación de Richards en su forma
unidimensional:
C
( )
t
=
z
K
( )
z
1
Ilustración 5. Relación entre la longitud de melga y el gasto de aporte óptimo para el
suelo franco de Montecillo para tres láminas de aplicación: 8, 10 y 12 cm. Ks en cm/s.
que debe resolverse sobre el dominio de
solución definido por una columna de sue-
lo. Como condición inicial para la solución
de la ecuación de Richards unidimensional
se debe especificar la distribución de las
presiones en el espacio:
( )
z
o
=
, en la su-
perficie del suelo se impone una condición
de frontera tipo Dirichlet, con el potencial
de presión igual al tirante de agua calcu-
lado mediante las ecuaciones de Saint-
Venant:
0 t,0 z,h
> = =
, en la frontera in-
ferior se tiene en cuenta una condición de
gradiente unitario:
0 t,P z,1 z/)z (
> = =
,
donde P es la profundidad de la columna
de suelo, que debe ser mayor que la máxi-
ma posición que el frente de humedeci-
miento alcanza durante el tiempo en que
transcurre el riego.
La solución de la ecuación de Richards hace
indispensable representar las propiedades
hidrodinámicas del suelo expresando el
potencial de presión ( ) como una función
del contenido volumétrico de agua ( ) y la
conductividad hidráulica
K
como una fun-
ción de . Como es señalado por Fuentes
