Enunciado

113 – Uma chaminé passando gases quentes de exaustão tem uma área de seção transversal quadrada de 300 mm de lado. As paredes são construídas de tijolos refratários de 150 mm de espessura com uma condutividade térmica de 0,85 W/m.K e difusividade térmica α = 5,5x10-7m2/s. Inicialmente sem gases na chaminé escoando, as paredes estão a uma temperatura uniforme de 25°C. A superfície interior é rapidamente exposta aos gases quentes a 350°C com um coeficiente convecção de 100W/m2.K, enquanto a superfície exterior experimenta uma convecção com ar a 25°C e com coeficiente convectivo de 5W/m2.K. Usando um método de diferenças finitas implícito com um espaçamento de 50 mm e incremento de tempo de 1 h, encontre a distribuição de temperatura na parede e a taxa de perda de calor por convecção do exterior da chaminé em 5, 10, 50 e 100 h depois da introdução dos gases na chaminé.

Passo 1

Vamos começar determinando a expressão para a temperatura que representará a distribuição de temperatura.

Vamos fazer um balanço de energia na parede da chaminé.

O balanço de energia global é:

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

Onde,

- E ˙ e n t r a e E ˙ s a i , são, respectivamente, a taxa de transferência de energia que entra no sistema e a que sai do sistema;

- E ˙ g , é a taxa de geração de energia dentro do sistema;

- E ˙ s i s t e m a , é taxa de energia acumulada no sistema, que varia com o tempo em um processo transiente.

A energia acumulada no sistema é:

E ˙ s i s t e m a = ρ c p ∂ T ∂ t

Desprezando a geração de energia:

E ˙ g = 0

A energia que entra no sistema é dada por:

E ˙ e n t r a = q c o n v , i + q c o n d

E, a energia que sai é:

E ˙ s a i = q c o n v , e

Passo 2

Assim, a Eq. 1 fica:

q c o n v , i + q c o n d - q c o n v , e = ρ c p ∂ T ∂ t

Onde,

q c o n v , i = h i A T ∞ , i - T i

q c o n v , e = h e A T ∞ , e - T e

q c o n d = k A d T d x

Assim,

h i A T ∞ , i - T i + k A d T d x - h e A T ∞ , e - T e = ρ c p ∂ T ∂ t

Passo 3

Aplicando diferenças finitas pele método implícito:

∂ T ∂ t = T m p + 1 - T m p ∆ t

d T d x = T m p + 1 - T m - 1 p + 1 ∆ x

Assim, para o interior da parede:

h i k A T ∞ , i - T m - 1 p + T m p + 1 - T m - 1 p + 1 ∆ x - h e k A T ∞ , e - T N p = 1 α T m p + 1 - T m p ∆ t

Esse é o balanço de energia global na chaminé discretizado com diferenças finitas pelo método implícito.

Passo 4

Na superfície interna:

h i k A T ∞ , i - T i p + T 1 p + 1 - T 0 p + 1 ∆ x = 1 α T 1 p + 1 - T 1 p ∆ t

Onde,

T i p = T 0 p

Assim,

h i k A T ∞ , i - T 0 p + T 1 p + 1 - T 0 p + 1 ∆ x = 1 α T 1 p + 1 - T 1 p ∆ t

T 1 p + 1 1 ∆ x - 1 α ∆ t - 1 ∆ x T 0 p + 1 = h i k A T ∞ , i - T 0 p - 1 α ∆ t T 1 p

Onde, para p = 0:

T 1 0 = T 0 0 = T i 0 = 25 ℃ = 298 K

Assim, aplicando os valores:

T 1 p + 1 1 0,05 - 1 ( 5,5 x 10 - 7 ) ( 3600 ) - 1 0,05 T 0 p + 1 = 100 0,85 ( 0,15 2 ) 623 - T 0 p - 1 ( 5,5 x 10 - 7 ) ( 3600 ) T 1 p

485 T 1 p + 1 + 20 T 0 p + 1 = 505 T 1 p - 2,65 62 3 - T 0 p

Passo 5

Na superfície externa:

T N p + 1 - T N - 1 p + 1 ∆ x - h e k A T ∞ , e - T e p = 1 α T N p + 1 - T N p ∆ t

