Tecnológico Nacional de México

Instituto Tecnológico de Cerro Azul

Nombre: Luis Enrique Pérez Ruiz

Numero de control: 19500375

Materia: Transferencia de calor.

Ing. Miguel Ángel Salas Bautista.

Carrera: Ing. Electromecánica.

Trabajo: Investigación de la U2

Fecha de entrega: 22/07/21

1
 Índice

2.1 ANÁLISIS POR PARAMETROS DEL TRANSITORIO---------------------------3

2.2 PARED PLANA-----------------------------------------------------------------------------7

2.3 sistemas radiales-------------------------------------------------------------------------10

2.4 aplicación de análisis numérico (diferencias finitas, volumen finito) ---------13

2
 2.1 ANÁLISIS POR PARAMETROS DEL TRANSITORIO.
Conducción transitoria
En la conducción transitoria la temperatura de un mismo punto del cuerpo varía en
el intervalo de tiempo, por lo que el flujo de calor no es constante, es decir la
conducción transitoria ocurre en los procesos de enfriamiento o calentamiento de
superficies sólidas. Ejemplo: calentamiento de un eje en un horno para darle
tratamiento térmico, enfriamiento o calentamiento de paredes de hornos, etc. En el
proceso de transferencia de calor en estado transitorio tiene una gran importancia
la propiedad física de la difusividadtérmica, su sentido físico y la forma de
determinarla fueron expuestas con anterioridad.

ASPECTOS GENERALES CON LA CONDUCCIÓN DE CALOR.

Cuando hablamos de conducción transitoria, estamos hablando de un cuerpo que se enfría o se

calienta en la medida que transcurre el tiempo ( ). Para ilustrar lo que explicaremos
posteriormente imaginemos una pieza de sección transversal plana (una tira de goma de sección
cuadrada) o cilíndrica (eje de un mecanismo cualquiera) que tiene una temperatura inicial (t0) y se
coloca en un horno con temperatura (tf) y que intercambiará calor con la pieza con una intensidad
( ) determinada. En la medida que transcurra el tiempo esta piezas calentará hasta alcanzar una
temperatura (tx) determinada. La velocidad con que ocurra el calentamiento dependerá, entre
otros factores de la difusividadtérmica (a), y tendrá una temperatura en su centro (tx=0) y otra
más alta en su superficie (tx=L). En este caso pudieran surgir algunas interrogantes.

ECUACIÓN DE ESTADO TRANSITORIO.

En el caso de la conducción en estado transitorio la variación de la temperatura con respecto al

tiempo es

• Diferente de cero, por tanto la solución de la ecuación diferencial para este caso es mucho más
engorroso que

• En la conducción estacionaria.

• Las condiciones de unicidad son las siguientes:

• 1. Los parámetros físicos: conductividad ( ), calor específico (c), la densidad ( ), difusividad

térmica (a), etc.

• 2. La forma geométrica y dimensiones del cuerpo: (L1, L2, L3...........LN.) que pueden ser la
longitud, el radio, si

• es una pared plana o cilíndrica.

• 3. La distribución de temperaturas en el estado inicial: =0: t=t0=f(x,y,z).

3
EJEMPLO.
• Ejemplo: Una tira de goma ordinaria dura con espesor 20mm tiene una
temperatura uniforme de 140 C y es
• introducida en agua a 15 C. El coeficiente de traspaso de calor de la tira al agua
es de 65W/m2
• C. Calcule:
• a) Tiempo en que la superficie de la tira alcanza los 25 C.
• b) Temperatura en el centro de la tira transcurrido ese tiempo...
• Qué conozco?
• Que una tira larga de goma con temperatura inicial (t0) de 140 C se introduce en
agua con temperatura
• (tf) de 15 ºC

4
• Intensidad con que el agua fría intercambia calor con la tira caliente (coeficiente
de transferencia de
• Calor) igual a 65W/m2
• C.
• Espesor de la tira ( ) 20 mm
• Análisis teórico del problema.
• El objeto de estudio es un sólido, por tanto el mecanismo predominante de la
transferencia de calo es la
• Conducción, y como éste se está enfriando la conducción es transitoria.
• Para la solución de este problema debo utilizar gráficos donde se relacionan tres
números
• Adimensionales: Fourier, Biot y temperatura adimensional.
SOLUCIÓN.
• Solución del problema
• a) Tiempo en que la superficie de la tira alcanza los 25 C.
• Como se está hablando de una superficie de sección transversal cuadrada y de
valores en su
• Superficie de ella se utiliza el gráfico 1.
• Según el planteamiento del problema debo determinar el tiempo, para ello debo
conocer el
• Valor del Número de Fourier y de él despejar el tiempo.
• Cálculo del número de Biot:

