Soluções Tabeladas para Condução Transiente

Ver Todos os Exercícios

Teoria

Introdução

Um dos principais objetivos quando estudamos a condução de calor é determinar a distribuição de temperaturas no interior de um corpo.

Para casos em que B i < 0,1 , podemos usar o Método dos Parâmetros Concentrados. Quando B i > 0,1 , temos a opção de usar o Método dos Efeitos Espaciais.

No entanto, quando estamos trabalhando com placas planas, podemos usar soluções tabeladas para resolver a equação de calor.

Vamos pensar na seguinte situação:

Uma placa de espessura L , confinada a uma região 0 ≤ x ≤ L , está inicialmente a uma temperatura uniforme T i . Em t > 0 , a superfície x = 0 é conservada isolada e a superfície em x = L é mantida a uma temperatura constante T ∞ .

O nosso objetivo aqui é desenvolver a expressão da distribuição transiente de temperatura na placa.

Formulação Matemática do Problema

A equação de calor é dada por:

∂ ∂ x k ∂ T ∂ x + ∂ ∂ y k ∂ T ∂ y + ∂ ∂ z k ∂ T ∂ z + q ˙ = ρ c p ∂ T ∂ t

Para representar o nosso problema, vamos fazer as seguintes considerações:

  • Transferência de calor unidimensional;
  • Transferência de calor sem geração térmica q ˙ = 0 ;
  • Propriedades térmicas constantes.

Dessa forma, a nossa equação de calor fica com essa cara:

∂ 2 T x , t ∂ x 2 = 1 α . ∂ T x , t ∂ t

Em que α = k / ρ c p é a difusividade térmica do material.

Para que as condições de contorno sejam homogêneas, vamos definir uma nova temperatura θ x , t , dada por:

θ x , t = T x , t - T ∞

Dessa forma, a equação de calor terá essa cara:

∂ 2 θ x , t ∂ x 2 = 1 α . ∂ θ x , t ∂ t

As condições de contorno serão de superfície isolada em x = 0 , e superfície a temperatura constante em x = L . Dessa forma:

∂ θ x , t ∂ x = 0 e m   x = 0

θ x , t = 0 e m   x = L

E a condição inicial será dada por:

θ x , t = T i - T ∞ = θ i e m   t = 0

Soluções Tabeladas

A equação de calor que encontramos é resolvida através do método de separação de variáveis. O mesmo pressupõe que a função θ x , t pode ser representada como o produto de duas funções que dependem de x e t , separadamente:

θ x , t = ψ x . Γ t

Aplicando esse método e resolvendo a equação de calor, podemos representar a temperatura θ x , t como:

Tá, eu sei que assusta, então vamos entender quem é quem! A função F x corresponde à distribuição inicial de temperatura.

No nosso caso, como a placa está inicialmente a uma temperatura uniforme T i :

F x = θ i

Além disso:

  • ψ λ n , x : Autofunções;
  • λ n : Autovalores;
  • 1 / N : Integral de normalização.

Esses três caras são justamente o que a nossa tabela irá fornecer. Assim, pra determinar esses caras, basta saber as condições de contorno, se liga:

Onde H = h / k .

Na tabela, temos as seguintes condições de contorno:

  • ψ = 0 : Temperatura constante;
  • d ψ / d x = 0 : Superfície isolada;
  • - k . d ψ / d x + h ψ = 0 : Convecção na superfície.

Voltando ao nosso problema...

No nosso caso, temos uma condição de superfície isolada em x = 0 d ψ / d x = 0 e temperatura constante em x = L ψ = 0 .

Isso corresponde ao caso 4 da tabela. Dessa forma, a autofunção será:

Onde λ n é o autovalor dado como solução da seguinte equação:

cos ⁡ λ n L = 0

O cosseno será nulo para quando o ângulo assumir algum dos seguintes valores:

λ n L = π 2 ,   3 π 2 ,   5 π 2 ,   7 π 2 ,   9 π 2 , …

Logo, podemos dizer que:

λ n L = 2 n - 1 π 2

n = 1 ,   2 ,   3 ,   …

E a integral de normalização será:

Dessa forma, substituindo tudo isso na expressão de θ x , t , teremos:

θ x , t = ∑ n = 1 ∞ 1 N . exp ⁡ - α . λ n 2 t . ψ λ n , x . ∫ 0 L F x . ψ λ n , x . d x

θ x , t = ∑ n = 1 ∞ 2 L . exp ⁡ - α . λ n 2 t . cos ⁡ λ n x . ∫ 0 L θ i . cos ⁡ λ n x . d x

θ x , t = 2 L . θ i . ∑ n = 1 ∞ exp ⁡ - α . λ n 2 t . cos ⁡ λ n x . ∫ 0 L cos ⁡ λ n x . d x

Vamos resolver a integral separada pra não bagunçar demais:

I = ∫ 0 L cos ⁡ λ n x . d x = sen ⁡ λ n x λ n 0 L

I = 1 λ n sen ⁡ λ n L = 1 λ n sen ⁡ 2 n - 1 π 2 L . L

I = 1 λ n sen ⁡ 2 n - 1 π 2

Repare que, quando n for par, o seno será igual a - 1 e, quando n for ímpar, o seno será igual a 1 . Assim, podemos reescrever o resultado acima como:

I = 1 λ n - 1 n + 1

Substituindo na expressão da temperatura, teremos:

Onde:

λ n = 2 n - 1 π 2 L

Agora estamos prontos pra cair nos exercícios!

Soluções Tabeladas

Voltando ao nosso problema...

Exercícios Resolvidos

Exercício Resolvido #1

Considere condução de calor transiente e unidimensional em uma placa plana com propriedades constante e espessura L .

A placa, que se encontra inicialmente na temperatura uniforme T i , troca calor através de sua superfície em x = 0 por convecção com um meio na temperatura T ∞ com coeficiente de transferência de calor h .

A superfície em x = L é mantida na temperatura constante T L . Sabendo-se que a geração de energia na placa é nula, forneça para este problema:

i Formulação matemática.

i i Solução analítica.

Passo 1

i Formulação matemática.

Para chegar à formulação matemática, vamos partir da equação de calor, dada por:

∂ ∂ x k ∂ T ∂ x + ∂ ∂ y k ∂ T ∂ y + ∂ ∂ z k ∂ T ∂ z + q ˙ = ρ c p ∂ T ∂ t

Vamos fazer as seguintes considerações:

  • Condução de calor unidimensional;
  • Condução transiente;
  • Propriedades constantes;
  • Não há geração térmica.

Assim, a equação de calor pode ser escrita da seguinte forma:

k ∂ 2 T ∂ x 2 = ρ c p ∂ T ∂ t

∂ 2 T ∂ x 2 = ρ c p k ∂ T ∂ t

Para condição inicial, o problema nos disse que:

Como condições de contorno, o problema nos disse que a superfície em x = 0 troca calor por convecção e a superfície em x = L é mantida à temperatura constante T L .

Assim:

Passo 2

i i Solução analítica.

Para encontrar a solução analítica, vamos partir do princípio que o problema possui uma solução homogênea (a qual vamos calcular a partir das soluções tabeladas) e uma solução particular (não-homogênea, que teremos que resolver), que depende somente de x .

Assim:

T x , t = T h x , t + T p x , t

Aplicando isso à nossa formulação matemática, teremos:

∂ 2 T h ∂ x 2 + ∂ 2 T p ∂ x 2 = 1 α ∂ T h ∂ t - k d T h d x + h T h - k d T p d x + h T p = h T ∞ p a r a   x = 0 T h + T p = T L p a r a   x = L T h = T i p a r a   t = 0

Separando os dois problemas, teremos:

∂ 2 T h ∂ x 2 = 1 α ∂ T h ∂ t - k d T h d x + h T h = 0 p a r a   x = 0 T h = 0 p a r a   x = L T h = T i p a r a   t = 0

∂ 2 T p ∂ x 2 = 0 - k d T p d x + h T p = h T ∞ p a r a   x = 0 T p = T L p a r a   x = L

Passo 3

Começando pela solução homogênea:

∂ 2 T h ∂ x 2 = 1 α ∂ T h ∂ t - k d T h d x + h T h = 0 p a r a   x = 0 T h = 0 p a r a   x = L T h = T i p a r a   t = 0

Pelas condições de contorno, este corresponde ao caso 7 da nossa tabela, onde:

ψ λ n , x = sen ⁡ λ n L - x

1 N = 2 λ n 2 + H 2 L λ n 2 + H 2 + H

