17/05/2023, 20:05 Lista 3 Transcal 2_Jonas_Costa - Jupyter Notebook

Nome: Jonas Costa dos Santos

In [1]:

#Bib para auxiliar na resolução do problema

import math
import CoolProp
from CoolProp.CoolProp import PropsSI
import numpy as np
import sympy
from sympy import *
x = sympy.Symbol("x")
import matplotlib.pyplot as plt

1)Óleo de motor escoa através de um tubo de 25 mm de diâmetro a uma vazão de 0,5 kg/s. O óleo entra no
tubo a uma temperatura de 25 °C, enquanto a temperatura da superfície do tubo é mantida a 100 °C. (0,5
pontos)

(a) Determine a temperatura do óleo na saída para tubos com comprimentos de 5 m e 100 m. Para cada
caso, compare a média logarítmica das diferenças de temperaturas com a média aritmética das diferenças
de temperaturas. (0,25 pontos)

(b) Para 5 ≤ L ≤ 100 m, calcule e represente gracamente o número de Nusselt médio NuD e a temperatura
de saída do óleo como uma função de L. (0,25 pontos)

localhost:8888/notebooks/Lista 3 Transcal 2_Jonas_Costa.ipynb 1/13

17/05/2023, 20:05 Lista 3 Transcal 2_Jonas_Costa - Jupyter Notebook

resolução

a)
In [2]:

#Dados do problema
D1=25*10**-3 #m
m1=0.5 #kg/s
Tme1=25 #c
Ts1=100 #c
L1=5
L_1=100
#Dados tabelados
cp1=2035
mi1=0.0836
pr1=1205
k1=0.141

localhost:8888/notebooks/Lista 3 Transcal 2_Jonas_Costa.ipynb 2/13

17/05/2023, 20:05 Lista 3 Transcal 2_Jonas_Costa - Jupyter Notebook

In [3]:

#Calculo de Reynolds
Rey1=(4*m1)/(math.pi*D1*mi1)

#Para l=5

#Calculo de nulselt
Gzd=(D1/L1)*Rey1*pr1
Nu1=3.661+((0.0668*Gzd)/(1+(0.04*Gzd**(2/3))))

#Calculo de h
h1=(Nu1*k1)/D1

#Calculo Tmsai

Tmsai1= Ts1 - (Ts1 - Tme1) *(math.exp(-(math.pi*D1*L1*h1)/(m1*cp1)))

#Para l=100

Gzd1=(D1/L_1)*Rey1*pr1
Nu_1=3.661+((0.0668*Gzd1)/(1+(0.04*Gzd1**(2/3))))

#Calculo de h
h_1=(Nu_1*k1)/D1

#Calculo Tmsai
Tmsai_1=Ts1-(Ts1-Tme1)*math.exp(-(math.pi*D1*L_1*h_1)/(m1*cp1))

#resltados dos calculos

print('Reynolds= ' "{:e}".format(Rey1))
print('Nusselt= ' "{0:.3f}".format(Nu1),"para L=5m")
print('Nusselt= ' "{0:.3f}".format(Nu_1),"para L=100m")
print('h= ' "{0:.3f}".format(h1),"para L=5m")
print('h= ' "{0:.3f}".format(h_1),"para L=100m")
print("Respostas alternativa a),")
print('A temperatura do óleo na saída tubo com comprimentos de 5m sera de ' "{0:.1f}".fo
print('A temperatura do óleo na saída tubo com comprimentos de 100m sera de ' "{0:.1f}".

Reynolds= 3.046028e+02
Nusselt= 21.185 para L=5m
Nusselt= 7.041 para L=100m
h= 119.481 para L=5m
h= 39.709 para L=100m
Respostas alternativa a),
A temperatura do óleo na saída tubo com comprimentos de 5m sera de 28.4 °
C.
A temperatura do óleo na saída tubo com comprimentos de 100m sera de 44.8
°C.

localhost:8888/notebooks/Lista 3 Transcal 2_Jonas_Costa.ipynb 3/13

17/05/2023, 20:05 Lista 3 Transcal 2_Jonas_Costa - Jupyter Notebook

b)
In [4]:

#gerando graficos
L = np.linspace(5, 100)
Nu = (3.66+((0.0668*Rey1*pr1*(D1/L))/(1+0.04*((Rey1*pr1*(D1/L))**(2/3)))))
hm = (k1/D1)*(3.66+((0.0668*Rey1*pr1*(D1/L))/(1+0.04*((Rey1*pr1*
(D1/L))**(2/3)))))
T = Ts1-(Ts1-Tme1)*np.exp((-math.pi*D1*L*hm)/(m1*cp1))

#grafico do numero de nul

plt.plot(L,Nu, color = "blue" )
plt.grid(axis='y', color='0.60')
plt.grid(axis='x', color='0.60')
plt.ylabel("Nusselt médio")
plt.xlabel("Comprimento (m)")
plt.title("Gráfico Nusselt médio com relação ao comprimento")

Out[4]:

Text(0.5, 1.0, 'Gráfico Nusselt médio com relação ao comprimento')

localhost:8888/notebooks/Lista 3 Transcal 2_Jonas_Costa.ipynb 4/13

17/05/2023, 20:05 Lista 3 Transcal 2_Jonas_Costa - Jupyter Notebook

In [5]:

#Grafico da Temperatura
plt.plot(L,T, color = "blue" )
plt.grid(axis='y', color='0.60')
plt.grid(axis='x', color='0.60')
plt.ylabel("Temperatura em °C")
plt.xlabel("Comprimento (m)")
plt.title("Gráfico Temperatura que sai com relação ao comprimento")

Out[5]:

Text(0.5, 1.0, 'Gráfico Temperatura que sai com relação ao comprimento')

2. Óleo de motor é aquecido através do escoamento em um tubo circular de diâmetro D = 50 mm e

comprimento L = 25 m e cuja superfície é mantida a 150 °C. (0,5 pontos)

(a) Se a vazão e a temperatura de entrada do óleo são 0,5 kg/s a 20 °C, qual é a temperatura de saída
Tm,s? Qual é a taxa total de transferência de calor q para o tubo? (0,25 pontos)

(b) Para vazões na faixa de 0, 5 ≤ m˙ ≤ 2, 0 kg/s, calcule e represente graficamente as variações de Tm,s e
q com m˙ . Para qual (is) vazão (ões) de q e Tm,s, são maximizadas? (0,25 pontos)

localhost:8888/notebooks/Lista 3 Transcal 2_Jonas_Costa.ipynb 5/13

17/05/2023, 20:05 Lista 3 Transcal 2_Jonas_Costa - Jupyter Notebook

In [15]:

# Dados do problema
m2= 0.5 #kg/s
Ti2= 20 #c
Ts2= 150 #c
L2= 25 #m
D2= 0.05#m
# Valor para o óleo 353K
rho = 852
mi2= 37.5*10**-6
cp2 = 2131
v2= 0.032
k2 = 0.138
pr2 = 490

a)
In [16]:

#Calculo de Reynolds
Rey2=(4*m2)/(math.pi*D2*v2)

#Calculo de nulselt
Gzd2=(D2/L2)*Rey2*pr2

Nu2=3.66+((0.0668*Gzd2)/(1+(0.04*Gzd2**(2/3))))

#Calculo de h
h2=(Nu2*k2)/D2

Tmsai2 = Ts2-(Ts2-Ti2)*np.exp((-np.pi*D2*L2*h2)/(m2*cp2))

q2= m2*cp2*(Tmsai2 - Ti2)

#resltados dos calculos

print('Reynolds= ' "{:e}".format(Rey2))
print('Nusselt= ' "{0:.3f}".format(Nu2))
print('h= ' "{0:.3f}".format(h2))
print("Respostas alternativa a)")
print('A temperatura de saída sera de ' "{0:.1f}".format(Tmsai2),"°C.")
print('A taxa total de transferência de calor ' "{0:.1f}".format(q2),"°W.")

Reynolds= 3.978874e+02
Nusselt= 11.969
h= 33.034
Respostas alternativa a)
A temperatura de saída sera de 34.9 °C.
A taxa total de transferência de calor 15877.8 °W.

localhost:8888/notebooks/Lista 3 Transcal 2_Jonas_Costa.ipynb 6/13

17/05/2023, 20:05 Lista 3 Transcal 2_Jonas_Costa - Jupyter Notebook

b)
In [24]:

#gerando graficos
m_2 = np.linspace(0.5, 2)
R_2 = (4*m_2)/(np.pi*D2*v2)
h_2 = (k2/D2)*(3.66+((0.0668*R_2*pr2*(D2/L2))/(1+0.04*((R_2*pr2*(D2/L2))**(2/3)))))
TM = Ts2-(Ts2-Ti2)*np.exp((-np.pi*D2*L2*h_2)/(m_2*cp2))
q_2 = m_2*cp2*(TM - Ti2)

#grafico do numero de nul

plt.plot(m_2,TM, color = "blue" )
plt.grid(axis='y', color='0.60')
plt.grid(axis='x', color='0.60')
plt.ylabel("Nusselt médio")
plt.xlabel("Comprimento (m)")
plt.title("Gráfico Nusselt médio com relação ao comprimento")

Out[24]:

Text(0.5, 1.0, 'Gráfico Nusselt médio com relação ao comprimento')

localhost:8888/notebooks/Lista 3 Transcal 2_Jonas_Costa.ipynb 7/13

17/05/2023, 20:05 Lista 3 Transcal 2_Jonas_Costa - Jupyter Notebook

In [27]:

#grafico do numero de nul

plt.plot(m_2,q_2, color = "blue" )
plt.grid(axis='y', color='0.60')
plt.grid(axis='x', color='0.60')
plt.ylabel("Transferenia de calor em (W)")
plt.xlabel("Fluxo de massa em (kg/s)")
plt.title("Gráfico da tranferencia de calor pela vazão massica")

Out[27]:

Text(0.5, 1.0, 'Gráfico da tranferencia de calor pela vazão massica')

3. O canal de ar para o resfriamento da pá de uma turbina a gás pode ser aproximado por um tubo de 3
mm de diâmetro e 75 mm de comprimento. A temperatura de operação da pá é de 650 °C e o ar entra
no tubo a 427 °C.(0,5 pontos)

(a) Para uma vazão de ar de 0,18 kg/h, calcule a temperatura de saída do ar e o calor removido da pá.(0,25
pontos)

(b) Gere um gráfico da temperatura de saída do ar em função da vazão no intervalo 0, 1 ≤ m˙ ≤ 0, 6 kg/h.

Compare esse resultado com aqueles obtidos para pás com canais com diâmetro de 2 e 4 mm,
considerando que as demais condições permaneçam inalteradas.(0,25 pontos)

localhost:8888/notebooks/Lista 3 Transcal 2_Jonas_Costa.ipynb 8/13

17/05/2023, 20:05 Lista 3 Transcal 2_Jonas_Costa - Jupyter Notebook

In [62]:

# Dados do problema
D3 = 0.003
L3 = 0.075
To3 = 650
Te3 = 427
VM3 = 0.18/3600
Patm = 101325
Tmed3 = ((To3+Te3)/2)+ 273.15
m3 = (0.18/3600)

#Dados coolpe

mi3=PropsSI('V','P',Patm,'T',Tmed3,'air')
rho3=PropsSI('D','P',Patm,'T',Tmed3,'air')
Cp3 = PropsSI('C','P',Patm,'T',Tmed3,'air')
Pr3 = PropsSI('Prandtl','P',Patm,'T',Tmed3,'air')
vc3 = mi3/rho3
k3=(vc3*rho3*Cp3)/Pr3
print('Densidade = ' "{0:.3f}".format(rho3),"kg/m3")
print('Calor Especifico a pressão constante(Cp) = ' "{0:.3f}".format(Cp3/1000),"J/Kg")
print('Visosidade dinâmica (u) = ' "{0:.7f}".format(mi3),"kg/m s")
print('Visosidade cinemática (V) = ' "{0:.8f}".format(vc3),"m2/ s")
print('Número de Prandtl (Pr) = ' "{0:.8f}".format(Pr3))
print('k = ',k3)

Densidade = 0.435 kg/m3

Calor Especifico a pressão constante(Cp) = 1.101 J/Kg
Visosidade dinâmica (u) = 0.0000377 kg/m s
Visosidade cinemática (V) = 0.00008679 m2/ s
Número de Prandtl (Pr) = 0.71801629
k = 0.05787501255095444

localhost:8888/notebooks/Lista 3 Transcal 2_Jonas_Costa.ipynb 9/13

17/05/2023, 20:05 Lista 3 Transcal 2_Jonas_Costa - Jupyter Notebook

a)
In [63]:

#calculo do reynolds
Rey3=(4*VM3)/(np.pi*(D3)*mi3)

#Calculo de nulselt
Gzd3=(D3/L3)*Rey3*Pr3

Nu3=3.66+((0.0668*Gzd2)/(1+(0.04*Gzd2**(2/3))))

