CARATULA

INFORME N° 1
MATERIA:

ESTUDIANTES:
- HERBAS CAHUAYA JOEL GUSTAVO
- MAMANI BLANCO JOSE LUIS
- MAMANI QUISPE SOLEDAD BETTY
- QUISPE COLQUE NELCY XIMENA
DOCENTE:
ING. JORGE VASQUEZ PEÑARANDA
 INDICE
OBJETIVOS: .......................................................................................................................................... 1
Objetivo General: ............................................................................................................................ 1
Objetivo Específico: ......................................................................................................................... 1
JUSTIFICACION: ................................................................................................................................... 1
FUNDAMENTO TEORICO ..................................................................................................................... 2
Análisis para un flujo Bidimensional ............................................................................................... 4
Análisis para un flujo Tridimensional .............................................................................................. 4
La Ley de Fourier ............................................................................................................................. 4
Flujo Multidimensional: .................................................................................................................. 5
Ley de Fourier para Flujo Multidimensional ............................................................................... 5
Ecuación de energía en conducción multidimensional (flujo estacionario) ............................... 6
PROCEDIMIENTO EXPERIMENTAL ....................................................................................................... 7
CALCULOS, RESULTADOS Y GRAFICOS ............................................................................................... 9
Simulación: Transferencia de calor para un tiempo igual a cero .................................................... 9
Simulación: Transferencia de calor para un tiempo igual a un minuto ....................................... 10
Gráfica: Tomando como puntos de referencia los extremos de nuestra figura para un tiempo
igual a 1 minuto............................................................................................................................. 11
Gráfica: Tomando como puntos de referencia uno de los extremos y el punto medio de nuestra
figura para un tiempo igual a 1 minuto........................................................................................ 12
Gráfica: Tomando como puntos de referencia uno de los extremos y un punto en las aletas de
nuestra figura para un tiempo igual a 1 minuto .......................................................................... 14
CONCLUSIONES, RECOMENSACIONES Y SUGERENCIAS .................................................................... 15
CUESTIONARIO .................................................................................................................................. 16
BIBLIOGRAFIA .................................................................................................................................... 40
 INDICE DE FIGURAS
Figura 1 ________________________________________________________________________________ 2
Figura 2 ________________________________________________________________________________ 3
Figura 3 ________________________________________________________________________________ 4
Figura 4 ________________________________________________________________________________ 6
Figura 5 ________________________________________________________________________________ 9
Figura 6 _______________________________________________________________________________ 10

INDICE DE GRAFICOS

Grafica 1 ______________________________________________________________________ 11
Grafica 2 ______________________________________________________________________ 11
Grafica 3 ______________________________________________________________________ 12
Grafica 4 ______________________________________________________________________ 13
Grafica 5 ______________________________________________________________________ 14
Grafica 6 ______________________________________________________________________ 14
 FLUJO MULTIDIMENSIONAL EN ENGRANAJES

ALEACION HIERRO- COBRE

PIÑON CILINDRO DE DIENTES HELICOIDALES 16 DIENTES Y 10 cm DE

DIAMETRO

OBJETIVOS:

Objetivo General:
Realizar el estudio del flujo multidimensional, mediante la simulación de
transferencia de calor en un engranaje con ayuda del programa Comsol.

Objetivo Específico:
Observar, analizar e interpretar los resultados obtenidos en la simulación.
Representar con la ayuda del programa comsol, graficas del flujo multidimensional
en engranajes.
Simular y analizar el cambio de temperatura que puede existir en un solido.

JUSTIFICACION:

Es de suma importancia el estudio del flujo multidimensional, este nos permite analizar e
interpretar las diferentes direcciones que puede tomar el flujo de calor en un cuerpo
estudiado, nos permite analizar de una manera más real el flujo de calor, que es la energía
que en forma de calor atraviesa una superficie, tomando más de una dirección.

La transferencia de calor es de suma importancia en la industria. El flujo de calor, que puede

tomar más de una dirección (flujo multidimensional) permite el cambio físico y químico de
muchos materiales, que luego serán usados para generar diferentes productos en ámbito de
la industria; Es por ello la importancia de su estudio.

pág. 1
También podemos hablar de la importancia de la transferencia de calor, flujo
multidimensional en la naturaleza, este permite el flujo de calor y proporciona diferentes
cambios físicos los cuales mantienen una estabilidad ecológica.

A lo largo de los años, se han visto gran cantidad de cambios y uno de los más trascendentales
ha sido el avance en el ámbito de la tecnología; al igual que todos los cambios, este ha tenido
sus pros y sus contras en cuanto a su aplicación.

Actualmente nos encontramos atravesando una dura situación con la pandemia, es así que los
programas informáticos como matlab, comsol, son verdaderas herramientas que aumentan
la productividad personal, académica, profesional y laboral. Para nuestro caso utilizaremos
el programa comsol, para realizar el estudio del flujo multidimensional en engranajes.

FUNDAMENTO TEORICO

Utilicemos de ejemplo un sólido el cual transfiere calor por conducción, si tomamos como
referencia el eje “x” y movemos el sólido en esa dirección, la temperatura variara en esa
dirección, originado de esta manera un flujo de calor conductivo en dirección “x” tendremos
, como se observa en la siguiente figura :

Y
qx
Figura 1

pág. 2
Ahora bien, puede suceder que si además de movernos dentro del solido en dirección-x
también nos movemos en dirección-y, además de variar la temperatura en dirección x,
también varié la temperatura en dirección “y” originando un flujo de calor conductivo en
dirección-y “ ”; también debemos analizar la posibilidad de movemos dentro del solido en
dirección “z”, también la temperatura además de varia en dirección “z” de esta manera se
originara el flujo de calor conductivo “ ”, como se representa en la siguiente figura:

X
qy

Y qx
Figura 2

El objetivo de los análisis de transferencia de calor suele ser determinar el flujo de calor o
bien la temperatura resultante de un flujo de calor. Si se conoce la distribución de
temperaturas en el material puede fácilmente determinarse el flujo de calor, este depende de
la dirección en la cual varia la temperatura, entonces obtendremos tres flujos conductivos en
las tres direcciones o coordenadas qx, qy, qz. Ello indica que el flujo de calor es una
magnitud vectorial que se encuentra como:

q= qx i+ qy j+ qz k