Sendo λ n as raízes da seguinte expressão:

λ n cot ⁡ λ n L = - H

Sendo H = h / k . Assim, a solução homogênea será dada por:

T h x , t = ∑ n = 1 ∞ 1 N . exp ⁡ - α . λ n 2 t . ψ λ n , x . ∫ 0 L F x . ψ λ n , x . d x

T h x , t = ∑ n = 1 ∞ 2 λ n 2 + H 2 L λ n 2 + H 2 + H exp ⁡ - α . λ n 2 t sen ⁡ λ n L - x ∫ 0 L T i sen ⁡ λ n L - x . d x

Agora é arrumar a casa:

T h x , t = 2 T i L ∑ n = 1 ∞ λ n 2 + H 2 λ n 2 + H 2 + H exp ⁡ - α . λ n 2 t sen ⁡ λ n L - x ∫ 0 L sen ⁡ λ n L - x . d x

Vamos resolver essa integral separadamente:

I = ∫ 0 L sen ⁡ λ n L - x . d x

I = cos ⁡ λ n L - x λ n 0 L

I = 1 - cos ⁡ λ n L λ n

Finalmente, a solução homogênea será:

T h x , t = 2 T i L ∑ n = 1 ∞ λ n 2 + H 2 λ n 2 + H 2 + H exp ⁡ - α . λ n 2 t sen ⁡ λ n L - x 1 - cos ⁡ λ n L λ n

Onde λ n é a raiz de:

λ n cot ⁡ λ n L = - H

A partir da expressão acima, podemos chegar à expressão de sen ⁡ λ n L da seguinte forma:

λ n cot ⁡ λ n L = - H → cot ⁡ λ n L = - H λ n

A partir dessa expressão da tangente, podemos pensar no seguinte triângulo retângulo:

Onde:

cos ⁡ λ n L = H H 2 + λ n 2

Finalmente:

Passo 4

Para a solução particular, temos que:

∂ 2 T p ∂ x 2 = 0 - k d T p d x + h T p = h T ∞ p a r a   x = 0 T p = T L p a r a   x = L

Integrando uma vez a equação duas vezes, temos que:

d T p d x = C 1

T p x = C 1 x + C 2

Aplicando à primeira condição de contorno, para x = 0 :

- k d T p d x + h T p = h T ∞

- k . C 1 + h . C 1 × 0 + C 2 = h T ∞

- k . C 1 + h . C 2 = h T ∞

C 2 = T ∞ + k h . C 1

Aplicando a segunda condição de contorno, para x = L :

C 1 L + C 2 = T L

Substituindo a expressão de C 2 que acabamos de encontrar:

C 1 L + T ∞ + k h . C 1 = T L

C 1 = T L - T ∞ h L + k

A constante C 2 será:

C 2 = T L - C 1 L

C 2 = T L - T L - T ∞ h L + k L

Dessa forma, a solução particular será:

T p x = C 1 x + C 2

T p x = T L - T ∞ h L + k x + T L - T L - T ∞ h L + k L

Passo 5

Finalmente, a solução analítica do problema será:

T x , t = T h x , t + T p x , t

Substituindo as soluções:

Resposta

i

∂ 2 T ∂ x 2 = 1 α ∂ T ∂ t

T x , t = T i p a r a   t = 0

- k d T d x + h T = h T ∞ p a r a   x = 0

T x , t = T L p a r a   x = L

i i

T x , t = 2 T i L ∑ n = 1 ∞ λ n 2 + H 2 λ n 2 + H 2 + H exp ⁡ - α . λ n 2 t sen ⁡ λ n L - x 1 λ n H H 2 + λ n 2 + T L + T L - T ∞ h L + k   x - L

Exercício Resolvido #2

(Adaptado) Um chip, com base quadrada de lado L = 5   m m e espessura t = 1   m m , é montado em um substrato. O substrato comporta-se como um isolante ideal. O chip é resfriado por convecção forçada de ar através da sua superfície superior, que é a única não isolada pelo substrato.

A temperatura do ar é 20 ° C e o coeficiente de transferência de calor é 150   W / m 2 K . O chip encontra-se inicialmente desligado e em equilíbrio com o ar. Quando o chip é ligado, energia é gerada uniformemente em seu interior com uma taxa volumétrica q ˙ = 9 × 10 6   W / m 3 .

