Enunciado

Seja o elemento combustível do Exemplo 5.11 . Inicialmente, o elemento está a uma temperatura uniforme de 250 ° C sem geração de energia. Subitamente, o elemento é inserido no núcleo do reator passando a ter uma taxa volumétrica de geração de energia uniforme d e   q ˙ =   108   W / m 3 . As superfícies são resfriadas por convecção com T ∞   =   250 ° C e h   =   1100   W / ( m 2   ·   K ) . Usando o método explícito com um incremento espacial de 2   m m , determine a distribuição de temperaturas 1,5   s após o elemento ser inserido no núcleo do reator.

Passo 1

Vamos lá! Nosso problema pede a distribuição de temperatura pra um incremento de 2 m m , mas antes vamos fazer um esboço para facilitar nossa vida.

Agora precisamos fazer nossas famosas considerações. Hehe (u.u)

I - C o n d u ç ã o   t r a n s i e n t e   u n i d i m e n s i o n a l

I I - P r o p r i e d a d e s   c o n s t a n t e s

I I I - t = 0 , q ˙ = 0 ,   p a r a   t > 0 , q ˙   é   u n i f o r m e

Passo 2

Vamos lá, com nossos dados em mãos vamos calcular a distribuição de temperatura pelo método explicito, como foi sugerido pelo problema, assim iremos usar as mesmas equações de diferença finita para nossa rede nodal do exemplo 5.11 . Então temos:

Nos nós interiores, m = 1 ,   2 ,   3 ,   4 .

T m p + 1 = F 0 T m - 1 p + T m + 1 p + q ˙ Δ x 2 k + 1 - 2 F 0 T m p

No nó central, m = 0 .

T 0 p + 1 = F 0 2 T 1 p + q ˙ Δ x 2 k + 1 - 2 F 0 T 0 p

E no nó de superfície, m = 5 .

T 5 p + 1 = 2 F 0 T 4 p + B i T ∞ + q ˙ Δ x 2 k + 1 - 2 F 0 - 2 B i F 0 T 5 p

Passo 3

Agora temos que o critério de estabilidade mais rígido se encontra nos nós de superfície em ambos os lados, pois 1 - 2 F 0 - 2 B i F 0 não pode ser menor que 0 , assim:

1 - 2 F 0 - 2 B i F 0 = 1 - 2 F 0 1 + B i ≥ 0

Organizando

F 0 1 + B i ≤ 1 2

Calculando o número de B i o t temos:

B i = h Δ x k = 1100 W m 2 ⋅ K ⋅ 0,002   m 30   W m ⋅ K = 0,0733

Assim:

F 0 ≤ 1 2 1 + B i → F 0 ≤ 1 2 1 + 0,0733

F 0 ≤ 0,466

Como F 0 = α Δ t Δ x 2 , vamos encontrar agora os intervalos de tempo, assim temos:

Δ t ≤ F 0 Δ x 2 α → Δ t ≤ 0,466 0,002 m 2 5   ⋅ 10 - 6 m 2 s = 0,373   s

Para não trabalharmos com números muito quebrados e ficarmos bem com nosso critério de estabilidade, vamos arredondar Δ t = 0,3   s , Assim:

F 0 = α Δ t Δ x 2 = 5   ⋅ 10 - 6 m 2 s 0,373   s 0,002 m 2 = 0,375

Por fim o tempo será dado por:

t = Δ t ⋅ p = 0,3 ⋅ p s

Passo 4

Agora vamos substituir todos os dados que encontramos nas nossas equações dos nós, com q ˙ = 100 M W m 3 , assim temos:

T 0 p + 1 = 0,375 2 T 1 p + 10 8   W m 3 0,002 2 30   W m ⋅ K + 1 - 2 ⋅ 0,375 T 0 p

T 0 p + 1 = 0,375 2 T 1 p + 13,333 + 0,25 T 0 p

Fazendo o mesmo processo para as demais temos:

T 1 p + 1 = 0,375 T 0 p + T 2 p + 13,333 + 0,25 T 1 p

T 2 p + 1 = 0,375 T 1 p + T 3 p + 13,333 + 0,25 T 2 p

T 3 p + 1 = 0,375 T 2 p + T 4 p + 13,333 + 0,25 T 3 p

T 4 p + 1 = 0,375 T 3 p + T 5 p + 13,333 + 0,25 T 4 p

T 5 p + 1 = 2 ⋅ 0,375 T 4 p + 0,0733 ⋅ 250 + 13,333 2 + 1 - 2 ⋅ 0,375 - 2 ⋅ 0,0733 ⋅ 0,375 T 5 p

T 5 p + 1 = 0,75 T 4 p + 24 , 99 + 0 , 19 5 T 5 p

Passo 5

Agora vamos fazer uma tabela com nossos dados para 5 intervalos, como a temperatura inicial é 250 ° C temos:

P.S: vale ressaltar que nos instantes iniciais do processo, o plano central não vai sentir o efeito do ar refrigerante.

Isso é tudo pessoal! \o/

Resposta

Ei, a resposta está no passo a passo :)

Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados

Seja a transferência de calor unidimensional (radial) em um sistema de coordenadas cilíndricas em regime estacionário com geração volumétrica de calor.(a) Deduza a equação de diferenças finitas para q
Seja a configuração nodal 4 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 superior do vértice externo está per
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 d
103 – Considere a nacela da turbina eólica no Exemplo 1.3. A nacela é formada de um material compósito termoplástico de espessura L=20 mm. A condutividade térmica, a densidade, e o calor especifico do
95 – Considere o coletor solar do Problema 3.89. Depois de alcançar uma operação em condição de estado estacionário com q”rad=800 W/m2, o fluxo de radiação liquido para a superfície é rapidamente redu