TECNLOGICO NACIONAL DE MEXICO.

INSTITUTO TECNOLOGICO DE CERRO AZUL, VER.

MATERIA: TRANFERECIA DE CALOR

INGENIERO. MIGUEL ANGEL SALAS BAUTISTA.

UNIDAD 2.

ALUMNA: SELINA MARTINEZ SALAS.

GRUPO :3
NO. CONTROL: 18500361
TEMA 2: Conducción en
estado transitorio.
2.1 Análisis por
parámetros del
transitorio.
2.2 Pared plana.
2.3 Sistemas radiales.
2.4 Aplicación de
análisis numérico
(diferencias finitas,
volumen finito).
TEMA 2. Conducción en estado transitorio.
2.1 Análisis por parámetros del transitorio.
En el análisis de la transferencia de calor, se observa que algunos cuerpos se
comportan como un “bulto” cuya temperatura interior permanece uniforme en
todo momento durante un proceso de transferencia de calor. La temperatura de
esos cuerpos se puede tomar sólo como una función del tiempo, T(t). El análisis
de la transferencia de calor que utiliza esta idealización se conoce como análisis
de sistemas concentrados.
es una cantidad positiva cuya dimensión es (tiempo)–1. El recíproco de b tiene
unidad de tiempo (por lo común s) y se llama constante de tiempo. En la figura
4-3 se tiene la gráfica de la ecuación 4-4 para diferentes valores de b. Se pueden
hacer dos observaciones a partir de esta figura y de la relación antes dada:

1. La ecuación 4-4 permite determinar la temperatura T(t) de un cuerpo en el

instante t= 0, de modo alternativo, el tiempo t requerido para alcanzar el valor
específico T(t).
2. La temperatura de un cuerpo se aproxima a la del medio ambiente, T, en
forma exponencial. Al principio, la temperatura del cuerpo cambia con rapidez,
pero, posteriormente, lo hace más bien con lentitud. Un valor grande de b indica
que el cuerpo tenderá a alcanzar la temperatura del medio ambiente en un
tiempo corto. Entre mayor sea el valor del exponente b, mayor será la velocidad
de decaimiento de la temperatura. Note que b es proporcional al área superficial,
pero inversamente proporcional a la masa y al calor específico del cuerpo. Esto
no es sorprendente, ya que tarda más tiempo en calentarse o enfriarse una masa
grande, en especial cuando tiene un calor específico grande.
Una vez que, con base en la ecuación 4-4, se cuenta con la temperatura T(t) en
el instante t, se puede determinar la razón de la transferencia de calor por
convección entre el cuerpo y su medio ambiente en ese tiempo a partir de la ley
de Newton del enfriamiento como

La cantidad total de transferencia de calor entre el cuerpo y el medio circundante

durante el intervalo desde tiempo de t= 0 hasta t es simplemente el cambio en
el contenido de energía de ese cuerpo:

La cantidad de transferencia de calor llega a su límite superior cuando el cuerpo

alcanza la temperatura T del medio circundante. Por lo tanto, la transferencia de
calor máxima entre el cuerpo y sus alrededores es (figura 4-4)
Criterios para el análisis de sistemas concentrados
El primer paso en el establecimiento de un criterio para la aplicabilidad del
análisis de sistemas concentrados es definir la longitud característica:

y un número de Biot, Bi

La longitud característica Lc que se utiliza en el número de Biot para formas

geométricas simples en las que la transferencia de calor es unidimensional,
como una pared plana larga de 2L de espesor, un cilindro largo de radio ro, y
una esfera de radio ro, se convierte en L (medio espesor), ro/2 y ro/3,
respectivamente. La ecuación 4-9 también se puede expresar como en la figura
4-6.

El número de Biot es la razón de la resistencia interna de un cuerpo a la

conducción de calor con respecto a su resistencia externa a la convección de
calor. Por lo tanto, un número pequeño de Biot representa poca resistencia a la
conducción del calor y, por lo tanto, gradientes pequeños de temperatura dentro
del cuerpo.
En el análisis de sistemas concentrados se supone una distribución uniforme de
temperatura en todo el cuerpo, el cual es el caso sólo cuando la resistencia
térmica de éste a la conducción de calor (la resistencia a la conducción) sea
cero. Por consiguiente, el análisis de sistemas concentrados es exacto cuando
Bi= 0 y aproximado cuando Bi>0.
En general se acepta que el análisis de sistemas concentrados es aplicable
si:

Cuando se satisface este criterio, las temperaturas dentro del cuerpo con
relación a la de los alrededores (es decir, T - T∞) permanecen dentro de un
margen de 5% entre sí, incluso para configuraciones geométricas bien
redondeadas como la de una esfera.
Como consecuencia, cuando Bi<0.1, la variación de la temperatura con la
ubicación dentro del cuerpo es ligera y, de manera razonable, se puede
considerar como si fuera uniforme.
2.2 Pared plana y 2.3. Sistemas radiales.
Los cuerpos relativamente pequeños de materiales intensamente conductores se
aproximan a este comportamiento. Sin embargo, en general, la temperatura
dentro de un cuerpo cambia de punto a punto, así como de tiempo en tiempo.
En esta sección se considera la variación de la temperatura con el tiempo y la
posición en problemas unidimensionales, como los asociados con una pared
plana grande, un cilindro largo y una esfera.
Problema de conducción transitoria unidimensional, en forma
adimensional
Procedimiento de resolución para el caso de una pared plana.
Considérese una pared plana de espesor 2L que, inicialmente, se encuentra a
una temperatura uniforme Ti , como se muestra en la figura 4-12a). En el
instante t= 0, la pared se sumerge en un fluido a la temperatura T∞ y se expone
a transferencia de calor por convección, desde ambos lados, con un coeficiente
de convección de h. La altura y el ancho de la pared son grandes en relación con
su espesor, de donde se puede considerar la conducción de calor en esa pared
como unidimensional. Asimismo, existe simetría térmica respecto al plano
medio que pasa por x= 0 y, como consecuencia, la distribución de temperaturas
debe ser simétrica respecto a ese plano medio. Por lo tanto, el valor de la
temperatura en cualquier punto -x del intervalo -L ≤ x ≤ 0 en el instante t debe
ser igual al valor en x del intervalo 0 ≤ x ≤ L, en el mismo instante. Esto significa
que se puede formular y resolver el problema de conducción de calor en la mitad
positiva del dominio, 0 ≤ x ≤ L, y después aplicar la solución a la otra mitad.
En las siguientes condiciones: propiedades termofísicas constantes, no hay
generación de calor, simetría térmica respecto al plano medio, temperatura
inicial uniforme y coeficiente constante de convección, el problema de
conducción transitoria unidimensional de calor en el semidominio 0 ≤ x ≤ L de
la pared plana se puede expresar como:
La ecuación de conducción de calor en coordenadas cilíndricas o esféricas se
puede expresar en forma adimensional de manera semejante. Obsérvese que la
expresión en forma adimensional reduce el número de variables independientes
y de parámetros, de ocho a tres: .
El resultado final es:
Solución exacta del problema de conducción transitoria unidimensional*

La ecuación diferencial en derivadas parciales en forma adimensional, dada en

las ecuaciones 4-12 junto con sus condiciones de la frontera e inicial, se puede
resolver con la aplicación de varias técnicas analíticas y numéricas, incluidos
los métodos de la transformada de Laplace u otra, el método de separación de
variables, el de diferencias finitas y el de elementos finitos.
Ahora se demostrará el uso del método de separación de variables, mediante su
aplicación al problema de la conducción transitoria unidimensional de calor,
dado en las ecuaciones 4-12. En primer lugar, se expresa la función de la
temperatura adimensional θ(X, τ) como un producto de una función sólo de X
y una función sólo de τ, como:

La parte izquierda de esta ecuación es función sólo de X y la parte derecha sólo

lo es de t. Si se considera que tanto X como τ pueden hacerse variar de manera
independiente, únicamente puede cumplirse la igualdad de la ecuación 4-15,
para cualquier valor de X y de τ, si esta ecuación es igual a una constante.
Además, debe ser una constante negativa, la cual se indicará como -λ² , ya que
una constante positiva hará que la función G(τ) crezca en forma ilimitada con
el tiempo (para hacerse infinita), lo cual carece de significado físico; un valor
de cero para esa constante significa que no hay dependencia respecto al tiempo,
lo cual una vez más no es coherente con el problema físico planteado. Al hacer
la ecuación 4-15 igual a -λ², da
Esta última ecuación se conoce como ecuación característica, y sus raíces
se llaman valores característicos o eigenvalores (o valores propios). En
este caso, la ecuación característica es implícita y, por lo tanto, se necesita
determinar numéricamente los valores característicos. Entonces, se
concluye que se tiene un número infinito de soluciones de la forma
, y la solución de este problema lineal de conducción de
calor es una combinación lineal de ellas
Con esto se completa el análisis para la resolución del problema de
conducción transitoria unidimensional de calor en una pared plana.
La aproximación del cilindro grande permite el supuesto de una
conducción unidimensional en la dirección radial. Es una aproximación
razonable para cilindros que tienen una razón L/r0 ≥ 10 de longitud (l)
con respecto al radio (r0). En la tabla 4-1 se resumen los resultados para
estas tres configuraciones geométricas. La solución para la pared plana
también es aplicable cuando se trata de una pared plana de espesor L cuya
superficie izquierda, en x ≤ 0, esté aislada y la derecha, en x ≤ L, esté
sujeta a convección, ya que éste es precisamente el problema matemático
que se resolvió.
Soluciones aproximadas, analíticas y
gráficas.
La solución analítica para la conducción transitoria unidimensional de
calor en una pared plana comprende series infinitas y ecuaciones
implícitas, las cuales son difíciles de evaluar. Por lo tanto, existen las
soluciones analíticas con el fin de presentar las soluciones en forma
tabular o gráfica, usando relaciones sencillas.
Las cantidades adimensionales definidas en los párrafos anteriores para
una pared plana también se pueden usar para un cilindro o una esfera, al
reemplazar la variable espacial x por r y el semiespesor L por el radio
exterior ro. Nótese que la longitud característica que se encuentra en la
definición del número de Biot se toma como el semiespesor L, para la
pared plana, y el radio ro, para el cilindro y la esfera, en lugar de V/A,
que se usa en el análisis de los sistemas concentrados. Con anterioridad
se mencionó que los términos de las soluciones en serie de la tabla 4-1
convergen con rapidez al aumentar el tiempo; para τ 0.2, si se conserva
el primer término de la serie y se desprecian todos los restantes, se tiene
como resultado un error por debajo de 2%. Suele haber interés en la
solución para tiempos con τ 0.2, por lo que resulta muy conveniente
expresar la solución usando esta aproximación de un término, dada como

Si se comparan los dos conjuntos de ecuaciones anteriores, se observa que en cualquier

parte de una pared plana, un cilindro o una esfera, las temperaturas adimensionales están
relacionadas con la temperatura en el centro por
lo cual muestra que la dependencia de la temperatura adimensional
respecto al tiempo, dentro de una configuración geométrica dada, es la
misma en toda la extensión. Es decir, si la temperatura adimensional en
el centro u0 disminuye 20% durante un tiempo especificado, del mismo
modo disminuye la tempera adimensional θ0 en cualquier otra parte del
medio, durante el mismo tiempo.