Forneça:

i Formulação matemática do problema.

i i Solução analítica.

Passo 1

i Formulação matemática do problema.

Para chegar à formulação matemática, vamos partir da equação de calor, dada por:

∂ ∂ x k ∂ T ∂ x + ∂ ∂ y k ∂ T ∂ y + ∂ ∂ z k ∂ T ∂ z + q ˙ = ρ c p ∂ T ∂ t

Vamos fazer as seguintes considerações:

  • Condução de calor unidimensional;
  • Condução transiente;
  • Propriedades constantes;

Assim, a equação de calor pode ser escrita da seguinte forma:

k ∂ 2 T ∂ x 2 + q ˙ = ρ c p ∂ T ∂ t

∂ 2 T ∂ x 2 + q ˙ k = ρ c p k ∂ T ∂ t

Para condição inicial, o problema nos disse que o chip encontra-se inicialmente em equilíbrio com o ar, assim:

Como condições de contorno, o problema nos disse que uma das superfícies está termicamente isolada x = 0 , e há troca de calor por convecção com o ar na outra superfície x = L .

Assim:

Passo 2

i i Solução analítica.

Nós vamos dividir a solução em uma solução homogênea (que vamos calcular a partir das soluções tabeladas) e uma solução particular (não-homogênea, que teremos que resolver), que depende somente de x .

Assim:

T x , t = T h x , t + T p x , t

Aplicando isso à nossa formulação matemática, teremos:

∂ 2 T h ∂ x 2 + ∂ 2 T p ∂ x 2 + q ˙ k = 1 α ∂ T h ∂ t d T h d x + d T p d x = 0 p a r a   x = 0 k d T h d x + h T h + k d T p d x + h T p = h T ∞ p a r a   x = L T h = T ∞ p a r a   t = 0

Separando os dois problemas, teremos:

∂ 2 T h ∂ x 2 = 1 α ∂ T h ∂ t d T h d x = 0 p a r a   x = 0 k d T h d x + h T h = 0 p a r a   x = L T h = T ∞ p a r a   t = 0

∂ 2 T p ∂ x 2 + q ˙ k = 0 d T p d x = 0 p a r a   x = 0 k d T p d x + h T p = h T ∞ p a r a   x = L

Passo 3

Começando pela solução homogênea:

∂ 2 T h ∂ x 2 = 1 α ∂ T h ∂ t d T h d x = 0 p a r a   x = 0 k d T h d x + h T h = 0 p a r a   x = L T h = T ∞ p a r a   t = 0

Pelas condições de contorno, este corresponde ao caso 6 da nossa tabela, onde:

ψ λ n , x = cos ⁡ λ n x

1 N = 2 λ n 2 + H 2 L λ n 2 + H 2 + H

Sendo λ n as raízes da seguinte expressão:

λ n tg ⁡ λ n L = H

Sendo H = h / k . Assim, a solução homogênea será dada por:

T h x , t = ∑ n = 1 ∞ 1 N . exp ⁡ - α . λ n 2 t . ψ λ n , x . ∫ 0 L F x . ψ λ n , x . d x

T h x , t = ∑ n = 1 ∞ 2 λ n 2 + H 2 L λ n 2 + H 2 + H exp ⁡ - α . λ n 2 t cos ⁡ λ n x ∫ 0 L T ∞ cos ⁡ λ n x . d x

Agora é arrumar a casa:

T h x , t = 2 T ∞ L ∑ n = 1 ∞ λ n 2 + H 2 λ n 2 + H 2 + H exp ⁡ - α . λ n 2 t cos ⁡ λ n x ∫ 0 L cos ⁡ λ n x . d x

Vamos resolver essa integral separadamente:

I = ∫ 0 L cos ⁡ λ n x . d x

I = sen ⁡ λ n x λ n 0 L

I = sen ⁡ λ n L λ n

Finalmente, a solução homogênea será:

T h x , t = 2 T ∞ L ∑ n = 1 ∞ λ n 2 + H 2 λ n 2 + H 2 + H L exp ⁡ - α . λ n 2 t cos ⁡ λ n x sen ⁡ λ n L λ n

Onde λ n é a raiz de:

λ n tg ⁡ λ n L = H