#Calculo de h
a = (3.66/(np.tanh(2.264*(Gzd3**(-1/3))+1.7*(Gzd3**(-2/3))))+0.0499*Gzd3*np.tanh(Gzd3**-1
b = np.tanh(2.432*Pr3**(1/6)*Gzd3**(-1/6))
h3 = (k3/D3)*(a/b)

Tmsai3 = (((To3-Te3)*np.exp(-(np.pi*D3*L3*h3)/(VM3*Cp3)))-650)*-1

#calculo do q
Q3 = m3*Cp3*(Tmsai3-Te3)

#resltados dos calculos

print('Reynolds= ' "{:e}".format(Rey3))
print('Nusselt= ' "{0:.3f}".format(Nu3))
print('h= ' "{0:.3f}".format(h3))
print('A temperatura de saída sera de ' "{0:.1f}".format(Tmsai3),"°C.")
print('O calor removido da pá sera de ' "{0:.1f}".format(Q3),"W.")

Reynolds= 5.624297e+02
Nusselt= 11.969
h= 97.077
A temperatura de saída sera de 585.9 °C.
O calor removido da pá sera de 8.7 W.

localhost:8888/notebooks/Lista 3 Transcal 2_Jonas_Costa.ipynb 10/13

17/05/2023, 20:05 Lista 3 Transcal 2_Jonas_Costa - Jupyter Notebook

b)
In [64]:

#Geração Grafico
vm3 = np.linspace(0.1, 0.6)
a3 = (3.66/(np.tanh(2.264*(Gzd3**(-1/3))+1.7*(Gzd3**(-2/3))))+0.0499*Gzd3*np.tanh(Gzd3**-
b3 = np.tanh(2.432*Pr3**(1/6)*Gzd3**(-1/6))
h3 = (k3/D3)*(a3/b3)
TsG3 = (((To3-Te3)*np.exp(-(np.pi*D3*C3*h3)/(vm3*Cp3)))-650)*-1

plt.plot(vm3,TsG3, color = "red")

plt.grid(axis='y', color='0.90')
plt.grid(axis='x', color='0.90')
plt.ylabel("Temperatura na saída (°C)")
plt.xlabel("Vazão mássica (Kg/s)")
plt.legend("")
plt.title("Temperatura na saída pela Vazão massica")

Out[64]:

Text(0.5, 1.0, 'Temperatura na saída pela Vazão massica')

4 Ar aquecido necessário para um processo de secagem de alimentos é gerado pela passagem de ar

ambiente a 20 °C através de longos tubos circulares (D = 50 mm, L = 5 m) posicionados no interior de um
condensador de vapor. Vapor de água saturado a pressão atmosférica condensa sobre a superfície externa
dos tubos, mantendo uma temperatura superficial uniforme de 100 °C. (0,5 pontos)

(a) Se uma vazão de ar de 0,01 kg/s for mantida em cada tubo, determine a temperatura de saída do ar
Tm,sai, e a taxa de transferência de calor total q no tubo. (0,25 pontos)

localhost:8888/notebooks/Lista 3 Transcal 2_Jonas_Costa.ipynb 11/13

17/05/2023, 20:05 Lista 3 Transcal 2_Jonas_Costa - Jupyter Notebook

(b) A temperatura de saída do ar pode ser controlada pelo ajuste da vazão mássica no tubo. Calcule e
represente graficamente Tm,sai, como uma função de m˙ para 0, 005 ≤ m˙ ≤ 0, 050 kg/s. Se um processo
específico de secagem necessita de aproximadamente 1 kg/s de ar a 75 °C, que 2 projeto e condições de
operação devem ser estabelecidos para o aquecedor de ar, sujeito à limitação de que o diâmetro e o
i t d t b i 50 5 ti t ? (0 25 t )

a)
In [78]:

#dados do problema
m4=0.01
D4=0.05
L4=5
Ts5=100
Te=20
#dados tabelados
mi5=200.17*10**-7
pr5=0.702
k5=0.0287
cp4=1008.32

In [81]:

#Calculo de Reynolds
Rey4=(4*m4)/(math.pi*D4*mi5)

#caulcuo de nul

Nu4=0.023*Rey4**(4/5)*pr5**0.4

h_4=(Nu4*k5)/D4

Tmsai4 = Ts5-(Ts5-Te)*np.exp((-np.pi*D4*L4*h_4)/(m4*cp4))

Q4 = m4*cp4*(Tmsai4-Te)

#resltados dos calculos

print('Reynolds= ' "{:e}".format(Rey4))
print('Nusselt= ' "{0:.3f}".format(Nu4))
print('h= ' "{0:.3f}".format(h_4))
print('A temperatura de saída sera de ' "{0:.1f}".format(Tmsai4),"°C.")
print('O calor removido da pá sera de ' "{0:.1f}".format(Q4),"W.")

Reynolds= 1.272158e+04
Nusselt= 38.362
h= 22.020
A temperatura de saída sera de 85.6 °C.
O calor removido da pá sera de 661.5 W.

localhost:8888/notebooks/Lista 3 Transcal 2_Jonas_Costa.ipynb 12/13

17/05/2023, 20:05 Lista 3 Transcal 2_Jonas_Costa - Jupyter Notebook

In [83]:

# Equaçoes para graficar:

Vm4 = np.linspace(0.005, 0.05)
Ra4 = (4*Vm4)/(np.pi*D4*mi5)
Nud4 = 0.023*(Ra4**(4/5))*pr5**0.4
h4 = Nud4*(k5/D4)
Tm_sai4 = (-np.exp((-np.pi*D4*L4*h4)/(Vm4*cp4))*(Ts5-Te))+Ts5
#Gráfico
plt.plot(Vm4,Tm_sai4, color = "blue")
plt.grid(axis='y', color='0.90')
plt.grid(axis='x', color='0.90')
plt.ylabel("Temperatura de saída (°C)")
plt.xlabel("Vazão mássica (Kg/s)")
plt.legend("")
plt.title("Gráfico da temperatura em função da vazão mássica")

Out[83]:

Text(0.5, 1.0, 'Gráfico da temperatura em função da vazão mássica')

localhost:8888/notebooks/Lista 3 Transcal 2_Jonas_Costa.ipynb 13/13