Una vez que se conoce el número Bi, se pueden usar estas relaciones para
determinar la temperatura en cualquier parte del medio. La determinación
de las constantes A1 y λ1 suele requerir interpolación. Para quienes
prefieren la lectura de gráficas en lugar de la interpolación, se han trazado
representaciones de estas relaciones y de las soluciones de aproximación
de un término, conocidas como gráficas de temperaturas transitorias.
Nótese que, a veces, las gráficas son difíciles de leer y, por ende, están sujetas
a errores de lectura. Por lo tanto, en su lugar debe preferirse las relaciones antes
dadas. Las gráficas de temperaturas transitorias de las figuras 4-16, 4-17 y 4-
18, para una pared plana grande, un cilindro largo y una esfera, fueron
presentadas por M. P. Heisler, en 1947, y se conocen como gráficas de Heisler.
La temperatura del cuerpo cambia de la temperatura inicial Ti a la de los
alrededores Ti al final del proceso transitorio de conducción de calor. Por lo
tanto, la cantidad máxima de calor que un cuerpo puede ganar (o perder si Ti
T∞) es sencillamente el cambio en el contenido de energía del cuerpo. Es decir,

Si se usan las relaciones apropiadas de temperatura adimensional basadas en la

aproximación de un término para la pared plana, el cilindro y la esfera, y se
realizan las integraciones indicadas, se obtienen las siguientes relaciones para
la fracción de transferencia de calor en esas configuraciones geométricas:
2.4 Aplicación de análisis numérico
(diferencias finitas, volumen finito).
DIFERENCIAS FINITAS
• Los métodos de resolución analítica se basan en la solución de la ecuación
diferencial que rige junto con las condiciones de frontera. Estos métodos
conducen a funciones soluciones para la temperatura en cada punto del medio
• Por otra parte, los métodos numéricos se basan en el reemplazo de la ecuación
diferencial por un conjunto de n ecuaciones algebraicas para las temperaturas
desconocidas en n puntos seleccionados en el medio y la solución simultánea
de estas ecuaciones conduce a valores de la temperatura en esos puntos
discretos.
• Existen varias maneras de obtener la formulación numérica de un problema
de conducción de calor, como los métodos de las diferencias finitas, de
elementos finitos, de elementos frontera y de balance de energía (o de volumen
de control).
Los métodos numéricos para resolver ecuaciones diferenciales se basan en el
reemplazo de las ecuaciones diferenciales por ecuaciones algebraicas. En
diferencias finitas, esto se realiza al reemplazar las derivadas por diferencias.
En los problemas en régimen transitorio, las temperaturas cambian con el
tiempo así como con la posición y, de este modo, la solución en diferencias
finitas de este tipo de problemas requiere la diferenciación en el tiempo y el
espacio.
Esto se realiza al seleccionar un intervalo apropiado de tiempo, t, y resolver para
las temperaturas nodales desconocidas varias veces para cada t hasta que se
obtiene la solución en el instante deseado.
Por ejemplo, considere un objeto metálico caliente que se extrae del horno a una
temperatura inicial de Ti, en el instante t =0, y se deja enfriar en el aire ambiente.
Si se elige un intervalo de tiempo de t 5 min, la determinación de la distribución
de temperatura en la pieza metálica después de 3 h requiere la determinación de
la temperatura 3X60/5 =36 veces, o sea, en 36 intervalos de tiempo. Por lo tanto,
el tiempo requerido de cálculo para este problema será 36 veces el
correspondiente a uno de estado estacionario. La elección de un t más pequeño
aumentará la precisión de la solución, pero también incrementará el tiempo de
cálculo.
La formulación en los problemas de conducción. Este término adicional aparece
como una primera derivada de la temperatura con respecto al tiempo en la
ecuación diferencial, y como un cambio en el contenido de energía interna
durante t en la formulación del balance de energía.
Los nodos y los elementos de volumen en los problemas en régimen transitorio
se seleccionan igual que en los de estado estacionario y, una vez más, si se
supone por conveniencia que toda la transferencia de calor es hacia el elemento,
el balance de energía sobre un elemento de volumen durante un intervalo de
tiempo t se puede expresar como
Las temperaturas nodales en los problemas en régimen transitorio por lo común
cambian durante cada intervalo de tiempo y el lector puede preguntarse si debe
usar las temperaturas en el previo intervalo de tiempo i o en el nuevo intervalo
de tiempo i +1 para los términos del primer miembro de la ecuación. Bien, los
dos procedimientos son razonables y ambos se aplican en la práctica. En el
primer caso, se dice que el procedimiento en diferencias finitas es el método
explícito y el segundo el método implícito, y se expresan en la forma general
como
Parece que la derivada con respecto al tiempo se expresa en la forma de
diferencia hacia adelante en el caso explícito y en la de diferencia hacia atrás en
el implícito.
El método explícito es fácil de poner en práctica, pero impone un límite sobre
el intervalo de tiempo admisible para evitar inestabilidades en la solución, y el
método implícito requiere que las temperaturas nodales se resuelvan en forma
simultánea pero no impone límite sobre la magnitud del intervalo de tiempo.
VOLUMENES FINITOS.
El método de los volúmenes de control finitos permite discretizar y resolver
numéricamente ecuaciones diferenciales. Es un método alternativo a los de
diferencias finitas y elementos finitos
• Consideremos una malla de discretización del espacio fluido. En torno a cada
punto de esta malla se construye un volumen de control que no se traslapa con
los de los puntos vecinos. De esta forma el volumen total de fluido resulta ser
igual a la suma de los volúmenes de control considerados.
•La ecuación diferencial a resolver se integra sobre cada volumen de control, lo
cual entrega como resultado una versión discretizada de dicha ecuación.
• Para realizar la integración se requiere especificar perfiles de variación de la
variable dependiente entre los puntos de la malla, de modo de poder evaluar las
integrales resultantes.
• La principal propiedad del sistema de ecuaciones discretizadas resultante, es
que la solución obtenida satisface en forma exacta las ecuaciones de
conservación consideradas, independientemente del tamaño de la malla.
De forma general envuelve los siguientes pasos:
• Descomponer el dominio en volúmenes de control.
• Formular las ecuaciones integrales de conservación para cada volumen de
control.
• Formular las ecuaciones integrales de conservación para cada volumen de
control.
• Aproximar los valores de las variables en las caras y las derivadas con la
información de las variables nodales.
• Ensamblar y resolver el sistema algebraico obtenido.
• El punto de partida del método de volúmenes finitos es la descomposición del
dominio en pequeños volúmenes de control (VCs) donde las variables son
almacenadas en los nodos.
• Los nodos, donde se guardan las variables a resolver, se ubican en los vértices
de esta malla o en los centros de los volúmenes (Fig. 2), y es la referencia con
la que se cuenta para aproximar la solución en el resto del dominio.

• Después de haber definido los VCs, las ecuaciones de conservación son

descritas en su forma integral para cada volumen.
• En este proceso el teorema de la divergencia de Gauss es usado para convertir
integrales de volumen sobre operadores de divergencia y gradiente a integrales
de superficie.
Por ejemplo la ecuación de transporte general:

QUEDA:

Donde, ϕ es la variable transportada a través de un medio de densidad ρ y

constante de difusión Γ que se mueve a una velocidad V con un término fuente
Sϕ.
• Para la obtención final del sistema algebraico de ecuaciones las integrales son
aproximadas, las volumétricas haciendo un cálculo simple basado en las
variables nodales, las superficiales aproximándolas por los balances a través de
las caras de los VCs.
• En dicho proceso aproximaciones de las variables en las caras y de los
gradientes son hechas en función de las variables nodales.

• Finalmente, el sistema de ecuaciones es resuelto a través de métodos

segregados o acoplados.