A partir da expressão acima, podemos chegar à expressão de sen ⁡ λ n L da seguinte forma:

λ n tg ⁡ λ n L = H → tg ⁡ λ n L = H λ n

A partir dessa expressão da tangente, podemos pensar no seguinte triângulo retângulo:

Onde:

sen ⁡ λ n L = H H 2 + λ n 2

Finalmente:

Passo 4

Para a solução particular, temos que:

∂ 2 T p ∂ x 2 + q ˙ k = 0 d T p d x = 0 p a r a   x = 0 k d T p d x + h T p = h T ∞ p a r a   x = L

Integrando uma vez a equação duas vezes, temos que:

d T p d x = - q ˙ k x + C 1

T p x = - q ˙ 2 k x 2 + C 1 x + C 2

Aplicando à primeira condição de contorno, para x = 0 :

d T p d x = 0

- q ˙ k × 0 + C 1 = 0

C 1 = 0

Assim:

d T p d x = - q ˙ k x

T p x = - q ˙ 2 k x 2 + C 2

Aplicando a segunda condição de contorno, para x = L :

k d T p d x + h T p = h T ∞

k - q ˙ k L + h - q ˙ 2 k L 2 + C 2 = h T ∞

C 2 = T ∞ + q ˙ L h + q ˙ L 2 2 k

Substituindo a expressão de C 2 que acabamos de encontrar:

T p x = - q ˙ 2 k x 2 + T ∞ + q ˙ L h + q ˙ L 2 2 k

Passo 5

Finalmente, a solução analítica do problema será:

T x , t = T h x , t + T p x , t

Substituindo as soluções:

Resposta

i

∂ 2 T ∂ x 2 + q ˙ k = 1 α ∂ T ∂ t T x , t = T ∞ t = 0 d T d x = 0 p a r a   x = 0 k d T d x + h T = h T ∞ p a r a   x = L

i i

T x , t = 2 T ∞ L ∑ n = 1 ∞ λ n 2 + H 2 λ n 2 + H 2 + H . exp ⁡ - α . λ n 2 t . cos ⁡ λ n x . 1 λ n . H H 2 + λ n 2 + q ˙ 2 k L 2 - x 2 + q ˙ L h + T ∞

Exercício Resolvido #3

Uma placa de aço ( α = 1,7 × 10 - 7   m 2 / s ) , com espessura L = 20   m m , encontra-se inicialmente na temperatura uniforme T i = 240 ° C . As duas superfícies da placa são então mantida na temperatura T s = 40 ° C . Para este problema, considerando propriedades físicas constantes, forneça:

i Formulação matemática;

i i Solução analítica;

i i i Temperatura do centro da placa 10 minutos depois de iniciado o resfriamento.

Passo 1

i Formulação matemática.

Para chegar à formulação matemática, vamos partir da equação de calor, dada por:

∂ ∂ x k ∂ T ∂ x + ∂ ∂ y k ∂ T ∂ y + ∂ ∂ z k ∂ T ∂ z + q ˙ = ρ c p ∂ T ∂ t

Vamos fazer as seguintes considerações:

  • Condução de calor unidimensional;
  • Condução transiente;
  • Propriedades constantes;
  • Não há geração térmica.

Assim, a equação de calor pode ser escrita da seguinte forma:

k ∂ 2 T ∂ x 2 = ρ c p ∂ T ∂ t

∂ 2 T ∂ x 2 = ρ c p k ∂ T ∂ t

Para condição inicial, o problema nos disse que:

Como condições de contorno, o problema nos disse que as superfícies são mantidas a 40 ° C . Assim:

Passo 2

i i Solução analítica;

Nesse caso, como as temperaturas em ambos os lados da placa são iguais, podemos transformar esse problema em homogêneo utilizando θ x , t , tal que:

θ x , t = T x , t - T s

Fazendo isso, a formulação matemática passa a ser dada por:

∂ 2 θ x , t ∂ x 2 = 1 α . ∂ θ x , t ∂ t θ x , t = 0 p a r a   x = 0 θ x , t = 0 p a r a   x = L = 20   m m θ x , t = θ i = 200 ° C p a r a   t = 0

Pelas condições de contorno, este corresponde ao caso 1 da nossa tabela, onde:

ψ λ n , x = sen ⁡ λ n x

1 N = 2 L

