Diferencia entre revisiones de «Modelo Térmico. (grupo 6C)»

Revisión del 16:46 11 mar 2015

Este artículo está en versión beta. El autor de este artículo no lo ha terminado todavía, por favor no lo edites hasta que elimine este mensaje.

1 INTRODUCCIÓN

El trabajo que a continuación desarrollaremos, se basa en la formulación de un modelo matemático que describa el comportamiento de la temperatura en el interior de un edificio a lo largo de 24 horas en función de:

• Temperatura exterior. (M(t))
• Calor que se genera dentro del edificio. (H(t))
• Enfriamiento y calentamiento producido por el sistema de aire acondicionado y calefacción respectivamente. (U(t))

Para formular el modelo hemos utilizado un análisis compartimental ya que consideramos el edificio como una sola habitación en el que la temperatura está definida como T(t)

2 PROBLEMA DE CAUCHY

La ley de enfriamiento de Newton establece que:

"Cuando la diferencia entre un cuerpo y su medio externo no es demasiado grande, el calor transferido en la unidad de tiempo hacia el cuerpo o desde el cuerpo por conducción, convección y radiación es aproximadamente proporcional a la diferencia de temperatura entre el cuerpo y el medio externo"

Teniendo esto en cuenta, esta ley establece que hay una razón de cambio de temperatura que es proporcional a la temperatura interior y exterior del edificio designadas en nuestro caso como T(t) y M(t) respectivamente. Dicha razón queda determinada según el P.V.I siguiente:

[math] \left \{ \begin{matrix} T'=k(M(t)-T)+H(t)+U(t), t\gt0 \\ T(to)=T_0 \end{matrix} \right . [/math]

k: es una constante real que depende del número de ventanas, aislamiento del edificio..

3 ANÁLISIS DE LA TEMPERATURA EN EL INTERIOR DEL EDIFICIO

3.1 EDIFICIO "CERRADO"

A la hora de analizar la temperatura del edificio vamos a ir suponiendo distintas situaciones,en está primera situación vamos a tratar de saber cuanto tiempo tarda en cambiar considerablemente la temperatura del edificio cuando esta cerrado, esto quiere decir que en el interior del edificio no hay actividad que afecte al calentamiento o enfriamiento del mismo ni tampoco hay calefacción ni aire acondicionado, por tanto las variables U(t) y H(t) son iguales a cero; por sólo vamos a tener en cuenta la temperatura exterior M(t).

Suponemos que la temperatura en el interior del edificio y a media noche es 14ºC, y la temperatura exterior es constante e igual 8ºC, (M(t)=8ºC), por otra parte, la constante de tiempo del edificio está definida como t'= [math]\dfrac{1}{k} [/math], nosotros suponemos que t'=3, por lo tanto la constante k=[math]\dfrac{1}{3} [/math].

Con las hipótesis expuestas en el párrafo anterior, obtenemos el siguiente P.V.I:

[math] \left \{ \begin{matrix} T'=\dfrac{1}{3}(8-T), t\gt0 \\ T(to)=14 \end{matrix} \right . [/math]

Utilizamos el método de Euler implícito para el PVI anterior, con una longitud de paso h=0,001 y compararemos este resultado con el resultado de la ecuación exacta obtenida analíticamente. Implementado el código que a continuación se muestra en Matlab, se ha llegado al gráfico que responde a nuestra cuestión inicial.

Ante todo cabe destacar que la solución exacta y la obtenida mediante el método numérico son prácticamente idénticas. Además en el gráfico se observa un brusco descenso de la temperatura durante las primeras horas del día, para luego estabilizarse cuando ya han pasado aproximadamente 10 horas, es decir el edificio trata de alcanzar el equilibrio térmico entre la temperatura exterior e interior perdiendo calor durante la noche.

3.2 EDIFICIO EN UN DÍA HÁBIL DURANTE LA PRIMAVERA U OTOÑO

En esta segunda situación vamos a analizar la variación de la temperatura en un día hábil, esto quiere decir que a diferencia del aparado anterior en el edificio se esta desarrollando una actividad que genera un calentamiento como por ejemplo las luces, máquinas... Dicho calentamiento está definido mediante la variable H(t). Al suponer estar en una estación del año intermedia no se precisa calefacción ni aire acondicionado por lo tanto, U(t)=0.