Page 20 - RELEVAMIENTO TOPOGRÁFICO VOLUMEN I
P. 20
Reescribiendo la ecuaci6n (11) de la siguiente forma:
F T?.t' + G T}+' + H Tjp++i' = I (35)
En cada intervalo de tiempo se resuelve un sistema de ecuaciones lineales, cuya matriz de
coeficientes es tridiagonal. El sistema de ecuaciones se soluciona mediante la aplicaci6n de un
M6todo de Doble Barrido aioache, 1972).
Condiciones de contorno
En superficie
KH£)superf=±(Qn+HI+H) (36)
=)`-o (37)
Condiciolles Iniciales
Debido a la derivada temporal de primer orden, se hace necesario contar con un perfil de
temperatura en la vertical, como condici6n ihicial:
T¢,0)=To V j=1,J (38)
Se adopt6 un perfil uniforme de temperaturas, compatible con una condici6n isotermal en el
embalse.
En la Figura 2.2 se muestra el esquema de discretizaci6n espacial. Las variables calculadas,
temperaturas y coeficientes de difusi6n, son asignadas al borde superior de cada capa.
A(zl T(I,I)
I I j-1 y
tLriIi.IEI I
1` Az J I _ I I
I,+ t i+1 al I
I
I I
Figura 2.2 ¢Ienderson-Sellers, 1984)