Diferencias Finitas en Transferencia de Calor

ESCUELA POLITÉCNICA NACIONAL FACULTAD DE INGENIERÍA MECÁNICA TRANSFERENCIA DE CALOR I Nombre: Francisco Tipan Zambrano Método de diferencias finitas para sistemas de estado estable

El método de diferencias finitas es una técnica numérica que permite la aplicación de geometrías y condiciones de frontera complejas en sistemas bidimensionales, y es la mejor alternativa para resolver problemas de este tipo. Pasos 1. Red nodal

La solución numérica permite determinar la temperatura en puntos discretos, que deben seleccionarse. Esto se logra dividiendo el medio de interés en pequeños elementos cuyo punto de referencia es su centro y se denomina punto nodal o nodo. En conjunto los elementos se llaman red nodal o malla. Los puntos se designan en un sistema bidimensional como m y n. Los nodos contienen información sobre las propiedades de los elementos y la posición y número de estos se escoge según la conveniencia geométrica y la precisión deseada.

2. Forma de diferencias finitas de la ecuación de calor

Un sistema bidimensional se caracteriza por la siguiente ecuación aproximada : de la segunda derivada 

 

  |  |                

Los gradientes de temperatura son:

|        |        Remplazando en la ecuación inicial:

        En términos de y sería:

       Con una red para la que y sustituyendo las ecuaciones anteriores en la ecuación de distribución de temperaturas se tiene:

       Esta es una forma aproximada de la ecuación de calor para nodos internos equidistantes. 3. Balance de energía

Se aplica la conservación de energía a un volumen de control alrededor de la región nodal. Se supone que el flujo de calor es hacia el nodo. Para sistemas estables se tiene:

 ̇  ̇   El intercambio de energía está influido por la conducción entre m,n y sus cuatro nodos contiguos, así como por la generación, y se obtiene:



∑ ̇   



Donde i se refiere a los nodos vecinos,  es la transferencia por conducción entre nodos, y se supone profundidad unitaria. Para la rapidez de conducción se supone que la transferencia ocurre exclusivamente en bandas en x y en y.

     ( )⁄

Donde es el área, es la aproximación por diferencias finitas al gradiente de temperaturas en la frontera entre nodos. Se puede realizar lo mismo para las demás velocidades de transferencia de calor.

Por conveniencia se resta la temperatura del nodo m,n de la del nodo contiguo, ya que se supone que el flujo de calor en m,n es como la figura 2. Si se tienen en el balance de energía:



 ̇            

Si no hay generación de energía energía.

̇ 

se tiene la ecuación inicial de balance de

La ecuación en diferencias finitas es necesaria para los nodos cuya temperatura es desconocida. Sin embargo no siempre se puede tomar todos los nodos como interiores y usar las ecuaciones dadas. Un ejemplo es una superficie aislada o una superficie que tiene condiciones convectivas, para las que se aplica el balance de energías.

Para ilustrar este método se analiza el nodo mostrado en la figura que intercambia energía por convección con un fluido a  . Las transferencias de calor por conducción son:



       

          En las condiciones en m.n también hay convección por lo que:

   ( ) ( ) Aquí se nota la suposición de que las superficies expuestas en la esquina se encuentran a una temperatura uniforme correspondiente a  , lo que es congruente con que todo nodo está caracterizado por una sola temperatura que es un promedio de temperaturas de la región.



Al final la ecuación queda:

    ( )       Donde la malla es tal que

.