Enunciado

99 – Um tubo fino de cobre de 20 mm de diâmetro, 1 m de comprimento, contém um material que gera energia térmica a uma taxa volumétrica de q ˙ = 0,2   M W / m 3 . A haste é resfriada por ar na pressão atmosférica em um escoamento cruzado com velocidade de 2 m/s. É proposto uma redução na temperatura do material por pulverização através de um revestimento fino e poroso de espessura t=2mm colocado no exterior do tubo de cobre. O revestimento é então saturado com água liquida, levando a um resfriamento evaporativo na superfície exposta do revestimento. A condutividade térmica do revestimento de água saturada é kC=3,5 W/m.K. Determine a temperatura do tubo de cobre com e sem revestimento. As condições ambientes são caracterizadas por uma temperatura de 25°C e uma umidade relativa de 45%.

Passo 1

Neste exercício temos duas situações. A primeira é um processo de troca térmica com um tubo de cobre trocando calor com ar atmosférico. Na segunda situação temos o tubo de cobre trocando calor com um revestimento saturado com água. Assim, a ideia é aplicar um balanço de energia nestas duas situações, com a primeira tendo o sistema como o tubo de cobre, e na segunda, tendo sistema composto pelo tubo de cobre e o revestimento.

Passo 2

Vamos aplicar um balanço de energia no tubo de cobre, sem o revestimento:

Lembrando que a expressão geral para um balanço de energia é a seguinte

E ˙ e n t r a - E ˙ s a i + E ˙ g = d E s i s t e m a d t         E q 1

Onde,

- E ˙ e n t r a , é a taxa de energia que entra no sistema, no caso o sistema é a lagoa;

- E ˙ s a i , é a taxa de energia que sai do sistema;

- E ˙ g , é a taxa de energia gerada no sistema;

- E s i s t e m a , é a quantidade de energia acumulada no sistema, que varia com o tempo.

Considerando uma condição de estado estacionário (sem variação com o tempo), assim

d E s i s t e m a d t = 0

Não há energia que entra no sistema, composto pelo tubo de cobre

E ˙ e n t r a = 0

A energia que sai é referente ao calor transferido por convecção, qconv, do tubo de cobre com o ar atmosférico

E ˙ s a i = q c o n v

E a energia gerada é referente a taxa volumétrica de energia térmica gerada dentro do tubo de cobre

E ˙ g = q ˙ = 0,2   M W / m 3

Assim, o balanço de energia no tubo de cobre, sem revestimento é:

- q c o n v + V q ˙ = 0             E q 2

Onde, V é o volume do cilindro que é dado por

V = π r 2 L = π 4 D 2 L

Passo 3

Agora um balanço de energia em um sistema composto pelo tubo de cobre e o revestimento:

Analogamente, pegando a Eq1

E ˙ e n t r a - E ˙ s a i + E ˙ g = d E s i s t e m a d t

Ainda em condição de estado estacionário (sem variação com o tempo)

  d E s i s t e m a d t = 0

Não há energia que entra no sistema

E ˙ e n t r a = 0

A energia que sai é referente ao calor transferido por condução, qcond, através do revestimento, também uma energia que é transferida para a água sobre o revestimento, que fará esta evaporar, qevap, e também um calor perdido por convecção, qconv, por causa do ar atmosférico, assim

E ˙ s a i = q c o n d + q e v a p + q c o n v

E a energia gerada é referente a taxa volumétrica de energia térmica gerada dentro do tubo de cobre

E ˙ g = q ˙ = 0,2   M W / m 3

Assim, o balanço de energia no tubo de cobre com revestimento é:

V q ˙ - ( q c o n d + q e v a p + q c o n v ) = 0             E q 3

Passo 4

Vamos agora equacionar os termos da Eq2 e Eq3

Vamos começar com o primeiro caso, da Eq2:

Para determinarmos qconv podemos usar a seguinte expressão

q c o n v = h A S T - T ∞

Onde,

-h é o coeficiente de transferência de calor;

-AS é a área superficial de troca térmica;

-T é a temperatura do tubo de cobre, que queremos encontrar;

- T ∞ , é a temperatura do ar.

A área de troca térmica é a área superficial do cilindro que é

A S = π D L

Onde,

-D é o diâmetro do cilindro;

-L é o comprimento do cilindro.

O coeficiente convectivo, h, nós podemos determinar usando a relação com o nº de Nusselt, Nu

N u = h L k

Onde, k á condutividade térmica do ar.

Isolando h

h = k N u L

Substituindo essas relações

q c o n v = k N u L π D L T - T ∞

q c o n v = π D k N u T - T ∞                       E q 4

Passo 5

Para determinarmos o valor de Nu, você pode lembrar que existem várias correlações empíricas o seu cálculo. Mas, você deve atentar que não é qualquer correlação empírica que é aplicada para este caso que se trata de fluxo cruzado em um cilindro. Assim temos que pegar uma correlação apropriada para este caso.

Uma correlação aplicada para fluxo cruzado em cilindros é a correlação de Hilpert que é mais prática e simples do que outras, então vamos usa-la. A correlação de Hilpert é a seguinte:

N u = h D k = C R e m P r 1 3

Onde,

-C e m são constantes encontradas em tabelas, e dependem do valor de Reynolds;

-Re e Pr são o número de Reynolds e número de Prandtl, respectivamente.

Passo 6

Vamos achar Reynolds (Re) e Prandtl (Pr):

A equação de Reynolds é a seguinte:

R e = V D v

Onde,

-D é o diâmetro do cilindro;

-V é a velocidade do fluido;

-ν é a viscosidade cinemática do fluido.

Vamos encontrar a viscosidade cinemática pela Tabela A.4

A T=25ºC=298K:

v = 15,7 x 10 - 6 m 2 s

Aplicando os valores

R e = 2 m s ( 0,02   m ) 15,7 x 10 - 6 m 2 s = 2548

Pela mesma Tabela A.4 achamos o número de Prandtl, Pr, e também o valor de k:

k = 2 6 , 1 x 10 - 3 W m K ; P r = 0,7 10

Passo 7

As constantes m e C, podem ser encontradas na Tabela 7.2 do Incropera, cap. 7, pegando o valor de Re encontrado.

Para Re=2548, os valores de C e m que estão na faixa de Re para esse valor é:

C = 0,683 ;   m = 0,466

Agora calculando Nusselt

N u = 0,683 ( 2548 ) 0,466 ( 0,710 ) 1 3

N u = 23,6

Passo 8

Pegando a Eq4 e a Eq2

- q c o n v + V q ˙ = 0

q c o n v = π D k N u T - T ∞

Chegamos que

V q ˙ = q c o n v = π D k N u T - T ∞

V q ˙ = π D k N u T - T ∞

Isolando T

T - T ∞ = V q ˙ π D k N u

T = T ∞ + V q ˙ π D k N u

Onde, para o volume V

V = π 4 D 2 L

Assim,

T = T ∞ + π 4 D 2 L q ˙ π D k N u = T ∞ + D L q ˙ 4 k N u

Aplicando os valores

T = 298   K + 0,02   m 1   m 0,2 x 10 6 W m 3 4 26,1 x 10 - 3 W m K 23,6

T = 1924   K

Passo 9

Vamos agora para o segundo caso:

A expressão que chegamos pelo balanço de energia foi

V q ˙ - ( q c o n d + q e v a p + q c o n v ) = 0

Para qconv fazendo de forma análoga ao primeiro caso chegaremos na expressão

q c o n v = π D R e v k N u T - T ∞

Onde, DRev é diâmetro do tubo com revestimento que é

D R e v = D + t = 20   m m + 2   m m = 22   m m = 0,022   m

O novo valor de Nusselt é

N u = C R e m ( P r ) 1 3

Precisamos do novo valor de Reynolds

R e = V D R e v v

R e = 2 m s ( 0,022   m ) 15,7 x 10 - 6 m 2 s = 2803

Os valores de C e m são:

C = 0,683 ; m = 0,466

Calculando Nusselt

N u = 0,683 ( 2 8 03 ) 0,466 ( 0,710 ) 1 3 = 24,6

Passo 10

Para qcond, podemos usar a Lei de Fourier, em que chegaremos na seguinte expressão

q c o n d = k c π D R e v ( T - T ∞ )

E para qevap que é a energia da água perdida por evaporação está relacionada diretamente com a entalpia de vaporização da água, Hvap, e também com o fluxo de massa que a água libera na evaporação. Dessa forma, a expressão para a taxa de evaporação é:

q e v a p = n á g u a ' ' A S H v a p

Onde,

-n”agua é o fluxo mássico de água evaporada;

-Hvap é a entalpia de vaporização da água, um dado tabelado que depende da temperatura

Onde,

A S = π D R e v L

Assim,

q e v a p = n á g u a ' ' π D R e v L H v a p

Passo 11

O fluxo mássico é dado pela seguinte expressão

n á g u a ' ' = h m ρ á g u a , S a t - ρ ∞

Onde,

-hm é o coeficiente de transferência de massa;

-ρágua,Sat é a densidade do vapor de água saturada na temperatura do processo;

- ρ ∞ , é a densidade da água saturada no ar.

As densidades podem ser determinadas em tabelas da água saturada pegando a temperatura de 25ºC, na qual para a densidade da água no ar temos que ela depende da umidade relativa da seguinte forma:

ρ ∞ = Φ ρ á g u a , S a t

Onde,

-φ é a umidade relativa.

Assim,

n á g u a ' ' = h m ρ á g u a , S a t - Φ ρ á g u a , S a t

n á g u a ' ' = h m ρ á g u a , S a t 1 - Φ

Passo 12

Agora vamos determinar uma forma de encontrar o coeficiente de transferência de massa. Os processos de transferência de massa e transferência de calor são muito similares, assim, você pode lembrar que podemos usar as mesmas correlações para o cálculo dos parâmetros da transferência de massa. Analogamente, a equação usada para o coeficiente de transferência de calor foi

N u = C R e m P r 1 3

Para o coeficiente de transferência de massa, basta trocarmos o Nusselt (Nu) por Sherwood (Sh) e o Prandtl (Pr) por Schmidt (Sc), assim

S h = C R e m S c 1 3

Lembrando que

S h = h m D R e v D A B ; S c = v D A B

Onde, DAB é o coeficiente difusivo de transferência de massa. Assim,

h m = S h D A B L

h m = D A B L C R e m S c 1 3

h m = D A B L C R e m v D A B 1 3

Passo 13

Pegando a expressão de hm e substituindo na equação do fluxo

n á g u a ' ' = h m ρ á g u a , S a t 1 - Φ

n á g u a ' ' = D A B L C R e m v D A B 1 3 ρ á g u a , S a t 1 - Φ

Agora vamos pegar essa e substituir na eq. de qevap

q e v a p = n á g u a ' ' π D R e v L H v a p

q e v a p = D A B L C R e m v D A B 1 3 ρ á g u a , S a t 1 - Φ π D R e v L H v a p

q e v a p = π D R e v H v a p D A B C R e m v D A B 1 3 ρ á g u a , S a t 1 - Φ

Passo 14

Pegando os dados para a água saturada, a T=25ºC=298 K, Tabela A.6:

  ρ á g u a , S a t = 1 v á g u a , S a t = 1 44,3 m 3 k g = 0,0226 k g m 3 ;   H v a p = 2442 x 10 3 J k g

-Dado para o coeficiente difusivo DAB, Tabela A.9, temperatura 20ºC=293K:

D A B = 0,26 x 10 - 4 m 2 s

Aplicando os valores

q e v a p = π 0,022   m 2442 x 10 3 J k g 0,26 x 10 - 4 m 2 s 0,683 2802 0,466 15,7 x 10 - 6 m 2 s 0,26 x 10 - 4 m 2 s 1 3 0,0226 k g m 3 1 - 0,45

q e v a p = 1,2 7   W

Passo 15

Pegando a expressão do balanço de energia

V q ˙ - ( q c o n d + q e v a p + q c o n v ) = 0

V q ˙ = ( q c o n d + q e v a p + q c o n v )

Substituindo as relações de qcond e qconv

V q ˙ = π D R e v k N u T - T ∞ + k c π D R e v T - T ∞ + q c o n v

V q ˙ = π D R e v T - T ∞ k N u + k c + q e v a p

Sendo V o volume do cilindro:

V = π 4 D 2 L

Isolando T

T - T ∞ = π 4 D 2 L q ˙ - q e v a p π D R e v k N u + k c

T = T ∞ + π 4 D 2 L q ˙ - q e v a p π D R e v k N u + k c

Aplicando os valores

T = 298   K + π 4 0,02   m 2 1,0   m 0,2 x 10 6 W m 3 - 1,27   W π ( 0,022   m ) 26,1 x 10 - 3 W m K 24,6 + 3,5 W m K

T = 513   K

Resposta

Sem Revestimento:

T = 1 924   K

Com Revestimento:

T = 513   K

Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados

O condensador de uma planta de potência a vapor é constituído por uma matriz quadrada alinhada com 625 tubos, cada um com diâmetro de 25 mm. Considere condições nas quais o vapor de água saturado a 0,
Determine a taxa de condensação total e a taxa de transferência de calor para o processo de condensação do Problema 10.46, estando o tubo inclinado com ângulos θ = 0, 30, 45 e 60° a partir da horizont
Vapor saturado oriundo de um processo químico condensa a uma taxa baixa sobre a superfície interna de um recipiente cilíndrico vertical com paredes finas, comprimento L e diâmetro D. A parede do recip
O condensador de uma planta de potência a vapor é construído com tubos de aço inoxidável AISI 302 ( k a = 15 W/(m · K)), com diâmetros externo e interno de D e = 30 mm e D i = 26 mm, respectivamente.
Determine a taxa de condensação sobre uma esfera com 100 mm de diâmetro e com uma temperatura superficial de 150°C, que se encontra imersa em vapor saturado de etilenoglicol a uma pressão de 1 atm. Ap