Si el flujo de calor tiene dos componentes que son cero y una diferente de cero se tiene flujo

unidimensional q = qx i y la temperatura depende solo de una variable en este caso

T=T(x)

pág. 3
Análisis para un flujo Bidimensional
Si el flujo de calor tiene una componente que es cero y dos diferentes de cero se tiene flujo
bidimensional el cual lo representaremos de la siguiente manera:

q= qx i+ qy j

La temperatura dependerá solo de estas dos variables, es así que obtendremos: T=T(x, y)

Análisis para un flujo Tridimensional

Si analizamos un cuerpo y si en este el flujo de calor tiene tres componentes diferentes de
cero se tiene flujo tridimensional, el cual se representa de la siguiente manera:

q= qx i+ qy j+ qz k
Z

qz

X
qy

Y qx
Figura 3

En este caso la temperatura dependerá de estas tres variables: T=T(x, y,z)

La Ley de Fourier
Para un flujo multidimensional tenemos la ley de Fourier nos indica

También lo podemos expresar:

pág. 4
Donde:

Flujo Multidimensional:
Ley de Fourier para Flujo Multidimensional
La ley de Fourier para un flujo multidimensional se expresa en términos de derivadas
parciales y no asi en forma de derivada total. Debemos analizar para cada eje de coordenada:

Para el flujo conductivo en dirección-x se tiene:

En flujo de calor se tiene:

Flujo de calor en y:

Flujo de calor en z:

Sumando los flujos de calor en las tres direcciones tenemos:

# # #
! " $% % ()
# #& #'

El término entre paréntesis es el vector gradiente de temperatura ∇ ; Obtendremos la

siguiente ecuación:
(q/A)=-K ∇T
Ley de Fourier para flujo multidimensional

pág. 5
Ecuación de energía en conducción multidimensional (flujo estacionario)
Utilizaremos el sistema de coordenadas cartesianas, para luego efectuar un balance de energía
sobre un volumen elemental

En el estado estacionario hay flujos de calor que entran y salen en las tres direcciones
coordenadas, como se muestra en la siguiente figura:

∆Z

∆X
∆Y

Figura 4

Realizando el balance de energía en estado estacionario obtendremos:

, - , ./ . ,/ , , - ,0 , , $ . ,/ ,
+ ! + !
$ $ . , . , $ $ . , . ,

∆ %∆ %∆

Utilizando la ley de Fourier

2
∆ 2
∆ ∆&∆'

2 2
∆ ∆
2 2

Sustituyendo obtenemos:

-∆&∆'∆ ( ! % ∆ ∆'∆ ( ! % ∆&∆ ∆ ( ! 0

Utilizando límites ∆ ≅ 0 ∆& ≅ 0, ∆' ≅ 0 obtendremos:

pág. 6
 # # # # # #
"( ) % "( ) % "( ) % 0
# # #& #& #' #'

Representa la ecuación la ecuación de conducción de calor para flujo multidimensional

estacionario en coordenadas cartesianas con conductividad variable con la temperatura.

PROCEDIMIENTO EXPERIMENTAL

Para el siguiente procedimiento se utilizó el programa COMSOL, en donde a la figura dada

que en este caso fue un piñón cilíndrico de dientes helicoidales se le pude crear una
simulación de cómo se comportaría en la vida real si este se llevara a distintas condiciones
de ambiente, para este caso en especifico la transferencia de calor con relación al tiempo para
ver los perfiles de temperatura que en este se desenvuelven.

Para esto necesitamos los datos del piñón (engrane) en un principio y el tipo de análisis el
cual llevaremos a cabo

Condiciones

- Flujo multidimensional

- Estado transitorio (variaciones en el tiempo)

Datos del Piñón

- Diámetro = 10 cm

- Numero de dientes = 16

Material

- Aleación de hierro y cobre (Acero)

Con todo esto también necesitaremos propiedades del acero a temperatura ambiente y presión
atmosférica que utilizamos para nuestra simulación en el programa

K = 50,2 w/m*k coeficiente de conductividad térmica del acero

D = 7800 Kg/m³ Densidad del acero

pág. 7
Cp = 460 J/Kg*K Calor especifico del acero

Finalmente, también utilizamos propiedades del aire a distintas temperaturas para verificar
el cambio de temperatura que existía en nuestro solido a analizar

H = 10 T = 293,15 K Coeficiente de conductividad térmica del aire a 20°C

H = 40 T = 263,15 K Coeficiente de conductividad térmica del aire a -10°C

Con todos estos datos vamos llenándolo las distintas interfaces que nos brinda el programa
COMSOL para luego realizar nuestro analizar en distintos puntos de nuestra figura. (esta
parte del informe se verá con más detalle en el video explicativo que se adjunta al presente
informe junto a las respectivas animaciones que nos dan un panorama más amplio de cómo
se desarrolla visualmente la transferencia de calor en nuestra figura).

INICIO

Abrir programa
COMSOL

Identifica r datos

Bus car en tabla

propiedades del
acero

Insertar datos

Realizar la
simulación

Analizar e
interpretar lo
observado

Fin

pág. 8
CALCULOS, RESULTADOS Y GRAFICOS

Simulación: Transferencia de calor para un tiempo igual a cero

Figura 5

En la figura observamos que para un tiempo igual a cero no existirá cambio en la

temperatura, ni flujo de calor.

pág. 9
Simulación: Transferencia de calor para un tiempo igual a un minuto

Figura 6

Para un tiempo igual a 1 minuto, observamos que existe una transferencia de calor en nuestra
figura, debido a la diferencia de temperaturas en los extremos del sólido.

pág. 10
Gráfica: Tomando como puntos de referencia los extremos de nuestra figura para un
tiempo igual a 1 minuto

Grafica 1

Grafica 2

pág. 11
Nuestra figura se encuentra sometido en sus extremos a temperaturas T= Tamb= 20 °C= 293
K y T= -10 °C= 263 K en un tiempo igual a cero. Para un tiempo igual a 1 minuto observamos
que la temperatura va disminuyendo conforme avanza en longitud del arco.

Gráfica: Tomando como puntos de referencia uno de los extremos y el punto medio
de nuestra figura para un tiempo igual a 1 minuto

Grafica 3

pág. 12
 Grafica 4

En la figura observamos el comportamiento de la temperatura, tomando como puntos de

referencia uno de los extremos y el punto medio de nuestra figura. En la gráfica vemos que
la temperatura va descendiendo de 293,1 a 291,7 esto puede ser debido al espesor que
presenta nuestra sólido.

pág. 13
Gráfica: Tomando como puntos de referencia uno de los extremos y un punto en las
aletas de nuestra figura para un tiempo igual a 1 minuto

Grafica 5

Grafica 6

pág. 14
Tomando como puntos de referencia uno de los extremos de la figura y un punto en las
aletas, encontramos que el comportamiento de la temperatura disminuye, presentando una
pendiente inclinada.

CONCLUSIONES, RECOMENSACIONES Y SUGERENCIAS

A través de la simulación en comsol logramos estudiar el flujo multidimensional en

engranajes, además de visualizar las direcciones que toma el flujo de calor dentro de un
cuerpo.

Se realizó la construcción de gráficas del flujo multidimensional en engranajes, gracias a los

datos brindados y otros encontrados en tablas.

Se tomaron diferentes tiempos, para así observar la trasferencia del flujo de calor en el solido.
Para un tiempo igual a 1 segundo, observamos que la transferencia calor se produce de
manera violenta, esto debido a la diferencia de temperaturas en los extremos del sólido.

Para un tiempo igual a 1 minuto observamos que la temperatura va disminuyendo conforme

avanza en longitud del arco, esto debido a que la temperatura en los extremos de la figura
permanece constante T= Tamb= 20 °C= 293 K y T= -10 °C= 263 K.

Tomando como otros puntos de estudio, uno de los extremos y el punto medio de nuestra
figura en tiempo igual a 1 minuto observamos que la temperatura va descendiendo de 293,1
a 291,7 esto puede ser debido al espesor pequeño que presenta nuestro sólido.

También se tomaron puntos de referencia a uno de los extremos de la figura y un punto en

las aletas, donde se observamos gracias a las gráficas que temperatura disminuye,
presentando una pendiente inclinada, donde una de las razones podría ser la forma elicoidal
que presenta nuestro piñón.

Se recomienda para la experimentación el manejo adecuado del programa, ya que al insertar

datos erróneos o al tanteo la construcción se dificulta.

pág. 15
CUESTIONARIO

1.-Se va a enfriar aire en la sección del evaporador de un refrigerador, al pasarlo sobre un

banco de tubos de 0.8 cm de diámetro exterior y 0.4 m de largo en el interior, de los cuales
se evapora el refrigerante a –20°C. El aire se aproxima al banco en la dirección perpendicular
a 0°C y 1 atm, con una velocidad media de 4 m/s. Los tubos se encuentran dispuestos en
forma alineada con pasos longitudinal y transversal de SL  ST  1.5 cm . Se tienen 30 filas
en la dirección del flujo con 15 tubos en cada una de ellas. Determine:
Propiedades: Asumiendo una temperatura de -5°C

k  0,02326 W m  K   1,317 kg m3
cp  d cos KJ kg  K Pr  0,737
5
  1,705  10 kg m  s Prs  PrTs 20 C  0,7408

La densidad a la entrada de temperatura de 0°C

Ts  20C
i  1, 292 kg m 3

V  4m s SL
Análisis: D  0,008m . SL  ST  1.5 cm T1  0C
V  4m s
ST
ST 0.015
Vmax  V  4 m s   8,571m s
ST  D 0,015  0,008

V D 1,317 kg m   8,571m s  0,008m 

3

Re D  max   5296
 1,705  10 5 kg m  s

 
0,25
NuD  0.27 Re0,63 Pr0,36 Pr Prs D
D

NuD  0, 27  5296   0,7375  0,7379 0,7408

0,63 0,36 0,25
 53,62
N L  16

2.- Un proceso de calentamiento de aire consiste en un tubo semicircular para el que la

superficie plana se mantiene a 1000 K y la otra superficie está bien aislada. El radio del tubo
es 20 mm, ambas superficies tienen una emisividad de 0.8. Si aire atmosférico fluye a través
del tubo a 0.01 Kg/s y Tm = 373 K, ¿Cuál es la rapidez a la que se debe suministrar calor por
10 m de longitud para mantener la superficie plana a 1000 K? ¿Cuál es la temperatura de la
superficie aislada?

pág. 16
 Aire
P=1 atm T2=?? Ɛ=0,8
Tm=373K

T1= 1000K Ɛ=0,8

Solución.
Es en estado estable, el flujo de aire es completamente desarrollado.
El ejercicio no nos dice a que temperatura entra el aire ni tampoco a la cual sale, pero si nos
dice que la temperatura media del aire dentro del tubo es de 373 K, en base a esta
temperatura se evaluara las propiedades del aire.

K  0.03095 [W / m*º C ] ;   2.181 * 105 [Kg / m * s ] ; CP  1009 [J / Kg *º C ] ; Pr  0.7111

Balance energético
Como la superficie semicircular está bien aislada y no hay adición externa de calor, un
balance de energía superficial da:

q2 rad  q2 rad

Dado que el recinto constituye un recinto de dos superficies, la transferencia neta de

radiación a la superficie 2 se puede evaluar a partir de la anterior ecuación mediante la
siguiente expresión:

pág. 17
 4 4
 (T 1  T 2)
 h * A2 (T2  Tm )
1  1 1 1 2
 
1 * A1 F12 * A1  2 * A2
F12  1
A1  0, 4[m 2 ]
A2  0, 62[m 2 ]
4* Asec cion transversal 2 *0, 02[m]
Dh    0, 0244[m]
P  2
 *V * Dh m * Dh (0, 01[kg / s ])* 0, 0244
Re   
 As ec c transv *   *(0, 02[m]) 2 * 2,181*105 [kg / m * s ]
2
Re  17835,1

Nu D  0, 023* Re0,8 Pr 0,4  0, 023*(17835,1)0,8 * (0, 7111)0,4

Nu D  50,52
h * Dh
Nu D   50, 52
K
50, 52* 0, 03095[W / m 2 * K ]
h  64,1[W / m * K ]
0, 0244m

5,67*108[W / m2 * K 4 ](10004  T24 )

 (64,1[W / m2 * K ])(0,62[m2 ])(T2  373)
1  0,8 1 1  0,8
 
0,8*0,4[m2 ] 1*0,4[m2 ] 0,8*0,62[m2 ]
1,60704*108 T24  39,742T2  30894,16  0

T2  777, 36[ K ]

De un balance de energia en la superficie caliente:

q1exter  q1rad  q1conv  q2 conv  q1conv
q1exter  q2 conv  q1conv  h * A1 *(T1  Tm )  h * A2 *(T2  Tm )
q1exter  64,1[W / m 2 * K ]*0, 62[m2 ]*(1000  373) K  64,1[W / m 2 * K ]*0, 62m2 *(777,36  373) K

q1exter  32,15[ KW ]

3.- Se tiene un tanque de agua caliente de 50 cm de diámetro y 95 cm de largo, el cual está

colocado sobre el techo de una casa. El agua que se encuentra en su interior se calienta
durante el día hasta 80 oC mediante un colector solar de placa plana. Entonces durante la
noche el tanque se expone al viento con una temperatura del aire de 18 oC y una velocidad
promedio de 40 km/h. Estime la temperatura del tanque después de un periodo de 45 min.
Suponga que la superficie del tanque está a la misma temperatura que el agua que se

pág. 18
encuentra en su interior y que el coeficiente de transferencia de calor sobre las superficies
superior e inferior es igual al correspondiente a la superficie lateral.

TANQUE DE AGUA
D=50

AIRE
V∞=40 [Km/h]

Propiedades Las propiedades del agua a 80 °C

  971,8[kg / m3 ]
Cp  4197[ J / kg C ]

Las propiedades del aire a 1 atm y suponiendo una temperatura de film de 50°C
k  0, 02735[W / mC ]
  1, 798*105 [m 2 / s ]
Pr  0, 7228

Solución
40*1000
[m / s]
V D 3600
Re    309, 015
 1, 798*105 [ m 2 / s ]

El número de Nusselt es:

4/5
0,62*Re0,5 *Pr1/3   Re  
5/8

Nu  0,3  1  
[1  (0, 4 / Pr)2/3 ]1/4   282,000  
5/8 4/5
0,62*309,0150,5 *0,72281/3   309,015  
Nu  0,3  1  
[1  (0, 4 / 0,7228)2/3 ]1/4   282,000  
Nu  484,9

pág. 19
El coeficiente convectivo de transferencia de calor
k 0, 02735[W / mC ]
h * Nu  * 484, 9  26, 53[W / m 2 C ]
D 0,5[m]

El área superficial del tanque

D2 (0,5)2
As   DL  2   (0,5)* (0, 95)  2  1,885[m 2 ]
4 4
La velocidad de transferencia de calor
 80  T2
Q  hAs (Ts  T )  (26,53[W / m 2 C ]) *(1,885[ m 2 ]) *(  18)
2
Donde T2 es la temperatura final del agua (80+T2)/2 se obtiene la temperatura promedio
del agua durante el proceso de enfriamiento

D2 (0,5) 2
m   *V   L  971,8[ kg / m3 ]*  * *0,95[ m]  181, 27[ kg ]
4 4
Q  mCp (T2  T1 )  (181, 27[ kg ]) *(4197[ J / kg C ]) *(80  T2 )C

La velocidad promedio de la transferencia de calor es:

 Q 181, 27[kg ]* 4197[ J / kg C ]*(80  T2 )C
Q 
t 45*(60[ s])

 80  T2  181, 27[kg ]* 4197[ J / kg C ]*(80  T2 )C
Q  (26,53[W / m 2 C ])(1,885m 2 ) *   18  C 
 2  45*(60[ s])

T2  69,9C

4.- Considere la placa base de una plancha domestica de 800 W con un espesor de L=0,6 cm
area de la base de A=1690 cm2 y conductividad termica de k= 200W/m oC. La superficie
interior de la placa base se sujeta a un flujo uniforme de calor generado por los calentadores
de resistencia interior. Cuando se alcanza la condciones estables de operación, la temperatura
de la superficie exterior de la placa es de 85 oC. Descartando cualquier perdida de calor a
traves de la pared superior de la plancha.
Determinar:
a) Exprese la ecuacion diferencial y las condiciones de frontera para la
condicion unidimensional de calor en estado estable a traves de la plancha.
b) Obtenga una relacion para la variacion de la temperatura en la placa
base, resolviendo la ecuacion diferencial.
c) Evalue la temperatura de la superficie interior