5
6
 2.2 PARED PLANA.
Conducción de calor en régimen transitorio en paredes planas grandes, cilindro largos y esferas
con efecto especiales

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.

Cuando la pared se expone por primera vez al medio circundante que está a T∞< Ti en t = 0, toda
la pared está a la temperatura inicial Ti. Pero la temperatura de la pared en las superficies y cerca
de éstas empieza a caer como resultado de la transferencia de calor de ella hacia el medio
circundante Éste crea un gradiente de temperatura en la pared y se inicia la conducción de calor
desde las partes internas de ella hacia sus superficies exteriores.

Problema de conducción transitoria unidimensional, en forma adimensional Considérese una

pared plana de espesor 2L que, inicialmente, se encuentra a una temperatura uniforme Ti, 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.

7
Ahora se intentará expresar en forma adimensional el problema mediante la
definición de una variable espacial adimensional

Donde:

El número de Fourier indica la medida de calor conducido a través de un cuerpo

en relación con el calor almacenado. Un valor grande de Fourier indica una
propagación más rápida de calor a través del cuerpo.
Ecuación generalizada en sistemas concentrados

8
Sin variables de espacio

Solución exacta del problema de conducción transitoria unidimensional

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.
Las constantes A1 y λ1 son funciones sólo del número Bi, y en la tabla 4-2 (del
libro de transferencia de calor y masa de Cengel) se da una lista de sus valores
con respecto al número Bi, para las tres configuraciones geométricas
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.

9
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 (Cengel)
La cantidad de transferencia de calor Q en un tiempo finito t es menor que este
máximo, y puede expresarse como la suma de los cambios de la energía interna
de toda la configuración geométrica. Un signo negativo para QMáx. Indica que el
cuerpo está rechazando calor

2.3 sistemas radiales.

Los sistemas cilíndricos y esféricos a menudo experimentan gradientes de
temperatura sólo
En la dirección radial, y por consiguiente se tratan como unidireccionales. Además
bajo
Condiciones de estado estacionario, sin generación de calor estos sistemas se
pueden
Analizar usando la expresión de la Ley de Fourier en las coordenadas adecuadas.

Paredes cilíndricas Considere el cilindro hueco de la figura 7, cuyas superficie

externa e interna se exponen a fluidos de diferente temperaturas. Para
condiciones de estado estacionario, sin generación interna de calor, la ley de
Fourier en coordenadas cilíndricas se expresa como

10
Siendo Qr una constante en la dirección radial. Si consideramos también la forma
del área de transferencia para esta geometría, nos queda

Donde A rL r = 2π es el área normal a la dirección de transferencia de calor.

Escribiendo la ecuación anterior en término de integrales con las condiciones de
frontera, Tr T ( 1 1 ) = s, y Tr T ( 2 2 ) = s, , obtenemos:

Si considerando k = constante y resolvemos, nos queda:

También es posible obtener la distribución de temperaturas en la dirección radial

en el cilindro, esto es:

En el caso de la pared cilíndrica, la distribución de temperaturas ya no es lineal,

sino logarítmica.
De este resultado, es evidente que la resistencia térmica para la conducción radial
es de la forma

11
 Paredes esféricas

Considere la esfera hueca de la figura 8. Para el estado estacionario y condiciones

unidimensionales, Qr es constante y la forma apropiada para la ley de Fourier es

Donde 2 A r r = 4π es el área normal a la dirección de transferencia de calor, luego

la ecuación anterior se puede expresar en forma integral como:

Resolviendo bajo la suposición que la conductividad térmica es constante,

obtenemos:

Y la distribución de temperaturas en la esfera vendrá dada por:

12
2.4 aplicación de análisis numérico (diferencias finitas, volumen finito)
Los métodos numéricos son aplicaciones de algoritmos por las cuales es posible
formular y solucionar problemas matemáticos usando operaciones aritméticas
menos complejas. Ellos también se conocen como métodos indirectos. El análisis
numérico idealiza y concibe métodos para «aprobar» de forma eficiente las
soluciones de problemas expresados matemáticamente. El objetivo principal del
análisis numérico es encontrar soluciones
«aproximadas» para problemas complejos.
Un algoritmo es un grupo finito de operaciones organizadas y ordenadas que
permite solucionar un cierto problema. Se trata de una serie de instrucciones o
reglas establecidas que, por medio de una sucesión de etapas, permiten acercar
el resultado real.
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 parala 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
por ejemplo, la geometría del objeto de estudio. Este dominio se considerará
continuo. Este continuo está divide en un número discreto de pequeños cuerpos
con formato específico denominados elementos finitos, e interconectados por
puntos comunes denominados puntos nodales o nudos.

13
El procedimiento de discretización (la división del continuo en partes más chicas),
ecuacionamiento y cálculo es adecuado para la programación y posterior uso en
computadoras, lo que hizo con que sea muy popular y útil para distintas industrias.
Resumiendo, la división de la geometría en elementos finitos permite solucionar un
problema complejo, subdividiéndolo en problemas más simples, lo que posibilita
que la computadora realice con eficiencia estas tareas.
El método de los elementos finitos resuelve distintos problemas, que son
ocasionados y reducidos a sistemas de ecuaciones diferenciales. Por ejemplo,
considerándose estructuras se llega a una ecuación, dicha ecuación del
movimiento, que resume el equilibrio entre esfuerzos internos de una estructura
(fuerza inercial, fuerza de amortiguación y fuerza elástica) y una fuerza externa. Si
el problema es estático o casi estático, se considera apenas el equilibrio entre
esfuerzos internos elásticos y los esfuerzos externos.
Examínese un problema simple estático, para entendimiento, cuyos fundamentos
se aplicarán para solución por el método de los elementos finitos.
En 1660, Robert Hooke observó y describió la llamada ley de la Elasticidad, que
tiene su nombre: la Ley de Hooke. Describe que la variación de la tensión con la
extensión en un resorte es lineal.

Figura 1: Relación (Ley de Hooke) entre esfuerzo externo (F), rigidez (K) y
desplazamiento (u). K.u representa el esfuerzo interno.
En un problema donde la geometría es más compleja que un resorte lineal, se
discretiza la geometría y, a partir de cada componente (elemento) de la estructura
14
discretizada y del conocimiento de las propriedades de materias que la
constituyen, se obtiene la rigidez del elemento. Los elementos se conectan por los
vértices, llamados nudos, formando la estructura continúa discretizada.

Matemáticamente, con la rigidez de cada elemento y conociéndose los

movimientos de cada nudo llamados grados de libertad (GDL) o, en inglés,
Degrees of Freedom (DOF), se forma una matriz de rigidez que representa la
rigidez de la estructura de geometría compleja.
Método de los Volúmenes Finitos (FVM)
El método de volúmenes finitos (Finite Volume Method – FVM en inglés) se
introdujo en la década de 1970 por McDonald, MacCormack y Paullay e,
históricamente, ha sido el método preferido por los científicos e ingenieros que
trabajan con la mecánica de fluidos, aunque él, el FVM, no se limite apenas a la
solución de problemas de mecánica de fluidos.
Comúnmente, el método de los volúmenes finitos se usa para solucionar
problemas de la mecánica de los fluidos, algunos de ellos considerados
complejos, como los que abarcan flujos multifásicos, reactivos o fuertemente
turbulentos. En la práctica, el FVM mostró que es el método más eficaz en el
cálculo y solución de distintos problemas de mecánica de los fluidos.
En el Métodos de los Volúmenes Finitos se usa la idea de observación de Euler,
es decir, material fluye por un volumen de control fijo. A partir de los valores
calculados en los nudos y centros de los volúmenes de control se obtiene una
solución que se transporta al resto del dominio.

Figura 4: Mallas superficiales coloridas por la razón de aspecto de la célula.

Los principios de conservación de la masa, momentum (cantidad de movimiento) y
energía, son la base de la modelación matemática en el Método de los Volúmenes
Finitos para la mecánica del continuo. Por definición, estos principios son
respetados por las ecuaciones montadas a partir de la discretización del continuo
que se realiza en este método. De forma general, el FVM abarca los siguientes
pasos:

15
• Descomponer el dominio en volúmenes de control;
• Formular las ecuaciones integrales de conservación para cada volumen de
control;
• Aproximar numéricamente las integrales;
• Aproximar los valores de las variables en las faces y las derivadas con la
información de las variables nodales;
• Montar y resolver el sistema algebraico obtenido;
El sistema de ecuaciones se resuelve y, como resultado, se obtienen respuestas
como presión, temperatura y velocidad. Obsérvese que es una solución numérica
aproximada. Un aspecto fundamental en los algoritmos de solución es el proceso
de descomposición de matrices donde se han desarrollado diversas técnicas para
aumentar la eficiencia del uso de recursos computacionales y la velocidad de la
obtención de resultados.

16