Onde,

T e p = T N - 1 p

Assim,

T N p + 1 - T N - 1 p + 1 ∆ x - h e k A T ∞ , e - T N - 1 p = 1 α T N p + 1 - T N p ∆ t

T N p + 1 1 ∆ x - 1 α ∆ t - 1 ∆ x T N - 1 p + 1 = h e k A T ∞ , e - T N - 1 p - 1 α ∆ t T N p

Onde, para p = 0:

T N 0 = T N - 1 0 = T e 0 = 25 ℃ = 298 K

O nodo N é o nodo de nº 4.

Assim, aplicando os valores:

T 4 p + 1 1 0,05 - 1 ( 5,5 x 10 - 7 ) ( 3600 ) - 1 0,05 T 3 p + 1 = 5 0,85 ( 0,15 2 ) 298 - T 3 p - 1 ( 5,5 x 10 - 7 ) ( 3600 ) T 4 p

Assim,

4 85 T 4 p + 1 + 20 T 3 p + 1 = 505 T 4 p - 0 , 13 298 - T 3 p

Passo 6

Agora para os nodos interiores, podemos usar a Eq. 5.97 do Incropera que é aplicada para nodos interiores com o método de diferenças finitas implícito. Assim, para os nodos de 1 até 3:

1 + 2 F O T m p + 1 - F O T m - 1 p + 1 + T m + 1 p + 1 = T m p

Onde, FO é o nº de Fourier:

F O = α ∆ t ∆ x 2

F O = ( 5,5 x 10 - 7 ) ( 3600 ) 0,05 = 0,04

Substituindo, vamos chegar que:

- 0,04 T m - 1 p + 1 + 1,08 T m p + 1 - 0,04 T m + 1 p + 1 = T m p

Essa é eq. para cada nodo m interior de 1 até 3,

Passo 7

Vamos agora montar o sistema linear que compõem as equações para cada nodo:

485 T 1 p + 1 + 20 T 0 p + 1 = 505 T 1 p - 2,65 62 3 - T 0 p

- 0,04 T 0 p + 1 + 1,08 T 1 p + 1 - 0,04 T 2 p + 1 = T 1 p

- 0,04 T 1 p + 1 + 1,08 T 2 p + 1 - 0,04 T 3 p + 1 = T 2 p

- 0,04 T 2 p + 1 + 1,08 T 3 p + 1 - 0,04 T 4 p + 1 = T 3 p

4 85 T 4 p + 1 + 20 T 3 p + 1 = 505 T 4 p - 0 , 13 298 - T 3 p

Onde, no passo de tempo de p = 0, isto é, no tempo inicial, todas as temperaturas são iguais a 25ºC:

T 0 0 = T 1 0 = T 2 0 = … = T 4 0 = 25 ℃ = 298 K

Passo 8

Resolvendo o sistema linear pelo EXCEL, obtemos a seguinte distribuição de temperatura na parede:

Para t = 5 h:

Para t = 10 h:

Resposta

Para t = 5 h:

Para t = 10 h:

Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados

Lembre-se da placa retangular bidimensional do Problema 4.2. Usando um método numérico apropriado, com Δx = Δy = 0,25 m, determine a temperatura no seu ponto central (1; 0,5).
Uma aleta plana com seção transversal uniforme é feita com um material de condutividade térmica igual a 50 W m K , tem espessura w = 6 mm e comprimento de L = 48 mm, e é muito grande na direção normal
Seja o sistema para resfriamento de turbina a gás mostrado no Exemplo 4.3. No Problema 3.23 são descritas vantagens associadas à aplicação de um revestimento de barreira térmica (RBT) sobre a superfíc
Uma placa de grande espessura, com difusividade térmica de 5,6×10 - 6m 2 / s e condutividade térmica de 20W / ( m·K ) , está inicialmente a uma temperatura uniforme de 325 ° C . De repente, a sua supe
Pilares de concreto reforçados com aço são usados em grandes construções. Falhas estruturais podem ocorrer em altas temperaturas em função de um incêndio por causa do enfraquecimento do núcleo metálic