pág. 20
 
 Q 800[W ]
q0   4 2
 50000[W / m 2 ]
Abase 160*10 [m ]
d 2T
0
dx 2
dT (0) 
k  q0  50000[W / m 2 ]
dx
T ( L)  T2  85C
dT
 C1
dx
T ( x)  C1 x  C2

C1 y C2 son constantes arbitrarias


 q
x0 kC1  q 0  C1   0
k
xL

q L
T ( L)  C1 L  C2  T2  C2  T2  C1 L  C2  T2  0
k
  
q q L q ( L  x)
T ( x)   0 x  T2  0  0  T2
k k k
(50000[W / m 2 ])(0, 006  x)[m]
T ( x)   85[C ]
20[W / mC ]
T ( x)  2500(0, 006  x)  85

x0
T (0)  2500(0, 006  0)  85  100C

T (0)  100C

Ejercicio de nivel intermedio

5.- Un conducto hueco de sección transversal cuadrada de dimensiones internas de 10*10
(cm) está construido con ladrillo de k = 0.21 [W/m°C] con un espesor de 10 [cm]. En
condiciones de equilibrio la temperatura interna y externa es de 300°C y 30°C
respectivamente. Estimar la perdida de calor a través del conducto.

pág. 21
 DATOS

Tw1  300C

Tw 2  30C

W
k L  0.21
mC
x  5cm  x  y

En el nodo 1 30  300  2T2  4T1  0

En el nodo 2 T1  30  T1  T3  4T2  0

En el nodo 3 2  30  2T2  4T3  0

La solución de las tres ecuaciones es:

T1  153.75C T2  142.5C T3  86.25C

El calor que entra en la pared

 T  T   W
QZe  8k L   w1 1  Tw1  T2   QZe  387.45
 2  m

El calor que sale de la pared

 T  T   W
QZs  8k L   1 w 2  T2  Tw 2   T3  Tw 2  QZs  387.45 
 2  m

