Vista/Abrir - RiuNet

Los resultados del modelo numérico, además, nos proporcionaron el valor de las
alturas piezométricas en cada una de las celdas del modelo.
Figura 3.23 Mapa de alturas piezométricas del modelo numérico (en cm)
Estas alturas piezométricas nos permiten el cálculo del campo de velocidades
locales en el medio. Para ello calculamos el valor de los flujos horizontal y vertical en
cada elemento y, puesto que conocemos el área de contacto entre celdas y la porosidad
efectiva de éstas, podemos obtener el valor de la velocidad de flujo.
Para el cálculo de las alturas piezométricas, MODFLOW-88 utiliza la media
geométrica de los valores de conductividad horizontal o vertical de las dos celdas entre
las que se calcula la diferencia de alturas. Es decir, si se calcula el flujo horizontal entre
la celda i,j y la celda i,j+1, MODFLOW lo calculará utilizando:
hi , j
hi , j1
·
Q K K A
x
h h h
i, j i, j i, j 1
hor
(3.4)
Donde K h i,j y K h i,j+1 son la conductividad horizontal de las celdas i,j e i,j+1,
respectivamente. Del mismo modo, para el cálculo del flujo en la dirección horizontal
entre las celdas i,j e i+1,j:
h, h
1,
·
Q K K A
y
v v v i j i j
i, j i, j i 1, j ver
(3.5)
Sin embargo, a partir de la versión MODFLOW-96, se permite al usuario escoger entre
utilizar la media armónica (más realista en el caso de flujo en columna) o la media
geométrica, más adecuada para flujo bidimensional y tridimensional.
97