Ecuación de Richards
Ecuación de Richards
Ecuación de Richards
FACULTAD DE CIENCIAS
INGENIERIA MATEMÁTICA
OPTATIVA 3
1 Introducción
2 Historia:
En la segunda mitad del siglo IXX, el ingeniero francés Henry Darcy, de-
sarrolló el primer estudio sistemático del movimiento del agua a través de un
medio poroso, del cual surgio la ley de Darcy, posteriormente se establecieron
los modelos matemáticos para este fenómeno, para ‡ujo saturado el modelo es
la ecuación de difusión y para el ‡ujo no saturado conocido como ecuación de
Richards (1931):
La ecuación de Richards, cuanti…ca los procesos hidrológicos que describe
de manera precisa la in…ltración del ‡ujo a través de la zona no saturada. La
ecuación de Richards es una combinación de la ecuación de ‡ujo no saturado de
Darcy-Buckinham y la ecuación de conservación de masa, cabe destacar que la
ecuación carece de soluciones analíticas explicitas, por lo que normalmente se
requiere de métodos numéricos para su resolución.
La ecuación de Richards es una ecuación no lineal tridimensional lo más
frecuente es el uso de un esquema numérico que combina el método de elementos
…nitos con el método de Newton
3 Objetivos Generales:
Deducir la ecuación de Richards, estudiar aplicación en los suelos, la existencia
y unicidad de la solución o soluciones de la ecuación de Richards para casos
especí…cos, aplicar el método de elementos …nitos para aproximar la solución y
implementar en el programa computacional FreeFemC++.
4 Objetivos especí…cos:
1. Deducir la ecuación de Richards.
1
2. Discretizar la ecuación de Richards bidimensional con el método de ele-
mentos …nitos para aproximar su solución mediante sistemas de ecuaciones
no lineales.
3. Validar el modelo.
4. Analizar e interpretar los resultados del modelado.
2
Q
2. q = A es decir: caudal que circula por m2 de sección.
3. K conductividad hidráulica.
dh
4. dl gradiente hidráulico expresado en incrementos in…nitesimales.
qz = Kz @H
@z = Kz @h
@z +1
(qx ; qy ; qz ) = !
q =K 5H
3
vacío de R3 limitado por una super…cie cerrada llamada super…cie de control.
Los volúmenes de control pueden ser …jos o móviles. El ‡uido puede entrar o
salir del volumen de control. El volumen de control puede sufrir deformaciones
con el tiempo.
Si F = 5g
R R
gdV = 5g !dS:
R R R
div(F )vdV + 5v F dV = (F !)vdS:
Si F = 5u
R R R
v udV + 5u 5vdV = (5u !)vdS:
f: [T1 ; T2 ] ! R
(x; t) 7 ! f (x; t)
4
donde ! v es la velocidad local de (t) y ! es el vector normal exterior a
(t). Se supone que la integral de línea existe.
5
5.4 Ecuación de Richards
La ley de Darcy-Buckingham ! q = K(h) 5 h en la ecuación de conservación de
masa para un sistema conservativo @% !
@t + div( q ) = 0 sobre T = ]T1 ; T2 [
entonces:
8 @u
< @t r ( (u) 5 u) = 0 sobre T
+ condiciones inciales
:
+condiciones de frontera
Tambien, por la ley de Darcy-Buckingham ! q = K(h) 5 H en la ecuación
de conservación de masa para un sistema conservativo @% + div(!
q ) = 0 sobre
@t
T = ]T1 ; T2 [ entonces:
@u
@t r ( (u) 5 H) = 0
@u @
@t @x K(u) @h
@x
@
@y K(u) @h
@y
@
@z K(u) @h
@z + K = 0
Así, remplazando estos datos se tiene que dado que se trabaja en un suelo
isotrópico
@u @ K(u) @u @ K(u) @u @ K(u) @u @K(u) @u
@t @x C(u) @x @y C(u) @y @z C(u) @z @u @z = 0:
Se ha obtenido
8
< @u@t r (D(u) 5 u) @K(u) @u
@u @z = 0 sobre T
+ condiciones inciales
:
+condiciones de frontera
6
5.6 Condiciones de frontera
7
3. el límite está lejos de la región de interés y se puede suponer que las
presiones y saturaciones de agua y aire son constantes e iguales a las
iniciales.