6.- Hallar el flujo de calor de una chimenea cuyo interior fluye gases de combustión de tal
manera que en el interior tiene 371°C y la superficie exterior está a 38°C las dimensiones de
la chimenea es de 60*30[cm] y el espesor es de 30[cm] está construido de ladrillo de
conductividad k = 1.2 [W/m°C].

pág. 22
DATOS:
W
T1w  371C T2 w  38C kC  1.2 y  15cm
mC

 T 
Q   k  A 
n  x 

Q
El área A  x  y    k  T 
Z n

Analizando en el nodo 1 2T2  38  371  4T1  0

Analizando en el nodo 2 T1  371  38  T3  4T2  0

Analizando en el nodo 3 T2  T4  2  38  4T3  0

Analizando en el nodo 4 T3  T5  371  38  4T4  0

Analizando en el nodo 5 T4  T6  371  38  4T5  0

Analizando en el nodo 6 2T5  371  38  4T6  0

Resolviendo estas ecuaciones se tiene:

T1  190.69C T2  176.88C T3  107.84C

T4  178.48C T5  197.06C T6  200.72C

El calor que entra al conducto

T  T T T 
QZi  4 kC   1w 1  T1w  T2   T1w  T4   T1w  T5   1w 6 
 2 2 

W
QZi  3532.2
m

El calor que sale del conducto

T  T T T 
QZe  4 kC   1 2 w  T2  T2 w   2 T3  T2 w   T4  T2 w    T5  T2 w   6 2 w 
 2 2 

pág. 23
 W
QZe  3531.864
m

7.- Una chimenea de sección cuadrada de 20*20[cm] está construida con ladrillo k=0.81 [W/
m°C] de 10 [cm] de espesor, los gases de la chimenea mantiene la temperatura interior de la
chimenea a 280°C el exterior está compuesto a un ambiente cuya temperatura es de 23°C y
un coeficiente de convección de 10[W/m°C]. Encontrar el flujo de calor a través de la
chimenea.

DATOS
W w
w  10cm k  0.81 h  10
m  C m  C
2

Tw1  280C T  23C

Nodo 1 2T2  280  T4  4T1  0

Nodo 2 T1  T5  T3  280  4T2  0

Nodo 3 2T6  2T2  4T3  0

2  h  w  h  w 
Nodo 4 2T1  2T5   T  2   2   T4  0
k  k 

2  h  w  h  w 
Nodo 5 2T2  T4  T6   T  2   2   T5  0
k  k 

pág. 24
 2  h  w  h  w 
Nodo 6 2T3  T5  T7   T  2   2   T6  0
k  k 

2  h  w  h  w 
Nodo 7 T6  T6   T  2   1  T7  0
k  k 

Resolviendo estas ecuaciones

T1  172.702C T2  161.566C T3  112.0C T4  87.497C

T5  81.919C T6  62.345C T7  40.607C

El calor que entra en la pared

 T  T   QZe  139.38723
W
QZe  k   w1 1  Tw1  T2   m
 2 
El calor que sale de la pared

 T  T   W
QZs  k   1 4  T2  T5   T3  T6   QZs  139.24265 
 2  m

8.- Determine el flujo de calor por unidad de profundidad en el segmento circular de la

figura siguiente. Supóngase que la conductividad térmica del material es de 0.7 [W/m°C].
Una de las superficies es isotérmica y las otras se encuentran a 100°C y 25°C.

DATOS
w r1  50cm r2  54cm
k  0.7
mC

T2  100C T1  25C Lw  1m

Por balance de energía

Eentrada  Egenerada  Esale  Ealmacenada

pág. 25
Condiciones
d
Eentrada  Q Egenerada  0 Esale  Q  Q d Ealmacenada  mCpT
d
 d  dT
Q   Q  Q d   mCpT T  r    r1
 d  d

dT dT
Q   kA  k  r  Lw  V  r  L  r1  d  m   V
dx r1  d 

d  dT   dT dT
    k  r  Lw    d    V  Cp      r  L  r1  d    Cp 
 d  r1  d    d d

1 d   d  dT
  k   T     Cp
r1 d   d  
2
d

Para estado estable

d   d 
T 
T2  T1     T
k  T  0
d    d     2

Su solución

C1 dT (T2  T1 )
T     C2 
k d 
Por condiciones de frontera

C2  T1 k
C1   T2  T1 

dT ( r  r )  Lw ( r1  r2 )  Lw
Q   k  A  k 1 2  T1  T2  Q  k   T1  T2 
r1  d    r1   r1

Q  1.3369  W

9. La distribución de temperaturas a través de una pared plana de 1 metro de espesor y área

frontal igual a 15 m2, es:
t(x) = a + b ∙ x + c x3

pág. 26
Siendo el valor de las constantes el siguiente: a = 500 ºC; b = -150 ºC/m; c = - 25 ºC/m3.
Existe una generación uniforme de energía igual a 500 W/m3. Las propiedades del material
son: densidad ρ = 1800 kg/m3; conductividad k = 30 W/mºC. Calor específico cp = 3 ⋅ 103
J/kg ºC.
Se pide determinar:
a. El calor que entra y sale del material.
b. El calor almacenado en el material.
c.La rapidez respecto al tiempo del cambio de temperatura del material en la coordenada x =
0.5.
SOL.
De acuerdo con el enunciado, se trata de una conducción unidimensional en la dirección x.
Según eso:
a)
t
Qx   k  A    k  A  (b  3  c  x 2 )
x

W  C 
Qentrada  Qx  0    k  A  b  30  15m2   150   67500W
mC  m

Qsalida  Qx 1  k  A  (b  3  c )  101250W

b) Qalmacenado  Qentrada  Qgenerado  Qsalida

Sustituyendo valores
Qalmacenad  Qentrada  Qgenerado  Qsalida  67500  500  15  1  101250  26250W

c) la tasa de cambio de la temperatura

t k  2t 1
  2 Q
   c p x   cp

Se realiza la derivada segunda

 2t 2
 (a  b  x  c  x3 )  6  c  x
x 2  x 2
Y se sustituyen valores:
t 30 1 C
  6   25   0.5   500  3.24  10 4
 1800  3  10 3
1800  3  10 3
s