Sendo λ n as raízes da seguinte expressão:

sen ⁡ λ n L = 0

Assim, a solução homogênea será dada por:

θ x , t = ∑ n = 1 ∞ 1 N . exp ⁡ - α . λ n 2 t . ψ λ n , x . ∫ 0 L F x . ψ λ n , x . d x

θ x , t = ∑ n = 1 ∞ 2 L exp ⁡ - α . λ n 2 t sen ⁡ λ n x ∫ 0 L θ i sen ⁡ λ n x . d x

Agora é arrumar a casa:

θ x , t = 2 θ i L ∑ n = 1 ∞ exp ⁡ - α . λ n 2 t sen ⁡ λ n x ∫ 0 L sen ⁡ λ n x . d x

Vamos resolver essa integral separadamente:

I = ∫ 0 L sen ⁡ λ n x . d x

I = - cos ⁡ λ n x λ n 0 L

I = 1 - cos ⁡ λ n L λ n

Finalmente, a solução homogênea será:

θ x , t = 2 θ i L ∑ n = 1 ∞ exp ⁡ - α . λ n 2 t sen ⁡ λ n x 1 - cos ⁡ λ n L λ n

Onde λ n é a raiz de:

sen ⁡ λ n L = 0

O seno é nulo para qualquer múltiplo de π . Assim:

λ n L = n π → λ n = n π L

Dessa forma, substituindo na expressão da solução homogênea:

θ x , t = 2 θ i L ∑ n = 1 ∞ exp ⁡ - α n 2 π 2 t L 2 sen ⁡ n π x L 1 - cos ⁡ n π . L n π

Repare que, para valores pares de n , 1 - cos ⁡ n π = 0 e, para valores ímpares, 1 - cos ⁡ n π = 2 . Assim, podemos levar somente em consideração os valores ímpares de n , substituindo o n na equação por 2 n - 1 . Dessa forma:

Substituindo de volta a temperatura T x , t :

Passo 3

i i i Temperatura do centro da placa 10 minutos depois de iniciado o resfriamento.

Como queremos a temperatura no centro da placa, vamos considerar x = 10   m m = 0,01   m . Assim, substituindo os valores, onde:

t = 10   m i n = 600   s

L = 20   m m

T i = 240 ° C

T s = 40 ° C

α = 1,7 × 10 - 7   m 2 / s

Considerando apenas o primeiro termo da série n = 1 , teremos:

T 0,01 ; 600 = 40 + 2 240 - 40 π × exp ⁡ - 1,7 × 10 - 7 × 2 × 1 - 1 2 π 2 × 600 0,02 2 × sen ⁡ 2 × 1 - 1 π × 0,01 0,02 × 2 2 × 1 - 1

Jogando isso tudo na calculadora:

Resposta

i

∂ 2 T ∂ x 2 = 1 α ∂ T ∂ t

T x , t = 240 ° C p a r a   t = 0

T x , t = 40 ° C p a r a   x = 0

T x , t = 40 ° C p a r a   x = L = 20   m m

i i

T x , t = T s + 2 T i - T s π ∑ n = 1 ∞ exp ⁡ - α 2 n - 1 2 π 2 t L 2 sen ⁡ 2 n - 1 π x L 2 2 n - 1  

i i i T 0,01 ; 600 = 60,55 ° C

Exercícios de Livros Relacionados

O Problema 4.9 tratou do armazenamento de lixo radioativo no subsolo em recipientes esféricos. Em função das incertezas nas propriedades térmicas do solo, deseja-se medir a temperatura em regime estac
Deduza uma expressão para a razão da energia total transferida a partir da superfície isotérmica de um cilindro infinito para o seu interior, Q / Q o , que seja válida para F o0,2 . Represente os seus
Os componentes estruturais de uma aeronave moderna são normalmente fabricados com materiais compósitos de alto desempenho. Esses materiais são fabricados pela impregnação de malhas de fibras extremame
Bombas de calor, com base no subsolo, operam usando o solo, em vez do ar ambiente, como a fonte de calor (ou sumidouro) para aquecimento (ou resfriamento) de construções. Um líquido troca energia com
86 – Durante os meses de inverno frio, um escritório é parcialmente aquecido extraindo energia térmica armazenado no solo abaixo, usando uma bomba de calor. Durante o verão, o prédio é parcialmente re