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

De MateWiki
Saltar a: navegación, buscar
(ANÁLISIS DE LA TEMPERATURA INTERIOR)
(EDIFICIO "CERRADO")
Línea 38: Línea 38:
 
<math> \left \{  
 
<math> \left \{  
 
       \begin{matrix}  
 
       \begin{matrix}  
         T'=)<math>\dfrac{1}{3} </math>(8-T),    t>0 \\
+
         T'=\dfrac{1}{3}(8-T),    t>0 \\
 
         T(to)=14
 
         T(to)=14
 
       \end{matrix}
 
       \end{matrix}

Revisión del 16:22 11 mar 2015

Warning.png 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.
Trabajo realizado por estudiantes
Título Modelo térmico. Grupo 22C
Asignatura Ecuaciones Diferenciales
Curso Curso 2014-15
Autores Antonio Diaz Margarit, Pablo Rodríguez Sandoval, Marta Serrano Grande, David Sánchez Carretero, Ricardo García Dengra
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

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. 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.