pág. 27
10.- Se considera una placa plana perfectamente aislada por una de sus superficies, mientras
que la otra se encuentra expuesta al ambiente, recibiendo una radiación solar igual a 1000
W/m2. La temperatura ambiente es igual a 25ºC, siendo la temperatura de cielo equivalente
igual a 15ºC. El coeficiente de transmisión de calor por convección es igual a 10 W/m2ºC. Se
pide calcular la temperatura de equilibrio de la placa.
Datos: la absortividad de la placa es igual a su emisividad e igual a 0.8.
Solución

En equilibrio, el balance de calor a la placa queda como sigue:

  G  h T  T f       T 4  Tc ,eq 4 

Donde T es la incógnita que queremos calcular.

Sustituyendo valores:


0.8  1000  10  T   25  273   0.8  5.67  108  T 4  15  273
4

T  345K  72C

11.- Una zanahoria cuya forma puede aproximarse con un cilindro de 8 pulgadas de longitud
y ¾ pulgadas de diámetro, inicialmente está a 70°F de temperatura y luego se deja caer en
agua hirviendo a la presión atmosférica. ¿Durante cuánto tiempo debe cocinarse para llegar
a 195°F?
Datos:
k=0,28 Btu/h pie°F; Cp=0,95 Btu/lb°F; ρ=64 lb/pie3; h=350 Btu/h pie2°F
Solución:
Agua hirviendo: T=195°F

T  212 F


 D2  l
L 4
 2 
2  D l    Dl
4 

pág. 28
 2
 3 1  8
   
4  4 12  12
L  0.015 pies
   3 1 2 8  3 1 8
2             
 4  4 12  12  4 12 12

hL 350  0,015
Bi    Bi  18,75
k 0,28

Como 18,75 > 0,1 se realizara método analítico grafico

k 0,28
     0,005
 C p 64  0.95

k 0,28 r 0
m   0,0256 n  0
rm h  0,75  rm rm
350   
 2  12 
T  T 195  212  0,005
Y   0,1197 X    5,12  
Tb  T 70  212 2
r
m  0,75 
2

 
 2  12 
Con los valores calculados hallamos el valor de X de grafica para cilindro:
X  0,46  0,46  5,12  

  0,1 Hrs   6  min 

12.- Determinar el tiempo de enfriamiento de 1000°F a 200°F, de una esfera de acero dulce
de 0,5 pulgadas de diámetro, sometida a un flujo de aire frio a 80°F que produce un h=20
Btu/Hr.pie2. °F.
La dimensión lineal característica es L=R/3. La conductividad térmica del acero dulce a
una temperatura promedio de 600°F = (1000°F+200°F)/2, es aproximadamente igual a 25
Btu/Hr pie °F.
Solución:

hL (20)  (0.25 /  3  12  pie )

Bi   0,0056  0,0056  0,1 analisis debloque
k 25
200  80
 0,1304  e  ( Bi )( Fo )   Bi  Fo   2,0369
1000  80

pág. 29
 2
R
363,73  
2,0369  3
( Fo )   363,73  2   
0,0056 L 
De tablas tenemos que:
Btu
25
Hr  pie   F pie2
  0,46
 lbm   Btu  Hr
 490 pie3   0,11 lbm   F 
  

Entonces:
2
 0,25 
(363,73)   pie2
  3  12   0,0381Hr  2,29 min
0,46 pie2
Hr
  2,29 min

13.- Para una placa de acero de gran superficie, de espesor de 30 cm. hallar la temperatura
en su plano medio después de haber transcurrido un tiempo de 14,85 mins.
Datos:

8 =371 °C, 9 =37,8 °C, : 7850 (>/-3, Cp =0,13 kcal/kg°C, k  37,2 kcal / m  h  C ,
considerar la resistencia despreciable h = ∞
Solución:
hL  L
Bi      0,1 temp. estratificada
k k

T  T 4   a1x  x 1 9 a1x 3 x 1 25a1x 5 rm 

 e sin  e sin  e sin 
T  Tb π  2rm 3 2rm 5 2rm 

37,8  T 4 r 1 3 rm 1 25a1x 5 rm 
  e a1x sin m  e9 a1x sin  e sin 
37,8  371   2rm 3 2rm 5 2rm 
2
  k 37,2 m2
a1    ;     0,036
2  C p 7850 0,13 hrs

pág. 30
  14,85 
0,036   
a  60   0,396
x 2 
rm 0,152

37,8  T 4  ( 2 )2 0,396  rm 1 9( 2 )2 0,396 3 rm 

 e sin  e sin 
37,8  371   2rm 3 2rm 

37,8  T
 0,48  T  197.7 C T  197.7 C
37,8  371

14.- Se tiene una placa de 10 [cm] de espesor a una temperatura uniforme de 20°C que se
introduce en un medio a 100 °C adquiriendo instantáneamente esta temperatura. Determinar
mediante técnicas numéricas el tiempo necesario para que el plano medio de la pieza alcance
una temperatura de 60 °C.

Datos
BCDE
T1 = 20 °C ∆A F

EM
T2 = 100 °C ∝ 6 ∙ 10KL
N

T = 60 °C
En la ecuación general de la conducción para coordenadas rectangulares y sin generación.

d2 1  d  d2 Ta  2T  Tb
T   T   T  T 1  Ta  Tb
dx 2
  d  dx 2
X 2

Ta  2T  Tb T 1  T 1  d  T 1  2 

X 2

X 2
 
  d
r 

 T1 T1  T
dT  
0 X 2
d

 T 2  T  2    X 2  T2 T 
ln        ln      144,40566 s
 T 2  T1  X 2
2  T 2  T1 

  144,40566 s

pág. 31
15.- A un tubo de 40 mm de diámetro exterior se le adosan aletas anulares de aluminio de 0.5
mm de espesor y 100 mm de radio exterior, separadas entre sí 5 mm. Las aletas están aisladas
térmicamente en su extremo. El coeficiente de película exterior vale 70 W/m2ºC.
Sabiendo que la temperatura de la superficie del tubo es de 70ºC y la del medio exterior de
20ºC, determinar:
a) El calor disipado por metro de tubería con aletas
b) La temperatura en el borde exterior de la aleta

Por ser aletas de aluminio, la conductividad vale sacado de (libro de tablas):

k = 236 W/m°C
a) Calor disipado por metro de tubería con aletas

Q W    A  At f  n f   ho   to  t f ,o 

