Um dos pontos fortes dos métodos numéricos é a sua capacidade de lidar com condições de contorno complexas. No esboço, a condição de contorno muda de fluxo térmico especificado, q s ' ' (para dentro do domínio), para convecção, na posição do nó (m,n). Escreva a equação de diferenças finitas bidimensional, em regime estacionário, neste nó.

Passo 1

Vamos começar desenhando as taxas de transferência de calor no volume de controle:

Vamos agora aplicar o balanço de energia a esse volume de controle em 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 3 + q 4 + q s ' ' A s

Passo 2

Vamos agora calcular separadamente cada taxa de transferência de calor:

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

q 1 = k A ∆ T L

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

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

q 2 = k A L ∆ T

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

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

q 3 = k A L ∆ T

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

q 4 = q c o n v , y

q 4 = h A ∆ T

q 4 = h ∆ x 2 ∙ 1 T ∞ - T m , n

q 5 = q s ' ' A s

q 5 = q s ' ' ∆ x 2 ∙ 1

Passo 3

Substituindo no balanço de energia:

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

Para ∆ x = ∆ y :

k 2 T m - 1 , n - T m , n + k T m , n - 1 - T m , n + k 2 T m + 1 , n - T m , n + h ∆ x 2 T ∞ - T m , n + q s ' ' ∆ x 2 = 0

T m - 1 , n - T m , n + 2 T m , n - 1 - T m , n + T m + 1 , n - T m , n + h ∆ x k T ∞ - T m , n + q s ' ' ∆ x k = 0

T m - 1 , n + 2 T m , n - 1 + T m + 1 , n + h ∆ x k T ∞ - 4 + h ∆ x k T m , n = - q s ' ' ∆ x k

Resposta

T m - 1 , n + 2 T m , n - 1 + T m + 1 , n + h ∆ x k T ∞ - 4 + h ∆ x k T m , n = - q s ' ' ∆ x k

Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados