Ecuación de Richards

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 8

UNIVERSIDAD CENTRAL DEL ECUADOR

FACULTAD DE CIENCIAS
INGENIERIA MATEMÁTICA
OPTATIVA 3

1 Introducción

La ecuación de Richards es usada para predecir el movimiento de agua en los


suelos que son variablemente saturados. Esta ecuación tiene aplicaciones im-
portantes en hidrología, agronomía, protección ambiental, y otros campos rela-
cionados con el suelo. En el campo de la protección ambiental, la predicción
del movimiento del agua permitirá estimar el modelo de transporte de contam-
inantes.

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.

5 Deducción de la ecución de Richards


5.1 Ley de Darcy
El experimento de Darcy:
Los aparatos de que usó Darcy se llaman permeámetros de carga constante
(…g. 3.1), el que consiste en un recipiente vertical se sección constantes A, el
cual se rellena de un medio poroso por ejemplo arena, se hace circular agua
conectando uno de sus extremos a un depósito elevado cuyo nivel se mantiene
constante de manera que la presión del ‡uido a la entrada es constante. Se mide
el volumen de agua que atraviesa el recipiente por undad de tiempo es decir el
caudal Q y las alturas piezométricas h1 y h2 medidas en dos puntos entre los
cuales existe una longitud 4l.

Darcy encontró que


4h
Q = KA
4l
donde K es la constante de proporcionalidad. La ley de Darcy se expresa de
esta forma:
dh
q= K
dl
donde:

1. h es el cambio de altura respecto a la horizontal.

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.

el signo menos se debe a que el caudal es una magnitud vectorial, cuya


dirección es hacia los 4h decrecientes; es decir, que 4h o dh es negativo y, por
tanto, el caudal será positivo.
Ley de Darcy en tres dimenciones
Considere tres direcciones x; y; z considerando suelos anisótropos (distintas
características en cada una de sus direcciones)
@h
qx = Kx
@x
@h
qy = Ky
@y
@H
qz = Kz
@z
Observese que para el eje Z se tiene la acción de la gravedad por lo cual se
tendria que la presión en condiciones normales de 1atm, despues la forma de
estado transitorio de la ecuación de Richards solo en la dirección vertical es

qz = Kz @H
@z = Kz @h
@z +1

Si consideramos un suelo isotropo (Término aplicado a los medios o materi-


ales cuyas propiedades en un punto dado son independientes de la dirección en
la que sean medidas) se tiene que Kx = Ky = Kz = K, se tiene que

(qx ; qy ; qz ) = !
q =K 5H

5.2 Ley de Darcy-Buckingham


En los suelos parcialmente saturados existen dos ‡uidos en los poros: agua y aire.
El coe…ciente de permeabilidad de suelos parcialmente saturados aumenta al
aumentar la presión del líquido, pues esto provoca un incremento en la cantidad
de gas disuelta y, por tanto, una disminución en el espacio ocupado por burbujas
gaseosas luego en este tipo de suelo K depende de h y la ley de darcy toma la
forma:
!
q = K(h) 5 H

5.3 Resultados Fundamentales del Análisis Vectorial


5.3.1 Volumen de control
Un volumen de control es una región …ja del plano o del espacio a través de la
cual circula un ‡uido. Se considera como volumen de control a un conjunto no

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.

5.3.2 Teorema de la divergencia de Gauss