Como nos pide por metro lineal de tubería

Q  W   At Af 
  
L  m   L L
 nf   ho   to  t f ,o 

El área del tubo es:
At
At  2  ro  1  n     Lt   2    ro  1  n   
Lt

El área de superficie aleteada es:

A f  2   re2  r02   n  Lt   2     re2  r02   n

At
Lt

No tenemos el dato del número de aletas por metro lineal de tubería (n). Sin embargo,
tenemos el dato de separación entre planos medios de aletas  1 

Se cumple que:

pág. 32
 1
n
1

Sustituyendo en las ecuaciones de las áreas, se tiene:

At
 0.113 m
Lt

El área de superficie aleteada es:

Af
 12.06 m
Lt

Calculo del rendimiento de la aleta anular:

El rendimiento de la aleta anular se calcula gráficamente a partir de la figura.
Para ello es necesario el cálculo de siguientes parámetros:

ro 2  ho
c   0.19  c  rec  n  rec   3,45
rec k 
n f  0.2

Sustituyendo en la fórmula del calor disipado por metro lineal de tubería con aletas, se
tiene:

Q W  W 
    0.113  12.06  0.2   70   70  20   8837.5  
Lm m
b) Temperatura en el borde extremo de la aleta:
Aplicando la ecuación de distribución de temperatura en una aleta anular:

 t  t f ,o I0  n  r   K1  n  re   K 0  n  r   I1  n  re 
 
 o to  t f ,o I 0  n  r0   K1  n  re   K0  n  r0   I1  n  re 

En este caso, la ecuación se calcula en r=rec (t=tc)

 te  t f ,o I 0  n  rec   K1  n  rec   K0  n  rec   I1  n  rec 

 
 o to  t f ,o I0  n  r0   K1  n  rec   K0  n  r0   I1  n  rec 

Empleando la tabla tabulada de la ecuación de Bessel.

Para ello necesitamos calcular:

pág. 33
 2  h0 2  70
n   34.44
k  236  0.5e  3
n  rec  3.45 ; n  ro  0.69

Sustituyendo e interpolando

I0  n  r0   I 0  0.69   1.12

I 0  n  rec   I 0  3.45  7.07

I1  n  rec   I1  3.45  5.93

K0  n  r0   K0  0.69   0.67

K0  n  rec   K0  3.45  0.021

K1  n  rec   K1  3.45  0.024

Sustituyendo en la ecuación
7.07  0.024  0.021 5.93
te  20  50   23.64C
1.12  0.024  0.67  5.93

16.- ¿Cuál es la carga de calor generada en un restaurante por una parrilla de 1m

por 0,8m que se mantiene a 134 ºC? La temperatura ambiente es 20 ºC.
Solución:
Los parámetros apropiados, evaluados a la temperatura media d película:
Ts  T
Tf   77C
2
Son:

m2
v   22,38 105  9,29  10 2   2,08  105
s
W
k   0,017351,729   0,0300
m  K
Pr  0,697

1 1
Para aire:     2,857 10 3   K 1 , y el número de grashof es:
T 350

g  (Ts  T ) L3
GrL 
v2

pág. 34
Donde L = 0,9m, es el valor promedio de las dos dimensiones de la parrilla.

GrL 
 9,8m / s  2,857 10 K  114K  0,9m
2 3 1 3

 5,38 109
4

 2,08 10  ms 5 2
2

Utilizando la ecuación correspondiente para una placa caliente mirando hacia arriba,
obtenemos:

1
hL
 (0,14)(GrL Pr) 3
k
 1

 0,14   5,38 109   0,697   3  7,25 2

0,0300 W
h
0,9 m  K

Y la rapidez de transferencia de calor es:


q  h A Ts  T    7,25W / m2  K  (1 0,8)m2  (114 K )  661,2W

q  661,2W
17.- Unos cojinetes de bolas de acero inoxidable (AISI 304), que se calientan de manera
uniforme a 850ºC, se endurecen al contemplar en un baño de aceite que se mantiene a 40ºC.
el diámetro de la bola es 20 mm. y el coeficiente de convección asociado con el baño de
aceite es 1000 W/m2K.
a) Si el templado no va a ocurrir sino hasta que la temperatura de la superficie de las bolas
alcance 100ºC, ¿Cuánto tiempo deben permanecer estas en el aceite? ¿Cuál es la
temperatura del centro al final del periodo de enfriamiento?
b) Si se templan 10000 bolas por hora ¿Cuál es la velocidad a la que el sistema de
enfriamiento de baño de aceite debe quitar energía a fin de mantener su temperatura a
40ºC?
solucion :

pág. 35
PROPIEDADES
T  500C : k  22, 2[W / mK ]
Cp  579[ J / kgK ]
  7900[kg / m3 ]
  4,85*10 6 [m 2 / s ]
PARA EL BIOT

h(r0 / 3) 1000[W / m 2 K ]*(0, 010 / 3)

Bi    0,15
k 22, 2[W / mK ]
k 22, 2[W / mK ]
Bi 1    2, 22
hr0 1000[W / m 2 K ]* 0, 01[ m]
 (r0 , t )
 0,80
0 (t )
 (r0 , t ) T (r0 , t )  T 100  40
   0, 074
0 (t ) Ti  T 850  40
0  (r0 , t ) / i 0, 074
   0, 093
i  (r0 , t ) /  0 0,80

Para:

0 / i  0, 093; Bi 1  2, 22
t*  Fo  2, 0
r0 2 Fo (0, 01m) 2 (2)
t   41[ s ]
 4,85*106 [m2 / s ]
0  T0  T  0, 093(Ti  T )  0, 093(850  40)
T0  115C

Para:

Bi 2 Fo  (1/ 2, 2) 2 * 2  0, 41
Bi  (hr0 / k )  0, 45
Q
 0, 93
Q0
Q  0,93 CpV (Ti  T )
Q  0,93*7900[kg / m3 ]*579[ J / kgK ]*  / 6* (0, 02[m])3 *810C
Q  1, 44*104 [ J ]

pág. 36
Para la transferencia de energia

q  104 Q / 3600s
q  4*104W

q  40[ KW ]

