CONVECCIÓN

PROBLEMA 1
Determinar el coeficiente de transmisión de calor por convección hacia la superficie
de una barra de acero 15% Cr 10% Ni, que se encuentra a 160ºC, y se calienta en
un horno con una corriente de aire a 600ºC y una velocidad de 10 m/s, perpendicular
al eje de la barra. La longitud de la barra es igual a 1.5 metros y el diámetro es igual
a 20 cm.
ts=160; %ºC
tf=600; %ºC
u=10; %m/s
l=1.5; %m
d=0.2; %m
%Propiedades aire a tm=380ºC
tm=(tf+ts)/2; %ºC
cp=1064; %J/kg/ºC
rho=0.5331; %kg/m3
my=32.61e-06; %Pa*s
k=49.1e-03; %W/m/ºC
Pr=0.7066;
Prf=0.7224;
Prs=0.6982;
Re=rho*u*d/my; %Re=3.2695e+04
C=0.26;
m=0.6;
n=0.37;
h=k*(C*Re^m*Pr^n*(Prf/Prs)^0.25)/d;

Solución
h = 28.9475 W/m2/ºC

PROBLEMA 2
Se quiere determinar el coeficiente de convección existente en torno a una placa
metálica de pequeño espesor, suspendida verticalmente, mediante la medida
experimental del cambio de temperatura durante un transitorio de enfriamiento.
Asumiendo que la placa no presenta en su interior un gradiente térmico significativo
y despreciando la transmisión de calor por radiación de la placa al recinto,
determinar el coeficiente de convección existente en el instante en que la placa tiene

10
una temperatura de 65°C y presenta una tasa de enfriamiento de 0,05°C/s por efecto
del aire a 25°C.
DATOS:
• Placa:
Dimensiones = 0.3 x 0.3 m
Masa = 3,75 kg
cp = 2770 J/kg°C

%% PROBLEMA
ts=65; %ºC
tf=25; %ºC
deltat_deltatau=0.05; %ºC/s
%Placa
A=0.3*0.3; %m2
m=3.75; %kg
cp=2770; %J/kg/ºC
%Calor por unidad de tiempo
Q_deltatau=m*cp*deltat_deltatau; %W
%La placa se enfría exclusivamente por convección, luego:
h=Q_deltatau/A/(ts-tf);

Solución
h = 144.2708 W/m2/ºC

PROBLEMA 3
Una tubería de acero al carbono ordinario, de 102 mm de diámetro exterior y
espesor normal (3.5 mm), se encuentra aislada mediante una capa de fibra de vidrio
de densidad 32 kg/m3 de 20 mm de espesor. Si por el interior de la tubería circula
agua sobrecalentada a 120°C a 1 m/s, calcular las pérdidas de calor por metro lineal,
cuando el viento sopla en dirección perpendicular a la tubería a una temperatura de
-10°C y una velocidad de 10 km/h.

r1=(102/2)*1e-03; %m
r2=r1+3.5e-03; %m
r3=r2+20e-03; %m

11
k12=43; %W/m/ºC. Conductividad del acero al carbono (pag. 77)
k23=0.034; %W/m/ºC. Conductividad de la fibra de vidrio (pag. 87)
%Interior: agua sobrecalentada
t_i=120; %ºC
v_i=1; %m/s
%Propiedades del agua líquida sobrecalentada a 120ºC (p.59) (Suponemos
que
%la temperatura media de masa es 120ºC)
cp_i=4246; %J/kg/ºC
rho_i=942.8; %kg/m3
nu_i=0.2321e-03; %kg/m/s
k_i=0.6833; %W/m/ºC
Pr_i=1.44;
%Exterior: aire
t_e=10; %ºC
v_e=10e03/3600; %m/s
%Propiedades del aire a 10ºC (p.21). Suponemos que la temperatura media
de
%película es 10ºC. Luego habrá que iterar
cp_e=1005.8; %J/kg/ºC
rho_e=1.2467; %kg/m3
nu_e=17.69e-06; %kg/m/s
k_e=24.87e-03; %W/m/ºC
Pr_e=0.716;
%%RESOLUCIÓN
%Calculo del calor disipado:
Q/l=2*pi*(ti-te)/(1/r1/hi+log(r2/r1)/k12+log(r3/r2)/k23+1/r3/he)

