Enunciado

Seja a configuração nodal 3 da tabela 4.2 . Deduza as equações de diferenças finitas, para condições de regime estacionário, nas seguintes situações.

(a) O contorno está isolado. Explique como a equação 4.42 pode ser modificada para concordar com o seu resultado.

(b) O contorno está sujeito a um fluxo de calor constante.

Passo 1

A configuração nodal 3 é

Desenhando as taxas de transferência de calor no volume de controle:

Antes de deduzir a equação de diferenças finitas vamos ressaltar que os contornos do vértice estão isolados, então q 3 = 0 .

Para deduzir vamos começar com o balanço de energia para regime estacionário:

d E d t = E ˙ e n t - E ˙ s a í d a

0 = E ˙ e n t

0 = q 1 + q 2 + q 4

Passo 2

Vamos agora calcular separadamente cada uma das taxas de transferência de calor. Todas as taxas de transferência são condutivas. Então, começando com q 1 :

q 1 = q ( m - 1 , n ) → ( m , n )

q 1 = k A L ∆ T

q 1 = k ∆ y ∙ 1 ∆ x T m - 1 , n - T m , n

Analogamente, para q 2 :

q 2 = q ( m , n - 1 ) → ( m , n )

q 2 = k A L ∆ T

q 2 = k ∆ x 2 ∙ 1 ∆ y T m , n - 1 - T m , n

A relação para q 4 é

q 4 = q ( m , n + 1 ) → ( m , n )

q 4 = k A L ∆ T

q 4 = k ∆ x 2 ∙ 1 ∆ y T m , n + 1 - T m , n

Então

k ∆ y ∙ 1 ∆ x T m - 1 , n - T m , n + k ∆ x 2 ∙ 1 ∆ y T m , n - 1 - T m , n + k ∆ x 2 ∙ 1 ∆ y T m , n + 1 - T m , n = 0

Para ∆ x = ∆ y :

k T m - 1 , n - T m , n + k 2 T m , n - 1 - T m , n + k 2 T m , n + 1 - T m , n = 0

2 T m - 1 , n - T m , n + T m , n - 1 - T m , n + T m , n + 1 - T m , n = 0

2 T m - 1 , n + T m , n - 1 + T m , n + 1 - 4 T m , n = 0

Comparando com a equação 4.42 podemos ver que chegaríamos na mesma relação se considerarmos que não há os termos de convecção ( h = 0 ) .

Passo 3

Vamos agora para o item (b). Neste item, devemos considerar que há um fluxo de calor constante no contorno do vértice. Então q 3 ≠ 0 .

Voltando para o balanço de energia:

q 1 + q 2 + q 3 + q 4 = 0

k ∆ y ∙ 1 ∆ x T m - 1 , n - T m , n + k ∆ x 2 ∙ 1 ∆ y T m , n - 1 - T m , n + k ∆ x 2 ∙ 1 ∆ y T m , n + 1 - T m , n + q 3 = 0

k ∆ y ∙ 1 ∆ x T m - 1 , n - T m , n + k ∆ x 2 ∙ 1 ∆ y T m , n - 1 - T m , n + k ∆ x 2 ∙ 1 ∆ y T m , n + 1 - T m , n + q ' ' ∆ y ∙ 1 = 0

Considerando ∆ x = ∆ y :

k T m - 1 , n - T m , n + k 2 T m , n - 1 - T m , n + k 2 T m , n + 1 - T m , n + q ' ' ∆ y ∙ 1 = 0

2 T m - 1 , n - T m , n + T m , n - 1 - T m , n + T m , n + 1 - T m , n + 2 q ' ' ∆ y k = 0

2 T m - 1 , n + T m , n - 1 + T m , n + 1 - 4 T m , n = - 2 q ' ' ∆ y k

Resposta

(a) 2 T m - 1 , n + T m , n - 1 + T m , n + 1 - 4 T m , n = 0 .

(b) 2 T m - 1 , n + T m , n - 1 + T m , n + 1 - 4 T m , n = - 2 q ' ' ∆ y k .

Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados

Um cilindro sólido redondo, feito com um material plástico ( α=6 × 10 - 7 m 2 / s ) , está inicialmente a uma temperatura uniforme de 20 ° C e se encontra isolado ao longo de sua superfície lateral e
Seja o Problema 5.9, exceto pelo fato de que agora o volume combinado do banho de óleo e da esfera é igual a V=1m 3 . O banho de óleo é bem misturado e isolado termicamente. (a) Considerando que as pr
Uma parede plana ( ρ=4000k g / m 3,c p=500J / ( k g·K ) ,k=10W / ( m·K ) ) , com espessura L=20m m , tem inicialmente uma distribuição de temperaturas estacionária linear, com as superfícies mantidas
A superfície superior de uma placa, incluindo os seus sulcos, é mantida a uma temperatura uniforme T 1=200 ° C . A superfície inferior se encontra a T 2=20 ° C . A condutividade térmica da placa é de
Uma fábula antiga descreve como um navio de madeira foi destruído por soldados que refletiram luz do sol nos seus escudos de bronze polido sobre o casco, fazendo o navio ficar em chamas. Para testar a