18.- Considere el sistema del problema 2. La superficie interior se expone a gases calientes
a 350 °C con un coeficiente de convección de 100 W/m2 K, mientras la superficie exterior
experimenta convección con aire a 25°C y un coeficiente de convección de 5 W/m2 K.
(a) Para un espaciado de malla de 75 mm, calcule el campo de temperaturas dentro del
sistema y determine la pérdida de calor por unidad de longitud por convección desde la
superficie externa del escape al aire. Compare este resultado con el calor ganado por
convección desde los gases calientes al aire.
(b)Determine el efecto del espaciado de la malla sobre el campo de temperaturas y la
pérdida de calor por unidad de longitud al aire. De manera específica, considere un
espaciado de malla de 25 mm. Para Δx = Δy = 25 mm, explore el efecto de cambios en los
coeficientes de convección sobre el campo de temperaturas, y la pérdida de calor.
SOLUCION.-

2hi x  h x 
 2T4  T2  T2   T ,i   i  2  T1  0
k  k 

2hi x  h 
 2T5  T3  T1   T , i  2  i x  2  T2  0
k  k 

2hi x  h x 
2 T6  T6   T2  T2   T ,i  2  3  i  T3  0
k  k 

pág. 37
 T8  T5  T1  T5   4T4  0
T9  T6  T2  T4   4T5  0
T10  T7  T3  T5   4T6  0

T11  T11  T6  T6   4T7  0

2ho x  h x 
 2T4  T9  T9   T , o  2  o  2  T8  0
k  k 

2ho x  h x 
 2T5  T70  T8   T , o  2  o  2  T9  0
k  k 

2ho x  h x 
 2T6  T11  T9   T, o  2  o  2  T10  0
k  k 

2ho x  h x 
 2T7  T12  T10   T, o  2  o  2  T11  0
k  k 

2ho x  h x 
T11  T11   T, o  2  o  1 T12  0
k  K 

T1k  0,09239T2k 1  0,09239T4k 1  285,3

T2k  0,04620T1k  0,04620T3k 1  0,09239T5k 1  285,3

T3k  0,8457T2k  0,1692T6k 1  261, 2

T4k  0, 25T1k  0,50T5k 1  0, 25T8k 1

T5k  0, 25T2k  0, 25T4k  0, 25T6k 1  0, 25T9k 1

T6k  0, 25T3k  0, 25T5k  0, 25T7k 1  0, 25T9k 1

T7k  0,50T6k  0,50T11k 1

T8k  0, 4096T4k  0, 4096T9k 1  4,52

T9k  0, 4096T6R  0, 2048T8k  0, 2048T10k 1  4,52

T10k  0, 4096T6k  0, 2048T9k  0, 2048T11k 1  4,52

pág. 38
 T11k  0, 4096T7k  0, 2048T10k  0, 2048T12k 1  4,52

T12k  0,6939T11k  7,65

k T1 T2 T3 T4 T5 T6 T7 T8 T9 T10 T11 T12

0 340 330 315 250 225 205 195 160 150 140 125 110
1 338, 336, 324, 237, 232, 225, 175, 163, 161, 155, 130, 98,3
9 3 3 2 1 4 2 1 7 6 7
2 338, 337, 328, 241, 241, 226, 178, 169, 170, 158, 130, 98,1
3 4 0 4 5 6 6 6 0 9 4
3 338, 338, 328, 247, 245, 230, 180, 175, 173, 161, 131, 98,9
8 4 2 7 7 6 5 6 7 2 6
4 339, 338, 328, 251, 248, 232, 182, 178, 176, 162, 132, 99,8
4 8 9 6 7 9 3 7 0 9 8
5 339, 339, 329, 254, 250, 234, 183, 180, 177, 164, 133, 100,
8 2 3 0 5 5 7 6 5 1 8 5
6 340, 339, 329, 255, 251, 235, 184, 181, 178, 164, 134, 101,
1 4 7 4 7 7 7 8 5 7 5 0
7 340, 339, 329, 256, 252, 236, 185, 182, 179, 165, 135, 101,
3 5 9 5 5 4 5 7 1 6 1 4

1 1 
q  h0  x  T8  T , o   T9  T, o   T10  T, o   T11  T , o   T12  T , o  
2 2 

W 1 1 
q  5  0, 075m  (182, 7  25)  (179,1  25)  (165, 6  25)  (135,1  25)  (101, 4  25)  C
m K
2
2 2 
q  195W m

1 1 

qtot  8  hi x  T ,i  T1   T ,i  T2   T ,i  T3    8  190,5W m  1,52 kW m
2 2 
b)

(0,0) (300,0)
x (mm)

y (mm) Δx + Δy =25mm
(mm)

pág. 39 (0,150) (150,150)

 y\x 0 25 50 75 100 125 150 175 200 225 250 275 30
0
0 180, 180, 178, 175, 171, 165, 158, 149, 140, 129, 119, 108, 98,
7 2 4 4 1 3 1 6 1 9 4 7 0
25 204, 203, 198, 198, 193, 186, 178, 168, 157, 145, 133, 121,
2 6 2 2 3 7 3 4 4 6 4 0
50 228, 228, 222, 222, 217, 209, 200, 188, 175, 161, 147,
9 3 6 6 2 7 1 4 4 6 5
75 255, 254, 252, 252, 243, 235, 223, 209, 194, 177,
0 4 4 4 1 0 9 8 1 8
10 282, 281, 280, 280, 271, 263, 250, 232, 213,
0 4 8 1 1 6 3 5 8 5
12 310, 310, 309, 309, 303, 296, 282, 257,
5 9 5 3 3 2 0 2 5
15 340, 340, 339, 339, 337, 335, 324,
0 0 0 6 6 9 3 7

q  8hi x T ,i  T1  2  T, i  T2   T ,i  T3   T, i  T4   T ,i  T5   T, i  T6   T ,i  T7  2 

q  1,52 kW m

BIBLIOGRAFIA

Ruge, A.(2018). Flujo Multimensional obtenido de:

https://www.studocu.com/co/document/universidad-nacional-de-
colombia/transferencia-de-calor/trabajo-tutorial/flujo-
multidimensional/9498757/view

Holman J. ( 1999) Transferencia de calor Octava edición . España: Concepción Fernández

Madrid

Tutorías con ingenio Universidad Nacional. (2016) método numérico flujo multidimensional
obtenido de: https://www.youtube.com/watch?v=bgGskM3K75g

pág. 40