Sea un dominio acotado de R3 con su frontera, = [ y F : ! R3
@f @g
con F = (f; g; h) un campo vectorial de…nido en . Supongamos que @x ; @y
y
@h
@z son integrables en , entonces se tiene:
R R
div F dV = F !dS:

Si F = 5g
R R
gdV = 5g !dS:

En términos de la derivada normal queda de la forma


@g
@!
= 5g !:

5.3.3 Teorema de Green - Gauss


Sea un dominio acotado de R3 con su frontera, = [ y F : ! R3
@v @v
con F = (f; g; h) un campo vectorial de…nido en . Supongamos que @x f; @y g
@v
y @z h son integrables en , entonces se tiene:

R R R
div(F )vdV + 5v F dV = (F !)vdS:

Si F = 5u
R R R
v udV + 5u 5vdV = (5u !)vdS:

5.3.4 Teorema de transporte de Reynolds

Sea R3 un dominio acotado con su frontera tipo C 1 a trozos, sean


T1 ; T2 2 R tales que 0 T1 < T2 y

f: [T1 ; T2 ] ! R
(x; t) 7 ! f (x; t)

Una función que para todo t 2 [T1 ; T2 ]; f ( ; t); @f


@t ( ; t) son integrables en
. Para t 2 [T1 ; T2 ], sea (t) un volumen de control y (t) su frontera.
Entonces
R R R
d
dt (t)
f (x; t)dV = (t) @f
@t (x; t)dV + (t)
f (x; t)(!v !)dS:

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.

La interpretación física del teorema del transporte de Reynolds es : la rapidez


de cambio de ‡ujo en el volumen de control es igual a la rapidez de ‡ujo que
entra en el volumen de control más el ‡ujo en la frontera.
Por el teorema de la divergencia de Gauss, tenemos
R R
(t)
f (x; t)(!
v !)dS = (t) div(f ! v )dV:

Así,de la ecuación del teorema del transporte de Reynolds se tiene


R R !
d
dt (t)
f (x; t)dV = (t) @f
@t (x; t) + div(f v )dV:

5.3.5 Ecuación de continuidad


Sea (t) un volumen de control que se mueve con una velocidad ! v , t 2 [T1 ; T2 ],
por lo que la densidad del ‡uido contenido en (t) es %(x; t) = dm
dx (t) , de donde
dm(t) = %(x; t)dV e integrando sobre (t) tenemos
R R
m(t) = (t) dm(t) = (t) %(x; t)dV:

Derivando la masa respecto al tiempo se tiene que


dm d
R
dx (t) = dt (t)
%(x; t)dV;

y por el teorema de transporte de Reynolds


R
dm
dx (t) =
@%
(t) @t
(x; t) + div(%!
v )dV:

Ahora, por la ley de la conservación de la masa dm dx (t) = 0, t 0. Esto es


valido para sistemas conservativos. Si exitiera algún factor externo que variara
la masa del ‡uido (perdida de material por reaciones químicas, evaporación,
entre otros) , se tendria que dmdx (t) 6= 0:
R
0 = dm
dx (t) =
@%
(t) @t
(x; t) + div(%!
v )dV:

Como el volumen de control (t), es …jo pero arbitrario para un t …jo


@%
@t (x; t) + div(%!
v ) = 0; casi todo punto de ]T1 ; T2 [:

la cual se la denomina ecuación de continuidad para sistemas conservativos.

Ahora, por de…nición de caudal


! dm d!
q = dm
dt = d!
x
x dt
= %!
v

Remplazando en la ecuación de continuidad para sitemas conservativos de


tiene
@% !
@t (x; t) + div q = 0; sobre ]T1 ; T2 [:

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

Introduciendo la capacidad especi…ca de un ‡uido C(u); dada por


@u
C(u) = @h ;

y la difusividad de un ‡uido en un suelo D(u); dada por


@h K(u)
D(u) = K(u) @u = C(u) :

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

5.5 Condiciones inciales


El estado inicial del sistema se puede de…nir en términos de las variables pri-
marias o en términos de otras variables, a partir de las cuales se pueden calcular
las primarias.
El contenido de agua o humedad es la cantidad de ‡uido (agua) contenida en
un material con la velocidad de ‡ujo, tal como el suelo (la humedad del suelo),
las rocas, la cerámica o la madera medida sobre la base de análisis volumétricos o
gravimétricos. Esta propiedad se utiliza en una amplia gama de áreas cientí…cas
y técnicas y se expresa como una proporción que puede ir de 0 (completamente
seca) hasta el valor de la porosidad de los materiales en el punto de saturación .
Así, para las condiciones iniciales t = T1 se tiene que ve cual es la humedad
presente en cada punto del volumen de control T jt=0 : Donde u(! x ; 0) = g(x; y; z).

6
5.6 Condiciones de frontera

1. Las condiciones de frontera de Dirichlet especi…can los valores de la solu-


ción en una parte determinada de la frontera;
2. Las condiciones de frontera de Neumann especi…can el valor de la derivada
espacial de la solución. En el caso de ‡ujo en medios porosos, este tipo
de condición generalmente se escribe en términos del ‡ujo de ‡uido en la
dirección normal al límite;
3. Condiciones de frontera de Robin, que especi…can una relación entre el
valor de la solución y su derivada, y pueden verse como una generalización
de las condiciones de frontera de Dirichlet y Neumann.
Se pueden especi…car varias condiciones de frontera en diferentes partes de
la frontera del dominio. Para un problema de ‡ujo constante, descrito por una
ecuación elíptica, es necesario proporcionar condiciones de frontera de Dirichlet
en al menos un punto de la frontera, de lo contrario el problema se vuelve mal
planteado. Para problemas inestables, descritos por ecuaciones parabólicas, no
existe tal restricción. En el caso de un ‡ujo de dos fases, es necesario especi-
…car las condiciones de contorno para cada fase, y se pueden utilizar diferentes
tipos de condiciones para cada ‡uido en la misma parte del límite. A contin-
uación se enumeran algunos ejemplos de condiciones de contorno que se utilizan
normalmente en el modelado de ‡ujo insaturado.
Las condiciones de contorno de Dirichlet proporcionan valores de las pre-
siones y / o saturaciones de las dos fases del ‡uido. Por lo general, se aplican
cuando:
1. se conoce la presión del agua, por ejemplo
1. (a) el límite corresponde al nivel freático,
(b) se conoce la distancia entre el límite y el nivel freático, y se puede
suponer una distribución hidrostática de la presión del agua entre
ellos,
(c) el límite está en contacto con una masa de agua libre que se
supone estática, por ejemplo en la pendiente hacia el agua de
un dique o en la super…cie del suelo estancada por una capa de
agua,
(d) se aplica un valor negativo especí…co de la presión del agua, por
ejemplo durante experimentos de in…ltración o drenaje en condi-
ciones controladas,
1. la presión del aire es conocida, por ejemplo. el límite está expuesto a la
atmósfera y no está completamente saturado de agua,
2. la saturación de agua es conocida, por ejemplo el límite está completa-
mente saturado y la presión del aire no se puede especi…car de manera
signi…cativa,

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.

Las condiciones de frontera de Neumann se usan típicamente cuando la fron-


tera es impermeable para una o ambas fases, %! v!n = 0 , donde ! n es un vector
unitario normal al límite del dominio . Esta condición también se usa para
representar los ‡ujos de in…ltración y evaporación en la super…cie del suelo.

También podría gustarte