%1. Cálculo de hi: convección forzada en el interior de tubería

lc_i=2*r1;
Re_i=rho_i*v_i*lc_i/nu_i; %(>2300: convección forzada)
hi=k_i*0.023*Re_i^0.8*Pr_i^0.3/lc_i; %Dittus Boelter. Suponemos
enfriamiento del fluido y que ts-tb < 6ºC (al ser agua, fluido poco
viscoso, se va a cumplir
%2. Cálculo de he: convección forzada, flujo externo perpendicular a
una Tubería

12
lc_e=2*r3;
Re_e=rho_e*v_e*lc_e/nu_e; %(>2300: convección forzada. Re_e>10^6:
aplicamos Churchill y Berstein
he=k_e*(0.3+0.62*Re_e^0.5*Pr_e^(1/3)*(1+(Re_e/2.82e05)^(5/8))^(4/5)/(1
+(0.4/Pr_e)^(2/3))^0.25)/lc_e;
%Calculo del calor total
Q_l=2*pi*(t_i-t_e)/(1/r1/hi+log(r2/r1)/k12+log(r3/r2)/k23+1/r3/he);
%W/m
%Comprobación temperatura t3: Q_l=2*pi*r3*he*(t3-te);
t3=t_e+Q_l/(2*pi*r3*he);
%Propiedades del aire a tm=(t_e+t3)/2=25ºC. Aproximamos a 30ºC
cp_e2=1006.4; %J/kg/ºC
rho_e2=1.1644; %kg/m3
nu_e2=18.65e-06; %kg/m/s
k_e2=26.01e-03; %W/m/ºC
Pr_e2=0.712;
%Recalculo
Re_e2=rho_e2*v_e*lc_e/nu_e2; %(>2300: convección forzada. Re_e2>10^6:
aplicamos Churchill y Berstein
he2=k_e2*(0.3+0.62*Re_e2^0.5*Pr_e2^(1/3)*(1+(Re_e2/2.82e05)^(5/8))^(4/
5)/(1+(0.4/Pr_e2)^(2/3))^0.25)/lc_e;
%Recalculo del calor total y temperatura pared exterior
Q_l2=2*pi*(t_i-t_e)/(1/r1/hi+log(r2/r1)/k12+log(r3/r2)/k23+1/r3/he2);
%W/m
t32=t_e+Q_l2/(2*pi*r3*he2);

Solución
hi = 5.3595e+03 W/m2/ºC
he = 6.5739 W/m2/ºC
Q_l = 69.0519 W/m

PROBLEMA 4
Un conducto de aire acondicionado de 0,70 m de ancho y 0,30 m de alto y de 20 m
de longitud, se encuentra a una temperatura superficial de 20°C en un local a 35°C.
Calcular el calor transmitido por convección desde el aire al conducto. Repetir el
ejercicio suponiendo 35°C en la superficie del conducto, que en este caso sería de
calefacción, en un local a 20°C. Suponer que β no varía.

13
%Datos
w=0.7; %m
h=0.3; %m
l=20; %m
dh=2*w*h/(w+h); %m
tf1=35; %ºC
ts1=20;
tf2=20; %ºC
ts2=35;
%Propiedades a tm=27.5ºC
cp_1=1006.4; %J/kg/ºC
rho_1=1.1644; %kg/m3
nu_1=18.65e-06; %kg/m/s
k_1=26.01e-03; %W/m/ºC
Pr_1=0.712;
g=9.81;
beta=1/(273+tf1);
alpha_1=k_1/rho_1/cp_1;
%Correlación de Churchill y Chu
Ra_1=rho_1*g*beta*abs(ts1-tf1)*dh^3/nu_1/alpha_1;
h_1=k_1*(0.6+0.387*Ra_1^(1/6)/(1+(0.559/Pr_1)^(9/16))^(8/27))^2/dh;
q=pi*dh*l*h_1*(tf1-ts1);

Solución
Ra_1 = 9.9567e+07
h_1 = 3.5004
q = 1.3856e+03 W

PROBLEMA 5
Un conductor eléctrico de 0,80 cm de diámetro colocado horizontalmente, disipa 10
kcal/h por metro lineal en un recinto en el que la temperatura del aire es 30°C. Se
desea conocer cuál es la temperatura superficial del conductor.

%Datos
convert_kcalh_W=4186/3600;
d=0.8e-02; %m
Q_l=10*convert_kcalh_W;
tf=30; %ºC

14
%Estimación de la temperatura superficial. Suponemos que h=10 W/m2/ºC.
tsi=Q_l/pi/d/10+tf;
ts_1=70;
%Suponemos ts_i=70: tm=50
cp_1=1007.4; %J/kg/ºC
rho_1=1.0924; %kg/m3
nu_1=19.57e-06; %kg/m/s
k_1=27.81e-03; %W/m/ºC
Pr_1=0.709;
g=9.81;
beta=1/(273+tf);
alpha_1=k_1/rho_1/cp_1;
%Correlación de Churchill y Chu
Ra_1=rho_1*g*beta*(ts_1-tf)*d^3/nu_1/alpha_1; %(Ra_1 < 10^9: régimen
laminar)
h_1=k_1*(0.6+0.387*Ra_1^(1/6)/(1+(0.559/Pr_1)^(9/16))^(8/27))^2/d;
%Recalculo
ts_2=Q_l/pi/d/h_1+tf;
tm_2=(ts_2+tf)/2;
%Mismas propiedades porque aproximamos a 50ºC. El problema converge,
luego ts=ts_2

Solución
h_1 = 9.8414 W/m2/ºC
ts_2 = 77.0113 ºC
tm_2 = 53.5056 ºC

15
CAMBIADORES DE CALOR
PROBLEMA 1
El aceite de un motor se enfría por aire en un intercambiador de calor de flujo
cruzado de un solo paso, donde ambos fluidos están no mezclados. Entra aire
atmosférico a 20ºC y 0.53 kg/s. Aceite a 0.026 kg/s entra a 80ºC y fluye por un tubo
de 10/12 mm de diámetro interior y exterior. En total hay 60 tubos, y el área total de
transferencia de calor es 1 m2. El coeficiente global de convección es 53 W/m2/ºC. Se
pide:
1. Calcular las temperaturas de salida de ambos fluidos. Dar por buenas las
iteraciones para diferencias de temperatura entre supuesto y calculado
inferiores a 5ºC.
2. Teniendo en cuenta que el aceite es un fluido altamente viscoso, calcular el
coeficiente de transmisión de calor por convección del lado del aceite.
Suponer que la temperatura interior de la pared del tubo es 10ºC superior a
la temperatura media de masa del aceite.
NOTA: No interpolar propiedades. Aproximar a la temperatura más cercana.

Propiedades del aceite del motor

cp µ k
ρ
temperatura (calor (viscosidad (conductividad
(densidad)
específico) dinámica) térmica)
ºC kg/m 3 J/kg/ºC kg/m/s (Pa·s) W/m/ºC
0 899.1 1796 3.848 0.147
20 888.2 1880 0.799 0.145
40 876.1 1964 0.21 0.144
60 864 2047 0.0725 0.14
65 861 2068 0.062375 0.1395
70 858 2089 0.05225 0.139
75 855 2110 0.042125 0.1385
80 852 2131 0.032 0.138
100 840 2219 0.0171 0.137
120 829 2307 0.0103 0.135
140 816.9 2395 0.00654 0.133
160 805.9 2483 0.00451 0.132

Solución
%Aire
mf=0.53; %kg/s
tfe=20; %ºC

16
%Aceite
mc=0.026; %kg/s
tce=80; %ºC
di=10e-03; %m
de=12e-03; %m
U=53; %W/m2/ºC
A=1; %m2
n=60; %número de tubos

%Longitud de los tubos

l=A/pi/de/n; %m.

%Suponemos:
tcs=60; %ºC
tfs=40; %ºC

%Propiedades del aire a tfm=30ºC;

cpf=1007; %J/kg/ºC
rhof=1.149; %kg/m3
myf=18.69e-06; %Pa*s
kf=26.62e-03; %W/m/ºC

%Propiedades del aceite a tcm=70ºC

rhoc=858; %kg/m3
cpc=2089; %J/kg/ºC
myc=0.05225; %Pa*s
kc=0.139; %W/m/ºC
Prc=cpc*myc/kc;

%1. Determinar las temperaturas de salida del aceite y del aire

Cf=mf*cpf; %Cf=533.7100;
Cc=mc*cpc; %Cc=54.3140:
Cmin=Cc;
NTU=U*A/Cmin;
CR=Cc/Cf;
%Por las gráficas:
epsilon=0.48;
Q=epsilon*Cmin*(tce-tfe);

17
tfs=tfe+Q/Cf;
tcs=tce-Q/Cc;

%Rehacemos cálculos con tcs=51.2; tfs=23;

%Propiedades del aceite a tfm=21.5ºC. Se aproxima a 20ºC
rhoc_1=888.2; %kg/m3
cpc_1=1880; %J/kg/ºC
myc_1=0.799; %Pa*s
kc_1=0.145; %W/m/ºC
Prc_1=cpc_1*myc_1/kc_1;

%Propidades del aire a tfm=65.6ºC. Aproximamos a 70ºC

cpf_1=1009; %J/kg/ºC
rhof_1=1.015; %kg/m3
myf_1=20.56e-06; %Pa*s
kf_1=29.52e-03; %W/m/ºC

%Recálculo de las propiedades

Cf_1=mf*cpf_1; %Cf=533.7100;
Cc_1=mc*cpc_1; %Cc=54.3140:
Cmin_1=Cc_1;
NTU_1=U*A/Cmin_1;
CR_1=Cc_1/Cf_1;

epsilon_1=0.5;
Q_1=epsilon_1*Cmin_1*(tce-tfe);

tfs_1=tfe+Q_1/Cf_1;
tcs_1=tce-Q_1/Cc_1;

%1.Cálculo del coeficiente de película interior

Ac=pi*di^2/4;
uc=mc/rhoc_1/Ac;
Rec=rhoc_1*uc*di/myc_1; %Rec=63.3574. Hay que utilizar correlación de
flujo laminar. Correlación de Sieder y Tate
%Viscosidad dinámica a la temperatura de la pared
mycs_1=32e-03;
hi=kc_1*1.86*((di/l)*Rec*Prc)^(1/3)*(myc_1/mycs_1)^0.14/di;

Solución
18
tfs_1 = 22.7421 ºC
tcs_1 = 50 ºC
hi = 177.3343 W/m2/ºC

PROBLEMA 2
Una caldera que se utiliza para generar vapor saturado tiene forma de un
intercambiador de calor en flujo cruzado, con agua que fluye por los tubos y un gas
a alta temperatura en flujo cruzado sobre los tubos. El gas, que tiene un calor
específico de 1140 J/kg/ºC y un flujo de masa de 12 kg/s, entra en el intercambiador
a 1200ºC. El agua, que tiene un flujo de 4 kg/s, entra como líquido saturado a 175ºC,
y sale como vapor saturado a la misma temperatura. Si el coeficiente global de
transferencia de calor es igual a 70 W/m2/ºC y hay 400 tubos, cada uno de diámetro
exterior igual a 0.03 m, ¿Cuál es la longitud del tubo que se requiere?.
Hipótesis y datos adicionales: Calor latente de evaporación del agua a 175ºC:
Δhlv=2.032·106 J/kg

%Datos
%Gas
cp_c=1140; %J/kg/ºC. Calor específico del gas
m_c=12; %kg/s. Caudal del gas
t_c_e=1200; %ºC. Temperatura de entrada caliente

C_c=m_c*cp_c;

%Agua
m_f=4; %kg/s
t_f_e=175; %ºC. Temperatura de entrada frío
t_f_s=175; %ºC
hlv=2.032e06; %J/kg

cp_f=4411; %J/kg/ºC
C_f=m_f*cp_f;

U=70; %W/m2/ºC
n=400; %Número de tubos
do=0.03; %m

Q=m_f*hlv;

19